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.

We use K-giant stars selected from the LAMOST DR5 to study the variation of the rotational velocity of the galactic halo at different space positions. Modelling the rotational velocity distribution with both the halo and disk components, we find that the rotational velocity of the halo population decreases almost linearly with increasing vertical distance to the galactic disk plane, $Z$, at fixed galactocentric radius, $R$. The samples are separated into two parts with $6<R<12$ kpc and $12<R<20$ kpc. We derive that the decreasing rates along $Z$ for the two subsamples are $-3.07\pm0.63$ and $-1.89\pm0.37$ km s$^{-1}$ kpc$^{-1}$, respectively. Compared with the TNG simulations, we suggest that this trend is probably caused by the interaction between the disk and halo. The results from the simulations show that only the oblate halo can provide a decreasing rotational velocity with an increasing $Z$. This indicates that the Galactic halo is oblate with galactocentric radius $R<20$ kpc. On the other hand, the flaring of the disk component (mainly the thick disk) is clearly traced by this study, with $R$ between 12 and 20 kpc, the disk can vertically extend to $6\sim10$ kpc above the disk plane. What is more interesting is that, we find the Gaia-Enceladus-Sausage (GES) component has a significant contribution only in the halo with $R<12$ kpc, i.e. a fraction of 23$-$47\%. While in the outer subsample, the contribution is too low to be well constrained.

QPOs are seen as peak features in the X-ray power spectral density of stellar mass black holes and neutron stars, of which frequencies can be measured with high precision. These QPO frequencies are believed to be related to fundamental frequencies of test particles, which are mainly determined by the background metric. We consider the metric introduced in Ghasemi-Nodehi and Bambi (Eur Phys J C 76:290, 2016). The fundamental frequencies in this metric are sensitive to some of the metric parameters but insensitive to other parameters, which means, the differences in fundamental frequencies in this metric and those in the Kerr ones can be significant for small changes of some but not all parameters around the Kerr value. By comparing with the QPO observations of GRO J1655-40, we find that only one parameter of the Ghasemi-Nodehi-Bambi metric can be strongly constrained, but other parameters cannot. We also use nested algorithm to investigate whether better constraints on the Ghasemi-Nodehi-Bambi metric parameters can be obtained from QPO observations of multiple objects by simulations. We find that four parameters can be strongly constrained while other parameters cannot . Our results suggest that QPOs may be important tools for testing the Kerr metric.

We report the discovery of a mid-infrared (MIR) flare using WISE data in the center of the nearby Seyfert 1.9 galaxy MCG-02-04-026. The MIR flare began in the first half of 2014, peaked around the end of 2015, and faded in 2017. During these years, energy more than $7\times10^{50}$ erg was released in the infrared, and the flare's MIR color was generally turning red. We detected neither optical nor ultraviolet (UV) variation corresponding to the MIR flare based on available data. We explained the MIR flare using a dust echo model in which the radiative transfer is involved. The MIR flare can be well explained as thermal reradiation from dust heated by UV-optical photons of a primary nuclear transient event. Although the transient event was not seen directly due to dust obscuration, we can infer that it may produce a total energy of at least $\sim10^{51}$ erg, most of which was released in less than $\sim$3 years. The nature of the transient event could be a stellar tidal disruption event by the central supermassive black hole (SMBH), or a sudden enhancement of the existing accretion flow onto the SMBH, or a supernova which was particularly bright.

(abridged) We present ALMA observations of hydrogen fluoride, HF J=1-0, H20 (220-211), and the 1.2 THz rest-frame continuum emission from the z=4.7 system BR1202-0725. BR1202-0725 is a galaxy group consisting of a QSO, a sub-millimeter galaxy (SMG), and two Ly-alpha emitters. We detected HF in emission in the QSO and possibly in absorption in the SMG, while water is detected in emission in both the QSO and SMG. The QSO is the most luminous HF emitter yet found and has the same ratio of HF emission line to infrared luminosity as a sample of local AGN and the Orion Bar. This consistency covers about 10 orders-of-magnitude in infrared luminosity, L_IR. Based on the conclusions of a study of HF emission in the Orion Bar and modeling, the HF emission in the QSO is either excited by collisions with electrons and H2 in molecular plasmas irradiated by the AGN and intense star formation or predominately by collisions with H2, with a modest contribution from electrons, in a relatively high temperature (~120 K), dense (~10^5 cm^-3) medium. Although HF should be an excellent tracer of molecular outflows, we do not find strong evidence for outflows in HF in either the QSO or the SMG. From a putative absorption feature in HF in the SMG, we estimate an upper limit on the outflow rate, dM/dt_outflow <~45 M_sun/yr. The ratio of the outflow rate to the star formation rate is <5% for the SMG. The broadness of the H2O line in the SMG, FWHM~1020 km/s, may suggest that either the gas on large scales (>4 kpc) is significantly more disturbed and turbulent due either to interactions and mass exchange with the other members of the group or to the dissipation of the energy of the intense star formation or both. The lack of significant molecular outflows in either source may imply that much of the energy from the intense star formation and AGN activity in this pair is being dissipated in their ISM.

Instruments flown on CubeSats are small. Meaningful applications of CubeSats in astronomical observations rely on the choice of a particular subject that is feasible for CubeSats. Here we report the result of a feasibility study for observing gamma-ray polarization from Cygnus X-1 using a small Compton polarimeter on board a 3U CubeSat. Silicon detectors and cerium bromide scintillators were employed in the instrument models that we discussed in this study. Through Monte Carlo simulations with Geant4-based MEGAlib package, we found that, with a 10-Ms on-axis, zenith-direction observation in a low inclination, low altitude earth orbit radiation background environment, the minimum detectable polarization degree can be down to about 10\% in 160 - 250 keV, 20\% in 250 - 400 keV, and 65\% in 400 - 2000 keV. A 3U CubeSat dedicated to observing Cygnus X-1 can therefore yield useful information on the polarization state of gamma-ray emissions from the brightest persistent X-ray black-hole binary in the sky.

We present 5.5 GHz observations with the VLA of a sample of nearby galaxies with energetic nuclear outbursts at mid-infrared (MIR) bands. These observations reach a uniform depth down to a median rms of ~10 uJy, representing one of most sensitive searches for radio emission associated with nuclear transients. We detect radio emission in 12 out of 16 galaxies at a level of >5sigma, corresponding to a detection rate of 75%. Such a high detection is remarkably different from previous similar searches in stellar tidal disruption events. The radio emission is compact and not resolved for the majority of sources on scales of ~<0.5" (<0.9 kpc at z<0.1). We find the possibility of the star-formation contributing to the radio emission is low, but an AGN origin remains a plausible scenario, especially for sources that show evidence of AGN activity in their optical spectra. If the detections could represent radio emission associated with nuclear transient phenomenon such as jet or outflow, we use the blast wave model by analogy with the GRB afterglows to describe the evolution of radio light curves. In this context, the observations are consistent with a decelerating jet with an energy of ~10^{51-52} erg viewed at 30\degree-60\degree off-axis at later times, suggesting that powerful jets may be ubiquitous among MIR-burst galaxies. Future continuous monitoring observations will be crucial to decipher the origin of radio emission through detections of potential flux and spectral evolution. Our results highlight the importance of radio observations to constrain the nature of nuclear MIR outbursts in galaxies.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2''$ to $1.5''$ ($1'' \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18'' (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius.

The Jiangmen Underground Neutrino Observatory is proposed to determine neutrino mass hierarchy using a 20~ktonne liquid scintillator detector. Strict radio-purity requirements have been put forward for all the components of the detector. According to the MC simulation results, the radon dissolved in the water Cherenkov detector should be below 200~mBq/m$^3$. Radium, the progenitor of radon, should also be taken seriously into account. In order to measure the radium concentration in water, a radium measurement system, which consists of a radium extraction system, a radon emanation chamber and a radon concentration measurement system, has been developed. In this paper, the updated radon concentration in gas measurement system with a one-day-measurement sensitivity of $\sim$5~mBq/m$^3$, the detail of the development of the radium concentration in water measurement system with a sensitivity of $\sim$23~mBq/m$^3$ as well as the measurement results of Daya Bay water samples will be presented.

We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $\beta\sim0.9$ to 0.7 for K giants and from $\beta\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<\beta<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system.

We study bottom flavor abundance in the early Universe near to a temperature $T_\mathrm{H}\simeq150\,\mathrm{MeV}$, the condition for hadronization of deconfined quark-gluon plasma (QGP). We show bottom flavor abundance nonequilibrium lasting microseconds. In our study we use that in both QGP, and the hadronic gas phase (HG) $b$ and $\bar b$ quarks near $T_\mathrm{H}$ are bound in B-mesons and antimesons subject to $CP$ violating weak decays. A coincident non-equilibrium abundance of bottom flavor can lead to matter genesis at required strength: a) The specific thermal yield per entropy is $n_b^{th}/\sigma=10^{-10}\sim10^{-13}$. b) Considering time scales, millions of cycles of B-meson decays, and $b\bar b$-pair recreation processes occur.

Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter search to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The "standard" EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition to it, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) range. The first is due to bremsstrahlung of electrons scattered on neutral atoms (so-called neutral bremsstrahlung, NBrS). The second, responsible for electron avalanche scintillations in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a SiPM matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillations effects. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter search. It should be emphasized that this is the first practical application of the NBrS effect in detection science.

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars.

The potential association between the blazar TXS 0506+056 and the neutrino event IceCube-170922A provides a unique opportunity to study the possible physical connection between the high-energy photons and neutrinos. We explore the correlated electromagnetic and neutrino emissions of blazar TXS 0506+056 by a self-consistent leptonic-hadronic model, taking into account particle stochastic acceleration and all relevant radiative processes self-consistently. The electromagnetic and neutrino spectra of blazar TXS 0506+056 are reproduced by the proton synchrotron and hybrid leptonic-hadronic models based on the proton-photon interactions. It is found that the hybrid leptonic-hadronic model can be used to better explain the observed X-ray and $\gamma$-ray spectra of blazar TXS 0506+056 than the proton synchrotron model. Moreover, the predicted neutrino spectrum of the hybrid leptonic-hadronic model is closer to the observed one compared to the proton synchrotron model. We suggest that the hybrid leptonic-hadronic model is more favored if the neutrino event IceCube-170922A is associated with the blazar TXS 0506+056.

We investigate the formation mechanism for the observed nearly polar aligned (perpendicular to the binary orbital plane) debris ring around the eccentric orbit binary 99 Herculis. An initially inclined nonpolar debris ring or disc will not remain flat and will not evolve to a polar configuration, due to the effects of differential nodal precession that alter its flat structure. However, a gas disc with embedded well coupled solids around the eccentric binary may evolve to a polar configuration as a result of pressure forces that maintain the disc flatness and as a result of viscous dissipation that allows the disc to increase its tilt. Once the gas disc disperses, the debris disc is in a polar aligned state in which there is little precession. We use three-dimensional hydrodynamical simulations, linear theory, and particle dynamics to study the evolution of a misaligned circumbinary gas disc and explore the effects of the initial disc tilt, mass, and size. We find that for a wide range of parameter space, the polar alignment timescale is shorter than the lifetime of the gas disc. Using the observed level of alignment of 3 deg. from polar, we place an upper limit on the mass of the gas disc of about 0.014 M_sun at the time of gas dispersal. We conclude that the polar debris disc around 99 Her can be explained as the result of an initially moderately inclined gas disc with embedded solids. Such a disc may provide an environment for the formation of polar planets.

We analyze the X-ray, optical, and mid-infrared data of a "changing-look" Seyfert galaxy \sdssj15 at $z\simeq0.086$. Over a period of one decade (2009 - 2018), its broad H$\alpha$ line intensity increased by a factor of $\sim$4. Meanwhile, the X-ray emission in 2014 as observed by \chandra\ was about five times brighter than that in 2010 by {\it Suzaku}, and the corresponding emissions in V-band, mid-infrared W1 band brighten by $\sim$ 0.18, 0.32 mag, respectively. Moreover, the absorption in X-rays is moderate and stable, i.e. ${\rm N_{H}}\sim 10^{21} {\rm cm^{-2}}$, but the X-ray spectrum becomes harder in the 2014 \chandra\ bright state (i.e. photon index $\Gamma = 1.52^{+0.06}_{-0.06}$) than that of the 2010 \suzaku\ low state ($\Gamma=2.03^{+0.22}_{-0.21}$). With an Eddington ratio being lower than a few percent, the inner region of the accretion disk in \sdssj15\ is likely a hot accretion flow. We then compile from literature the X-ray data of "changing-look" AGNs, and find that they generally follow the well-established "V"-shaped correlation in AGNs, that is, above a critical turn-over luminosity the X-ray spectra soften with the increasing luminosity, and below that luminosity the trend is reversed in a way of "harder when brighter". This presents a direct evidence that CL-AGNs have distinctive changes in not only the optical spectral type, but also the X-ray spectral shape. The similarity in the X-ray spectral evolution between CL-AGNs and black hole X-ray binaries indicates that the observed CL-AGNs phenomena may relate to the state transition in accretion physics.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $\beta$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $\gamma$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector.

Submillimeter rotational lines of H2O are a powerful probe in warm gas regions of the ISM, tracing scales and structures ranging from kpc disks to the most compact and dust-obscured regions of galactic nuclei. The ortho-H2O(423-330) line at 448 GHz, which was recently detected in a local luminous infrared galaxy (Pereira-Santaella et al. 2017), offers a unique constraint on the excitation conditions and ISM properties in deeply buried galaxy nuclei since the line requires high far-IR optical depths to be excited. In this letter, we report the first high-redshift detection of the 448 GHz H2O(423-330) line using ALMA, in a strongly lensed submillimeter galaxy (SMG) at z=3.63. After correcting for magnification, the luminosity of the 448 GHz H2O line is ~10^6 L_sun. In combination with three other previously detected H2O lines, we build a model that "resolves" the dusty ISM structure of the SMG, and find that it is composed of a ~1 kpc optically thin (optical depth at 100{\mu}m {\tau}_{100}~0.3) disk component with dust temperature T_{dust} \approx 50 K emitting a total infrared power of 5e12 L_sun with surface density \Sigma_{IR}=4e11 L_sun kpc^{-2}, and a very compact (0.1 kpc) heavily dust-obscured ({\tau}_{100} \gtrsim 1) nuclear core with very warm dust (100 K) and \Sigma_{IR}=8e12 L_sun kpc^{-2}. The H2O abundance in the core component, X_{H2O}~(0.3-5)e{-5}, is at least one order of magnitude higher than in the disk component. The optically thick core has the characteristic properties of an Eddington-limited starburst, providing evidence that radiation pressure on dust is capable of supporting the ISM in buried nuclei at high redshifts. The multi-component ISM structure revealed by our models illustrates that dust and molecules such as H2O are present in regions characterized by highly differing conditions and scales, extending from the nucleus to more extended regions of SMGs.

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

Impact craters, as "lunar fossils", are the most dominant lunar surface features and occupy most of the Moon's surface. Their formation and evolution record the history of the Solar System. Sixty years of triumphs in the lunar exploration projects accumulated a large amount of lunar data. Currently, there are 9137 existing recognized craters. However, only 1675 of them have been determined age, which is obviously not satisfactory to reveal the evolution of the Moon. Identifying craters is a challenging task due to their enormous difference in size, large variations in shape and vast presence. Furthermore, estimating the age of craters is extraordinarily difficult due to their complex and different morphologies. Here, in order to effectively identify craters and estimate their age, we convert the crater identification problem into a target detection task and crater age estimation into a taxonomy structure. From an initial small number of available craters, we progressively identify craters and estimate their age from Chang'E data by transfer learning (TL) using deep neural networks. For comprehensive identification of multi-scale craters, a two-stage craters detection approach is developed. Thus 117240 unrecognized lunar craters that range in diameter from 532 km to 1 km are identified. Then, a two-stage classification approach is developed to estimate the age of craters by simultaneously extracting their morphological features and stratigraphic information. The age of 79243 craters larger than 3 km in diameter is estimated. These identified and aged craters throughout the mid and low-latitude regions of the Moon are crucial for reconstructing the dynamic evolution process of the Solar System.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

Recent observations have revealed a gallery of substructures in the dust component of nearby protoplanetary discs, including rings, gaps, spiral arms, and lopsided concentrations. One interpretation of these substructures is the existence of embedded planets. Not until recently, however, most of the modelling effort to interpret these observations ignored the dust back reaction to the gas. In this work, we conduct local-shearing-sheet simulations for an isothermal, inviscid, non-self-gravitating, razor-thin dusty disc with a planet on a fixed circular orbit. We systematically examine the parameter space spanned by planet mass ($0.1M_\text{th} \le M_\text{p} \le 1M_\text{th}$, where $M_\text{th}$ is the thermal mass), dimensionless stopping time ($10^{-3} \le \tau_\text{s} \le 1$), and solid abundance ($0 < Z \le 1$). We find that when the dust particles are tightly coupled to the gas ($\tau_\text{s} < 0.1$), the spiral arms are less open and the gap driven by the planet becomes deeper with increasing $Z$, consistent with a reduced speed of sound in the approximation of a single dust-gas mixture. By contrast, when the dust particles are marginally coupled ($0.1 \lesssim \tau_\text{s} \lesssim 1$), the spiral structure is insensitive to $Z$ and the gap structure in the gas can become significantly skewed and unidentifiable. When the latter occurs, the pressure maximum radially outside of the planet is weakened or even extinguished, and hence dust filtration by a low-mass ($M_\text{p} < M_\text{th}$) planet could be reduced or eliminated. Finally, we find that the gap edges where the dust particles are accumulated as well as the lopsided large-scale vortices driven by a massive planet, if any, are unstable, and they are broken into numerous small-scale dust-gas vortices.

Application of cubesats in astronomical observations has been getting more and more mature in recent years. Here we report a concept study of a small Compton polarimeter to fly on a cubesat for observing polarization of soft gamma-rays from a black-hole X-ray binary, Cygnus X-1. Polarization states provide very useful diagnostics on the emission mechanism and the origin of those gamma rays. In our study, we conducted Monte Carlo simulations to decide the basic design of this small polarimeter. Silicon detectors and cerium bromide scintillators were employed in this study. We estimated its on-axis Compton efficiency at different energies and its data telemetry requirement when flying in a low earth orbit. Our results indicate that it is feasible to achieve high signal-to-noise ratio for observing Cyg X-1 with such a small instrument. Based on this study, we will proceed to have a more realistic design and look for opportunities of a cubesat space mission.

We report measurements of the charged daughter fraction of $^{218}$Po as a result of the $^{222}$Rn alpha decay, and the mobility of $^{218}$Po$^+$ ions, using radon-polonium coincidences from the $^{238}$U chain identified in 532 live-days of DarkSide-50 WIMP-search data. The fraction of $^{218}$Po that is charged is found to be 0.37$\pm$0.03 and the mobility of $^{218}$Po$^+$ is (8.6$\pm$0.1)$\times$10$^{-4}$$\frac{\text{cm}^2}{\text{Vs}}$.

We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$\alpha$ and H$\beta$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$ km s$^{-1}$, fulfilling the conventional definition of a narrow line Seyfert 1 (NLS1) galaxy. However, its NELs are unprecedentedly strong when compared to type 1AGNs. This, together with its large MIR excess ($g-W_4 = 13.172$ mag), implies that the observer and the narrow emission line region (NELR) might see a different ionization continuum. Our optical spectropolarimetry confirms its type 2 nature by detecting a polarized blue continuum and Balmer BELs ($FWHM_{\rm Polarized~BEL} \approx 1,183$ km s$^{-1}$), with a high polarization degree of $> 20$% in the blue wing. The mass and Eddington rate of the central black hole is estimated based on both transmitted and scattered AGN radiation, which is $M_{\bullet} < 2.9 \times 10^7 M_\odot$ and $L_{\rm bol}/L_{\rm Edd} > 1.5$. Severe extinction of the AGN emission also enables a clear view of the compact host galaxy. Discovery of J1203+1624 suggests that NLS1 follows AGN unfication scheme, studying of its analogs could blaze a new trail for exploring the connection between black hole growth and star formation in the host galaxy. Interesting features of J1203+1624 like the peculiar NELs, inflowing scattering clouds within sublimation radius are worth detailed follow-ups in the future.

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population.

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $\alpha$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$\alpha$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $\alpha$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

We study the dust properties of 192 nearby galaxies from the JINGLE survey using photometric data in the 22-850micron range. We derive the total dust mass, temperature T and emissivity index beta of the galaxies through the fitting of their spectral energy distribution (SED) using a single modified black-body model (SMBB). We apply a hierarchical Bayesian approach that reduces the known degeneracy between T and beta. Applying the hierarchical approach, the strength of the T-beta anti-correlation is reduced from a Pearson correlation coefficient R=-0.79 to R=-0.52. For the JINGLE galaxies we measure dust temperatures in the range 17-30 K and dust emissivity indices beta in the range 0.6-2.2. We compare the SMBB model with the broken emissivity modified black-body (BMBB) and the two modified black-bodies (TMBB) models. The results derived with the SMBB and TMBB are in good agreement, thus applying the SMBB, which comes with fewer free parameters, does not penalize the measurement of the cold dust properties in the JINGLE sample. We investigate the relation between T and beta and other global galaxy properties in the JINGLE and Herschel Reference Survey (HRS) sample. We find that beta correlates with the stellar mass surface density (R=0.62) and anti-correlates with the HI mass fraction (M(HI)/M*, R=-0.65), whereas the dust temperature correlates strongly with the SFR normalized by the dust mass (R=0.73). These relations can be used to estimate T and beta in galaxies with insufficient photometric data available to measure them directly through SED fitting.

We construct a large halo K-giant sample by combining the positions, distances, radial velocities, and metallicities of over 13,000 LAMOST DR5 halo K giants with the Gaia DR2 proper motions, which covers a Galactocentric distance range of 5-120 kpc. Using a position-velocity clustering estimator (the 6Distance), we statistically quantify the presence of position-velocity substructure at high significance: K giants have more close pairs in position-velocity space than a smooth stellar halo. We find that the amount of substructure in the halo increases with increasing distance and metallicity. With a percolation algorithm named friends-of-friends (FoF) to identify groups, we identify members belonging to Sagittarius (Sgr) Streams, Monoceros Ring, Virgo overdensity, Hercules-Aquila Cloud, Orphan Streams and other unknown substructures and find that the Sgr streams account for a large part of grouped stars beyond 20 kpc and enhance the increase of substructure with distance and metallicity. For the first time, we identify spectroscopic members of Monoceros Ring in the south and north Galactic hemisphere, which presents a rotation of about 185 km s^{-1} and mean metallicity is -0.66 dex.

Changing look active-galatic-nuclei (CL AGNs) can yield considerable insight into accretion physics as well as the co-evolution of black holes and their host galaxies. A large sample of these CL AGNs is essential to achieve the latter goal. We propose an effective method to select CL candidates from spectroscopic quasar catalogs using the mid-infrared (MIR) variability information derived from ALL-WISE/NEOWISE data releases. Our primary selection criteria include both a large amplitude flux variation and a transition of MIR color from an AGN to a normal galaxy. A pilot spectroscopic follow-up of 7 candidates among about 300 candidates selected from Sloan Digital Sky Survey low-redshift (z<0.5) AGN sample results in 6 new turn-off CL AGNs. We expect to obtain hundreds of CL AGNs once full spectroscopic follow-up of the sample is carried out.

We present the SCUBA-2 850 ${\mu}m$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 ${\mu}m$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $\beta$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample.

Follow-up observations of (sub-)mm-selected gravitationally-lensed systems have allowed a more detailed study of the dust-enshrouded phase of star-formation up to very early cosmic times. Here, the case of the gravitationally lensed merger in HATLAS J142935.3-002836 (also known as H1429-0028; z_lens=0.218, z_bkg=1.027) is revisited following recent developments in the literature and new APEX observations targeting two carbon monoxide (CO) rotational transitions J_up=3 and 6. We show that the line-profiles comprise three distinct velocity components, where the fainter high-velocity one is less magnified and more compact. The modelling of the observed spectral line energy distribution of CO J_up=2 to 6 and [CI]3P_1-3P_0 assumes a large velocity gradient scenario, where the analysis is based on four statistical approaches. Since the detected gas and dust emission comes exclusively from only one of the two merging components (the one oriented North-South, NS), we are only able to determine upper-limits for the companion. The molecular gas in the NS component in H1429-0028 is found to have a temperature of ~70K, a volume density of log(n[/cm3])~3.7, to be expanding at ~10km/s/pc, and amounts to M_H2=4(-2,+3)*1e9 Msun. The CO to H2 conversion factor is estimated to be alpha_CO=0.4(-0.2,+0.3) Msun/(K.km/s.pc2). The NS galaxy is expected to have a factor of >10x more gas than its companion (M_H2<3e8 Msun). Nevertheless, the total amount of molecular gas in the system comprises only up to 15 per cent (1sigma upper-limit) of the total (dynamical) mass.

Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$\mu$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $\mu$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged)

There are still some important unanswered questions about the detailed particle acceleration and escape occurring during the quiescent epoches. As a result, the particle distribution that is adopted in the blazar quiescent spectral model have numerous unconstrained shapes. To help remedy this problem, we introduce a analytical particle transport model to reproduce quiescent broadband spectral energy distribution of blazar. In this model, the exact electron distribution is solved from a generalized transport equation that contains the terms describing first-order and secondary-order \emph{Fermi} acceleration, escape of particle due to both the advection and spatial diffusion, energy losses due to synchrotron emission and inverse-Compton scattering of an assumed soft photon field. We suggest that the advection is a significant escape mechanism in blazar jet. We find that in our model the advection process tends to harden the particle distribution, which enhances the high energy components of resulting synchrotron and synchrotron self-Comptom spectrum from jet. Our model is able to roughly reproduce the observed spectra of extreme BL Lac object 1ES 0414+009 with reasonable assumptions about the physical parameters.

We present a detailed study of a transient in the center of SDSS1115+0544 based on the extensive UV, optical, mid-IR light curves (LC) and spectra over 1200 days. The host galaxy is a quiescent early type galaxy at $z$ = 0.0899 with a blackhole mass of $2\times10^7M_\odot$. The transient underwent a 2.5 magnitude brightening over $\sim120$ days, reaching a peak $V$-band luminosity (extinction corrected) of $-20.9$ magnitude, then fading 0.5 magnitude over 200 days, settling into a plateau of $>600$ days. Following the optical brightening are the significant mid-IR flares at $3.4$ and $4.5\mu$m, with a peak time delay of $\sim180$ days. The mid-IR LCs are explained as the echo of UV photons by a dust medium with a radius of $5\times10^{17}$ cm, consistent with $\rm E(B-V)$ of 0.58 inferred from the spectra. This event is very energetic with an extinction corrected $L_{bol} \sim 4\times10^{44}$ erg s$^{-1}$. Optical spectra over 400 days in the plateau phase revealed newly formed broad H$\alpha, \beta$ emission with a FWHM of $\sim3750$ km s$^{-1}$ and narrow coronal lines such as [Fe VII], [Ne V]. This flare also has a steeply rising UV continuum, detected by multi-epoch $Swift$ data at $+700$ to $+900$ days post optical peak. The broad Balmer lines and the UV continuum do not show significant temporal variations. The slow evolving LCs over 1200 days, the constant Balmer lines and UV continuum at late-times rule out TDE and SN IIn as the physical model for this event. We propose that this event is a `turn-on' AGN, transitioning from a quiescent state to a type 1 AGN with a sub-Eddington accretion rate of $0.017M_\odot$/yr. This change occurred on a very short time scale of $\sim 120- 200$ days. The discovery of such a rapid `turn-on' AGN poses challenges to accretion disk theories and may indicate such event is not extremely rare.

Based on the [Fe/H] versus [Mg/Fe] diagram and distances from APOGEE data release 14, we compare the spatial distributions, the l-Vlos diagram and the abundance gradients between high-[Mg/Fe] and low-[Mg/Fe] sequences. The two sequences are clearly shown at 5<|Z|<10 kpc in the metallicity range of -1.6 <[Fe/H] <-0.7, where the halo at |Z| > 10 kpc consists of low-[Mg/Fe] stars only. In the intermediate-metallicity range of -1.1 <[Fe/H]<-0.7, a [Mg/Fe] gradient is detected for stars at |Z|=10-30 kpc and it flattens out at |Z|>30 kpc. The l-Vlos diagram is adopted to separate halo stars from the disk by defining the transition metallicity, which is of [Fe/H]~ -1.1 dex for the high-[Mg/Fe] sequence and of [Fe/H]~-0.7 dex for the low-[Mg/Fe] sequence. The R and |Z| distributions for the high-[Mg/Fe] sequence, the thick disk at -1.1<[Fe/H]<-0.7 and the in situ halo at -1.6<[Fe/H]<-1.1, have a cutoff at R~15 kpc and |Z|~10 kpc, beyond which low-[Mg/Fe] halo stars are the main contributions. In the metallicity range of -1.6<[Fe/H]<-0.7, there is a negative metallicity gradient for the high-[Mg/Fe] halo at |Z|<8-10 kpc, while only a marginal or no slope in the [Fe/H] versus |Z| diagram for the low-[Mg/Fe] halo at |Z|<8-10 kpc, beyond which both the high-[Mg/Fe] halo and low-[Mg/Fe] halo flatten out toward |Z| > 20 kpc. These results indicate a complicated formation history of the Galaxy and we may see a hint of a three-section halo, i.e. the inner in situ halo within $|Z|~8-10$ kpc, the intermediately outer dual-mode halo at |Z|~10-30 kpc, and the extremely outer accreted halo with |Z|>30 kpc.

We show that the non-integer effective number of neutrinos $N^{\mathrm{eff}}_\nu$ can be understood as an effect of lepton $L$ asymmetry in the early Universe carried by the Dirac neutrino cosmic background. We show that $N_\nu^{\mathrm{eff}}=3.36\pm0.34$ (CMB only) and $N_\nu^{\mathrm{eff}}= 3.62\pm0.25$ (CMB and $H_0$) require a ratio between baryon number $B$ and lepton number to be $1.16 \times 10^{-9}\leqslant B/|L|\leqslant 1.51 \times 10^{-9}$. These values are close to the baryon-to-photon ratio $0.57\times 10^{-9}\leqslant B/N_\gamma \leqslant 0.67\times10^{-9}$. Thus instead of the usual $|L|\ll N_\gamma$ and $B\simeq |L|$, we propose to use $0.4 \leqslant |L|/N_\gamma\leqslant 0.52$ and $B\ll|L|$ as another natural choice, which resolves the tension between Planck-CMB and $H_0$ and leads to a non-integer value of $N_\nu^{\mathrm{eff}}>3$.

The bright [CI] 1-0 and [CI] 2-1 lines of atomic carbon are becoming more and more widely employed tracers of the cold neutral gas in high-redshift galaxies. Here we present observations of these lines in the 11 galaxies of the set of Planck's Dusty GEMS, the brightest gravitationally lensed galaxies on the extragalactic submillimeter sky probed by the Planck satellite. We have [CI] 1-0 measurements for seven, and [CI] 2-1 measurements for eight galaxies, including four galaxies where both lines are measured. We use our observations to constrain the gas excitation mechanism, excitation temperatures, optical depths, atomic carbon and molecular gas masses, and carbon abundances. Ratios of L_CI/L_ FIR are similar to those found in the local Universe, and suggest that the total cooling budget through atomic carbon has not strongly changed in the last 12 Gyr. Both lines are optically thin and trace 1 - 6 x 10^7 M_sun of atomic carbon. Carbon abundance ratios with H_2, X_CI, are between 2.5 and 4 x 10^-5, for a "ULIRG" CO-to-H_2 conversion factor of alpha_CO=0.8 M_sun/ [K km s^-1 pc^2]. Ratios of molecular gas masses derived from [CI] 1-0 and CO agree within the measurement uncertainties for five galaxies, and to better than a factor of 2 for another two with [CI] 1-0 measurements, after taking CO excitation carefully into account. This does not support the idea that intense, high-redshift starburst galaxies host large quantities of "CO-dark" gas. These results also support the common assumptions underlying most molecular gas mass estimates made for massive, dusty, high-redshift starburst galaxies, although the good agreement between the masses obtained with both tracers cannot be taken as an independent confirmation of either alpha_CO or X_CI.

We present an extensive CO emission-line survey of the Planck's dusty GEMS, a small set of 11 strongly lensed dusty star-forming galaxies at z = 2-4 discovered with Planck and Herschel satellites, using EMIR on the IRAM 30-m telescope. We detected a total of 45 CO rotational lines from Jup=3 to Jup=11, and up to eight transitions per source, allowing a detailed analysis of the gas excitation and interstellar medium conditions within these extremely bright, vigorous starbursts. We applied radiative transfer models using the large velocity gradient approach to infer the spatially-averaged molecular gas densities, $n_{H_2}$~10$^{2.6-4.1}$ cm$^{-3}$, and kinetic temperatures, $T_k$~30-1000 K. In five sources, we find evidence of two distinct gas phases with different properties and model their CO ladder with two excitation components. The warm (70-320 K) and dense gas reservoirs in these galaxies are highly excited, while the cooler (15-60 K) and more extended low-excitation components cover a range of gas densities. In two sources, the latter is associated with diffuse Milky Way-like gas phases, which provides evidence that a significant fraction of the total gas masses of dusty starburst galaxies can be embedded in cool, low-density reservoirs. Finally, we show that the CO line luminosity ratios are consistent with those predicted by models of photon-dominated regions and disfavor scenarios of gas clouds irradiated by intense X-ray fields from active galactic nuclei. By combining CO, [CI] and [CII] line diagnostics, we obtain average PDR gas densities significantly higher than in normal star-forming galaxies at low-redshift, and far-ultraviolet radiation fields at least 100 times more intense than in the Milky Way. These spatially-averaged conditions are consistent with those in high-redshift SMGs and in a range of low-redshift environments.

The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

We report the discovery of an active intermediate-mass black hole (IMBH) candidate in the center of nearby barred bulgeless galaxy $\rm NGC~3319$. The point X-ray source revealed by archival Chandra and XMM-Newton observations is spatially coincident with the optical and UV galactic nuclei from Hubble Space Telescope observations. The spectral energy distribution derived from the unresolved X-ray and UV-optical flux is comparable with active galactic nuclei (AGNs) rather than ultra-luminous X-ray sources, although its bolometric luminosity is only $3.6\times10^{40}~\rm erg~s^{-1}$. Assuming an Eddington ratio range between 0.001 and 1, the black hole mass ($M_\rm{BH}$) will be located at $3\times10^2-3\times10^5~M_{\odot}$, placing it in the so-called IMBH regime and could be the one of the lowest reported so far. Estimates from other approaches (e.g., fundamental plane, X-ray variability) also suggest $M_\rm{BH}\lesssim10^5~M_{\odot}$.Similar to other BHs in bulgeless galaxies, the discovered IMBH resides in a nuclear star cluster with mass of $\sim6\times10^6~M_{\odot}$. The detection of such a low-mass BH offers us an ideal chance to study the formation and early growth of SMBH seeds, which may result from the bar-driven inflow in late-type galaxies with a prominent bar such as $\rm NGC~3319$.

The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MHD) effects meets this requirement. However, simulations have shown that the midplane of a dead zone is not completely quiescent. In order to examine the concentration of solids in such an environment, we use the local-shearing-box approximation to simulate a particle-gas system with an Ohmic dead zone including mutual drag force between the gas and the solids. We systematically compare the evolution of the system with ideal or non-ideal MHD, with or without back-reaction drag force from particles on gas, and with varying solid abundances. Similar to previous investigations of dead zone dynamics, we find that particles of dimensionless stopping time $\tau_s=0.1$ do not sediment appreciably more than those in ideal magneto-rotational turbulence, resulting in a vertical scale height an order of magnitude larger than in a laminar disk. Contrary to the expectation that this should curb the formation of planetesimals, we nevertheless find that strong clumping of solids still occurs in the dead zone when solid abundances are similar to the critical value for a laminar environment. This can be explained by the weak radial diffusion of particles near the mid-plane. The results imply that the sedimentation of particles to the mid-plane is not a necessary criterion for the formation of planetesimals by the streaming instability.

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valpara\'iso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/\sigma_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

We study the radial metallicity gradient $\Delta[M/H]/\Delta R_g$ as a function of [Mg/Fe] and $|Z|$ with the help of a guiding radius based on the Apache Point Observatory Galactic Evolution Experiment and Gaia and then analyze the radial migration effect on the radial metallicity gradient and metallicity-rotation gradient between the Galactic thin and thick disks. The derived trend of gradient $\Delta[M/H]/\Delta R_g$ versus [Mg/Fe] shows a transition at [Mg/Fe] $\sim 0.18$ dex, below which the gradient is negative and varies a little as [Mg/Fe] increases; however, it changes sharply in [Mg/Fe] ranges of 0.16-0.18, above which the gradient increases linearly with increasing [Mg/Fe], being a positive value at [Mg/Fe]$\gtrsim 0.22$ dex. These positive gradients in the high-[Mg/Fe] populations are found at $|Z| < 0.8$ kpc, and there are nearly no gradients toward higher $|Z|$. By comparing the metallicity distributions, the radial metallicity gradients $\Delta[M/H]/\Delta R$ and the metallicity-rotation gradients between the total sample and $|R-R_g|<2$ kpc subsample (or $|R-R_g|>2$ kpc subsample), we find that, for the thick disk, blurring flattens the gradient $\Delta[M/H]/\Delta R$ and favors metal-poor high-eccentricity stars. These stars are responsible for the measured positive metallicity-rotation gradient of the thick disk.

The radial drift and diffusion of dust particles in protoplanetary disks affect both the opacity and temperature of such disks as well as the location and timing of planetesimal formation. In this paper, we present results of numerical simulations of particle-gas dynamics in protoplanetary disks that include dust grains with various size distributions. We consider three scenarios in terms of particle size ranges, one where the Stokes number $\tau_{\rm{s}} = 10^{-1} - 10^0$, one where $\tau_{\rm{s}} = 10^{-4} - 10^{-1}$ and finally one where $\tau_{\rm{s}} = 10^{-3} - 10^{0}$. Moreover, we consider both discrete and continuous distributions in particle size. In accordance with previous works we find in our multi-species simulations that different particle sizes interact via the gas and as a result their dynamics changes compared to the single-species case. The larger species trigger the streaming instability and create turbulence that drives the diffusion of the solid materials. We measure the radial equilibrium velocity of the system and find that the radial drift velocity of the large particles is reduced in the multi-species simulations and that the small particle species move on average outwards. We also vary the steepness of the size distribution, such that the exponent of the solid number density distribution, $\rm{d}\it{N}/\rm{d}\it{a} \propto \it{a^{-q}}$, is either $q = 3$ or $q = 4$. We measure the scale height of the particles and observe that small grains are stirred up well above the sedimented midplane layer where the large particles reside. Our measured diffusion and drift parameters can be used in coagulation models for planet formation as well as to understand relative mixing of the components of primitive meteorites (matrix, chondrules and CAIs) prior to inclusion in their parent bodies.

The most intensively star-forming galaxies are extremely luminous at far-infrared (FIR) wavelengths, highly obscured at optical and ultraviolet wavelengths, and lie at $z\ge 1-3$. We present a programme of ${\it Herschel}$ FIR spectroscopic observations with the SPIRE FTS and photometric observations with PACS, both on board ${\it Herschel}$, towards a sample of 45 gravitationally lensed, dusty starbursts across $z\sim 1-3.6$. In total, we detected 27 individual lines down to 3-$\sigma$, including nine $[\rm C{\small II}]$ 158-$\mu$m lines with confirmed spectroscopic redshifts, five possible $[\rm C{\small II}]$ lines consistent with their far-infrared photometric redshifts, and in some individual sources a few $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 145-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, and OH 119-$\mu$m (in absorption) lines. To derive the typical physical properties of the gas in the sample, we stack all spectra weighted by their intrinsic luminosity and by their 500-$\mu$m flux densities, with the spectra scaled to a common redshift. In the stacked spectra, we detect emission lines of $[\rm C{\small II}]$ 158-$\mu$m, $[\rm N{\small II}]$ 122-$\mu$m, $[\rm O{\small III}]$ 88-$\mu$m, $[\rm O{\small III}]$ 52-$\mu$m, $[\rm O{\small I}]$ 63-$\mu$m, and the absorption doublet of OH at 119-$\mu$m, at high fidelity. We find that the average electron densities traced by the $[\rm N{\small II}]$ and $[\rm O{\small III}]$ lines are higher than the average values in local star-forming galaxies and ULIRGs, using the same tracers. From the $[\rm N{\small II}]/[\rm C{\small II}]$ and $[\rm O{\small I}]/[\rm C{\small II}]$ ratios, we find that the $[\rm C{\small II}]$ emission is likely dominated by the photo-dominated regions (PDR), instead of by ionised gas or large-scale shocks.

We describe the theory of Fermi-type acceleration, including first-order Fermi acceleration at a parallel shock front and second-order Fermi acceleration in a test particle limit. Including the theory of the turbulent acceleration and the derivation of the general Fokker-Planck equation, we take into account the basic particle-transport equation and the Fokker-Planck equation with spatially homogeneous isotropic distribution. In the cases of some special physical processes, we construct the particle-transport equations, and compile the analytical or semi-analytical solutions in 11 different cases. Even though, traditionally, one of the electron energy distributions for these different cases can be used to reproduce the multi-wavelength emission of a blazar, due to very long integration times of high-energy observations, a time-integrated electron energy distribution should be introduced in approaches to modeling the multi-wavelength spectral energy distribution of blazars, especially in modeling the flaring state.

The present work considers a plane shock front propagating along a cylindrical jet. Electrons experience the diffusive shock acceleration around the shock front, and subsequently drift away into the downstream flow in which they emit most of their energy. Assuming a proper boundary condition at the interface between the shock zone and the downstream zone, we solve the transport equation for the electrons in the downstream flow zone, where the combined effects of escape, synchrotron and IC cooling in the Thomson regime are taken into account. Using the electron spectrum obtained in this manner we calculate the multi-wavelength spectral energy distribution of Mrk 501 in the synchrotron self-Compton scenario. We check numerically if the Klein-Nishina cross-section could be approximated to the Thomsom regime. We consider whether the model results yield physically reasonable parameters, and further discuss some of implications of the model results. It suggests that the process of diffusive shock acceleration operates in the outflow of Mrk 501.

Current precision big bang nucleosynthesis (BBN) studies motivate us to revisit the neutron lifespan in the plasma medium of the early universe. The mechanism we explore is the Fermi-blocking of decay electrons and neutrinos by plasma. As result, neutrons live longer and we find a significant 6.4\% modification of neutron abundance in the BBN era arising from in plasma modification of the neutron lifetime. This effect can influence the final abundances of the light elements in BBN.

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.

The anisotropy parameter $\beta$ characterizes the extent to which orbits in stellar systems are predominantly radial or tangential, and is likely to constrain, for the stellar halo of the Milky Way, scenarios for its formation and evolution. We have measured the anisotropy $\beta$ as a function of Galactocentric radius from $5-100$ kpc for over 8600 metal poor ([Fe/H] $<-1.3$) halo K giants from the LAMOST catalog with line-of-sight velocities and distances, matched to proper motions from the second $Gaia$ data release. We construct full 6-D positions and velocities for the K giants to directly measure the 3 components of the velocity dispersion $(\sigma_r, \sigma_\theta, \sigma_\phi)$ (in spherical coordinates). We find that the orbits in the halo are radial over our full Galactocentric distance range reaching over 100 kpc. The anisotropy remains remarkably unchanged with Galactocentric radius from approximately 5 to 25 kpc, with an amplitude that depends on the metallicity of the stars, dropping from $\beta \approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $\beta \approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $\beta\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $\beta<0.3$ past 100 kpc. Within 8 kpc, $\beta<0.8$. The measurement of anisotropy is affected by substructure and streams, particularly beyond a Galactocentric distance of approximately 25 kpc, where the Sagittarius stream is prominent in the data. These results are complimentary to recent analysis of simulations by Loebman et al. and of SDSS/$Gaia$ DR1 data by Belokurov et al.

Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atomic Carbon in high-z galaxies, at previously inaccessible frequencies from the ground. We have targeted three gravitationally lensed galaxies at redshift of about 3 and conducted a comparative study of the observed high-J CO/CI ~ratios with well-studied nearby galaxies. Atomic Carbon (CI(2-1)) was detected in one of the three targets and marginally in a second, while in all three targets the $J=7\to6$ CO line is detected. The CO(7-6)/CI(2-1), CO(7-6)/CO(1-0) line ratios and the CO(7-6)/(far-IR continuum) luminosity ratio are compared to those of nearby objects. A large excitation status in the ISM of these high-z objects is seen, unless differential lensing unevenly boosts the CO line fluxes from the warm and dense gas more than the CO(1-0), CI(2-1), tracing a more widely distributed cold gas phase. We provide estimates of total molecular gas masses derived from the atomic Carbon and the Carbon monoxide CO(1-0), which within the uncertainties turn out to be equal.

We report the discovery of a sample of 19 low redshift (z<0.22) spectroscopically non-Seyfert galaxies that show slow declining mid-infrared (MIR) light-curves (LCs), similar to those of tidal disruption event (TDE) candidates with extreme coronal lines. Two sources also showed a relatively fast rising MIR LCs. They consist of 61% sample of the WISE MIR variable non-Seyfert galaxies with SDSS spectra. In a comparison sample of optically selected Seyfert galaxies, the fraction of sources with such a LC is only 15%. After rejecting 5 plausible obscured Seyfert galaxies with red MIR colours, remaining 14 objects are studied in detail in this paper. We fit the declining part of LC with an exponential law, and the decay time is typically one year. The observed peak MIR luminosities ($\nu L_\nu$) after subtracting host galaxies are in the range of a few 10^42 to 10^44 erg~s^-1 with a median of 5x10^43 erg~s^-1 in the W2 band. The black hole masses distribute in a wide range with more than half in between 10^7 to 10^8 ~M_sun, but significantly different from that of optical/UV selected TDEs. Furthermore, MIR luminosities are correlated with black hole masses, the stellar mass or luminosity of their host bulges. Most galaxies in the sample are red and luminous with an absolute magnitude of r between -20 to -23. We estimate the rate of event about 10^-4 gal^-1~yr^-1 among luminous red galaxies. We discuss several possibilities for the variable infrared sources, and conclude that most likely, they are caused by short sporadic fueling to the supermassive black holes via either the instability of accretion flows or tidal disruption of stars.

In this work we present new APEX/SEPIA Band-5 observations targeting the CO ($J=2\text{-}1$) emission line of 24 Herschel-detected galaxies at $z=0.1-0.2$. Combining this sample {with} our recent new Valpara\'iso ALMA Line Emission Survey (VALES), we investigate the star formation efficiencies (SFEs = SFR/$M_{\rm H_{2}}$) of galaxies at low redshift. We find the SFE of our sample bridges the gap between normal star-forming galaxies and Ultra-Luminous Infrared Galaxies (ULIRGs), which are thought to be triggered by different star formation modes. Considering the $\rm SFE'$ as the SFR and the $L'_{\rm CO}$ ratio, our data show a continuous and smooth increment as a function of infrared luminosity (or star formation rate) with a scatter about 0.5 dex, instead of a steep jump with a bimodal behaviour. This result is due to the use of a sample with a much larger range of sSFR/sSFR$_{\rm ms}$ using LIRGs, with luminosities covering the range between normal and ULIRGs. We conclude that the main parameters controlling the scatter of the SFE in star-forming galaxies are the systematic uncertainty of the $\alpha_{\rm CO}$ conversion factor, the gas fraction and physical size.

We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.

We present X-ray observations of PWN G16.73+0.08/SNR G16.7+0.1 using archival data of {\it Chandra} ACIS. The X-ray emission peak location of this pulsar wind nebula is found to be offset by 24 arcsec from the centre of the 1.4-GHz emission of this nebula. The X-ray nebula is elongated in the direction from the X-ray peak to the 1.4-GHz emission centre. This offset suggests that G16.73+0.08 is an evolved pulsar wind nebula interacting with the supernova remnant reverse shock. We identify a point source, CXO J182058.16-142001.5, near the location of the X-ray peak. The spectrum of the X-ray nebula can be described by an absorbed power law of photon index $0.98^{+0.79}_{-0.71}$ and hydrogen column density $N_{\rm H}=4.99^{+2.75}_{-2.28}\times 10^{22}$ cm$^{-2}$. CXO J182058.16-142001.5 is likely a pulsar. We estimate its spin-down power to be about $2.6\times 10 ^{36}$ erg s$^{-1}$. Assuming its age at 3000 and 10,000 years, its dipole magnetic field strength at the polar surface is estimated to be about $4.2 \times 10^{13}$ G and $1.3 \times 10^{13}$ G, respectively.

PS16dtm was classified as a candidate tidal disruption event (TDE) in a dwarf Seyfert 1 galaxy with low-mass black hole ($\sim10^6M\odot$) and has presented various intriguing photometric and spectra characteristics. Using the archival WISE and the newly released NEOWISE data, we found PS16dtm is experiencing a mid-infrared (MIR) flare which started $\sim11$ days before the first optical detection. Interpreting the MIR flare as a dust echo requires close pre-existing dust with a high covering factor, and suggests the optical flare may have brightened slowly for some time before it became bright detectable from the ground. More evidence is given at the later epochs. At the peak of the optical light curve, the new inner radius of the dust torus has grown to much larger size, a factor of 7 of the initial radius due to strong radiation field. At $\sim150$ days after the first optical detection, the dust temperature has dropped well below the sublimation temperature. Other peculiar spectral features shown by PS16dtm are the transient, prominent FeII emission lines and outflows indicated by broad absorption lines detected during the optical flare. Our model explains the enhanced FeII emission from iron newly released from the evaporated dust. The observed broad absorption line outflow could be explained by accelerated gas in the dust torus due to the radiation pressure.

(abridged) We present the IRAM-30m observations of multiple-J CO and CI line emission in a sample of redshift ~2-4 Herschel-ATLAS SMGs. A non-negligible effect of differential lensing is found for the CO emission lines, which could have caused significant underestimations of the linewidths, hence of the dynamical masses. The CO SLEDs are found to be similar to those of the local starburst-dominated ULIRGs and of the previously studied SMGs. After correcting for lensing amplification, we derived the global properties of the bulk of molecular gas in the SMGs using non-LTE radiative transfer modelling. The gas thermal pressure is found to be correlated with star formation efficiency. Further decomposing the CO SLEDs into two excitation components, we find a low-excitation component, which is less correlated with star formation, and a high-excitation one which is tightly related to the on-going star-forming activity. Additionally, tight linear correlations between the FIR and CO line luminosities have been confirmed for the $J \ge 5$ CO lines, implying that these CO lines are good tracers of star formation. The [CI](2-1) lines follow the tight linear correlation between the luminosities of the [CI](2-1) and the CO(1-0) line found in local starbursts, indicating that CI lines could serve as good total molecular gas mass tracers for high-redshift SMGs. The total mass of the molecular gas reservoir, $(1-30) \times 10^{10} M_\odot$, suggests a typical molecular gas depletion time ~20-100 Myr and a gas to dust mass ratio ${\delta}_{\rm GDR}$~30-100. The ratio between CO line luminosity and the dust mass appears to be slowly increasing with redshift for the SMGs, which need to be further confirmed. Finally, through comparing the linewidth of CO and H2O lines, we find that they agree well in almost all our SMGs, confirming that the emitting regions are co-spatially located.

Astrophysical polarization measurements in the soft gamma-ray band are becoming more feasible as detectors with high position and energy resolution are deployed. Previous work has shown that the minimum detectable polarization (MDP) of an ideal Compton polarimeter can be improved by $\sim 21\%$ when an unbinned, maximum likelihood method is used instead of the standard approach of fitting a sinusoid to a histogram of azimuthal scattering angles. Here we outline a procedure for implementing this maximum likelihood approach for real, non-ideal polarimeters. As an example, we use the recent observation of GRB 160530A with the Compton Spectrometer and Imager. We find that the MDP for this observation is reduced by $20\%$ when the maximum likelihood method is used instead of the standard method.

A long duration gamma-ray burst, GRB 160530A, was detected by the Compton Spectrometer and Imager (COSI) during the 2016 COSI Super Pressure Balloon campaign. As a Compton telescope, COSI is inherently sensitive to the polarization of gamma-ray sources in the energy range 0.2-5.0 MeV. We measured the polarization of GRB 160530A using 1) a standard method (SM) based on fitting the distribution of azimuthal scattering angles with a modulation curve, and 2) an unbinned, maximum likelihood method (MLM). In both cases, the measured polarization level was below the $99\%$ confidence minimum detectable polarization levels of $72.3 \pm 0.8\%$ (SM) and $57.5 \pm 0.8\%$ (MLM). Therefore, COSI did not detect polarized gamma-ray emission from this burst. Our most constraining $90\%$ confidence upper limit on the polarization level was $46\%$ (MLM).

We present a model-dependent method to estimate the redshifts of three TeV BL Lac objects (BL Lacs) through fltting their (quasi-) simultaneous multi-waveband spectral energy distributions (SEDs) by one-zone leptonic synchrotron self-Compton (SSC) model. Considering the impact of electron energy distributions (EEDs) on the results, we use three types of EEDs, such as the power-law EED with exponential cut-ofi (PLC), the log-parabola (PLLP) and the broken power-law (BPL) EEDs, to flt the SEDs. We also use a parameter fi to describe the uncertainties of the extragalactic background light (EBL) models, as in (Abdo et al. 2010b). We then use Markov Chain Monte Carlo (MCMC) method to explore multi-dimensional parameter space and obtain the uncertainties of the model parameters based on the observational data. We apply our method to obtain the redshifts of three TeV BL Lac objects in the marginalized 68% confldence, and flnd that the PLC EED does not flt the SEDs. For 3C66A, the redshift is 0.14 - 0.31 and 0.16 - 0.32 in the BPL and PLLP EEDs; for PKS1424+240, the redshift is 0.55 - 0.68 and 0.55 - 0.67 in the BPL and PLLP EEDs; for PG1553+113, the redshift is 0.22 - 0.48 and 0.22 - 0.39 in the BPL and PLLP EEDs. We also estimate the redshift of PKS1424+240 in the high stage to be 0.46 - 0.67 in the PLLP EED, roughly consistent with that in the low stage.

We present the Herschel Bright Sources (HerBS) sample, a sample of bright, high-redshift Herschel sources detected in the 616.4 square degree H-ATLAS survey. The HerBS sample contains 209 galaxies, selected with a 500 {\mu}m flux density greater than 80 mJy and an estimated redshift greater than 2. The sample consists of a combination of HyLIRGs and lensed ULIRGs during the epoch of peak cosmic star formation. In this paper, we present SCUBA-2 observations at 850 ${\mu}$m of 189 galaxies of the HerBS sample, 152 of these sources were detected. We fit a spectral template to the Herschel-SPIRE and 850 ${\mu}$m SCUBA-2 flux densities of 22 sources with spectroscopically determined redshifts, using a two-component modified blackbody spectrum as a template. We find a cold- and hot-dust temperature of 21.29 K and 45.80 K, a cold-to-hot dust mass ratio of 26.62 and a $\beta$ of 1.83. The poor quality of the fit suggests that the sample of galaxies is too diverse to be explained by our simple model. Comparison of our sample to a galaxy evolution model indicates that the fraction of lenses is high. Out of the 152 SCUBA-2 detected galaxies, the model predicts 128.4 $\pm$ 2.1 of those galaxies to be lensed (84.5%). The SPIRE 500 ${\mu}$m flux suggests that out of all 209 HerBS sources, we expect 158.1 $\pm$ 1.7 lensed sources, giving a total lensing fraction of 76 per cent.

It is known that some active galactic nuclei (AGNs) transited from type 1 to type 2 or vice versa. There are two explanations for the so-called changing look AGNs: one is the dramatic change of the obscuration along the line-of-sight, the other is the variation of accretion rate. In this paper, we report the detection of large amplitude variations in the mid-infrared luminosity during the transitions in 10 changing look AGNs using WISE and newly released NEOWISE-R data. The mid-infrared light curves of 10 objects echoes the variability in the optical band with a time lag expected for dust reprocessing. The large variability amplitude is inconsistent with the scenario of varying obscuration, rather supports the scheme of dramatic change in the accretion rate.

The rest-frame UV spectra of three recent tidal disruption events (TDEs), ASASSN-14li, PTF15af and iPTF16fnl display strong nitrogen emission lines but weak or undetectable carbon lines. In these three objects, the upper limits of the C III] $\lambda 1908$/N III] $\lambda 1750$ ratio are about two orders of magnitude lower than those of quasars, suggesting a high abundance ratio of [N/C]. With detailed photoionization simulations, we demonstrate that $C^{2+}$ and $N^{2+}$ are formed in the same zone, so the C III]/N III] ratio depends only moderately on the physical conditions in the gas and weakly on the shape of the ionizing continuum. There are smaller than $0.5$ dex variations in the line ratio over wide ranges of gas densities and ionization parameters at a given metallicity. This allows a robust estimate of the relative abundance ratio nitrogen to carbon. We derive a relative abundance ratio of [N/C]$>1.5$ for ASASSN-14li, and even higher for PTF15af and iPTF16fnl. This suggests that the broad line region in those TDE sources is made of nitrogen-enhanced core material that falls back at later times. Based on stellar evolution models, the lower limit of the disrupted star should be larger than 0.6M$_{\odot}$. The chemical abundance of the line emitting gas gives a convincing evidence that the flares origin from stellar tidal disruptions. The coincidence of the weakness of the X-ray emission with the strong broad absorption lines in PTF15af, iPTF16fnl and the strong X-ray emission without such lines in ASASSN-li14 are in analogy to quasars with and without broad absorption lines.

It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon excess appeared at a time coincident with the solar flare of 8 of November 2000. This muon excess corresponds to above 40 GeV primary protons which came from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation is estimated to be about 0.1%. It has been convinced that solar protons could be accelerated to tens of GeV in those Class X solar flares which usually arose solar cosmic ray ground level enhancement (GLE) events. However, whether a Class M solar flare like the non-GLE event of 8 November 2000 may also accelerate solar protons to such high energies? It is interesting and noteworthy.

Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time $\tau_\mathrm{s} = 10^{-2}$ and $10^{-3}$ -- representing mm- and cm-sized particles interior of the water ice line -- concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of $\tau_\mathrm{s} \ll 1$. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability.

The origin of narrow line region (NLR) outflows remains unknown. In this paper, we explore the scenario in which these outflows are circumnuclear clouds driven by energetic accretion disk winds. We choose the well-studied nearby Seyfert galaxy NGC 4151 as an example. By performing 3D hydrodynamical simulations, we are able to reproduce the radial distributions of velocity, mass outflow rate and kinetic luminosity of NLR outflows in the inner 100 pc deduced from spatial resolved spectroscopic observations. The demanded kinetic luminosity of disk winds is about two orders of magnitude higher than that inferred from the NLR outflows, but is close to the ultrafast outflows (UFO) detected in X-ray spectrum and a few times lower than the bolometric luminosity of the Seyfert. Our simulations imply that the scenario is viable for NGC 4151. The existence of the underlying disk winds can be confirmed by their impacts on higher density ISM, e.g., shock excitation signs, and the pressure in NLR.

The Einstein's equivalence principle can be tested by using parameterized post-Newtonian parameters, of which the parameter $\gamma$ has been constrained by comparing the arrival times of photons with different energies. It has been constrained by a variety of astronomical transient events, such as gamma-ray bursts (GRBs), fast radio bursts as well as pulses of pulsars, with the most stringent constraint of $\Delta \gamma \lesssim 10^{-15}$. In this letter, we consider the arrival times of lights with different circular polarization. For a linearly polarized light, it is the combination of two circularly polarized lights. If the arrival time difference between the two circularly polarized lights is too large, their combination may lose the linear polarization. We constrain the value of $\Delta \gamma_{\rm p} < 1.6 \times 10^{-27}$ by the measurement of the polarization of GRB 110721A, which is the most stringent constraint ever achieved.

We present the mid-infrared (MIR) light curves (LCs) of a tidal disruption event (TDE) candidate in the center of a nearby ultraluminous infrared galaxy (ULIRG) F01004-2237 using archival {\it WISE} and {\it NEOWISE} data from 2010 to 2016. At the peak of the optical flare, F01004-2237 was IR quiescent. About three years later, its MIR fluxes have shown a steady increase, rising by 1.34 and 1.04 mag in $3.4$ and $4.6\mu$m up to the end of 2016. The host-subtracted MIR peak luminosity is $2-3\times10^{44}$\,erg\,s$^{-1}$. We interpret the MIR LCs as an infrared echo, i.e. dust reprocessed emission of the optical flare. Fitting the MIR LCs using our dust model, we infer a dust torus of the size of a few parsecs at some inclined angle. The derived dust temperatures range from $590-850$\,K, and the warm dust mass is $\sim7\,M_{\odot}$. Such a large mass implies that the dust cannot be newly formed. We also derive the UV luminosity of $4-11\times10^{44}$\,erg\,s$^{-1}$. The inferred total IR energy is $1-2\times10^{52}$\,erg, suggesting a large dust covering factor. Finally, our dust model suggests that the long tail of the optical flare could be due to dust scattering.

The galactocentric distance of quasar absorption outflows are conventionally determined using absorption troughs from excited states, a method hindered by severely saturated or self-blended absorption troughs. We propose a novel method to estimate the size of a broad absorption line (BAL) region which partly obscures an emission line region by assuming virialized gas in the emission region surrounding a supermassive black hole with known mass. When a spiky Lya1216 line emission is present at the flat bottom of the deep Nv1240 absorption trough, the size of BAL region can be estimated. We have found 3 BAL quasars in the SDSS database showing such Lya lines. The scale of their BAL outflows are found to be 3-26 pc, moderately larger than the theoretical scale (0.01-0.1 pc) of trough forming regions for winds originating from accretion discs, but significantly smaller than most outflow sizes derived using the absorption troughs of the excited states of ions. For these three outflows, the lower limits of ratio of kinetic luminosity to Eddington luminosity are 0.02%-0.07%. These lower limits are substantially smaller than that is required to have significant feedback effect on their host galaxies.

In this paper we present an analysis of absorption line variability in mini-BAL quasar LBQS 1206+1052. The SDSS spectrum demonstrates that the absorption troughs can be divided into two components of blueshift velocities of $\sim$700 km s$^{-1}$ and $\sim$1400 km s$^{-1}$ relative to the quasar rest-frame. The former component shows rare Balmer absorption, which is an indicator of high density absorbing gas, thus the quasar is worth follow-up spectroscopic observations. Our follow-up optical and near-infrared spectra using MMT, YFOSC, TripleSpec and DBSP reveal that the strengths of the absorption lines vary for both of the two components, while the velocities do not change. We reproduce all of the spectral data by assuming that only the ionization state of the absorbing gas is variable and that all other physical properties are invariable. The variation of ionization is consistent with the variation of optical continuum from the V-band light-curve. Additionally, we can not interpret the data by assuming that the variability is due to a movement of the absorbing gas. Therefore, our analysis strongly indicates that the absorption line variability in LBQS 1206+1052 is photoionization-driven. As shown from photo-ionization simulations, the absorbing gas with blueshift velocity of $\sim$700 km s$^{-1}$ has a density in the range of $10^9$ to $10^{10}$ cm$^{-3}$ and a distance of $\sim$1 pc, and the gas with blueshift velocity of $\sim$1400 km s$^{-1}$ has a density of $10^3$ cm$^{-3}$ and a distance of $\sim$1 kpc.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

We present a statistical analysis of the variability of broad absorption lines (BALs) in quasars using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 12 (SDSS DR12). We divide the sample into two groups according to the pattern of the variation of C iv BAL with respect to that of continuum: the equivalent widths (EW) of the BAL decreases (increases) when the continuum brightens (dims) as group T1; and the variation of EW and continuum in the opposite relation as group T2. We find that T2 has significantly (P_T<10-6 , Students T Test) higher EW ratios (R) of Si iv to C iv BAL than T1. Our result agrees with the prediction of photoionization models that C +3 column density increases (decreases) if there is a (or no) C +3 ionization front while R decreases with the incident continuum. We show that BAL variabilities in at least 80% quasars are driven by the variation of ionizing continuum while other models that predict uncorrelated BAL and continuum variability contribute less than 20%. Considering large uncertainty in the continuum flux calibration, the latter fraction may be much smaller. When the sample is binned into different time interval between the two observations, we find significant difference in the distribution of R between T1 and T2 in all time-bins down to a deltaT < 6 days, suggesting that BAL outflow in a fraction of quasars has a recombination time scale of only a few days.

We present ALMA and VLA detections of the dense molecular gas tracers HCN, HCO$^+$ and HNC in two lensed, high-redshift starbursts selected from the {\it Herschel}-ATLAS survey: {\it H}-ATLAS\,J090740.0$-$004200 (SDP.9, $z \sim 1.6$) and {\it H}-ATLAS\,J091043.1$-$000321 (SDP.11, $z \sim 1.8$). ALMA observed the $J = 3-2$ transitions in both sources, while the VLA observed the $J = 1-0$ transitions in SDP.9. We have detected all observed HCN and HCO$^+$ lines in SDP.9 and SDP.11, and also HNC(3--2) in SDP.9. The amplification factors for both galaxies have been determined from sub-arcsec resolution CO and dust emission observations carried out with NOEMA and the SMA. The HNC(1--0)/HCN(1--0) line ratio in SDP.9 suggests the presence of photon-dominated regions, as it happens to most local (U)LIRGs. The CO, HCN and HCO$^+$ SLEDs of SDP.9 are compatible to those found for many local, infrared (IR) bright galaxies, indicating that the molecular gas in local and high-redshift dusty starbursts can have similar excitation conditions. We obtain that the correlation between total IR ($L_{\rm IR}$) and dense line ($L_{\rm dense}$) luminosity in SDP.9 and SDP.11 and local star-forming galaxies can be represented by a single relation. The scatter of the $L_{\rm IR} - L_{\rm dense}$ correlation, together with the lack of sensitive dense molecular gas tracer observations for a homogeneous sample of high-redshift galaxies, prevents us from distinguishing differential trends with redshift. Our results suggest that the intense star formation found in some high-redshift dusty, luminous starbursts is associated with more massive dense molecular gas reservoirs and higher dense molecular gas fractions.

The Compton Spectrometer and Imager (COSI) is a balloon-borne gamma-ray (0.2-5 MeV) telescope with inherent sensitivity to polarization. COSI's main goal is to study astrophysical sources such as $\gamma$-ray bursts, positron annihilation, Galactic nucleosynthesis, and compact objects. COSI employs a compact Compton telescope design utilizing 12 high-purity cross strip germanium detectors (size: $8\times8\times1.5$ cm$^3$, 2 mm strip pitch). We require well-benchmarked simulations to simulate the full instrument response used for data analysis, to optimize our analysis algorithms, and to better understand our instrument and the in-flight performance. In order to achieve a reasonable agreement, we have built a comprehensive mass model of the instrument and developed a detailed detector effects engine, which takes into account the individual performance of each strip as well as the characteristics of the overall detector system. We performed detailed Monte-Carlo simulations with Cosima/Geant4, applied the detector effects engine, and benchmarked the results with pre-flight calibrations using radioactive sources. After applying the detector effects engine, the simulations closely resemble the measurements, and the standard calibration, event reconstruction, and imaging pipeline used for measurements can also be applied to the simulations. In this manuscript, we will describe the detector effects engine, the benchmarking tests with calibrations, and the application to preliminary results from COSI's 46-day balloon flight in 2016.

The Compton Spectrometer and Imager (COSI) is a balloon-borne, soft-gamma ray imager, spectrometer, and polarimeter with sensitivity from 0.2 to 5 MeV. Utilizing a compact Compton telescope design with twelve cross-strip, high-purity germanium detectors, COSI has three main science goals: study the 511 keV positron annihilation line from the Galactic plane, image diffuse emission from stellar nuclear lines, and perform polarization studies of gamma-ray bursts and other extreme astrophysical environments. COSI has just completed a successful 46-day flight on NASA's new Super Pressure Balloon, launched from Wanaka, New Zealand, in May 2016. We present an overview of the instrument and the 2016 flight, and discuss COSI's main science goals, predicted performance, and preliminary results.

A model-dependent method is proposed to determine the location of the $\gamma$-ray emitting region for a given flat spectrum radio quasar (FSRQ). In the model, the extra-relativistic electrons are injected at the base of the jet and non-thermal photons are produced by both synchrotron radiation and inverse-Comtpon (IC) scattering in the energy dissipation region. The target photons dominating inverse-Comtpon scattering originate from both synchrotron photons and external ambient photon fields, and the energy density of external radiation field is a function of the distance between the position of dissipation region and a central super-massive black hole, and their spectra are seen in the comoving frame. Moreover, the energy dissipation region could be determined by the model parameter through reproducing the $\gamma$-ray spectra. Such a model is applied to reproduce the quasi-simultaneous multi-wavelength observed data for 36 FSRQs. In order to define the width of the broad-line region shell and dusty molecular torus shell, a simple numerical constraint is used to determine the outer boundary of the broad-line region and dusty molecular torus. Our results show that 1) the $\gamma$-ray emitting regions are located at the range from 0.1 pc to 10 pc; 2) the $\gamma$-ray emitting regions are located outside the broad-line regions and within the dusty molecular tori; and 3) the $\gamma$-ray emitting region are located closer to the dusty molecular torus ranges than the broad-line regions. Therefore, it may concluded that a direct evidence for the \emph{far site} scenario could be obtained on the basis of the model results.

The sporadic accretion following the tidal disruption of a star by a super-massive black hole (TDE) leads to a bright UV and soft X-ray flare in the galactic nucleus. The gas and dust surrounding the black hole responses to such a flare with an echo in emission lines and infrared emission. In this paper, we report the detection of long fading mid-IR emission lasting up to 14 years after the flare in four TDE candidates with transient coronal lines using the WISE public data release. We estimate that the reprocessed mid-IR luminosities are in the range between $4\times 10^{42}$ and $2\times 10^{43}$ erg~s$^{-1}$ and dust temperature in the range of 570-800K when WISE first detected these sources three to five years after the flare. Both luminosity and dust temperature decreases with time. We interpret the mid-IR emission as the infrared echo of the tidal disruption flare. We estimate the UV luminosity at the peak flare to be 1 to 30 times $10^{44}$ erg s$^{-1}$ and for warm dust masses to be in the range of 0.05-1.3 Msun within a few parsecs. Our results suggest that the mid-infrared echo is a general signature of TDE in the gas-rich environment.

The streaming instability is a mechanism to concentrate solid particles into overdense filaments that undergo gravitational collapse and form planetesimals. However, it remains unclear how the initial mass function of these planetesimals depends on the box dimensions of numerical simulations. To resolve this, we perform simulations of planetesimal formation with the largest box dimensions to date, allowing planetesimals to form simultaneously in multiple filaments that can only emerge within such large simulation boxes. In our simulations, planetesimals with sizes between 80 km and several hundred kilometers form. We find that a power law with a rather shallow exponential cutoff at the high-mass end represents the cumulative birth mass function better than an integrated power law. The steepness of the exponential cutoff is largely independent of box dimensions and resolution, while the exponent of the power law is not constrained at the resolutions we employ. Moreover, we find that the characteristic mass scale of the exponential cutoff correlates with the mass budget in each filament. Together with previous studies of high-resolution simulations with small box domains, our results therefore imply that the cumulative birth mass function of planetesimals is consistent with an exponentially tapered power law with a power-law exponent of approximately -1.6 and a steepness of the exponential cutoff in the range of 0.3-0.4.

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

We report the detection of a significant infrared variability of the nearest tidal disruption event (TDE) ASASSN-14li using Wide-field Infrared Survey Explorer} and newly released Near-Earth Object WISE Reactivation} data. In comparison with the quiescent state, the infrared flux is brightened by 0.12 and 0.16 magnitude in the W1 ($3.4\mu$m) and W2 ($4.6\mu$m) bands at 36 days after the optical discovery (or $\sim110$ days after the peak disruption date). The flux excess is still detectable $\sim170$ more days later. Assuming that the flare-like infrared emission is from the dust around the black hole, its blackbody temperature is estimated to be $\sim2.1\times10^3$~K, slightly higher than the dust sublimation temperature, indicating that the dust is likely located close to the dust sublimation radius. The equilibrium between the heating and radiation of the dust claims a bolometric luminosity of $\sim10^{43}-10^{45}$~\lum, comparable with the observed peak luminosity. This result has for the first time confirmed the detection of infrared emission from the dust echoes of TDEs.

Recent studies have shown that outflows in at least some broad absorption line (BAL) quasars are extended well beyond the putative dusty torus. Such outflows should be detectable in obscured quasars. We present four WISE selected infrared red quasars with very strong and peculiar ultraviolet Fe ii emission lines: strong UV Fe II UV arising from transitions to ground/low excitation levels, and very weak Fe II at wavelengths longer than 2800 {\AA}. The spectra of these quasars display strong resonant emission lines, such as C IV, Al III and Mg II but sometimes, a lack of non-resonant lines such as C III], S III and He II. We interpret the Fe II lines as resonantly scattered light from the extended outflows that are viewed nearly edge-on, so that the accretion disk and broad line region are obscured by the dusty torus, while the extended outflows are not. We show that dust free gas exposed to strong radiation longward of 912 {\AA} produces Fe II emission very similar to that observed. The gas is too cool to collisionally excite Fe II lines, accounting for the lack of optical emission. The spectral energy distribution from the UV to the mid-infrared can be modeled as emission from a clumpy dusty torus, with UV emission being reflected/scattered light either by the dusty torus or the outflow. Within this scenario, we estimate a minimum covering factor of the outflows from a few to 20% for the Fe II scattering region, suggesting that Fe II BAL quasars are at a special stage of quasar evolution.

We report on simultaneous NuSTAR and XMM-Newton observations of the magnetar 1E 1048.1$-$5937 obtained in July 2013, along with Rossi X-ray Timing Explorer (RXTE) data for the same source obtained from December 2002 through March 2012. The NuSTAR data provide a clear detection of this magnetar's persistent emission up to 20 keV. We detect a previously unreported small secondary peak in the average pulse profile in the 7-10 keV band, which grows to an amplitude comparable to that of the main peak in the 10--20 keV band. We show using RXTE data that this secondary peak is likely transient. We find that the pulsed fraction increases with energy from a value of $\sim$0.55 at $\sim$2 keV to a value of $\sim$0.75 near 8 keV but shows evidence for decreasing at higher energies. After filtering out multiple bright X-ray bursts during the observation, we find that the phase-averaged spectrum from combined NuSTAR and XMM data is well described by an absorbed double blackbody plus power-law model, with no evidence for the spectral turn-up near $\sim$10 keV as has been seen in some other magnetars. Our data allow us to rule out a spectral turn-up similar to those seen in magnetars 4U 0142+61 and 1E 2259+586 of $\Delta\Gamma >= 2$, where $\Delta\Gamma$ is the difference between the soft-band and hard-band photon indexes. The absence of a significant spectral turn-up is consistent with what has been observed from a particularly active subset of magnetars having high spin-inferred magnetic fields, and with previously reported trends suggesting the degree of spectral turn-up is correlated with spin-down rate and/or spin-inferred magnetic field strength.

We employ a length-dependent conical jet model for the jet structure and emission properties of flat spectrum radio quasar 3C 279 in the steady state. In the model, ultra-relativistic leptons are injected at the base of jet and propagate along the jet structure, non-thermal photons are produced by both synchrotron emission and inverse Comtpon scattering off synchrotron photons and external soft photons at each segment of jet. We derive the total energy spectra contribution through integrating every segment. We apply the model to the quasi-simultaneous multi-wavelength observed data of two quiescent epoches. Using the observed radio data of source, we determine the length of jet $L\sim 100$ pc, and magnetic field $B_{0}\sim 0.1-1$ G at the base of jet. Assuming a steady geometry of the jet structure and suitable physical parameters, we reproduce the multi-wavelength spectra during two quiescent observed epoches, respectively. Our results show that the initial $\gamma$-ray emission site takes place at $\sim 0.5$ pc from black hole.

Context. Very high-energy (VHE) $\gamma$-ray measurements of distant TeV blazars can be nicely explained by TeV spectra induced by ultra high-energy cosmic rays. Aims. We develop a model for a plausible origin of hard spectra in distant TeV blazars. Methods. In the model, the TeV emission in distant TeV blazars is dominated by two mixed components. The first is the internal component with the photon energy around 1 TeV produced by inverse Compton scattering of the relativistic electrons on the synchrotron photons (SSC) with a correction for extragalactic background light absorbtion and the other is the external component with the photon energy more than 1 TeV produced by the cascade emission from high-energy protons propagating through intergalactic space. Results. Assuming suitable model parameters, we apply the model to observed spectra of distant TeV blazars of 1ES 0229+200. Our results show that 1) the observed spectrum properties of 1ES 0229+200, especially the TeV $\gamma$-ray tail of the observed spectra, could be reproduced in our model and 2) an expected TeV $\gamma$-ray spectrum with photon energy $>$1 TeV of 1ES 0229+200 should be comparable with the 50-hour sensitivity goal of the Cherenkov Telescope Array (CTA) and the differential sensitivity curve for the one-year observation with the Large High Altitude Air Shower Observatory (LHAASO). Conclusions. We argue that strong evidence for the Bethe-Heitler cascades along the line of sight as a plausible origin of hard spectra in distant TeV blazars could be obtained from VHE observations with CTA, LHAASO, HAWC, and HiSCORE.

The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. Mathematically the Graph- BLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of the GraphBLAS. Graphs represent connections between vertices with edges. Matrices can represent a wide range of graphs using adjacency matrices or incidence matrices. Adjacency matrices are often easier to analyze while incidence matrices are often better for representing data. Fortunately, the two are easily connected by matrix mul- tiplication. A key feature of matrix mathematics is that a very small number of matrix operations can be used to manipulate a very wide range of graphs. This composability of small number of operations is the foundation of the GraphBLAS. A standard such as the GraphBLAS can only be effective if it has low performance overhead. Performance measurements of prototype GraphBLAS implementations indicate that the overhead is low.

Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff, particularly because of small particles and/or strong local solid concentration, that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can for example be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

In recent decades, inorganic crystals have been widely used in dark matter direct search experiments. To contribute to the understanding of the capabilities of $CsI(Na)$ and $CaF_{2}(Eu)$ crystals, a mono-energetic neutron beam is utilized to study the properties of nuclear recoils, which are expected to be similar to signals of dark matter direct detection. The quenching factor of nuclear recoils in $CsI(Na)$ and $CaF_{2}(Eu)$, as well as an improved discrimination factor between nuclear recoils and ${\gamma}$ backgrounds in $CsI(Na)$, are reported.

We present a qualitative analysis of the variability of quasar broad absorption lines using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 10. We confirm that variations of absorption lines are highly coordinated among different components of the same ion or the same absorption component of different ions for C IV, Si IV and N V. Furthermore, we show that the equivalent widths of the lines decrease or increase statistically when the continuum brightens or dims. This is further supported by the synchronized variations of emission and absorption line equivalent width, when the well established intrinsic Baldwin effect for emission lines is taken into account. We find that the emergence of an absorption component is usually accompanying with dimming of the continuum while the disappearance of an absorption line component with brightening of the continuum. This suggests that the emergence or disappearance of a C IV absorption component is only the extreme case, when the ionic column density is very sensitive to continuum variations or the continuum variability amplitude is larger. These results support the idea that absorption line variability is driven mainly by changes in the gas ionization in response to continuum variations, that the line-absorbing gas is highly ionized, and in some extreme cases, too highly ionized to be detected in UV absorption lines. Due to uncertainties in the spectroscopic flux calibration, we cannot quantify the fraction of quasars with asynchronized continuum and absorption line variations.

We present correlations between 9 CO transition ($J=4-3$ to $12-11$) and beam-matched far-infrared (Far-IR) luminosities ($L_{\mathrm{FIR},\,b}$) among 167 local galaxies, using {\it{Herschel}} Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE; FTS) spectroscopic data and Photoconductor Array Camera and Spectrometer (PACS) photometry data. We adopt entire-galaxy FIR luminosities ($L_{\mathrm{FIR},\,e}$) from the {\it{IRAS}} Revised Bright Galaxy Sample and correct to $L_{\mathrm{FIR},\,b}$ using PACS images to match the varying FTS beam sizes. All 9 correlations between $L'_{\mathrm{CO}}$ and $L_{\mathrm{FIR},\,b}$ are essentially linear and tight ($\sigma$=0.2-0.3 dex dispersion), even for the highest transition, $J=12-11$. This supports the notion that the star formation rate (SFR) is linearly correlated with the dense molecular gas ($n_{\mathrm{H}_2}\gtrsim10^{4-6}\,cm^{-3}$). We divide the entire sample into three subsamples and find that smaller sample sizes can induce large differences in the correlation slopes. We also derive an average CO spectral line energy distribution (SLED) for the entire sample and discuss the implied average molecular gas properties for these local galaxies. We further extend our sample to high-{\it{z}} galaxies with CO $J=5-4$ data from the literature as an example, including submillimeter galaxies (SMGs) and "normal" star-forming BzKs. BzKs have similar FIR/CO(5-4) ratios as that of local galaxies, an follow well the locally-determined correlation, whereas SMG ratios fall around or slightly above the local correlation with large uncertainties. Finally, by including Galactic CO($J=10-9$) data as well as very limited high-{\it{z}} CO $J=10-9$ data, we verify that the CO(10-9) -- FIR correlation successfully extends to Galactic young stellar objects, suggesting that linear correlations are valid over 15 orders of magnitude.

The bright TeV source HESS J1640-465 is positionally coincident with the young SNR G338.3-0.0, and the nearby HESS J1641-463 with TeV gamma-ray emission seems to be closely associated with it. Based on the nonlinear diffusion shock acceleration (NLDSA) model, we explore the emission from these two TeV sources, the particle diffusion is assumed to be different inside and outside the absorbing boundary of the particles accelerated in the SNR shock. The results indicate that (1) the GeV to TeV emission from the region of the HESS J1640-465 is produced as a result of the particle acceleration inside the SNR G338.3-0.0; and (2) the runaway cosmic-ray particles outside the SNR are interacting with nearby dense molecular cloud (MC) at the region of the HESS J1641-463, corresponding $\pi^0$ decay gamma-ray in proton-proton collision contribute to the TeV emission from the HESS J1641-463. Also we investigate the possible X-ray emission in molecular cloud from synchrotron procedure by secondary $e^\pm$ produced through escaped protons interaction with the MC.

We present temporal and spectral analysis of simultaneous 0.5-79 keV Swift-XRT and NuSTAR observations of the magnetar 4U 0142+61. The pulse profile changes significantly with photon energy between 3 and 35 keV. The pulse fraction increases with energy, reaching a value of ~20%, similar to that observed in 1E 1841-045 and much lower than the ~80% pulse fraction observed in 1E 2259+586. We do not detect the 55-ks phase modulation reported in previous Suzaku-HXD observations. The phase-averaged spectrum of 4U 0142+61 above 20 keV is dominated by a hard power law with a photon index, $\Gamma$ ~ 0.65, and the spectrum below 20 keV can be described by two blackbodies, a blackbody plus a soft power law, or by a Comptonized blackbody model. We study the full phase-resolved spectra using the electron-positron outflow model of Beloborodov (2013). Our results are consistent with the parameters of the active j-bundle derived from INTEGRAL data by Hascoet et al. (2014). We find that a significant degeneracy appears in the inferred parameters if the footprint of the j-bundle is allowed to be a thin ring instead of a polar cap. The degeneracy is reduced when the footprint is required to be the hot spot inferred from the soft X-ray data.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

Analysis of cosmic microwave background radiation fluctuations favors an effective number of neutrinos, $N_\nu>3$. This motivates a reinvestigation of the neutrino freeze-out process. Here we characterize the dependence of $N_\nu$ on the Standard Model (SM) parameters that govern neutrino freeze-out. We show that $N_\nu$ depends on a combination $\eta$ of several natural constants characterizing the relative strength of weak interaction processes in the early Universe and on the Weinberg angle $\sin^2\theta_W$. We determine numerically the dependence $N_\nu(\eta,\sin^2\theta_W)$ and discuss these results. The extensive numerical computations are made possible by two novel numerical procedures: a spectral method Boltzmann equation solver adapted to allow emerging chemical non-equilibrium, and a method to evaluate Boltzmann equation collision integrals that generates a smooth integrand.

We report the identification of an unusual absorption line system in the quasar SDSS J080248.18$+$551328.9 and present a detailed study of the system, incorporating follow-up optical and NIR spectroscopy. A few tens of absorption lines are detected, including He I*, Fe II* and Ni II* that arise from metastable or excited levels, as well as resonant lines in Mg I, Mg II, Fe II, Mn II, and Ca II. All of the isolated absorption lines show the same profile of width $\Delta v\sim 1,500$km s$^{-1}$ centered at a common redshift as that of the quasar emission lines, such as [O II], [S II], and hydrogen Paschen and Balmer series. With narrow Balmer lines, strong optical Fe II multiplets, and weak [O III] doublets, its emission line spectrum is typical for that of a narrow-line Seyfert 1 galaxy (NLS1). We have derived reliable measurements of the gas-phase column densities of the absorbing ions/levels. Photoionization modeling indicates that the absorber has a density of $n_{\rm H} \sim (1.0-2.5)\times 10^5~ {\rm cm}^{-3}$ and a column density of $N_{\rm H} \sim (1.0-3.2)\times 10^{21} \sim {\rm cm}^{-2}$, and is located at $R\sim100-250$ pc from the central super-massive black hole. The location of the absorber, the symmetric profile of the absorption lines, and the coincidence of the absorption and emission line centroid jointly suggest that the absorption gas is originated from the host galaxy and is plausibly accelerated by stellar processes, such as stellar winds \zhy{and/or} supernova explosions. The implications for the detection of such a peculiar absorption line system in an NLS1 are discussed in the context of co-evolution between super-massive black hole growth and host galaxy build-up.

Aims. Two-dimensional MHD simulations are used to model the emission properties of TeV-bright shell-type supernova remnants (SNRs) and to explore their nature. Methods. In the leptonic scenario for the TeV emission, the $\gamma$-ray emission is produced via Inverse Compton scattering of background soft photons by high-energy electrons accelerated by the shocks of the SNRs. The TeV emissivity is proportional to the magnetic field energy density and MHD simulations can be used to model the TeV structure of such remnants directly. 2D MHD simulations for SNRs are then performed under the assumption that the ambient interstellar medium is turbulent with the magnetic field and density fluctuations following a Kolmogorov-like power-law spectrum. Results. (1) As expected, these simulations confirm early 1D and 2D modelings of these sources, namely the hydrodynamical evolution of the shock waves and amplification of magnetic field by Rayleigh-Taylor convective flows and by shocks propagating in a turbulent medium; (2) We reproduce rather complex morphological structure for $\gamma$-rays, suggesting intrinsic variations of the source morphology not related to the structure of the progenitor and environment; (3)The observed radial profile of several remnants are well reproduced with an ambient medium density of $0.1-1$ cm$^{-3}$. An even lower ambient density leads to a sharper drop of the TeV brightness with radius than what is observed near the outer edge of these remnants. Conclusions. In a turbulent background medium, we can reproduce the observed characteristics of several shell-type TeV SNRs with reasonable parameters except for a higher ambient density than that inferred from X-ray observations.

We present the results of a Palomar Transient Factory (PTF) archival search for blue transients which lie in the magnitude range between "normal" core-collapse and superluminous supernovae (i.e. with $-21\,{\leq}M_{R\,(peak)}\,{\leq}-19$). Of the six events found after excluding all interacting Type~IIn and Ia-CSM supernovae, three (PTF09ge, 09axc and 09djl) are coincident with the centers of their hosts, one (10iam) is offset from the center, and for two (10nuj and 11glr) a precise offset can not be determined. All the central events have similar rise times to the He-rich tidal disruption candidate PS1-10jh, and the event with the best-sampled light curve also has similar colors and power-law decay. Spectroscopically, PTF09ge is He-rich, while PTF09axc and 09djl display broad hydrogen features around peak magnitude. All three central events are in low star-formation hosts, two of which are E+A galaxies. Our spectrum of the host of PS1-10jh displays similar properties. PTF10iam, the one offset event, is different photometrically and spectroscopically from the central events and its host displays a higher star formation rate. Finding no obvious evidence for ongoing galactic nuclei activity or recent star formation, we conclude that the three central transients likely arise from the tidal disruption of a star by a super-massive black hole. We compare the spectra of these events to tidal disruption candidates from the literature and find that all of these objects can be unified on a continuous scale of spectral properties. The accumulated evidence of this expanded sample strongly supports a tidal disruption origin for this class of nuclear transients.

The inorganic CsI(Tl) crystal scintillator is a candidate served as an anti-compton detector for China Dark matter EXperiment (CDEX). Studying the intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The timing, energy and spatial correlations, as well as the capability of pulse shape discrimination provide powerful methods for the measurement of intrinsic radiopurities. The experimental design, detector performance and event-selection algorithms are described. A total of 359*3 kg-days data from three prototypes of CsI(Tl) crystals were taken at China Jinping underground laboratory where offers a good shielding environment. The contamination levels of internal isotopes from Cs-137, Th-232 and U-238 series, as well as the upper bounds of U-235 series are reported. Identification of the whole alpha8 peaks from U/Th decay chains and derivation of those corresponding quenching factors are achieved.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

The streaming instability is a promising mechanism to overcome the barriers in direct dust growth and lead to the formation of planetesimals. Most previous studies of the streaming instability, however, were focused on a local region of a protoplanetary disk with a limited simulation domain such that only one filamentary concentration of solids has been observed. The characteristic separation between filaments is therefore not known. To address this, we conduct the largest-scale simulations of the streaming instability to date, with computational domains up to 1.6 gas scale heights both horizontally and vertically. The large dynamical range allows the effect of vertical gas stratification to become prominent. We observe more frequent merging and splitting of filaments in simulation boxes of high vertical extent. We find multiple filamentary concentrations of solids with an average separation of about 0.2 local gas scale heights, much higher than the most unstable wavelength from linear stability analysis. This measures the characteristic separation of planetesimal forming events driven by the streaming instability and thus the initial feeding zone of planetesimals.

The magnetorotational instability (MRI) drives magnetized turbulence in sufficiently ionized regions of protoplanetary disks, leading to mass accretion. The dissipation of the potential energy associated with this accretion determines the thermal structure of accreting regions. Until recently, the heating from the turbulence has only been treated in an azimuthally averaged sense, neglecting local fluctuations. However, magnetized turbulence dissipates its energy intermittently in current sheet structures. We study this intermittent energy dissipation using high resolution numerical models including a treatment of radiative thermal diffusion in an optically thick regime. Our models predict that these turbulent current sheets drive order unity temperature variations even where the MRI is damped strongly by Ohmic resistivity. This implies that the current sheet structures where energy dissipation occurs must be well resolved to correctly capture the flow structure in numerical models. Higher resolutions are required to resolve energy dissipation than to resolve the magnetic field strength or accretion stresses. The temperature variations are large enough to have major consequences for mineral formation in disks, including melting chondrules, remelting calcium-aluminum rich inclusions, and annealing silicates; and may drive hysteresis: current sheets in MRI active regions could be significantly more conductive than the remainder of the disk.

We report NuSTAR observations of the millisecond pulsar - low mass X-ray binary (LMXB) transition system PSR J1023+0038 from June and October 2013, before and after the formation of an accretion disk around the neutron star. Between June 10-12, a few days to two weeks before the radio disappearance of the pulsar, the 3-79 keV X-ray spectrum was well fit by a simple power law with a photon index of Gamma=1.17 +/-0.08 (at 90% confidence) with a 3-79 keV luminosity of 7.4+/-0.4 x 10^32 erg/s. Significant orbital modulation was observed with a modulation fraction of 36+/-10%. During the October 19-21 observation, the spectrum is described by a softer power law (Gamma=1.66+/-0.06) with an average luminosity of 5.8+/-0.2 x 10^33 erg/s and a peak luminosity of ~1.2 x 10^34 erg/s observed during a flare. No significant orbital modulation was detected. The spectral observations are consistent with previous and current multi-wavelength observations and show the hard X-ray power law extending to 79 keV without a spectral break. Sharp edged, flat bottomed `dips' are observed with widths between 30-1000 s and ingress and egress time-scales of 30-60 s. No change in hardness ratio was observed during the dips. Consecutive dip separations are log-normal in distribution with a typical separation of approximately 400 s. These dips are distinct from dipping activity observed in LMXBs. We compare and contrast these dips to observations of dips and state changes in the similar transition systems PSR J1824-2452I and XSS J1227.0-4859 and discuss possible interpretations based on the transitions in the inner disk.

We report the timing and spatial resolution from the Muon Telescope Detector (MTD) installed in the STAR experiment at RHIC. Cosmic ray muons traversing the STAR detector have an average transverse momentum of 6 GeV/c. Due to their very small multiple scattering, these cosmic muons provide an ideal tool to calibrate the detectors and measure their timing and spatial resolution. The values obtained were ~100 ps and ~1-2 cm, respectively. These values are comparable to those obtained from cosmic-ray bench tests and test beams.

We perform a model independent study of the neutrino momentum distribution at freeze-out, treating the freeze-out temperature as a free parameter. Our results imply that measurement of neutrino reheating, as characterized by the measurement of the effective number of neutrinos $N_\nu$, amounts to the determination of the neutrino kinetic freeze-out temperature within the context of the standard model of particle physics where the number of neutrino flavors is fixed and no other massless (fractional) particles arise. At temperatures on the order of the neutrino mass, we show how cosmic background neutrino properties i.e. energy density, pressure, particle density, are modified in a physically consistent way as a function of neutrino mass and $N_\nu$.

As part of the DarkSide program of direct dark matter searches using liquid argon TPCs, a prototype detector with an active volume containing 10 kg of liquid argon, DarkSide-10, was built and operated underground in the Gran Sasso National Laboratory in Italy. A critically important parameter for such devices is the scintillation light yield, as photon statistics limits the rejection of electron-recoil backgrounds by pulse shape discrimination. We have measured the light yield of DarkSide-10 using the readily-identifiable full-absorption peaks from gamma ray sources combined with single-photoelectron calibrations using low-occupancy laser pulses. For gamma lines of energies in the range 122-1275 keV, we get consistent light yields averaging 8.887+-0.003(stat)+-0.444(sys) p.e./keVee. With additional purification, the light yield measured at 511 keV increased to 9.142+-0.006(stat) p.e./keVee.

We present results of follow-up optical spectroscopic observations of seven rare, extreme coronal line emitting galaxies reported by Wang et al. (2012) with Multi-Mirror Telescope (MMT). Large variations in coronal lines are found in four objects, making them strong candidates of tidal disruption events (TDE). For the four TDE candidates, all the coronal lines with ionization status higher than [Fe VII] disappear within 5-9 years. The [Fe VII] faded by a factor of about five in one object (J0952+2143) within 4 years, whereas emerged in other two without them previously. A strong increment in the [O III] flux is observed, shifting the line ratios towards the loci of active galactic nucleus on the BPT diagrams. Surprisingly, we detect a non-canonical [O III]5007/[O III]4959 2 in two objects, indicating a large column density of O$^{2+}$ and thus probably optical thick gas. This also requires a very large ionization parameter and relatively soft ionizing spectral energy distribution (e.g. blackbody with $T < 5\times 10^4$ K). Our observations can be explained as echoing of a strong ultraviolet to soft X-ray flare caused by tidal disruption events, on molecular clouds in the inner parsecs of the galactic nuclei. Re-analyzing the SDSS spectra reveals double-peaked or strongly blue-shouldered broad lines in three of the objects, which disappeared in the MMT spectra in two objects, and faded by a factor of ten in 8 years in the remaining object with a decrease in both the line width and centroid offset. We interpret these broad lines as arising from decelerating biconical outflows. Our results demonstrate that the signatures of echoing can persist for as long as ten years, and can be used to probe the gas environment in the quiescent galactic nuclei.

Current observations of supernova remnant (SNR) RX J1713.7-3946 favor the leptonic scenario for the TeV emission, where the radio to X-ray emission is produced via the synchrotron process and the $\gamma$-ray emission is produced via the inverse Comptonization of soft background photons, and the electron distribution can be inferred from the observed $\gamma$-ray spectrum with a spectral inversion method. It is shown that the observed correlation between the X-ray and $\gamma$-ray brightness of SNR RX J1713.7-3946 can be readily explained with the assumption that the energy density of energetic electrons is proportional to that of the magnetic field in such a scenario. A 2D magnetohydrodynamic simulation is then carried out to model the overall emission spectrum. It is found that the total energy of electrons above $\sim 1$ GeV is equal to that of the magnetic field. This is the first piece of observational evidence for energy equipartition between energetic electrons and magnetic field in the downstream of strong collision-less astrophysical shocks of SNRs.

We report the first systematic study of the submillimeter water vapor rotational emission lines in infrared (IR) galaxies based on the Fourier Transform Spectrometer (FTS) data of {\it Herschel} SPIRE. Among the 176 galaxies with publicly available FTS data, 45 have at least one H$_2$O emission line detected. The H$_2$O line luminosities range from $\sim 1 \times 10^5$ $L_{\odot}$ to $\sim 5 \times 10^7 L_{\odot}$ while the total IR luminosities ($L_\mathrm{IR}$) have a similar spread ($\sim 1-300 \times 10^{10} L_{\odot}$). In addition, emission lines of H$_2$O$^+$ and H$_2^{18}$O are also detected. H$_2$O is found, for most galaxies, to be the strongest molecular emitter after CO in FTS spectra. The luminosity of the five most important H$_2$O lines is near-linearly correlated with $L_\mathrm{IR}$, no matter whether strong active galactic nucleus signature is present or not. However, the luminosity of H$_2$O($2_{11}-2_{02}$) and H$_2$O($2_{20}-2_{11}$) appears to increase slightly faster than linear with $L_\mathrm{IR}$. Although the slope turns out to be slightly steeper when $z\sim 2-4$ ULIRGs are included, the correlation is still closely linear. We find that $L_\mathrm{H_2O}/L_\mathrm{IR}$ decreases with increasing $f_{25}/f_{60}$, but see no dependence on $f_{60}/f_{100}$, possibly indicating that very warm dust contributes little to the excitation of the submillimeter H$_2$O lines. The average spectral line energy distribution (SLED) of the entire sample is consistent with individual SLEDs and the IR pumping plus collisional excitation model, showing that the strongest lines are H$_2$O($2_{02}-1_{11}$) and H$_2$O($3_{21}-3_{12}$).

CRL 618 is a well-studied pre-planetary nebula. We have mapped its central region in continuum and molecular lines with the Submillimeter Array at 350 GHz at ~ 0.3" to 0.5" resolutions. Two components are seen in 350 GHz continuum: (1) a compact emission at the center tracing the dense inner part of the H II region previously detected in 23 GHz continuum and it may trace a fast ionized wind at the base, and (2) an extended thermal dust emission surrounding the H II region, tracing the dense core previously detected in HC3N at the center of the circumstellar envelope. The dense core is dusty and may contain mm-sized dust grains. It may have a density enhancement in the equatorial plane. It is also detected in carbon chain molecules HC3N and HCN, and their isotopologues, with higher excitation lines tracing closer to the central star. It is also detected in C2H3CN toward the innermost part. Most of the emission detected here arises within ~ 630 AU (~0.7") from the central star. A simple radiative transfer model is used to derive the kinematics, physical conditions, and the chemical abundances in the dense core. The dense core is expanding and accelerating, with the velocity increasing roughly linearly from ~ 3 km/s in the innermost part to ~ 16 km/s at 630 AU. The mass-loss rate in the dense core is extremely high with a value of ~ 1.15x10^{-3} solar mass per year. The dense core has a mass of ~ 0.47 solar mass and a dynamical age of ~ 400 yrs. It could result from a recent enhanced heavy mass-loss episode that ends the AGB phase. The isotopic ratios of 12C/13C and 14N/15N are 9$\pm$4 and 150$\pm50$, respectively, both lower than the solar values.

NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in \citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_\sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.

We clarify in a quantitative way the impact that distinct chemical $T_c$ and kinetic $T_k$ freeze-out temperatures have on the reduction of the neutrino fugacity $\Upsilon_\nu$ below equilibrium, i.e. $\Upsilon_\nu<1$, and the increase of the neutrino temperature $T_\nu$ via partial reheating. We establish the connection between $\Upsilon_\nu$ and $T_k$ via the modified reheating relation $T_\nu(\Upsilon_\nu)/T_\gamma$, where $T_\gamma$ is the temperature of the background radiation. Our results demonstrate that one must introduce the chemical nonequilibrium parameter, i.e., the fugacity, $\Upsilon_\nu$, as an additional standard cosmological model parameter in the evaluation of CMB fluctuations as its value allows measurement of $T_k$.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

Observations show that radial metallicity gradients in disk galaxies are relatively shallow, if not flat, especially at large galactocentric distances and for galaxies in the high-redshift universe. Given that star formation and metal production are centrally concentrated, this requires a mechanism to redistribute metals. However, the nature of this mechanism is poorly understood, let alone quantified. To address this problem, we conduct magnetohydrodynamical simulations of a local shearing sheet of a thin, thermally unstable, gaseous disk driven by a background stellar spiral potential, including metals modeled as passive scalar fields. Contrary to what a simple \alpha\ prescription for the gas disk would suggest, we find that turbulence driven by thermal instability is very efficient at mixing metals, regardless of the presence or absence of stellar spiral potentials or magnetic fields. The timescale for homogenizing randomly distributed metals is comparable to or less than the local orbital time in the disk. This implies that turbulent mixing of metals is a significant process in the history of chemical evolution of disk galaxies.

Quasars with redshifts greater than 4 are rare, and can be used to probe the structure and evolution of the early universe. Here we report the discovery of six new quasars with $i$-band magnitudes brighter than 19.5 and redshifts between 2.4 and 4.6 from the YFOSC spectroscopy of the Lijiang 2.4m telescope in February, 2012. These quasars are in the list of $z>3.6$ quasar candidates selected by using our proposed $J-K/i-Y$ criterion and the photometric redshift estimations from the SDSS optical and UKIDSS near-IR photometric data. Nine candidates were observed by YFOSC, and five among six new quasars were identified as $z>3.6$ quasars. One of the other three objects was identified as a star and the other two were unidentified due to the lower signal-to-noise ratio of their spectra. This is the first time that $z>4$ quasars have been discovered using a telescope in China. Thanks to the Chinese Telescope Access Program (TAP), the redshift of 4.6 for one of these quasars was confirmed by the Multiple Mirror Telescope (MMT) Red Channel spectroscopy. The continuum and emission line properties of these six quasars, as well as their central black hole masses and Eddington ratios, were obtained.

Tidal disruption of stars by supermassive black holes at the centers of galaxies is expected to produce unique emission line signatures, which have not yet been explored adequately. Here we report the discovery of extremely strong coronal lines from [Fe X] up to [Fe XIV] in a sample of seven galaxies (including two recently reported cases), that we interpret as such signatures. This is the first systematic search for objects of this kind, by making use of the immense database of the Sloan Digital Sky Survey. The galaxies, which are non-active as evidenced by the narrow-line ratios, show broad emission lines of complex profiles in more than half of the sample. Both the high ionization coronal lines and the broad lines turn out to be fading on time scales of years in objects observed with spectroscopic follow-ups, suggesting their transient nature. Variations of inferred non-stellar continua, which have absolute magnitudes of at least -16 to -18 mag in the g band, are also detected in more than half of the sample. These extreme coronal line emitters reside in sub-L_* disk galaxies (-21.3 < M_i < -18.5) with small stellar velocity dispersions. The sample seems to form two distinct types based on the presence or absence of the [Fe VII] lines, with the latter having relatively low luminosities of [O III], [Fe XI], and the host galaxies. These characteristics can most naturally be understood in the context of transient accretion onto intermediate mass black holes at galactic centers following tidal disruption of stars in a gas-rich environment. We estimate the incidence of such events to be around 10^-5 per year for a galaxy with -21.5 < M_i < -18.5.

Due to the gravitational influence of density fluctuations driven by magneto-rotational instability in the gas disk, planetesimals and protoplanets undergo diffusive radial migration as well as changes in other orbital properties. The magnitude of the effect on particle orbits can have important consequences for planet formation scenarios. We use the local-shearing-box approximation to simulate an ideal, isothermal, magnetized gas disk with vertical density stratification and simultaneously evolve numerous massless particles moving under the gravitational field of the gas and the host star. We measure the evolution of the particle orbital properties, including mean radius, eccentricity, inclination, and velocity dispersion, and its dependence on the disk properties and the particle initial conditions. Although the results converge with resolution for fixed box dimensions, we find the response of the particles to the gravity of the turbulent gas correlates with the horizontal box size, up to 16 disk scale heights. This correlation indicates that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence. Based on heuristic arguments, nevertheless, the criterion L_h / R ~ O(1), where L_h is the horizontal box size and R is the distance to the host star, is proposed to possibly circumvent this conundrum. If this criterion holds, we can still conclude that magneto-rotational turbulence seems likely to be ineffective at driving either diffusive migration or collisional erosion under most circumstances.

Recently, the CO J=3-2 observational result of the envelope of the 21 micrometer PPN IRAS 07134+1005 has been reported. Assuming that the CO J=3-2 line was optically thin, the mass-loss rate of the superwind in this PPN was found to be at least 2 orders of magnitude lower than the typical range. In order to obtain a more accurate mass-loss rate, we reexamine this data and construct a radiative transfer model to compare with the data. Also, in order to better resolve the superwind, we adopt a different weighting on the data to obtain maps at higher resolution. Our result shows that the CO J=3-2 emission is located slightly further away from the central source than the mid-IR emission, probably because that the material is cooler in the outer part and thus better traced by the CO emission. At lower resolution, however, the CO emission appeared to be spatially coincident with the mid-IR emission. Our model has two components, an inner ellipsoidal shell-like superwind with an equatorial density enhancement and an outer spheroidal AGB wind. The thick torus in previous model could be considered as the dense equatorial part of our ellipsoidal superwind. With radiative transfer, our model reproduces more observed features than previous model and obtains an averaged superwind mass-loss rate of ~ 1.8 x 10^-5 solar mass per year, which is typical for a superwind. The mass-loss rate in the equatorial plane is 3 x 10^-5 solar mass per year, also the same as that derived before from modeling CO J=1-0 emission.

State-to-state differential cross sections (DCSs) for rotationally inelastic scattering of H2O by H2 have been measured at 71.2 meV (574 cm-1) and 44.8 meV (361 cm-1) collision energy using crossed molecular beams combined with velocity map imaging. A molecular beam containing variable compositions of the (J = 0, 1, 2) rotational states of hydrogen collides with a molecular beam of argon seeded with water vapor that is cooled by supersonic expansion to its lowest para or ortho rotational levels (JKaKc= 000 and 101, respectively). Angular speed distributions of fully specified rotationally excited final states are obtained using velocity map imaging. Relative integral cross sections are obtained by integrating the DCSs taken with the same experimental conditions. Experimental state-specific DCSs are compared with predictions from fully quantum scattering calculations on the most complete H2O-H2 potential energy surface. Comparison of relative total cross sections and state-specific DCSs show excellent agreement with theory in almost all details

A wideband analog correlator has been constructed for the Yuan-Tseh Lee Array for Microwave Background Anisotropy. Lag correlators using analog multipliers provide large bandwidth and moderate frequency resolution. Broadband IF distribution, backend signal processing and control are described. Operating conditions for optimum sensitivity and linearity are discussed. From observations, a large effective bandwidth of around 10 GHz has been shown to provide sufficient sensitivity for detecting cosmic microwave background variations.

Sharp density features in protoplanetary discs, for instance at the edge of a magnetically dead zone, have recently been proposed as effective barriers to slow down or even stop the problematically fast migration of planetary cores into their central star. Density features on a radial scale approaching the disc vertical scale height might not exist, however, since they could be Rayleigh (or more generally Solberg-Hoiland) unstable. Stability must be checked explicitly in one-dimensional viscous accretion disc models because these instabilities are artificially eliminated in the process of reducing the full set of axisymmetric equations. The disc thermodynamics, via the entropy stratification, and its vertical structure also influence stability when sharp density features are present. We propose the concept of Rayleigh adjustment for viscous disc models: any density feature that violates Rayleigh stability (or its generalization) should be diffused radially by hydrodynamical turbulence on a dynamical time-scale, approaching marginal stability in a quasi-continuous manner.

The dynamics of planetesimals and planetary cores may be strongly influenced by density perturbations driven by magneto-rotational turbulence in their natal protoplanetary gas disks. Using the local shearing box approximation, we perform numerical simulations of planetesimals moving as massless particles in a turbulent, magnetized, unstratified gas disk. Our fiducial disk model shows turbulent accretion characterized by a Shakura-Sunyaev viscosity parameter of $\alpha \sim 10^{-2}$, with root-mean-square density perturbations of $\sim$10%. We measure the statistical evolution of particle orbital properties in our simulations including mean radius, eccentricity, and velocity dispersion. We confirm random walk growth in time of all three properties, the first time that this has been done with direct orbital integration in a local model. We find that the growth rate increases with the box size used at least up to boxes of eight scale heights in horizontal size. However, even our largest boxes show velocity dispersions sufficiently low that collisional destruction of planetesimals should be unimportant in the inner disk throughout its lifetime. Our direct integrations agree with earlier torque measurements showing that type I migration dominates over diffusive migration by stochastic torques for most objects in the planetary core and terrestrial planet mass range. Diffusive migration remains important for objects in the mass range of kilometer-sized planetesimals. Discrepancies in the derived magnitude of turbulence between local and global simulations of magneto-rotationally unstable disks remains an open issue, with important consequences for planet formation scenarios.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

Utilizing deep Hubble Space Telescope imaging from the two largest field galaxy surveys, the Extended Groth Strip (EGS) and the COSMOS survey, we examine the structural properties, and derive the merger history for 21,902 galaxies with M_*>10^{10} M_0 at z<1.2. We examine the structural CAS parameters of these galaxies, deriving merger fractions, at 0.2<z<1.2, based on the asymmetry and clumpiness values of these systems. We find that the merger fraction between z=0.2 and z=1.2 increases from roughly f_m=0.04+/-0.01 to f_m=0.13+/-0.01. We explore several fitting formalisms for parameterising the merger fraction, and compare our results to other structural studies and pair methods within the DEEP2, VVDS, and COSMOS fields. We also re-examine our method for selecting mergers, and the inherent error budget and systematics associated with identifying mergers using structure. For galaxies selected by M_*>10^{10} M_0, the merger fraction can be parameterised by f_m = f_0*(1+z)^m with the power-law slope m=2.3+/-0.4. By using the best available z = 0 prior the slope increases to m=3.8+/-0.2, showing how critical the measurement of local merger properties are for deriving the evolution of the merger fraction. We furthermore show that the merger fraction derived through structure is roughly a factor of 3-6 higher than pair fractions. Based on the latest cosmological simulations of mergers we show that this ratio is predicted, and that both methods are likely tracing the merger fraction and rate properly. We calculate, utilising merger time scales from simulations, and previously published merger fractions that the merger rate of galaxies with M_*>10^{10} M_0 increases linearly between z = 0.7 and z = 3, and that a typical M_*>10^{10} M_0 galaxy undergoes between 1-2 major mergers at z<1.2.

Large-scale star formation in disk galaxies is hypothesized to be driven by global gravitational instability. The observed gas surface density is commonly used to compute the strength of gravitational instability, but according to this criterion star formation often appears to occur in gravitationally stable regions. One possible reason is that the stellar contribution to the instability has been neglected. We have examined the gravitational instability of the Large Magellanic Cloud (LMC) considering the gas alone, and considering the combination of collisional gas and collisionless stars. We compare the gravitationally unstable regions with the on-going star formation revealed by Spitzer observations of young stellar objects. Although only 62% of the massive young stellar object candidates are in regions where the gas alone is unstable, some 85% lie in regions unstable due to the combination of gas and stars. The combined stability analysis better describes where star formation occurs. In agreement with other observations and numerical models, a small fraction of the star formation occurs in regions with gravitational stability parameter Q > 1. We further measure the dependence of the star formation timescale on the strength of gravitational instability, and quantitatively compare it to the exponential dependence expected from numerical simulations.

Evidence for triggered star formation is difficult to establish because energy feedback from massive stars tend to erase the interstellar conditions that led to the star formation. Young stellar objects (YSOs) mark sites of {\it current} star formation whose ambient conditions have not been significantly altered. Spitzer observations of the Large Magellanic Cloud (LMC) effectively reveal massive YSOs. The inventory of massive YSOs, in conjunction with surveys of interstellar medium, allows us to examine the conditions for star formation: spontaneous or triggered. We examine the relationship between star formation and gravitational instability on a global scale, and we present evidence of triggered star formation on local scales in the LMC.

We obtain the equations that define the equilibrium of a homogeneous relativistic gas of neutrons, protons and electrons in a constant magnetic field as applied to the conditions that probably occur near the center of neutron stars. We compute the relative densities of the particles at equilibrium and the Fermi momentum of electrons in the strong magnetic field as function of the density of neutrons and the magnetic field induction. Novel features are revealed as to the ratio of the number of protons to the number of neutrons at equilibrium in the presence of large magnetic fields.

We have applied a new restoration technique to archival [O~I], H$\alpha$, and continuum HST images of DG~Tau. The restored [O~I] and H$\alpha$ images show that DG~Tau has a jet with a projected length of 25~AU and width $\leq$10~AU, and is already collimated at a projected distance of $\sim$~40~AU (0\farcs25) from the star. Such a narrow width and short collimation distance for a stellar jet places important constraints on theoretical models of jet formation.