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 present Giant Metrewave Radio Telescope (GMRT) 323 MHz radio continuum observations toward 13 radio-loud quasars at $z>5$, sampling the low-frequency synchrotron emission from these objects. Among the 12 targets successfully observed, we detected 10 above $4\sigma$ significance, while 2 remain undetected. All of the detected sources appear as point sources. Combined with previous radio continuum detections from the literature, 9 quasars have power-law spectral energy distributions throughout the radio range; for some the flux density drops with increasing frequency while it increases for others. Two of these sources appear to have spectral turnover. For the power-law-like sources, the power-law indices have a positive range between 0.18 and 0.67 and a negative values between $-0.90$ and $-0.27$. For the turnover sources, the radio peaks around $\sim1$ and $\sim10$ GHz in the rest frame, the optically thin indices are $-0.58$ and $-0.90$, and the optically thick indices are 0.50 and 1.20. A magnetic field and spectral age analysis of SDSS J114657.59+403708.6 at $z=5.01$ may indicate that the turnover is not caused by synchrotron self-absorption, but rather by free-free absorption by the high-density medium in the nuclear region. Alternatively, the apparent turnover may be an artifact of source variability. Finally, we calculated the radio loudness $R_{2500\rm\, \AA}$ for our sample, which spans a very wide range from 12$^{+13}_{-13}$ to 4982$^{+279}_{-254}$.

We report sensitive [\ion{C}{1}]~$^3P_1$--$^3P_0$ and $^{12}$CO~$J$=4--3 observations of the nearby merging galaxy NGC 6052 using the Morita (Atacama Compact) Array of ALMA. We detect $^{12}$CO~$J$=4--3 toward the northern part of NGC 6052, but [\ion{C}{1}]~$^3P_1$--$^3P_0$ is not detected with a [\ion{C}{1}]~$^3P_1$--$^3P_0$ to $^{12}$CO~$J$=4--3 line luminosity ratio of$~\lesssim0.07$. According to models of photodissociation regions, the unusual weakness of [\ion{C}{1}]~$^3P_1$--$^3P_0$ relative to $^{12}$CO~$J$=4--3 can be explained if the interstellar medium has a hydrogen density larger than $10^5\,{\rm cm}^{-3}$, conditions that might arise naturally in the ongoing merging process in NGC 6052. Its [\ion{C}{1}]~$^3P_1$--$^3P_0$ emission is also weaker than expected given the molecular gas mass inferred from previous measurements of $^{12}$CO~$J$=1--0 and $^{12}$CO~$J$=2--1. This suggests that [\ion{C}{1}]~$^3P_1$--$^3P_0$ may not be a reliable tracer of molecular gas mass in this galaxy. NGC 6052 is a unique laboratory to investigate how the merger process impacts the molecular gas distribution.

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

Using a high luminosity ($L_{\rm bol} \sim 10^{47.5} - 10^{48.3}$ erg s$^{-1}$), high redshift ($3.2 < z < 3.8$) quasar sample of 19 quasars with optical and near-infrared spectroscopy, we investigate the reliability of the CIV-based black hole mass estimates ($M_{\rm BH}$). The median logarithm of the CIV- and H$\beta$-based $M_{\rm BH}$ ratios is 0.110 dex with the scatter of 0.647 dex. The CIV-to-H$\beta$ BH mass differences are significantly correlated with the CIV FWHMs, blueshifts and asymmetries. Corrections of the CIV FWHM using the blueshift and asymmetry reduce the scatter of the mass differences by $\sim$ 0.04-0.2 dex. Quasars in our sample accrete at the Eddington ratio $R_{\rm EDD}>0.3$ and cover a considerable range of blueshifts, with 18/19 of the quasars showing CIV blueshifts (with the median value of 1126 km s$^{-1}$) and 14/19 of the quasars showing CIV blueshifts larger than 500 km s$^{-1}$. It suggests that not all quasars with high Eddington ratios show large blueshifts. The Baldwin effect between the CIV rest-frame equivalent width (REW) and the continuum luminosity at 1350 $\overset{\lower.5em\circ}{\mathrm{A}}$ is not seen, likely due to the limited luminosity range of our sample. We find a lack of flux in the red wing of the composite spectrum with larger CIV blueshift, and detect a higher ratio of [OIII] quasars with REW$_{\rm [OIII]}>5$ $\overset{\lower.5em\circ}{\mathrm{A}}$ in the subsample with lower CIV blueshift. It is more likely that they are caused by the combination of the Eddington ratio and the orientation effect.

We present an ALMA study of the ~180 brightest sources in the SCUBA-2 map of the COSMOS field from the S2COSMOS survey, as a pilot study for AS2COSMOS - a full survey of the ~1,000 sources in this field. In this pilot we have obtained 870-um continuum maps of an essentially complete sample of the brightest 182 sub-millimetre sources (S_850um=6.2mJy) in COSMOS. Our ALMA maps detect 260 sub-millimetre galaxies (SMGs) spanning a range in flux density of S_870um=0.7-19.2mJy. We detect more than one SMG counterpart in 34+/-2 per cent of sub-millimetre sources, increasing to 53+/-8 per cent for SCUBA-2 sources brighter than S_850um>12mJy. We estimate that approximately one-third of these SMG-SMG pairs are physically associated (with a higher rate for the brighter secondary SMGs, S_870um>3mJy), and illustrate this with the serendipitous detection of bright [CII] 157.74um line emission in two SMGs, AS2COS0001.1 & 0001.2 at z=4.63, associated with the highest significance single-dish source. Using our source catalogue we construct the interferometric 870um number counts at S_870um>6.2mJy. We use the extensive archival data of this field to construct the multiwavelength spectral energy distribution of each AS2COSMOS SMG, and subsequently model this emission with MAGPHYS to estimate their photometric redshifts. We find a median photometric redshift for the S_870um>6.2mJy AS2COSMOS sample of z=2.87+/-0.08, and clear evidence for an increase in the median redshift with 870-um flux density suggesting strong evolution in the bright-end of the 870um luminosity function.

Hardening of cosmic proton- and nuclei-spectra is explained by using the gluon condensation (GC) model, which states that a large amount of gluons in proton may condense near the high energy threshold. The results present the universality of the GC-effect, which may be a common origin of a series of anomalous astrophysical phenomena including the broken power-law in gamma-ray spectra, the excess in positron- and electron-spectra and hardening of proton- and nuclei-spectra.

Over the last two decades, the Submillimetre Common-User Bolometer Array (SCUBA) and SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) achieved gread success in discovering the population of dusty starburst galaxies in the early universe. The SCUBA-2 surveys at 450 micron and 850 micron set important constraints on the obscured star formation over cosmic time, and in combination of deep optical and near-IR data, allows the study of protoclusters and structure formation. However, the current submillimeter (submm) surveys by JCMT are still limited by area of sky coverage (confusion limit mapping of only a few deg^2), which prevent a systematic study of large samples of the obscured galaxy population. In this white paper, we review the studies of the submm galaxies with current submillimeter/millimeter (submm/mm) observations, and discuss the important science with the new submm instruments in the next decade. In particular, with a 10 times faster mapping speed of the new camera, we will expect deep 850 micron surveys over 10 to 100 times larger sky area to i) largely increase the sample size of submm detections toward the highest redshift, ii) improve our knowledge of galaxy and structure formation in the early universe.

The fundamental stellar atmospheric parameters T_eff and log g and 13 chemical abundances are derived for medium-resolution spectroscopy from LAMOST Medium-Resolution Survey (MRS) data sets with a deep-learning method. The neural networks we designed, named as SPCANet, precisely map LAMOST MRS spectra to stellar parameters and chemical abundances. The stellar labels derived by SPCANet are with precisions of 119 K for T_eff and 0.17 dex for log g. The abundance precision of 11 elements including [C/H], [N/H], [O/H], [Mg/H], [Al/H], [Si/H], [S/H], [Ca/H], [Ti/H], [Cr/H], [Fe/H], and [Ni/H] are 0.06~0.12 dex, while of [Cu/H] is 0.19 dex. These precisions can be reached even for spectra with signal-to-noise as low as 10. The results of SPCANet are consistent with those from other surveys such as APOGEE, GALAH and RAVE, and are also validated with the previous literature values including clusters and field stars. The catalog of the estimated parameters is available at \url{http://paperdata.china-vo.org/LAMOST/MRS_parameters_elements.csv}.

Based on statistical analysis of synchrotron polarization intensity, we study the anisotropic properties of compressible magnetohydrodynamic (MHD) turbulence. The second-order normalized structure function, quadrupole ratio modulus and anisotropic coefficient are synergistically used to characterize the anisotropy of the polarization intensity. On the basis of pre-decomposition data cubes, we first explore the anisotropy of the polarization intensity in different turbulence regimes and find that the most significant anisotropy occurs in the sub-Alfv\'enic regime. Using post-decomposition data cubes in this regime, we then study the anisotropy of the polarization intensity from Alfv\'en, slow and fast modes. Statistics of polarization intensity from Alfv\'en and slow modes demonstrate the significant anisotropy while statistics of polarization intensity from fast modes show isotropic structures, which is consistent with the earlier results provided in Cho & Lazarian (2002). As a result, both quadrupole ratio modulus and anisotropic coefficient for polarization intensities can quantitatively recover the anisotropy of underlying compressible MHD turbulence. The synergistic use of the two methods helps enhance the reliability of the magnetic field measurement.

Recently, a full-scale data processing workflow of the Square Kilometre Array (SKA) Phase 1 was successfully executed on the world's fastest supercomputer Summit, proving that scientists have the expertise, software tools and computing resources to process the SKA data. The SKA-Summit experiment shows the importance of multidisciplinary cooperation between astronomy, computer science and others communities. The SKA science cannot be achieved without the joint efforts of talents from multiple fields.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO $(8-7)$, $(9-8)$, $\rm H_{2}O (2_{0,2}-1_{1,1})$ and $\rm OH^{+} (1_{1}-0_{1})$ and NOrthern Extended Millimeter Array (NOEMA) observations of CO $(5-4)$, $(6-5)$, $(12-11)$ and $(13-12)$ towards the $z = 6.003$ quasar SDSS J231038.88+185519.7, aiming to probe the physical conditions of the molecular gas content of this source. We present the best sampled CO spectral line energy distribution (SLED) at $z = 6.003$, and analyzed it with the radiative transfer code MOLPOP-CEP. Fitting the CO SLED to a one-component model indicates a kinetic temperature $T_{\rm kin} = 228 \ \rm K$, molecular gas density $log (n(\rm H_{2})/\rm cm^{-3}$ )=4.75, and CO column density $log(N(\rm CO)/\rm cm^{-2}) =17.5$, although a two-component model better fits the data. In either case, the CO SLED is dominated by a "warm" and "dense" component. Compared to samples of local (Ultra) Luminous Infrared Galaxies ((U)LIRGs), starburst galaxies and high redshift Submillimeter Galaxies (SMGs), J2310+1855 exhibits higher CO excitation at ($J \geq 8$), like other high redshift quasars. The high CO excitation, together with the enhanced $L_{\rm H_{2}O}/ L_{IR} $, $L_{\rm H_{2}O}/ L_{CO} $ and $L_{OH^{+}}/L_{\rm H_{2}O} $ ratios, suggests that besides the UV radiation from young massive stars, other mechanisms such as shocks, cosmic rays and X-rays might also be responsible for the heating and ionization of the molecular gas. In the nuclear region probed by the molecular emissions lines, any of these mechanisms might be present due to the powerful quasar and the starburst activity.

We present Atacama Large Millimiter/submillimiter Array (ALMA) observations of eight highly excited CO (J$_{\rm up}>8$) lines and continuum emission in two $z\sim6$ quasars: SDSS J231038.88+185519.7 (hereafter J2310), for which CO(8-7), CO(9-8), and CO(17-16) lines have been observed, and ULAS J131911.29+095951.4 (J1319), observed in the CO(14-13), CO(17-16) and CO(19-18) lines. The continuum emission of both quasars arises from a compact region ($< 0.9$ kpc). By assuming a modified black-body law, we estimate dust masses of Log$(M_{\rm dust}/M_{\odot})=8.75\pm0.07$ and Log$(M_{\rm dust}/M_{\odot})=8.8\pm0.2$ and dust temperatures of $T_{\rm dust}=76\pm3~{\rm K}$ and $T_{\rm dust}=66^{+15}_{-10}~{\rm K}$, respectively for J2310 and J1319. Only CO(8-7) and CO(9-8) in J2310 are detected, while $3\sigma$ upper limits on luminosities are reported for the other lines of both quasars. The CO line luminosities and upper limits measured in J2310 and J1319 are consistent with those observed in local AGN and starburst galaxies, and other $z\sim 6$ quasars, except for SDSS J1148+5251 (J1148), the only quasar at $z=6.4$ with a previous CO(17-16) line detection. By computing the CO SLEDs normalised to the CO(6-5) line and FIR luminosities for J2310, J1319, and J1149, we conclude that different gas heating mechanisms (X-ray radiation and/or shocks) may explain the different CO luminosities observed in these $z\sim6$ quasar. Future J$_{\rm up}>8$ CO observations will be crucial to understand the processes responsible for molecular gas excitation in luminous high-$z$ quasars.

We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

We present analysis of more than one hundred large-amplitude bipolar electrostatic structures in a quasi-perpendicular supercritical Earth's bow shock crossing, measured by the Magnetospheric Multiscale spacecraft. The occurrence of the bipolar structures is shown to be tightly correlated with magnetic field gradients in the shock transition region. The bipolar structures have negative electrostatic potentials and spatial scales of a few Debye lengths. The bipolar structures propagate highly oblique to the shock normal with velocities (in the plasma rest frame) of the order of the ion-acoustic velocity. We argue that the bipolar structures are ion phase space holes produced by the two-stream instability between incoming and reflected ions. This is the first identification of the ion two-stream instability in collisionless shocks. The implications for electron acceleration are discussed.

The properties of the molecular gas can shed light on the physical conditions of quasar host galaxies and the effect of feedback from accreting supermassive black holes. We present a new CO(2-1) survey of 23 z<0.1 Palomar-Green quasars conducted with the Atacama Large Millimeter/submillimeter Array. CO emission was successfully detected in 91% (21/23) of the objects, from which we derive CO luminosities, molecular gas masses, and velocity line widths. Together with CO(1-0) measurements in the literature for 32 quasars (detection rate 53%), there are 15 quasars with both CO(1-0) and CO(2-1) measurements and in total 40 sources with CO measurements. We find that the line ratio R_21 = L'_CO(2-1)/L'_CO(1-0) is subthermal, broadly consistent with nearby galaxies and other quasars previously studied. No clear correlation is found between R_21 and the intensity of the interstellar radiation field or the luminosity of the active nucleus. As with the general galaxy population, quasar host galaxies exhibit a strong, tight, linear L_IR-L_CO relation, with a normalization consistent with that of starburst systems. We investigate the molecular-to-total gas mass fraction with the aid of total gas masses inferred from dust masses previously derived from infrared observations. Although the scatter is considerable, the current data do not suggest that the CO-to-H_2 conversion factor of quasar host galaxies significantly differs from that of normal star-forming galaxies.

We propose a revision of the system developed by L'epine et al. (2007) for spectroscopic M subdwarf classification. Based on an analysis of subdwarf spectra and templates from Savcheva et al. (2014), we show thatthe CaH1 feature originally proposed by Gizis (1997) is important in selecting reliable cool subdwarf spectra. This index should be used in combination with the [TiO5, CaH2+CaH3] relation provided by L\'epine et al. (2007) to avoid misclassification results. In the new system, the dwarf-subdwarf separators are first derived from a sample of more than 80,000 M dwarfs and a "labeled" subdwarf subsample, these objects being all visually identified from their optical spectra. Based on these two samples, we re-fit the initial [TiO5, CaH1] relation, and propose a new [CaOH, CaH1] relation supplementing the [TiO5, CaH1] relation to reduce the impact of uncertainty in flux calibration on classification accuracy. In addition, we recalibrate the $\zeta_{TiO/CaH}$ parameter defined in L'epine et al. (2007) to enable its successful application to LAMOST spectra. Using this new system, we select candidates from LAMOST Data Release 4 and finally identify a set of 2791 new M subdwarf stars, covering the spectral sequence from type M0 to M7. This sample contains a large number of objects located at low Galactic latitudes, especially in the Galactic anti-center direction, expanding beyond previously published halo- and thick disk-dominated samples. Besides, we detect magnetic activity in 141 objects. We present a catalog for this M subdwarf sample, including radial velocities, spectral indices and errors, activity flags, with a compilation of external data (photometric and GAIA DR2 astrometric parameters). The catalog is provided on-line, and the spectra can be retrieved from the LAMOST Data Release web portal.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

In order to find clues to the origin of the "winged" or "X-shaped" radio galaxies (XRGs) we investigate here the parent galaxies of a large sample of 106 XRGs for optical-radio axes alignment, interstellar medium, black hole mass, and large-scale environment. For 41 of the XRGs it was possible to determine the optical major axis and the primary radio axis and the strong tendency for the two axes to be fairly close is confirmed. However, several counter-examples were also found and these could challenge the widely discussed backflow diversion model for the origin of the radio wings. Comparison with a well-defined large sample of normal FR II radio galaxies has revealed that: (i) XRGs possess slightly less massive central black holes than the normal radio galaxies (average masses being log$M_{\rm BH} \sim$ 8.81 $M_{\odot}$ and 9.07 $M_{\odot}$, respectively); (ii) a much higher fraction of XRGs ($\sim$ 80%) exhibits red mid-IR colors ($W2 - W3 > 1.5$), indicating a population of young stars and/or an enhanced dust mass, probably due to relatively recent galaxy merger(s). A comparison of the large-scale environment (i.e., within $\sim$ 1 Mpc) shows that both XRGs and FRII radio galaxies inhabit similarly poor galaxy clustering environments (medium richness being 8.94 and 11.87, respectively). Overall, the origin of XRGs seems difficult to reconcile with a single dominant physical mechanism and competing mechanisms seem prevalent.

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

An unexpected strong geomagnetic storm occurred on 2018 August 26, which was caused by a slow coronal mass ejection (CME) from a gradual eruption of a large quiet-region filament. We investigate the eruption and propagation characteristics of this CME in relation to the strong geomagnetic storm with remote sensing and in situ observations. Coronal magnetic fields around the filament are extrapolated and compared with EUV observations. We determine the propagation direction and tilt angle of the CME flux rope near the Sun using a graduated cylindrical shell (GCS) model and the Sun-to-Earth kinematics of the CME with wide-angle imaging observations from STEREO A. We reconstruct the flux-rope structure using a Grad-Shafranov technique based on the in situ measurements at the Earth and compare it with those from solar observations and the GCS results. Our conclusions are as follows: (1) the eruption of the filament was unusually slow and occurred in the regions with relatively low critical heights of the coronal field decay index; (2) the axis of the CME flux rope rotated in the corona as well as in interplanetary space, which tended to be aligned with the local heliospheric current sheet; (3) the CME was bracketed between slow and fast solar winds, which enhanced the magnetic field inside the CME at 1 AU; (4) the geomagnetic storm was caused by the enhanced magnetic field and a southward orientation of the flux rope at 1 AU from the rotation of the flux rope.

The future of exoplanet detection lies in the mid-infrared (MIR). The MIR region contains the blackbody peak of both hot and habitable zone exoplanets, making the contrast between starlight and planet light less extreme. It is also the region where prominent chemical signatures indicative of life exist, such as ozone at 9.7 microns. At a wavelength of 4 microns the difference in emission between an Earth-like planet and a star like our own is 80 dB. However a jovian planet, at the same separation exhibits 60 dB of contrast, or only 20 dB if it is hot due to its formation energy or being close to its host star. A two dimensional nulling interferometer, made with chalcogenide glass, has been measured to produce a null of 20 dB, limited by scattered light. Measures to increase the null depth to the theoretical limit of 60 dB are discussed.

In order to have a comprehensive view of the propagation and evolution of coronal mass ejections (CMEs) from the Sun to deep interplanetary space beyond 1 au, we carry out a kinematic analysis of 7 CMEs in solar cycle 23. The events are required to have coordinated coronagraph observations, interplanetary type II radio bursts, and multi-point in-situ measurements at the Earth and Ulysses. A graduated cylindrical shell model, an analytical model without free parameters and a magnetohydrodynamic model are used to derive CME kinematics near the Sun, to quantify the CME/shock propagation in the Sun-Earth space, and to connect in-situ signatures at the Earth and Ulysses, respectively. We find that each of the 7 CME-driven shocks experienced a major deceleration before reaching 1 au and thereafter propagated with a gradual deceleration from the Earth to larger distances. The resulting CME/shock propagation profile for each case is roughly consistent with all the data, which verifies the usefulness of the simple analytical model for CME/shock propagation in the heliosphere. The statistical analysis of CME kinematics indicates a tendency that the faster the CME, the larger the deceleration, and the shorter the deceleration time period within 1 au. For several of these events, the associated geomagnetic storms were mainly caused by the southward magnetic fields in the sheath region. In particular, the interaction between a CME-driven shock and a preceding ejecta significantly enhanced the preexisting southward magnetic fields and gave rise to a severe complex geomagnetic storm.

Massive, quiescent galaxies are already observed at redshift $z\sim4$, i.e. $\sim$1.5 Gyr after the Big Bang. Current models predict them to be formed via massive, gas-rich mergers at $z>6$. Recent ALMA observations of the cool gas and dust in $z\gtrsim$6 quasars have discovered [CII]- and far infrared-bright galaxies adjacent to several quasars. In this work, we present sensitive imaging and spectroscopic follow-up observations, with HST/WFC3, Spitzer/IRAC, VLT/MUSE, Magellan/FIRE and LBT/LUCI-MODS, of ALMA-detected, dust-rich companion galaxies of four quasars at $z\gtrsim 6$, specifically acquired to probe their stellar content and unobscured star formation rate. Three companion galaxies do not show significant emission in the observed optical/IR wavelength range. The photometric limits suggest that these galaxies are highly dust-enshrouded, with unobscured star formation rates SFR$_{\rm UV}<$few M$_{\odot}$ yr$^{-1}$, and a stellar content of $M_{*}<$10$^{10}$ M$_{\odot}$ yr$^{-1}$. However, the companion to PJ167-13 shows bright rest-frame UV emission (F140W AB = 25.48). Its SED resembles that of a star-forming galaxy with a total SFR$\sim$50 M$_{\odot}$ yr$^{-1}$ and $M_{*}\sim 9 \times 10^{9}$ M$_{\odot}$. All the companion sources are consistent with residing on the galaxy main sequence at $z\sim$6. Additional, deeper data from future facilities, such as JWST, are needed in order to characterize these gas-rich sources in the first Gyr of cosmic history.

We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of the galaxies. However, the measured [CII] deficits are fully consistent with those found at lower redshifts. Four of the galaxies show [CII] velocity fields that are consistent with ordered rotation, while the remaining six galaxies show no clear velocity gradient. All galaxies have high (~80-200 km/s) velocity dispersions, consistent with the interpretation that the interstellar medium (ISM) of these high redshift galaxies is turbulent. By fitting the galaxies with kinematic models, we estimate the dynamical mass of these systems, which range between (0.3 -> 5.4) x 1E10 Msun. For the three closest separation galaxy pairs, we observe dust and [CII] emission from gas in between and surrounding the galaxies, which is an indication that tidal interactions are disturbing the gas in these systems. Although gas exchange in these tidal interactions could power luminous quasars, the existence of quasars in host galaxies without nearby companions suggests that tidal interactions are not the only viable method for fueling their active centers. These observations corroborate the assertion that accreting supermassive black holes do not substantially contribute to the [CII] and dust continuum emission of the quasar host galaxies, and showcase the diverse ISM properties of galaxies when the universe was less than one billion years old.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

The Laser Ranging Interferometer (LRI) instrument on the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has provided the first laser interferometric range measurements between remote spacecraft, separated by approximately 220 km. Autonomous controls that lock the laser frequency to a cavity reference and establish the 5 degree of freedom two-way laser link between remote spacecraft succeeded on the first attempt. Active beam pointing based on differential wavefront sensing compensates spacecraft attitude fluctuations. The LRI has operated continuously without breaks in phase tracking for more than 50 days, and has shown biased range measurements similar to the primary ranging instrument based on microwaves, but with much less noise at a level of $1\,{\rm nm}/\sqrt{\rm Hz}$ at Fourier frequencies above 100 mHz.

We investigate a global extreme-ultraviolet (EUV) wave associated with a coronal mass ejection (CME)-driven shock on 2017 September 10. The EUV wave is transmitted by north- and south-polar coronal holes (CHs), which is observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory A (STEREO-A) from opposite sides of the Sun. We obtain key findings on how the EUV wave interacts with multiple coronal structures, and on its connection with the CME-driven shock: (1) the transmitted EUV wave is still connected with the shock that is incurvated to the Sun, after the shock has reached the opposite side of the eruption; (2) the south CH transmitted EUV wave is accelerated inside an on-disk, low-density region with closed magnetic fields, which implies that an EUV wave can be accelerated in both open and closed magnetic field regions; (3) part of the primary EUV wavefront turns around a bright point (BP) with a bipolar magnetic structure when it approaches a dim, low-density filament channel near the BP; (4) the primary EUV wave is diffused and apparently halted near the boundaries of remote active regions (ARs) that are far from the eruption, and no obvious AR related secondary waves are detected; (5) the EUV wave extends to an unprecedented scale of ~360{\deg} in latitudes, which is attributed to the polar CH transmission. These results provide insights into the effects of coronal density and magnetic field distributions on the evolution of an EUV wave, and into the connection between the EUV wave and the associated CME-driven shock.

Among more than 200 quasars known at $z\gtrsim6$, only one object, J0100+2802 (z=6.327), was found hosting a $>10^{10}M_\odot$ super-massive black hole (SMBH). In order to investigate the host galaxy properties of J0100+2802, we performed multi-band ALMA observations, aiming at mapping the dust continuum, [CII] and CO(6-5) emission lines with sub-kiloparsec scale resolution, as well as detecting high-J CO lines in CO(11-10), CO(10-9), and CO(7-6). The galaxy size is measured to be $R_{\rm major}=3.6\pm0.2$ kpc from the high resolution continuum observations. No ordered motion on kilo-parsec scales was found in both [CII] and CO(6-5) emissions. The velocity dispersion is measured to be 161$\pm$7 km s$^{-1}$, $\sim$3 times smaller than that estimated from the local M-$\sigma$ relation. In addition, we found that the CO emission is more concentrate (a factor of 1.8$\pm$0.4) than the [CII] emission. Together with CO(2-1) detected by VLA, we measured the CO Spectral Line Energy Distribution (SLED), which is best fit by a two-components model, including a cool component at $\sim24$ K with a density of $n_{\rm (H_2)}$=10$^{4.5}$ cm$^{-3}$, and a warm component at $\sim224$ K with a density of $n_{\rm (H_2)}$=10$^{3.6}$ cm$^{-3}$. We also fitted the dust continuum with a graybody model, which indicates that it has either a high dust emissivity $\beta\gtrsim2$ or a hot dust temperature $T_{\rm dust}\gtrsim60$ K, or a combination of both factors. The highly excited CO emission and hot dust temperature suggest that the powerful AGN in J0100+2802 could contribute to the gas and dust heating although future observations are needed to confirm this.

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $\mu$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn.

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

The lattice Hamiltonian method is developed for solving the Vlasov equation with nuclear mean-field based on the Skyrme pseudopotential up to next-to-next-to-next-to leading order. The ground states of nuclei are obtained through varying the total energy with respect to the density distribution of nucleons. Owing to the self-consistent treatment of initial nuclear ground state and the exact energy conservation in the lattice Hamiltonian method, the present framework of solving the Vlasov equation exhibits very stable nuclear ground state evolution. As a first application of the new lattice Hamiltonian Vlasov method, we explore the iso-scalar giant monopole and iso-vector giant dipole modes of finite nuclei. The obtained results are shown to be comparable to that from random-phase approximation and consistent with the experimental data, indicating the capability of the present method in dealing with the long-time near-equilibrium nuclear dynamics.

LAMOST DR5 released more than 200,000 low resolution spectra of early-type stars with S/N>50. Searching for metallic-line (Am) stars in such a large database and study of their statistical properties are presented in this paper. Six machine learning algorithms were experimented with using known Am spectra, and both the empirical criteria method(Hou et al. 2015) and the MKCLASS package(Gray et al. 2016) were also investigated. Comparing their performance, the random forest (RF) algorithm won, not only because RF has high successful rate but also it can derives and ranks features. Then the RF was applied to the early type stars of DR5, and 15,269 Am candidates were picked out. Manual identification was conducted based on the spectral features derived from the RF algorithm and verified by experts. After manual identification, 9,372 Am stars and 1,131 Ap candidates were compiled into a catalog. Statistical studies were conducted including temperature distribution, space distribution, and infrared photometry. The spectral types of Am stars are mainly between F0 and A4 with a peak around A7, which is similar to previous works. With the Gaia distances, we calculated the vertical height Z from the Galactic plane for each Am star. The distribution of Z suggests that the incidence rate of Am stars shows a descending gradient with increasing jZj. On the other hand, Am stars do not show a noteworthy pattern in the infrared band. As wavelength gets longer, the infrared excess of Am stars decreases, until little or no excess in W1 and W2 bands.

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

We present Karl G. Jansky Very Large Array (VLA) observations of the CO (2$-$1) line emission towards three far-infrared luminous quasars at $z\sim6$: SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7 with $\sim0\farcs6$ resolution and SDSS J205406.42$-$000514.8 with $\sim2\farcs1$ resolution. All three sources are detected in the CO (2$-$1) line emission -- one source is marginally resolved, and the other two appear as point sources. Measurements of the CO (2$-$1) line emission allow us to calculate the molecular gas mass even without a CO excitation model. The inferred molecular gas masses are (0.8$-$4.3) $\times$ 10$^{10}$ $M_{\odot}$. The widths and redshifts derived from the CO (2$-$1) line are consistent with previous CO (6$-$5) and [\ion{C}{2}] measurements. We also report continuum measurements using the Herschel for SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7, and for SDSS J231038.88+185519.7, data obtained at $\sim140$ and $\sim300$ GHz using the Atacama Large Millimeter/submillimeter Array (ALMA). In the case of SDSS J231038.88+185519.7, we present a detailed analysis of the spectral energy distribution and derive the dust temperature ($\sim40$ K), the dust mass ($\sim10^{9}$ $M_{\odot}$), the far-infrared luminosity (8$-$1000 $\mu$m; $\sim10^{13}$ $ L_{\odot}$) and the star formation rate (2400$-$2700 $M_{\odot}$ yr$^{-1}$). Finally, an analysis of the photo-dissociation regions associated with the three high redshift quasars indicates that the interstellar medium in these sources has similar properties to local starburst galaxies.

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

Enormous Ly$\alpha$ Nebulae (ELANe), unique tracers of galaxy density peaks, are predicted to lie at the nodes and intersections of cosmic filamentary structures. Previous successful searches for ELANe have focused on wide-field narrowband surveys, or have targeted known sources such as ultraluminous quasi-stellar-objects (QSOs) or radio galaxies. Utilizing groups of coherently strong Ly$\alpha$ absorptions (CoSLAs), we have developed a new method to identify high-redshift galaxy overdensities and have identified an extremely massive overdensity, BOSS1441, at $z=2-3$ (Cai et al. 2016a). In its density peak, we discover an ELAN that is associated with a relatively faint continuum. To date, this object has the highest diffuse Ly$\alpha$ nebular luminosity of $L_{\rm{nebula}}=5.1\pm0.1\times10^{44}$ erg s$^{-1}$. Above the 2$\sigma$ surface brightness limit of SB$_{\rm{Ly\alpha}}= 4.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, this nebula has an end-to-end spatial extent of 442 kpc. This radio-quiet source also has extended \civ\ $\lambda1549$ and \heii\ $\lambda1640$ emission on $\gtrsim30$ kpc scales. Note that the Ly$\alpha$, \heii\ and \civ\ emission all have double-peaked line profiles. Each velocity component has a full-width-half-maximum (FWHM) of $\approx700 - 1000$ km s$^{-1}$. We argue that this Ly$\alpha$ nebula could be powered by shocks due to an AGN-driven outflow or/and photoionization by a strongly obscured source.

We present ALMA band 8 observations of the [OIII] 88um line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: ~60 kpc). Each galaxy hosts dust-obscured star formation at rates > 100 M_sun/yr, but only the quasar shows evidence for an accreting 10^9 M_sun black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, ~1 Gyr after the Big Bang. Bright [OIII] 88um emission from ionized gas is detected in both systems; the positions and line-widths are consistent with earlier [CII] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [OIII] 88um/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [OIII] 5008A emission) does not dominantly determine the quasar's FIR [OIII] 88um luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources.

Based on an extended Skyrme interaction that includes the terms in relative momenta up to sixth order, corresponding to the so-called Skyrme pseudopotential up to next-to-next-to-next-to leading order (N3LO), we derive the expressions of Hamiltonian density and single nucleon potential under general non-equilibrium conditions which can be applied in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei. While the conventional Skyrme interactions, which include the terms in relative momenta up to second order, predict an incorrect behavior as a function of energy for nucleon optical potential in nuclear matter, the present extended N3LO Skyrme interaction can give a nice description for the empirical nucleon optical potential. We also construct three interaction sets with different high-density behaviors of the symmetry energy, by fitting both the empirical nucleon optical potential up to energy of $1$ GeV and the empirical properties of isospin asymmetric nuclear matter. These extended N3LO Skyrme interactions will be useful in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei at intermediate and high energies, and they can also be useful in nuclear structure studies within the mean-field model.

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

We study the magnetic field evolution in the active region (AR) 12673 that produced the largest solar flare in the past decade on 2017 September 6. Fast flux emergence is one of the most prominent features of this AR. We calculate the magnetic helicity from photospheric tangential flows that shear and braid field lines (shear-helicity), and from normal flows that advect twisted magnetic flux into the corona (emergence-helicity), respectively. Our results show that the emergence-helicity accumulated in the corona is $-1.6\times10^{43}~Mx^2$ before the major eruption, while the shear-helicity accumulated in the corona is $-6\times10^{43}~Mx^2$, which contributes about 79\% of the total helicity. The shear-helicity flux is dominant throughout the overall investigated emergence phase. Our results imply that the emerged fields initially contain relatively low helicity. Much more helicity is built up by shearing and converging flows acting on preexisted and emerging flux. Shearing motions are getting stronger with the flux emergence, and especially on both sides of the polarity inversion line of the core field region. The evolution of the vertical currents shows that most of the intense currents do not appear initially with the emergence of the flux, which implies that most of the emerging flux is probably not strongly current-carrying. The helical magnetic fields (flux rope) in the core field region are probably formed by long-term photospheric motions. The shearing and converging motions are continuously generated driven by the flux emergence. AR 12673 is a representative as photospheric motions contribute most of the nonpotentiality in the AR with vigorous flux emergence.

We present Atacama Large Millimeter Array 1mm observations of the rest-frame far-infrared (FIR) dust continuum in 27 quasars at redshifts 6.0 < z < 6.7. We detect FIR emission at >3sigma in all quasar host galaxies with flux densities at ~1900GHz in the rest-frame of 0.12 < S_rest,1900GHz < 5.9mJy, with a median (mean) flux density of 0.88mJy (1.59mJy). The implied FIR luminosities range from L_FIR = (0.27-13)x10^12 L_sun, with 74% of our quasar hosts having L_FIR > 10^12 L_sun. The estimated dust masses are M_dust = 10^7-10^9 M_sun. If the dust is heated only by star formation, then the star formation rates in the quasar host galaxies are between 50 and 2700 M_sun/yr. In the framework of the host galaxy-black hole coevolution model a correlation between ongoing black hole growth and star formation in the quasar host galaxy would be expected. However, combined with results from the literature to create a luminosity-limited quasar sample, we do not find a strong correlation between quasar UV luminosity (a proxy for ongoing black hole growth) and FIR luminosity (star formation in the host galaxy). The absence of such a correlation in our data does not necessarily rule out the coevolution model, and could be due to a variety of effects (including different timescales for black hole accretion and FIR emission).

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

Galaxy clusters trace the largest structures of the Universe and provide ideal laboratories for studying galaxy evolution and cosmology. Clusters with extended X-ray emission have been discovered at redshifts up to z ~ 2.5. Meanwhile, there has been growing interest in hunting for protoclusters, the progenitors of clusters, at higher redshifts. It is, however, very challenging to find the largest protoclusters at early times when they start to assemble. Here we report a giant protocluster of galaxies at redshift z = 5.7, when the Universe was only one billion years old. This protocluster occupies a volume of about 35x35x35 cubic co-moving megaparsecs. It is embedded in an even larger overdense region with at least 41 spectroscopically confirmed, luminous Lyman-alpha emitting galaxies (Lyman-alpha Emitters, or LAEs), including several previously reported LAEs. Its LAE density is 6.6 times the average density at z ~ 5.7. It is the only one of its kind in an LAE survey in four square degrees on the sky. Such a large structure is also rarely seen in current cosmological simulations. This protocluster will collapse into a galaxy cluster with a mass of (3.6+/-0.9) x 10^{15} solar masses, comparable to those of the most massive clusters or protoclusters known to date.

Bright high-redshift quasars ($z>6$), hosting supermassive black holes ($\rm M_{BH} > 10^8 M_{\odot}$), are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2\,mm ALMA dust continuum maps of 35 bright quasars at $6<z<7$ and search the primary beam for excess dust continuum emission from sources with L$_{\rm IR} \gtrsim{ 10^{12}}$\,L$_{\odot}$ as evidence for early protoclusters. We compare the detection rates of continuum sources at $\geq5\sigma$ significance in the fields surrounding the quasars (A$_{\rm eff}$ = 4.3\,arcmin$^2$) with millimeter number counts in blank field surveys. We discover 15 millimeter sources in the fields excluding the quasars themselves, corresponding to an overdensity $\delta_{\rm gal} \equiv (N_{\rm gal} - N_{\rm exp})/N_{\rm exp} = -0.07\pm0.56$, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [CII] overdensities found on the same scales to similarly-selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1\,cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line of sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in dust continuum.

The solar eruption on 2012 January 27 resulted in a wide-spread solar energetic particle (SEP) event observed by STEREO A and the near-Earth spacecraft (separated by 108\degree). The event was accompanied by an X-class flare, extreme-ultraviolet (EUV) wave and fast coronal mass ejection (CME). We investigate the particle release by comparing the release times of particles at the spacecraft and the times when magnetic connectivity between the source and the spacecraft was established. The EUV wave propagating to the magnetic footpoint of the spacecraft in the lower corona and the shock expanding to the open field line connecting the spacecraft in the upper corona are thought to be responsible for the particle release. We track the evolution of the EUV wave and model the propagation of the shock using EUV and white-light observations. No obvious evidence indicates that the EUV wave reached the magnetic footpoint of either STEREO A or L1-observers. Our shock modeling shows that the release time of the particles observed at L1 was consistent with the time when the shock first established contact with the magnetic field line connecting L1-observers. The release of the particles observed by STEREO A was delayed relative to the time when the shock was initially connected to STEREO A via the magnetic field line. We suggest that the particle acceleration efficiency of the portion of the shock connected to the spacecraft determines the release of energetic particles at the spacecraft.

We report the observations of a moderate but relatively intense geo-effective solar eruption on 2015 November 4 from the peripheral diffusive polarities of active region 12443. We use space-borne Solar Dynamics Observatory and ACE observations. EUV images identified helical pattern along a filament channel and we regard this channel as flux-rope structure. Flow velocity derived from tracked magnetograms infers converging motion along the polarity inversion line beneath the filament channel. An associated magnetic cancellation process was detected in the converging region. Further, the pre-eruptive EUV brightening was observed in the converging region, the most intense part of which appeared in the magnetic cancellation region. These observations imply that the converging and cancelling flux probably contributed to the formation of the helical magnetic fields associated with the flux rope. A filament-height estimation method suggests that the middle part of the filament probably lies at a low altitude and was consistent with the initial place of the eruption. A thick current channel associated with the flux rope is also determined. For a expanding thick current channel, the critical height of the decay index for torus instability lies in the range of 37 - 47 Mm. Southward magnetic fields in the sheath and the ejecta induced a geomagnetic storm with a Dst global minimum of ~-90 nT.

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

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

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

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

We present analysis of the magnetic field in seven solar flare regions accompanied by the pulsations of hard X-ray (HXR) emission. These flares were studied by Kuznetsov et al. (2016) (Paper~I), and chosen here because of the availability of the vector magnetograms for their parent active regions (ARs) obtained with the SDO/HMI data. In Paper~I, based on the observations only, it was suggested that a magnetic flux rope (MFR) might play an important role in the process of generation of the HXR pulsations. The goal of the present paper is to test this hypothesis by using the extrapolation of magnetic field with the non-linear force-free field (NLFFF) method. Having done this, we found that before each flare indeed there was an MFR elongated along and above a magnetic polarity inversion line (MPIL) on the photosphere. In two flare regions the sources of the HXR pulsations were located at the footpoints of different magnetic field lines wrapping around the central axis, and constituting an MFR by themselves. In five other flares the parent field lines of the HXR pulsations were not a part of an MFR, but surrounded it in the form of an arcade of magnetic loops. These results show that, at least in the analyzed cases, the "single flare loop" models do not satisfy the observations and magnetic field modeling, while are consistent with the concept that the HXR pulsations are a consequence of successive episodes of energy release and electron acceleration in different magnetic flux tubes (loops) of a complex AR. An MFR could generate HXR pulsations by triggering episodes of magnetic reconnection in different loops in the course of its non-uniform evolution along an MPIL. However, since three events studied here were confined flares, actual eruptions may not be required to trigger sequential particle acceleration episodes in the magnetic systems containing an MFR.

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

We use high spatial and temporal resolution observations, simultaneously obtained with the New Vacuum Solar Telescope and Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, to investigate the high-frequency oscillations above a sunspot umbra. A novel time--frequency analysis method, namely the synchrosqueezing transform (SST), is employed to represent their power spectra and to reconstruct the high-frequency signals at different solar atmospheric layers. A validation study with synthetic signals demonstrates that SST is capable to resolving weak signals even when their strength is comparable with the high-frequency noise. The power spectra, obtained from both SST and the Fourier transform, of the entire umbral region indicate that there are significant enhancements between 10 and 14 mHz (labeled as 12 mHz) at different atmospheric layers. Analyzing the spectrum of a photospheric region far away from the umbra demonstrates that this 12~mHz component exists only inside the umbra. The animation based on the reconstructed 12 mHz component in AIA 171 \AA\ illustrates that an intermittently propagating wave first emerges near the footpoints of coronal fan structures, and then propagates outward along the structures. A time--distance diagram, coupled with a subsonic wave speed ($\sim$ 49 km s$^{-1}$), highlights the fact that these coronal perturbations are best described as upwardly propagating magnetoacoustic slow waves. Thus, we first reveal the high-frequency oscillations with a period around one minute in imaging observations at different height above an umbra, and these oscillations seem to be related to the umbral perturbations in the photosphere.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

We present an investigation of the rotation and non-radial motion of a coronal mass ejection (CME) from AR 12468 on 2015 December 16 using observations from SDO, SOHO, STEREO A and Wind. The EUV and HMI observations of the source region show that the associated magnetic flux rope (MFR) axis pointed to the east before the eruption. We use a nonlinear fore-free field (NLFFF) extrapolation to determine the configuration of the coronal magnetic field and calculate the magnetic energy density distributions at different heights. The distribution of the magnetic energy density shows a strong gradient toward the northeast. The propagation direction of the CME from a Graduated Cylindrical Shell (GCS) modeling deviates from the radial direction of the source region by about 45 degr in longitude and about 30 degr in latitude, which is consistent with the gradient of the magnetic energy distribution around the AR.The MFR axis determined by the GCS modeling points southward, which has rotated counterclockwise by about 95 degr compared with the orientation of the MFR in the low corona.The MFR reconstructed by a Grad-Shafranov (GS) method at 1 AU has almost the same orientation as the MFR from the GCS modeling, which indicates that the MFR rotation occurred in the low corona. It is the rotation of the MFR that caused the intense geomagnetic storm with the minimum Dst of -155 nT. These results suggest that the coronal magnetic field surrounding the MFR plays a crucial role in the MFR rotation and propagation direction.

In this contribution, we have made use of the contemporary equation of states (EOSs) for the complete neutron star structure, and confronted them with one particular glitch constrain for the crustal moment of inertia (MOI). We find that with these EOSs, the radii of three millisecond pulsars selected by NICER: PSR J0437-4715, PSR J1614-2230, PSR J0751+1807, are all around 12.5 km. Also, a star with $M \lesssim 1.55 M_{\odot}$ would fulfill the MOI calculation for glitch constrain with the latest neutron superfluidity density, and the glitch crisis might not be present.

We present a survey of the [CII] 158 $\mu$m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 z>6 quasars using the Atacama Large Millimeter Array (ALMA) at ~1" resolution. The [CII] line was significantly detected (at >5-sigma) in 23 sources (85%). We find typical line luminosities of $L_{\rm [CII]}=10^{9-10}$ L$_\odot$, and an average line width of ~385 km/s. The [CII]-to-far-infrared luminosity ratio ([CII]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [CII] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between $2\times10^{10}$ and $2\times10^{11}$ M$_\odot$, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ($>3\times10^8$ M$_\odot$, assuming Eddington-limited accretion). In stacked ALMA [CII] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [CII] luminosity or equivalent width. This survey (with typical on-source integration times of 8 min) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

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

Until recently, determining the rotational properties of galaxies in the early universe (z>4, Universe age <1.5Gyr) was impractical, with the exception of a few strongly lensed systems. Combining the high resolution and sensitivity of ALMA at (sub-) millimeter wavelengths with the typically high strength of the [CII] 158um emission line from galaxies and long-developed dynamical modeling tools raises the possibility of characterizing the gas dynamics in both extreme starburst galaxies and normal star forming disk galaxies at z~4-7. Using a procedure centered around GIPSY's ROTCUR task, we have fit tilted ring models to some of the best available ALMA [CII] data of a small set of galaxies: the MS galaxies HZ9 & HZ10, the Damped Lyman-alpha Absorber (DLA) host galaxy ALMA J0817+1351, the submm galaxies AzTEC/C159 and COSMOS J1000+0234, and the quasar host galaxy ULAS J1319+0950. This procedure directly derives rotation curves and dynamical masses as functions of radius for each object. In one case, we present evidence for a dark matter halo of O(10^11) solar masses. We present an analysis of the possible velocity dispersions of AzTEC/C159 and ULAS J1319+0950 based on matching simulated observations to the integrated [CII] line profiles. Finally, we test the effects of observation resolution and sensitivity on our results. While the conclusions remain limited at the resolution and signal-to-noise ratios of these observations, the results demonstrate the viability of the modeling tools at high redshift, and the exciting potential for detailed dynamical analysis of the earliest galaxies, as ALMA achieves full observational capabilities.

We present a spectroscopic survey of high-redshift, luminous galaxies over four square degrees on the sky, aiming to build a large and homogeneous sample of Ly$\alpha$ emitters (LAEs) at $z\approx5.7$ and 6.5, and Lyman-break galaxies (LBGs) at $5.5<z<6.8$. The fields that we choose to observe are well-studied, such as SXDS and COSMOS. They have deep optical imaging data in a series of broad and narrow bands, allowing efficient selection of galaxy candidates. Spectroscopic observations are being carried out using the multi-object spectrograph M2FS on the Magellan Clay telescope. M2FS is efficient to identify high-redshift galaxies, owing to its 256 optical fibers deployed over a circular field-of-view 30 arcmin in diameter. We have observed $\sim2.5$ square degrees. When the program is completed, we expect to identify more than 400 bright LAEs at $z\approx5.7$ and 6.5, and a substantial number of LBGs at $z\ge6$. This unique sample will be used to study a variety of galaxy properties and to search for large protoclusters. Furthermore, the statistical properties of these galaxies will be used to probe cosmic reionization. We describe the motivation, program design, target selection, and M2FS observations. We also outline our science goals, and present a sample of the brightest LAEs at $z\approx5.7$ and 6.5. This sample contains 32 LAEs with Ly$\alpha$ luminosities higher than 10$^{43}$ erg s$^{-1}$. A few of them reach $\ge3\times10^{43}$ erg s$^{-1}$, comparable to the two most luminous LAEs known at $z\ge6$, `CR7' and `COLA1'. These LAEs provide ideal targets to study extreme galaxies in the distant universe.

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.

We establish a relation between the equation of state (EOS) of nuclear matter and the fourth-order symmetry energy $a_{\rm{sym,4}}(A)$ of finite nuclei in a semi-empirical nuclear mass formula by self-consistently considering the bulk, surface and Coulomb contributions to the nuclear mass. Such a relation allows us to extract information on nuclear matter fourth-order symmetry energy $E_{\rm{sym,4}}(\rho_0)$ at normal nuclear density $\rho_0$ from analyzing nuclear mass data. Based on the recent precise extraction of $a_{\rm{sym,4}}(A)$ via the double difference of the "experimental" symmetry energy extracted from nuclear masses, for the first time, we estimate a value of $E_{\rm{sym,4}}(\rho_0) = 20.0\pm4.6$ MeV. Such a value of $E_{\rm{sym,4}}(\rho_0)$ is significantly larger than the predictions from mean-field models and thus suggests the importance of considering the effects of beyond the mean-field approximation in nuclear matter calculations.

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

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum and [C II] 158 $\mu$m fine structure line emission towards a far-infrared-luminous quasar, ULAS J131911.29$+$095051.4 at $z=6.13$, and combine the new Cycle 1 data with ALMA Cycle 0 data. The combined data have an angular resolution $\sim$ $0.3$, and resolve both the dust continuum and the [C II] line emission on few kpc scales. The [C II] line emission is more irregular than the dust continuum emission which suggests different distributions between the dust and [C II]-emitting gas. The combined data confirm the [C II] velocity gradient that we previously detected in lower resolution ALMA image from Cycle 0 data alone. We apply a tilted ring model to the [C II] velocity map to obtain a rotation curve, and constrain the circular velocity to be 427 $\pm$ 55 km s$^{-1}$ at a radius of 3.2 kpc with an inclination angle of 34$^\circ$. We measure the dynamical mass within the 3.2 kpc region to be 13.4$_{-5.3}^{+7.8}$ $\times 10^{10}\,M_{\odot}$. This yields a black hole and host galaxy mass ratio of 0.020$_{-0.007}^{+0.013}$, which is about 4$_{-2}^{+3}$ times higher than the present-day $M_{\rm BH}$/$M_{\rm bulge}$ ratio. This suggests that the supermassive black hole grows the bulk of its mass before the formation of the most of stellar mass in this quasar host galaxy in the early universe.

We present ALMA observations of the dust continuum and [C II] 158um line emission from the z=6.0695 Lyman Break Galaxy WMH5. These observations at 0.3" spatial resolution show a compact (~3kpc) main galaxy in dust and [C II] emission, with a 'tail' of emission extending to the east by about 5kpc (in projection). The [C II] tail is comprised predominantly of two distinct sub-components in velocity, separated from the core by ~100 and 250km/s, with narrow intrinsic widths of about 80km/s, which we call 'sub-galaxies'. The sub-galaxies themselves are extended east-west by about 3kpc in individual channel images. The [C II] tail joins smoothly into the main galaxy velocity field. The [C II] line to continuum ratios are comparable for the main and sub-galaxy positions, within a factor 2. In addition, these ratios are comparable to z~5.5 LBGs. We conjecture that the WMH5 system represents the early formation of a galaxy through the accretion of smaller satellite galaxies, embedded in a smoother gas distribution, along a possibly filamentary structure. The results are consistent with current cosmological simulations of early galaxy formation, and support the idea of very early enrichment with dust and heavy elements of the accreting material.

The existence of massive ($10^{11}$ solar masses) elliptical galaxies by redshift z~4 (when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation rates exceeding 100 solar masses per year at z>6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z>6 are, with only one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than $10^9$ solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon ([CII] at a wavelength of 158 micrometres) in four galaxies at z>6 that are companions of quasars, with velocity offsets of less than 600 kilometers per second and linear offsets of less than 600 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the [CII] measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in [CII] brightness, linewidth and implied dynamical masses, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of twenty-five z>6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the [CII] luminosity function, then they can account for the population of massive elliptical galaxies at z~4 in terms of cosmic space density.

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

We investigate the coronal and interplanetary evolution of a coronal mass ejection (CME) launched on 2010 September 4 from a source region linking two active regions (ARs) 11101 and 11103, using extreme ultraviolet imaging, magnetogram, white-light and in situ observations from SDO, STEREO, SOHO, VEX and Wind. A potential-field source-surface model is employed to examine the configuration of the coronal magnetic field surrounding the source region. The graduated cylindrical shell model and a triangulation method are applied to determine the kinematics of the CME in the corona and interplanetary space. From the remote sensing and in situ observations we obtain some key results: (1) the CME was deflected in both the eastward and southward directions in the low corona by the magnetic pressure from the two ARs and possibly interacted with another ejection, which caused that the CME arrived at VEX that was longitudinally distant from the source region; (2) although VEX was closer to the Sun, the observed and derived CME arrival times at VEX are not earlier than those at Wind, which suggests the importance of determining both the frontal shape and propagation direction of the CME in interplanetary space; (3) the ICME was compressed in the radial direction while the longitudinal transverse size was extended.

In this paper, we investigate the Casacore Table Data System (CTDS) used in the casacore and CASA libraries, and methods to parallelize it. CTDS provides a storage manager plugin mechanism for third-party devel- opers to design and implement their own CTDS storage managers. Hav- ing this in mind, we looked into various storage backend techniques that can possibly enable parallel I/O for CTDS by implementing new storage managers. After carrying on benchmarks showing the excellent parallel I/O throughput of the Adaptive IO System (ADIOS), we implemented an ADIOS based parallel CTDS storage manager. We then applied the CASA MSTransform frequency split task to verify the ADIOS Storage Manager. We also ran a series of performance tests to examine the I/O throughput in a massively parallel scenario.

In this Technical Design Report (TDR) we describe the LZ detector to be built at the Sanford Underground Research Facility (SURF). The LZ dark matter experiment is designed to achieve sensitivity to a WIMP-nucleon spin-independent cross section of three times ten to the negative forty-eighth square centimeters.

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

Propagation of coronal mass ejections (CMEs) from the Sun far into interplanetary space is not well understood due to limited observations. In this study we examine the propagation characteristics of two geo-effective CMEs, which occurred on 2005 May 6 and 13, respectively. Significant heliospheric consequences associated with the two CMEs are observed, including interplanetary CMEs (ICMEs) at the Earth and Ulysses, interplanetary shocks, a long-duration type II radio burst, and intense geomagnetic storms. We use coronagraph observations from SOHO/LASCO, frequency drift of the long-duration type II burst, in situ measurements at the Earth and Ulysses, and magnetohydrodynamic (MHD) propagation of the observed solar wind disturbances at 1 AU to track the CMEs from the Sun far into interplanetary space. We find that both of the two CMEs underwent a major deceleration within 1 AU and thereafter a gradual deceleration when they propagated from the Earth to deep interplanetary space due to interactions with the ambient solar wind. The results also reveal that the two CMEs interacted with each other in the distant interplanetary space even though their launch times on the Sun were well separated. The intense geomagnetic storm for each case was caused by the southward magnetic fields ahead of the CME, stressing the critical role of the sheath region in geomagnetic storm generation, although for the first case there is a corotating interaction region involved.

We report Very Long Baseline Array (VLBA) observations of the 1.5 GHz radio continuum emission of the {\it z=6.326} quasar SDSS J010013.02+280225.8 (hereafter J0100+2802). J0100+2802 is, by far the most optically luminous, and radio-quiet quasar with the most massive black hole known at z>6. The VLBA observations have a synthesized beam size of 12.10 mas $\times$5.36 mas (FWHM), and detected the radio continuum emission from this object with a peak surface brightness of 64.6+/-9.0 micro-Jy/beam and a total flux density of 88+/-19 micro-Jy. The position of the radio peak is consistent with that from SDSS in the optical and Chandra in the X-ray. The radio source is marginally resolved by the VLBA observations. A 2-D Gaussian fit to the image constrains the source size to $\rm (7.1+/-3.5) mas x (3.1+/-1.7) mas. This corresponds to a physical scale of (40+/-20) pc x (18+/-10) pc. We estimate the intrinsic brightness temperature of the VLBA source to be T_{B}=(1.6 +/- 1.2) x 10^{7} K. This is significantly higher than the maximum value in normal star forming galaxies, indicating an AGN origin for the radio continuum emission. However, it is also significantly lower than the brightness temperatures found in highest redshift radio-loud quasars. J0100+2802 provides a unique example to study the radio activity in optically luminous and radio quiet active galactic nuclei in the early universe. Further observations at multiple radio frequencies will accurately measure the spectral index and address the dominant radiation mechanism of the radio emission.

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

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

We present the discovery of nine quasars at $z\sim6$ identified in the Sloan Digital Sky Survey (SDSS) imaging data. This completes our survey of $z\sim6$ quasars in the SDSS footprint. Our final sample consists of 52 quasars at $5.7<z\le6.4$, including 29 quasars with $z_{\rm AB}\le20$ mag selected from 11,240 deg$^2$ of the SDSS single-epoch imaging survey (the main survey), 10 quasars with $20\le z_{\rm AB}\le20.5$ selected from 4223 deg$^2$ of the SDSS overlap regions (regions with two or more imaging scans), and 13 quasars down to $z_{\rm AB}\approx22$ mag from the 277 deg$^2$ in Stripe 82. They span a wide luminosity range of $-29.0\le M_{1450}\le-24.5$. This well-defined sample is used to derive the quasar luminosity function (QLF) at $z\sim6$. After combining our SDSS sample with two faint ($M_{1450}\ge-23$ mag) quasars from the literature, we obtain the parameters for a double power-law fit to the QLF. The bright-end slope $\beta$ of the QLF is well constrained to be $\beta=-2.8\pm0.2$. Due to the small number of low-luminosity quasars, the faint-end slope $\alpha$ and the characteristic magnitude $M_{1450}^{\ast}$ are less well constrained, with $\alpha=-1.90_{-0.44}^{+0.58}$ and $M^{\ast}=-25.2_{-3.8}^{+1.2}$ mag. The spatial density of luminous quasars, parametrized as $\rho(M_{1450}<-26,z)=\rho(z=6)\,10^{k(z-6)}$, drops rapidly from $z\sim5$ to 6, with $k=-0.72\pm0.11$. Based on our fitted QLF and assuming an IGM clumping factor of $C=3$, we find that the observed quasar population cannot provide enough photons to ionize the $z\sim6$ IGM at $\sim90$\% confidence. Quasars may still provide a significant fraction of the required photons, although much larger samples of faint quasars are needed for more stringent constraints on the quasar contribution to reionization.

Cosmological simulations suggest a strong correlation between high optical-depth Ly$\alpha$ absorbers, which arise from the intergalactic medium (IGM), and 3-D mass overdensities on scales of $10-30$ $h^{-1}$ comoving Mpc. By examining the absorption spectra of $\sim$ 80,000 QSO sight-lines over a volume of 0.1 Gpc$^3$ in the Sloan Digital Sky Survey III (SDSS-III), we have identified an extreme overdensity, BOSS1441, which contains a rare group of strong Ly$\alpha$ absorbers at $z=2.32\pm 0.02$. This absorber group is associated with six QSOs at the same redshift on a 30 comoving Mpc scale. Using Mayall/MOSAIC narrowband and broadband imaging, we detect Ly$\alpha$ emitters (LAEs) down to $0.7\times L_{\rm{Ly\alpha}}^*$, and reveal a large-scale structure of Ly$\alpha$ emitters (LAEs) in this field. Our follow-up Large Binocular Telescope (LBT) observations have spectroscopically confirmed 19 galaxies in the density peak. We show that BOSS1441 has an LAE overdensity of $10.8\pm 2.6$ on a 15 comoving Mpc scale which could collapse to a massive cluster with $M\gtrsim10^{15}$ M$_\odot$ at $z\sim0$. This overdensity is among the most massive large-scale structures at $z\sim2$ discovered to date.

The 2015 March 15 coronal mass ejection as one of the two that together drove the largest geomagnetic storm of solar cycle 24 so far was associated with sympathetic filament eruptions. We investigate the relations between the different filaments involved in the eruption. A surge-like small-scale filament motion is confirmed as the trigger that initiated the erupting filament with multi-wavelength observations and using a forced magnetic field extrapolation method. When the erupting filament moved to an open magnetic field region, it experienced an obvious acceleration process and was accompanied by a C-class flare and the rise of another larger filament that eventually failed to erupt. We measure the decay index of the background magnetic field, which presents a critical height of 118 Mm. Combining with a potential field source surface extrapolation method, we analyze the distributions of the large-scale magnetic field, which indicates that the open magnetic field region may provide a favorable condition for F2 rapid acceleration and have some relation with the largest solar storm. The comparison between the successful and failed filament eruptions suggests that the confining magnetic field plays an important role in the preconditions for an eruption.

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

We analyze multi-spacecraft observations associated with the 2012 July 12 Coronal Mass Ejection (CME), covering the source region on the Sun from SDO, stereoscopic imaging observations from STEREO, magnetic field characteristics at MESSENGER, and type II radio burst and in situ measurements from Wind. A triangulation method based on STEREO stereoscopic observations is employed to determine the kinematics of the CME, and the outcome is compared with the result derived from the type II radio burst with a solar wind electron density model. A Grad-Shafranov technique is applied to Wind in situ data to reconstruct the flux-rope structure and compare it with the observation of the solar source region, which helps understand the geo-effectiveness associated with the CME structure. Conclusions are as follows: (1) the CME undergoes an impulsive acceleration, a rapid deceleration before reaching MESSENGER, and then a gradual deceleration out to 1 AU, which should be noticed in CME kinematics models; (2) the type II radio burst was probably produced from a high-density interaction region between the CME-driven shock and a nearby streamer or from the shock flank with lower heights, which implies uncertainties in the determination of CME kinematics using solely type II radio bursts; (3) the flux-rope orientation and chirality deduced from in situ reconstruction at Wind agree with those obtained from solar source observations; (4) the prolonged southward magnetic field near the Earth is mainly from the axial component of the largely southward inclined flux rope, which indicates the importance of predicting both the flux-rope orientation and magnetic field components in geomagnetic activity forecasting.

We carry out two-dimensional global particle-in-cell simulations of the interaction between the solar wind and a dipole field to study the formation of the bow shock and magnetosphere. A self-reforming bow shock ahead of a dipole field is presented by using relatively high temporal-spatial resolutions. We find that (1) the bow shock and the magnetosphere are formed and reach a quasi-stable state after several ion cyclotron periods, and (2) under the Bz southward solar wind condition the bow shock undergoes a self-reformation for low \b{eta}i and high MA. Simultaneously, a magnetic reconnection in the magnetotail is found. For high \b{eta}i and low MA, the shock becomes quasi-stationary, and the magnetotail reconnection disappears. In addition, (3) the magnetopause deflects the magnetosheath plasmas. The sheath particles injected at the quasi-perpendicular region of the bow shock can be convected to downstream of an oblique shock region. A fraction of these sheath particles can leak out from the magnetosheath at the wings of the bow shock. Hence, the downstream situation is more complicated than that for a planar shock produced in local simulations.

This is the second paper in a series on a new luminous z ~ 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z ~ 5 quasar sample to date, with 99 quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution (LEDE) model. Using the best fit QLF, the contribution of quasars to the ionizing background at z ~ 5 is found to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.

We study the solar energetic particle (SEP) event associated with the 2012 July 23 extreme solar storm, for which STEREO and the spacecraft at L1 provide multi-point remote sensing and in situ observations. The extreme solar storm, with a superfast shock and extremely enhanced ejecta magnetic fields observed near 1 AU at STEREO A, was caused by the combination of successive coronal mass ejections (CMEs). Meanwhile, energetic particles were observed by STEREO and near-Earth spacecraft such as ACE and SOHO, suggestive of a wide longitudinal spread of the particles at 1 AU. Combining the SEP observations with in situ plasma and magnetic field measurements we investigate the longitudinal distribution of the SEP event in connection with the associated shock and CMEs. Our results underscore the complex magnetic configuration of the inner heliosphere formed by solar eruptions. The examinations of particle intensities, proton anisotropy distributions, element abundance ratios, magnetic connectivity and spectra also give important clues for the particle acceleration, transport and distribution.

High-redshift quasars are important tracers of structure and evolution in the early universe. However, they are very rare and difficult to find when using color selection because of contamination from late-type dwarfs. High-redshift quasar surveys based on only optical colors suffer from incompleteness and low identification efficiency, especially at $z\gtrsim4.5$. We have developed a new method to select $4.7\lesssim z \lesssim 5.4$ quasars with both high efficiency and completeness by combining optical and mid-IR Wide-field Infrared Survey Explorer (WISE) photometric data, and are conducting a luminous $z\sim5$ quasar survey in the whole Sloan Digital Sky Survey (SDSS) footprint. We have spectroscopically observed 99 out of 110 candidates with $z$-band magnitudes brighter than 19.5 and 64 (64.6\%) of them are quasars with redshifts of $4.4\lesssim z \lesssim 5.5$ and absolute magnitudes of $-29\lesssim M_{1450} \lesssim -26.4$. In addition, we also observed 14 fainter candidates selected with the same criteria and identified 8 (57.1\%) of them as quasars with $4.7<z<5.4$ . Among 72 newly identified quasars, 12 of them are at $5.2 < z < 5.7$, which leads to an increase of $\sim$36\% of the number of known quasars at this redshift range. More importantly, our identifications doubled the number of quasars with $M_{1450}<-27.5$ at $z>4.5$, which will set strong constraints on the bright end of the quasar luminosity function. We also expand our method to select quasars at $z\gtrsim5.7$. In this paper we report the discovery of four new luminous $z\gtrsim5.7$ quasars based on SDSS-WISE selection.

As a follow-up study on Sun-to-Earth propagation of fast coronal mass ejections (CMEs), we examine the Sun-to-Earth characteristics of slow CMEs combining heliospheric imaging and in situ observations. Three events of particular interest, the 2010 June 16, 2011 March 25 and 2012 September 25 CMEs, are selected for this study. We compare slow CMEs with fast and intermediate-speed events, and obtain key results complementing the attempt of \citet{liu13} to create a general picture of CME Sun-to-Earth propagation: (1) the Sun-to-Earth propagation of a typical slow CME can be approximately described by two phases, a gradual acceleration out to about 20-30 solar radii, followed by a nearly invariant speed around the average solar wind level, (2) comparison between different types of CMEs indicates that faster CMEs tend to accelerate and decelerate more rapidly and have shorter cessation distances for the acceleration and deceleration, (3) both intermediate-speed and slow CMEs would have a speed comparable to the average solar wind level before reaching 1 AU, (4) slow CMEs have a high potential to interact with other solar wind structures in the Sun-Earth space due to their slow motion, providing critical ingredients to enhance space weather, and (5) the slow CMEs studied here lack strong magnetic fields at the Earth but tend to preserve a flux-rope structure with axis generally perpendicular to the radial direction from the Sun. We also suggest a "best" strategy for the application of a triangulation concept in determining CME Sun-to-Earth kinematics, which helps to clarify confusions about CME geometry assumptions in the triangulation and to improve CME analysis and observations.

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

We examine two successive flare eruptions (X5.4 and X1.3) on 2012 March 7 in the NOAA active region 11429 and investigate the magnetic field reconfiguration associated with the two eruptions. Using an advanced non-linear force-free field (NLFFF) extrapolation method based on the SDO/HMI vector magnetograms, we obtain a stepwise decrease in the magnetic free energy during the eruptions, which is roughly 20%-30% of the energy of the pre-flare phase. We also calculate the magnetic helicity, and suggest that the changes of the sign of the helicity injection rate might be associated with the eruptions. Through the investigation of the magnetic field evolution, we find that the appearance of the "implosion" phenomenon has a strong relationship with the occurrence of the first X-class flare. Meanwhile, the magnetic field changes of the successive eruptions with implosion and without implosion were well observed.

We study the role of the coronal magnetic field configuration of an active region in determining the propagation direction of a coronal mass ejection (CME). The CME occurred in the active region 11944 (S09W01) near the disk center on 2014 January 7 and was associated with an X1.2 flare. A new CME reconstruction procedure based on a polarimetric technique is adopted, which shows that the CME changed its propagation direction by around 28$^\circ$ in latitude within 2.5 R$_\odot$ and 43$^\circ$ in longitude within 6.5 R$_\odot$ with respect to the CME source region. This significant non-radial motion is consistent with the finding of M$\ddot{o}$stl et al. (2015). We use nonlinear force-free field (NLFFF) and potential field source surface (PFSS) extrapolation methods to determine the configurations of the coronal magnetic field. We also calculate the magnetic energy density distributions at different heights based on the extrapolations. Our results show that the active region coronal magnetic field has a strong influence on the CME propagation direction. This is consistent with the "channelling" by the active region coronal magnetic field itself, rather than deflection by nearby structures. These results indicate that the active region coronal magnetic field configuration has to be taken into account in order to determine CME propagation direction correctly.

Exploring nucleon drip lines and astrophysical rapid neutron capture process (r-process) paths in the nuclear landscape is extremely challenging in nuclear physics and astrophysics. While various models predict similar proton drip line, their predictions for neutron drip line and the r-process paths involving heavy neutron-rich nuclei exhibit a significant variation which hampers our accurate understanding of the r-process nucleosynthesis mechanism. Using microscopic density functional theory with a representative set of non-relativistic and relativistic interactions, we demonstrate for the first time that this variation is mainly due to the uncertainty of nuclear matter symmetry energy $E_{\rm{sym}}(\rho_{\rm{sc}})$ at the subsaturation cross density $\rho_{\rm{sc}}=0.11/0.16\times\rho_0$ ($\rho_0$ is saturation density), which reflects the symmetry energy of heavy nuclei. Using the recent accurate constraint on $E_{\rm{sym}}(\rho_{\rm{sc}})$ from the binding energy difference of heavy isotope pairs, we obtain quite precise predictions for the location of the neutron drip line, the r-process paths and the number of bound nuclei in the nuclear landscape. Our results have important implications on extrapolating the properties of unknown neutron-rich rare isotopes from the data on known nuclei.

The largest geomagnetic storms of solar cycle 24 so far occurred on 2015 March 17 and June 22 with $D_{\rm st}$ minima of $-223$ and $-195$ nT, respectively. Both of the geomagnetic storms show a multi-step development. We examine the plasma and magnetic field characteristics of the driving coronal mass ejections (CMEs) in connection with the development of the geomagnetic storms. A particular effort is to reconstruct the in situ structure using a Grad-Shafranov technique and compare the reconstruction results with solar observations, which gives a larger spatial perspective of the source conditions than one-dimensional in situ measurements. Key results are obtained concerning how the plasma and magnetic field characteristics of CMEs control the geomagnetic storm intensity and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed for the multi-step development of the 2015 March 17 and June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the CME flux-rope characteristics generate intense geomagnetic storms are found, which indicates that a southward flux-rope orientation is not a necessity for a strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geomagnetic storm resulted from the interaction between two successive CMEs plus the compression by a high-speed stream from behind, which is essentially the "perfect storm" scenario proposed by \citet[][i.e., a combination of circumstances results in an event of unusual magnitude]{liu14a}, so the "perfect storm" scenario may not be as rare as the phrase implies.

The transition between the supersonic solar wind and the subsonic heliosheath, the termination shock (TS), was observed by Voyager 2 (V2) on 2007 August 31-September 1 at a distance of 84 AU from the Sun. The data reveal multiple crossings of a complex, quasi-perpendicular supercritical shock. These experimental data are the starting point for a more sophisticated analysis that includes computer modeling of a shock in the presence of pickup ions (PUIs). here, we present two-dimensional (2-D) particle-in-cell (PIC) simulations of the TS including PUIs self-consistently. We also report the ion velocity distribution across the TS using the Faraday cup data from V2. A relatively complete plasma and magnetic field data set from V2 gives us the opportunity to do a full comparison between the experimental data and PIC simulation results. Our results show that: (1) The nonstationarity of the shock front is mainly caused by the ripples along the shock front and these ripples from even if the percentage of PUIs is high. (2) PUIs play a key role in the energy dissipation of the TS, and most of the incident ion dynamic energy is transferred to the thermal energy of PUIs instead of solar wind ions (SWIs). (3) The simulated composite heliosheath ion velocity distribution function is a superposition of a cold core formed by transmitted SWIs, the shoulders contributed by the hot reflected SWIs and directly transmitted PUIs, and the wings of the distribution dominated by the very hot reflected PUIs. (4) The V2 Faraday cups observed the cool core of the distribution, so they saw only a tip of the iceberg. For the evolution of the cool core distribution function across the TS, the computed results agree reasonably well with the V2experimental results.

We report the discovery of an ultra-luminous quasar J030642.51+185315.8 (hereafter J0306+1853) at redshift 5.363, which hosts a super-massive black hole (SMBH) with $M_{BH} = (1.07 \pm 0.27) \times10^{10}~M_\odot$. With an absolute magnitude $M_{1450}=-28.92$ and bolometric luminosity $L_{bol}\sim3.4\times10^{14} L_{\odot}$, J0306+1853 is one of the most luminous objects in the early Universe. It is not likely to be a beamed source based on its small flux variability, low radio loudness and normal broad emission lines. In addition, a $z=4.986$ Damped Ly$\alpha$ system (DLA) with $\rm [M/H]=-1.3\pm0.1$, among the most metal rich DLAs at $z \gtrsim 5$, is detected in the absorption spectrum of this quasar. This ultra-luminous quasar puts strong constraint on the bright-end of quasar luminosity function and massive-end of black hole mass function. It will provide a unique laboratory to the study of BH growth and the co-evolution between BH and host galaxy with multi-wavelength follow-up observations. The future high resolution spectra will give more insights to the DLA and other absorption systems along the line-of-sight of J0306+1853.

We present Karl G. Jansky Very Large Array (VLA) observations of the CO ($J = 2 \rightarrow 1$) line emission towards the $z = 6.419$ quasar SDSS J$114816.64+525150.3$ (J$1148+5251$). The molecular gas is found to be marginally resolved with a major axis of $0.9"$ (consistent with previous size measurements of the CO ($J = 7 \rightarrow 6$) emission). We observe tentative evidence for extended line emission towards the south west on a scale of ~$1.4"$, but this is only detected at $3.3\sigma$ significance and should be confirmed. The position of the molecular emission region is in excellent agreement with previous detections of low frequency radio continuum emission as well as [C ii] line and thermal dust continuum emission. These CO ($J = 2 \rightarrow 1$) observations provide an anchor for the low excitation part of the molecular line SED. We find no evidence for extended low excitation component, neither in the spectral line energy distribution nor the image. We fit a single kinetic gas temperature model of 50 K. We revisit the gas and dynamical masses in light of this new detection of a low order transition of CO, and confirm previous findings that there is no extended reservoir of cold molecular gas in J$1148+5251$, and that the source departs substantially from the low $z$ relationship between black hole mass and bulge mass. Hence, the characteristics of J$1148+5251$ at $z = 6.419$ are very similar to $z$~$2$ quasars, in the lack of a diffuse cold gas reservoir and kpc-size compactness of the star forming region.

We present Atacama Large Millimeter Array (ALMA) detections of atomic carbon line and dust continuum emission in two UV-luminous galaxies at redshift 6. The far-infrared (FIR) luminosities of these galaxies are substantially lower than similar starbursts at later cosmic epochs, indicating an evolution in the dust properties with redshift, in agreement with the evolution seen in ultraviolet (UV) attenuation by dust. The [CII] to FIR ratios are found to be higher than at low redshift showing that [CII] should be readily detectable by ALMA within the reionization epoch. One of the two galaxies shows a complex merger nature with the less massive component dominating the UV emission and the more massive component dominating the FIR line and continuum. Using the interstellar atomic carbon line to derive the systemic redshifts we investigate the velocity of Lyman alpha emission emerging from high-z galaxies. In contrast to previous work, we find no evidence for decreasing Lyman alpha velocity shifts at high-redshift. We observe an increase in velocity shifts from z$\sim$2 to z$\sim$6, consistent with the effects of increased IGM absorption.

In the present paper the results obtained in the investigation of possible diurnal effects for low-energy single-hit scintillation events of DAMA/LIBRA-phase1 (1.04 ton $\times$ yr exposure) have been analysed in terms of an effect expected in case of Dark Matter (DM) candidates inducing nuclear recoils and having high cross-section with ordinary matter, which implies low DM local density in order to fulfill the DAMA/LIBRA DM annual modulation results. This effect is due to the different Earth depths crossed by those DM candidates during the sidereal day.

We present the discovery of eight quasars at z~6 identified in the Sloan Digital Sky Survey (SDSS) overlap regions. Individual SDSS imaging runs have some overlap with each other, leading to repeat observations over an area spanning >4000 deg^2 (more than 1/4 of the total footprint). These overlap regions provide a unique dataset that allows us to select high-redshift quasars more than 0.5 mag fainter in the z band than those found with the SDSS single-epoch data. Our quasar candidates were first selected as i-band dropout objects in the SDSS imaging database. We then carried out a series of follow-up observations in the optical and near-IR to improve photometry, remove contaminants, and identify quasars. The eight quasars reported here were discovered in a pilot study utilizing the overlap regions at high galactic latitude (|b|>30 deg). These quasars span a redshift range of 5.86<z<6.06 and a flux range of 19.3<z_AB<20.6 mag. Five of them are fainter than z_AB=20 mag, the typical magnitude limit of z~6 quasars used for the SDSS single-epoch images. In addition, we recover eight previously known quasars at z~6 that are located in the overlap regions. These results validate our procedure for selecting quasar candidates from the overlap regions and confirming them with follow-up observations, and provide guidance to a future systematic survey over all SDSS imaging regions with repeat observations.

So far, roughly 40 quasars with redshifts greater than z=6 have been discovered. Each quasar contains a black hole with a mass of about one billion solar masses ($10^9 M_\odot$). The existence of such black holes when the Universe was less than 1 billion years old presents substantial challenges to theories of the formation and growth of black holes and the coevolution of black holes and galaxies. Here we report the discovery of an ultra-luminous quasar, SDSS J010013.02+280225.8, at redshift z=6.30. It has an optical and near-infrared luminosity a few times greater than those of previously known z>6 quasars. On the basis of the deep absorption trough on the blue side of the Ly $\alpha$ emission line in the spectrum, we estimate the proper size of the ionized proximity zone associated with the quasar to be 26 million light years, larger than found with other z>6.1 quasars with lower luminosities. We estimate (on the basis of a near-infrared spectrum) that the black hole has a mass of $\sim 1.2 \times 10^{10} M_\odot$, which is consistent with the $1.3 \times 10^{10} M_\odot$ derived by assuming an Eddington-limited accretion rate.

The Square Kilometre Array will be a revolutionary instrument for the study of gas in the distant Universe. SKA1 will have sufficient sensitivity to detect and image atomic 21 cm HI in individual galaxies at significant cosmological distances, complementing ongoing ALMA imaging of redshifted high-J CO line emission and far-infrared interstellar medium lines such as [CII] 157.7 um. At frequencies below ~50 GHz, observations of redshifted emission from low-J transitions of CO, HCN, HCO+, HNC, H2O and CS provide insight into the kinematics and mass budget of the cold, dense star-forming gas in galaxies. In advance of ALMA band 1 deployment (35 to 52 GHz), the most sensitive facility for high-redshift studies of molecular gas operating below 50~GHz is the Karl G. Jansky Very Large Array (VLA). Here, we present an overview of the role that the SKA could play in molecular emission line studies during SKA1 and SKA2, with an emphasis on studies of the dense gas tracers directly probing regions of active star-formation.

We present new near-infrared (IR) observations of the H$\beta\ \lambda4861$ and MgII $\lambda2798$ lines for 32 luminous quasars with $3.2<z<3.9$ using the Palomar Hale 200 inch telescope and the Large Binocular Telescope. We find that the MgII Full Width at Half Maximum (FWHM) is well correlated with the H$\beta$ FWHM, confirming itself as a good substitute for the H$\beta$ FWHM in the black hole mass estimates. The continuum luminosity at 5100 \AA\ well correlates with the continuum luminosity at 3000 \AA\ and the broad emission line luminosities (H$\beta$ and MgII). With simultaneous near-IR spectroscopy of the H$\beta$ and MgII lines to exclude the influences of flux variability, we are able to evaluate the reliability of estimating black hole masses based on the MgII line for high redshift quasars. With the reliable H$\beta$ line based black hole mass and Eddington ratio estimates, we find that the $z\sim3.5$ quasars in our sample have black hole masses $1.90\times10^{9} M_{\odot} \lesssim M_{\rm BH} \lesssim 1.37\times10^{10} M_{\odot}$, with a median of $\sim 5.14\times10^{9} M_{\odot}$ and are accreting at Eddington ratios between 0.30 and 3.05, with a median of $\sim1.12$. Assuming a duty cycle of 1 and a seed black hole mass of $10^{4} M_{\odot}$, we show that the $z\sim3.5$ quasars in this sample can grow to their estimated black hole masses within the age of the Universe at their redshifts.

Only very few z>5 quasars discovered to date are radio-loud, with a radio-to-optical flux ratio (radio-loudness parameter) higher than 10. Here we report the discovery of an optically luminous radio-loud quasar, SDSS J013127.34-032100.1 (J0131-0321 in short), at z=5.18+-0.01 using the Lijiang 2.4m and Magellan telescopes. J0131-0321 has a spectral energy distribution consistent with that of radio-loud quasars. With an i-band magnitude of 18.47 and radio flux density of 33 mJy, its radio-loudness parameter is ~100. The optical and near-infrared spectra taken by Magellan enable us to estimate its bolometric luminosity to be L_bol ~ 1.1E48 erg/s, approximately 4.5 times greater than that of the most distant quasar known to date. The black hole mass of J0131-0321 is estimated to be 2.7E9 solar masses, with an uncertainty up to 0.4 dex. Detailed physical properties of this high-redshift, radio-loud, potentially super-Eddington quasar can be probed in the future with more dedicated and intensive follow-up observations using multi-wavelength facilities.

This paper summarizes in a simple and intuitive way why the neutrons, the muons and the solar neutrinos cannot give any significant contribution to the DAMA annual modulation results. A number of these elements have already been presented in individual papers; they are recalled here. Afterwards, few simple considerations are summarized which already demonstrate the incorrectness of the claim reported in PRL 113 (2014) 081302.

On 2012 September 30 - October 1 the Earth underwent a two-step geomagnetic storm. We examine the Sun-to-Earth characteristics of the coronal mass ejections (CMEs) responsible for the geomagnetic storm with combined heliospheric imaging and in situ observations. The first CME, which occurred on 2012 September 25, is a slow event and shows an acceleration followed by a nearly invariant speed in the whole Sun-Earth space. The second event, launched from the Sun on 2012 September 27, exhibits a quick acceleration, then a rapid deceleration and finally a nearly constant speed, a typical Sun-to-Earth propagation profile for fast CMEs \citep{liu13}. These two CMEs interacted near 1 AU as predicted by the heliospheric imaging observations and formed a complex ejecta observed at Wind, with a shock inside that enhanced the pre-existing southward magnetic field. Reconstruction of the complex ejecta with the in situ data indicates an overall left-handed flux rope-like configuration, with an embedded concave-outward shock front, a maximum magnetic field strength deviating from the flux rope axis and convex-outward field lines ahead of the shock. While the reconstruction results are consistent with the picture of CME-CME interactions, a magnetic cloud-like structure without clear signs of CME interactions \citep{lugaz14} is anticipated when the merging process is finished.

In merger-driven models of massive galaxy evolution, the luminous quasar phase is expected to be accompanied by vigorous star formation in quasar host galaxies. In this paper, we use high column density Damped Lyman Alpha (DLA) systems along quasar sight lines as natural coronagraphs to directly study the far-UV (FUV) radiation from the host galaxies of luminous background quasars. We have stacked the spectra of $\sim$2,000 DLA systems (N_HI>10^{20.6} cm^{-2}) with a median absorption redshift < z > = 2.6 selected from quasars observed in the SDSS-III Baryon Oscillation Spectroscopic Survey. We detect residual flux in the dark troughs of the composite DLA spectra. The level of this residual flux significantly exceeds systematic errors in the SDSS fiber sky subtraction; furthermore, the residual flux is strongly correlated with the continuum luminosity of the background quasar, while uncorrelated with DLA column density or metallicity. We conclude that the flux could be associated with the average FUV radiation from the background quasar host galaxies (with medium redshift < z > = 3.1) that is not blocked by the intervening DLA. Assuming all of the detected flux originates from quasar hosts, for the highest quasar luminosity bin (< L >= 2.5x 10^{13} L_sun), the host galaxy has a FUV intensity of 1.5 +/- 0.2 x 10^{40} erg s^{-1} A^{-1}; this corresponds to an unobscured UV star formation rate of 9 M_sun/yr.

The results obtained in the search for possible diurnal effect in the single-hit low energy data collected by DAMA/LIBRA-phase1 (total exposure: 1.04 ton x yr) deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. are presented. At the present level of sensitivity the presence of any significant diurnal variation and of diurnal time structures in the data can be excluded for both the cases of solar and sidereal time. In particular, the diurnal modulation amplitude expected, because of the Earth diurnal motion, on the basis of the DAMA Dark Matter annual modulation results is below the present sensitivity.

The results obtained with the total exposure of 1.04 ton x yr collected by DAMA/LIBRA-phase1 deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. during 7 annual cycles (i.e. adding a further 0.17 ton x yr exposure) are presented. The DAMA/LIBRA-phase1 data give evidence for the presence of Dark Matter (DM) particles in the galactic halo, on the basis of the exploited model independent DM annual modulation signature by using highly radio-pure NaI(Tl) target, at 7.5 sigma C.L.. Including also the first generation DAMA/NaI experiment (cumulative exposure 1.33 ton x yr, corresponding to 14 annual cycles), the C.L. is 9.3 sigma and the modulation amplitude of the single-hit events in the (2-6) keV energy interval is: (0.0112 \pm 0.0012) cpd/kg/keV; the measured phase is (144 \pm 7) days and the measured period is (0.998 \pm 0.002) yr, values well in agreement with those expected for DM particles. No systematic or side reaction able to mimic the exploited DM signature has been found or suggested by anyone over more than a decade.

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars.

Experimental observations and theoretical arguments at Galaxy and larger scales have suggested that a large fraction of the Universe is composed by Dark Matter particles. This has motivated the DAMA experimental efforts to investigate the presence of such particles in the galactic halo by exploiting a model independent signature and very highly radiopure set-ups deep underground. Few introductory arguments are summarized before presenting a review of the present model independent positive results obtained by the DAMA/NaI and DAMA/LIBRA set-ups at the Gran Sasso National Laboratory of the INFN. Implications and model dependent comparisons with other different kinds of results will be shortly addressed. Some arguments put forward in literature will be confuted.

A 3-D coronal magnetic field is reconstructed for NOAA 11158 on Feb 14, 2011. A GPU-accelerated direct boundary integral equation (DBIE) method is implemented. This is about 1000 times faster than the original DBIE used on solar NLFFF modeling. Using the SDO/HMI vector magnetogram as the bottom boundary condition, the reconstructed magnetic field lines are compared with the projected EUV loop structures from different views three-dimensionally by SDO/AIA and STEREO A/B spacecraft simultaneously for the first time. They show very good agreement so that the topological configurations of the magnetic fields can be analyzed, thus its role in the flare process of the active region can be better understood. A quantitative comparison with some stereoscopically reconstructed coronal loops shows that the present averaged misalignment angles are at the same order as the state-of-the-art results obtained with reconstructed coronal loops as prescribed conditions and better than other NLFFF methods. It is found that the observed coronal loop structures can be grouped into bundles of closed and open loops with some central bright coronal loops around the polarity inversion line (PIL). The reconstructed highly-shearing magnetic field lines agree very well with the low-lying S-shaped filament channel along PIL. They are in a pivot position to all other surrounding coronal structures, and a group of electric current lines co-aligned with the central bright EUV loops overlying the filament channel is also obtained. This central lower-lying magnetic field loop system must have played a key role in powering the flare. It should be noted that while a strand-like coronal feature along PIL may be related to the filament, one cannot simply attribute all the coronal bright features along PIL to manifestation of the filament without any stereoscopically information. It shows that DBIE is rigorous and effective.

Next generation radio telescopes will require orders of magnitude more computing power to provide a view of the universe with greater sensitivity. In the initial stages of the signal processing flow of a radio telescope, signal correlation is one of the largest challenges in terms of handling huge data throughput and intensive computations. We implemented a GPU cluster based software correlator with various data distribution models and give a systematic comparison based on testing results obtained using the Fornax supercomputer. By analyzing the scalability and throughput of each model, optimal approaches are identified across a wide range of problem sizes, covering the scale of next generation telescopes.

The DAMA/LIBRA experiment, running at the Gran Sasso National Laboratory of the I.N.F.N. in Italy, has a sensitive mass of about 250 kg highly radiopure NaI(Tl). It is mainly devoted to the investigation of Dark Matter (DM) particles in the Galactic halo by exploiting the model independent DM annual modulation signature. The present DAMA/LIBRA experiment and the former DAMA/NaI one (the first generation experiment having an exposed mass of about 100 kg) have released so far results corresponding to a total exposure of 1.17 ton yr over 13 annual cycles. They provide a model independent evidence of the presence of DM particles in the galactic halo at 8.9 sigma C.L.. A short summary of the obtained results is presented and future perspectives of the experiment mentioned.

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles.

A few final comments on arXiv:1210.7548 are given to confute incorrect arguments claimed there.

We have been carrying out a systematic survey of the star formation and ISM properties in the host galaxies of z~6 quasars. Our 250 GHz observations, together with available data from the literature, yield a sample of 14 z~6 quasars that are bright in millimeter dust continuum emission with estimated FIR luminosities of a few 10^12 to 10^13 Lsun. Most of these millimeter-detected z~6 quasars have also been detected in molecular CO line emission, indicating molecular gas masses on order of 10^10 Msun. We have searched for [C II] 158 micron fine structure line emission toward four of the millimeter bright z~6 quasars with ALMA and all of them have been detected. All these results suggest massive star formation at rates of about 600 to 2000 Msun/yr over the central few kpc region of these quasar host galaxies.

We briefly remind references and arguments, already discussed in the past, which confute erroneous claims in arXiv:1210.5501.

We present one of the most ultraviolet (UV) luminous Lyman Break Galaxies (LBGs) (J1432+3358) at z=2.78, discovered in the NOAO Deep Wide-Field Survey (NDWFS) Bootes field. The R-band magnitude of J1432+3358 is 22.29 AB, more than two magnitudes brighter than typical L* LBGs at this redshift. The deep z-band image reveals two components of J1432+3358 separated by 1.0" with flux ratio of 3:1. The high signal-to-noise ratio (S/N) rest-frame UV spectrum shows Lya emission line and interstellar medium absorption lines. The absence of NV and CIV emission lines, the non-detection in X-ray and radio wavelengths and mid-infrared (MIR) colors indicate no or weak active galactic nuclei (AGN) (<10%) in this galaxy. The galaxy shows broader line profile with the full width half maximum (FWHM) of about 1000 km/s and larger outflow velocity (~500 km/s) than those of typical z~3 LBGs. The physical properties are derived by fitting the spectral energy distribution (SED) with stellar synthesis models. The dust extinction, E(B-V)=0.12, is similar to that in normal LBGs. The star formation rates (SFRs) derived from the SED fitting and the dust-corrected UV flux are consistent with each other, ~300 Msun/yr, and the stellar mass is 1.3e11 Msun. The SFR and stellar mass in J1432+3358 are about an order of magnitude higher than those in normal LBGs. The SED-fitting results support that J1432+3358 has a continuous star formation history with the star formation episode of 630 Myr. The morphology of J1432+3358 and its physical properties suggest that J1432+3358 is in an early phase of 3:1 merger process. The unique properties and the low space number density (~1e-7 Mpc^{-3})are consistent with the interpretation that such galaxies are either found in a short unobscured phase of the star formation or that small fraction of intensive star-forming galaxies are unobscured.

This paper gathers arguments and reasons why muons surviving the Gran Sasso mountain cannot mimic the Dark Matter annual modulation signature exploited by the DAMA/NaI and DAMA/LIBRA experiments. A number of these items have already been presented in individual papers. Further arguments have been addressed here in order to present a comprehensive collection and to enable a wider community to correctly approach this point.

Some theoretical and experimental aspects regarding the direct dark matter field are mentioned. In particular some arguments, which play a relevant role in the evaluation of model dependent interpretations of experimental results and in comparisons, are shortly addressed.

Based on the analysis of the microwave observations at frequency of 2.60 -- 3.80 GHz in a solar X1.3 flare event observed at Solar Broadband RadioSpectrometer in Huairou (SBRS/Huairou) on 2005 July 30, an interesting reversed drifting quasi-periodic pulsating structure (R-DPS) is confirmed. The R-DPS is mainly composed of two drifting pulsating components: one is a relatively slow very short-period pulsation (VSP) with period of about 130 -- 170 ms, the other is a relatively fast VSP with period of about 70 -- 80 ms. The R-DPS has a weak left-handed circular polarization. Based on the synthetic investigations of Reuven Ramaty High Energy Solar Spectroscopic Imaging (RHESSI) hard X-ray, Geostationary Operational Environmental Satellite (GOES) soft X-ray observation, and magnetic field extrapolation, we suggest the R-DPS possibly reflects flaring dynamic processes of the emission source regions.

We present new millimeter and radio observations of nine z~6 quasars discovered in deep optical and near-infrared surveys. We observed the 250 GHz continuum in eight of the nine objects and detected three of them. New 1.4 GHz radio continuum data have been obtained for four sources, and one has been detected. We searched for molecular CO (6-5) line emission in the three 250 GHz detections and detected two of them. We study the FIR and radio emission and quasar-host galaxy evolution with a sample of 18 z~6 quasars that are faint at UV/optical wavelengths (rest-frame 1450A magnitudes of m_1450\ge20.2). The average FIR-to-AGN UV luminosity ratio of this faint quasar sample is about two times higher than that of the bright quasars at z~6 (m_1450<20.2). A fit to the average FIR and AGN bolometric luminosities of both the UV/optically faint and bright z~6 quasars, and the average luminosities of samples of submillimeter /millimeter-observed quasars at z~2 to 5, yields a relationship of L_{FIR} {L_{bol}}^{0.62}. Five of the 18 faint z~6 quasars have been detected at 250 GHz. These 250 GHz detections, as well as most of the millimeter-detected optically bright z~6 quasars, follow a shallower trend of L_{FIR} {L_{bol}}^{0.45} defined by the starburst-AGN systems in local and high-z universe. The millimeter continuum detections in the five objects and molecular CO detections in three of them reveal a few x10^8 M_sun of FIR-emitting warm dust and 10^10 M_sun of molecular gas in the quasar host galaxies. All these results argue for massive star formation in the quasar host galaxies, with estimated star formation rates of a few hundred M_sun yr^{-1}. Additionally, the higher FIR-to-AGN luminosity ratio found in these 250 GHz-detected faint quasars also suggests a higher ratio between star formation rate and supermassive black hole accretion rate than the UV/optically most luminous quasars at z~6.

The identifications of quasars in the redshift range 2.2<z<3 are known to be very inefficient as their optical colors are indistinguishable from those of stars. Recent studies have proposed to use optical variability or near-IR colors to improve the identifications of the missing quasars in this redshift range. Here we present a case study by combining both factors. We select a sample of 70 quasar candidates from variables in SDSS Stripe 82, which are non-UV excess sources and have UKIDSS near-IR public data. They are clearly separated into two parts on the Y-K/g-z color-color diagram, and 59 of them meet or lie close to a newly proposed Y-K/g-z selection criterion for z<4 quasars. 44 of these 59 sources have been previously identified as quasars in SDSS DR7, and 35 among them are quasars at 2.2<z<3. We present spectroscopic observations of 14 of 15 remaining quasar candidates using the Bok 2.3m telescope and the MMT 6.5m telescope, and successfully identify all of them as new quasars at z=2.36 to 2.88. We also apply this method to a sample of 643 variable quasar candidates with SDSS-UKIDSS nine-band photometric data selected from 1875 new quasar candidates in SDSS Stripe 82 given by Butler & Bloom based on the time-series selections, and find that 188 of them are probably new quasars with photometric redshifts at 2.2<z<3. Our results indicate that the combination of optical variability and optical/near-IR colors is probably the most efficient way in finding 2.2<z<3 quasars and very helpful for constructing a complete quasar sample. We discuss its implications to the ongoing and upcoming large optical and near-IR sky surveys.

We report new observations of CO (2-1) line emission toward five z~6 quasars using the Ka-band receiver system on the Expanded Very Large Array (EVLA). Strong detections were obtained in two of them, SDSS J092721.82+200123.7 and CFHQS J142952.17+544717.6, and a marginal detection was obtained in another source, SDSS J084035.09+562419.9. Upper limits of the CO (2-1) line emission have been obtained for the other two objects. The CO (2-1) line detection in J0927+2001, together with previous measurements of the CO (6-5) and (5-4) lines, reveals important constraints on the CO excitation in the central ~10 kpc region of the quasar host galaxy. The CO (2-1) line emission from J1429+5447 is resolved into two distinct peaks separated by 1.2" (~6.9 kpc), indicating a possible gas-rich, major merging system, and the optical quasar position is consistent with the west peak. This result is in good agreement with the picture in which intense host galaxy star formation is coeval with rapid supermassive black hole accretion in the most distant universe. The two EVLA detections are ideal targets for further high-resolution imaging (e.g., with ALMA or EVLA observations) to study the gas distribution, dynamics, and SMBH-bulge mass relation in these earliest quasar-host galaxy systems.

We analytically investigate the condition for a particle accelerator to be active in the outer magnetosphere of a rotation-powered pulsar. Within the accelerator (or the gap), magnetic-field-aligned electric field accelerates electrons and positrons, which emit copious gamma-rays via curvature process. If one of the gamma-rays emitted by a single pair materializes as a new pair on average, the gap is self-sustained. However, if the neutron-star spin-down rate decreases below a certain limit, the gap becomes no longer self-sustained and the gamma-ray emission ceases. We explicitly compute the multiplicity of cascading pairs and find that the obtained limit corresponds to a modification of previously derived outer-gap death line. In addition to this traditional death line, we find another death line, which becomes important for millisecond pulsars, by separately considering the threshold of photon-photon pair production. Combining these traditional and new death lines, we give predictions on the detectability of gamma-ray pulsars with Fermi and AGILE. An implication on the X-ray observations of heated polar-cap emission is also discussed.

Line and continuum studies at centimeter through submillimeter wavelengths address probe deep into the earliest, most active and dust obscured phases of galaxy formation, and reveal the molecular and cool atomic gas. We summarize the techniques of radio astronomy to perform these studies, then review the progress on radio studies of galaxy formation. The dominant work over the last decade has focused on massive, luminous starburst galaxies (submm galaxies and AGN host galaxies). The far infrared luminosities are ~ 1e13 Lsun, implying star formation rates, SFR > 1e3 Msun/year. Molecular gas reservoirs are found with masses: M(H_2) > 1e10 (alpha/0.8}) Msun. The CO excitation in these luminous systems is much higher than in low redshift spiral galaxies. Imaging of the gas distribution and dynamics suggests strongly interacting and merging galaxies, indicating gravitationally induced, short duration (~ 1e7 year) starbursts. These systems correspond to a major star formation episode in massive galaxies in proto-clusters at intermediate to high redshift. Recently, radio observations have probed the more typical star forming galaxy population (SFR ~ 100 Msun/year), during the peak epoch of Universal star formation (z ~ 1.5 to 2.5). These observations reveal massive gas reservoirs without hyper-starbursts, and show that active star formation occurs over a wide range in galaxy stellar mass. The conditions in this gas are comparable to those found in the Milky Way disk. A key result is that the peak epoch of star formation in the Universe also corresponds to an epoch when the baryon content of star forming galaxies was dominated by molecular gas, not stars. We consider the possibility of tracing out the dense gas history of the Universe, and perform initial, admittedly gross, calculations. ABRIDGED

Estimating black hole masses of blazars is still a big challenge. Because of the contamination of jets, using the previously suggested size -- continuum luminosity relation can overestimate the broad line region (BLR) size and black hole mass for radio-loud AGNs, including blazars. We propose a new relation between the BLR size and $H_{\beta}$ emission line luminosity and present evidences for using it to get more accurate black hole masses of radio-loud AGNs. For extremely radio-loud AGNs such as blazars with weak/absent emission lines, we suggest to use the fundamental plane relation of their elliptical host galaxies to estimate the central velocity dispersions and black hole masses, if their velocity dispersions are not known but the host galaxies can be mapped. The black hole masses of some well-known blazars, such as OJ 287, AO 0235+164 and 3C 66B, are obtained using these two methods and the M - $\sigma$ relation. The implications of their black hole masses on other related studies are also discussed.

The present DAMA/LIBRA experiment and the former DAMA/NaI have cumulatively released so far the results obtained with the data collected over 13 annual cycles (total exposure: 1.17 ton $\times$ yr). They give a model independent evidence of the presence of DM particles in the galactic halo on the basis of the DM annual modulation signature at 8.9 $\sigma$ C.L. for the cumulative exposure.

We summarize cm through submm observations of the host galaxies of z ~ 6 quasars. These observations reveal the cool molecular gas (the fuel for star formation), the warm dust (heated by star formation), the fine structure line emission (tracing the CNM and PDRs), and the synchrotron emission. Our results imply active star formation in ~ 30% of the host galaxies, with star formation rates ~ 10^3 M_sun/year, and molecular gas masses ~ 10^10 M_sun. Imaging of the [CII] emission from the most distant quasar reveals a 'maximal starburst disk' on a scale ~ 1.5 kpc. Gas dynamical studies suggest a departure of these galaxies from the low-z M_{BH} -- M_{bulge} relation, with the black holes being, on average, 15 times more massive than expected. Overall, we are witnessing the co-eval formation of massive galaxies and supermassive black holes within 1 Gyr of the Big Bang.

We analyze the size evolution of HII regions around 27 quasars between z=5.7 to 6.4 ('quasar near-zones' or NZ). We include more sources than previous studies, and we use more accurate redshifts for the host galaxies, with 8 CO molecular line redshifts and 9 MgII redshifts. We confirm the trend for an increase in NZ size with decreasing redshift, with the luminosity normalized proper size evolving as: R_{NZ,corrected} = (7.4 \pm 0.3) - (8.0 \pm 1.1) \times (z-6) Mpc. While derivation of the absolute neutral fraction remains difficult with this technique, the evolution of the NZ sizes suggests a decrease in the neutral fraction of intergalactic hydrogen by a factor ~ 9.4 from z=6.4 to 5.7, in its simplest interpretation. Alternatively, recent numerical simulations suggest that this rapid increase in near-zone size from z=6.4 to 5.7 is due to the rapid increase in the background photo-ionization rate at the end of the percolation or overlap phase, when the average mean free path of ionizing photons increases dramatically. In either case, the results are consistent with the idea that z ~ 6 to 7 corresponds to the tail end of cosmic reionization. The scatter in the normalized NZ sizes is larger than expected simply from measurement errors, and likely reflects intrinsic differences in the quasars or their environments. We find that the near-zone sizes increase with quasar UV luminosity, as expected for photo-ionization dominated by quasar radiation.

We report our new observations of redshifted carbon monoxide emission from six z~6 quasars, using the PdBI. CO (6-5) or (5-4) line emission was detected in all six sources. Together with two other previous CO detections, these observations provide unique constraints on the molecular gas emission properties in these quasar systems close to the end of the cosmic reionization. Complementary results are also presented for low-J CO lines observed at the GBT and the VLA, and dust continuum from five of these sources with the SHARC-II bolometer camera at the CSO. We then present a study of the molecular gas properties in our combined sample of eight CO-detected quasars at z~6. The detections of high-order CO line emission in these objects indicates the presence of highly excited molecular gas, with estimated masses on the order of 10^10 M_sun within the quasar host galaxies. No significant difference is found in the gas mass and CO line width distributions between our z~6 quasars and samples of CO-detected $1.4\leq z\leq5$ quasars and submillimeter galaxies. Most of the CO-detected quasars at z~6 follow the far infrared-CO luminosity relationship defined by actively star-forming galaxies at low and high redshifts. This suggests that ongoing star formation in their hosts contributes significantly to the dust heating at FIR wavelengths. The result is consistent with the picture of galaxy formation co-eval with supermassive black hole (SMBH) accretion in the earliest quasar-host systems. We investigate the black hole--bulge relationships of our quasar sample, using the CO dynamics as a tracer for the dynamical mass of the quasar host. The results place important constraints on the formation and evolution of the most massive SMBH-spheroidal host systems at the highest redshift.

DAMA/LIBRA is running at the Gran Sasso National Laboratory of the I.N.F.N.. Here the results obtained with a further exposure of 0.34 ton x yr are presented. They refer to two further annual cycles collected one before and one after the first DAMA/LIBRA upgrade occurred on September/October 2008. The cumulative exposure with those previously released by the former DAMA/NaI and by DAMA/LIBRA is now 1.17 ton x yr, corresponding to 13 annual cycles. The data further confirm the model independent evidence of the presence of Dark Matter (DM) particles in the galactic halo on the basis of the DM annual modulation signature (8.9 sigma C.L. for the cumulative exposure). In particular, with the cumulative exposure the modulation amplitude of the single-hit events in the (2 -- 6) keV energy interval measured in NaI(Tl) target is (0.0116 +- 0.0013) cpd/kg/keV; the measured phase is (146 +- 7) days and the measured period is (0.999 +- 0.002) yr, values well in agreement with those expected for the DM particles.

The DAMA/LIBRA set-up (about 250 kg highly radiopure NaI(Tl) sensitive mass) is running at the Gran Sasso National Laboratory of the I.N.F.N.. The first DAMA/LIBRA results confirm the evidence for the presence of a Dark Matter particle component in the galactic halo, as pointed out by the former DAMA/NaI set-up; cumulatively the data support such evidence at 8.2 sigma C.L. and satisfy all the many peculiarities of the Dark Matter annual modulation signature. The main aspects and prospects of this model independent experimental approach will be outlined.

We report new continuum observations of fourteen z~6 quasars at 250 GHz and fourteen quasars at 1.4 GHz. We summarize all recent millimeter and radio observations of the sample of the thirty-three quasars known with 5.71<=z<=6.43, and present a study of the rest frame far-infrared (FIR) properties of this sample. These quasars were observed with the Max Plank Millimeter Bolometer Array (MAMBO) at 250 GHz with mJy sensitivity, and 30% of them were detected. We also recover the average 250 GHz flux density of the MAMBO undetected sources at 4 sigma, by stacking the on-source measurements. The derived mean radio-to-UV spectral energy distributions (SEDs) of the full sample and the 250 GHz non-detections show no significant difference from that of lower-redshift optical quasars. Obvious FIR excesses are seen in the individual SEDs of the strong 250 GHz detections, with FIR-to-radio emission ratios consistent with that of typical star forming galaxies. Most 250 GHz-detected sources follow the L_{FIR}--L_{bol} relationship derived from a sample of local IR luminous quasars (L_{IR}>10^{12}L_{\odot}), while the average L_{FIR}/L_{bol} ratio of the non-detections is consistent with that of the optically-selected PG quasars. The MAMBO detections also tend to have weaker Ly\alpha emission than the non-detected sources. We discuss possible FIR dust heating sources, and critically assess the possibility of active star formation in the host galaxies of the z~6 quasars. The average star formation rate of the MAMBO non-detections is likely to be less than a few hundred M_{\odot} yr^{-1}, but in the strong detections, the host galaxy star formation is probably at a rate of \gtrsim10^{3} M_{\odot} yr^{-1}, which dominates the FIR dust heating.

We use a recently released SDSS catalog of X-ray emitting AGNs in conjunction with the FIRST 20cm radio survey to investigate the black hole fundamental plane relationship between the 1.4GHz radio luminosity ($L_r$), 0.1-2.4 keV X-ray luminosity ($L_X$), and the black hole mass ($M$), namely, $log L_r=\xi_{RX}log L_X+\xi_{RM}log M +const$. For this purpose, we have compiled a large sample of 725 broad-line AGNs, which consists of 498 radio-loud sources and 227 radio-quiet sources. Our results are generally consistent with those in our previous work based on a smaller sample of 115 SDSS AGNs. We confirm that radio-loud objects have a steeper slope ($\xi_{RX}$) in the radio-X-ray relationship with respect to radio-quiet objects, and the dependence of the black hole fundamental plane on the black hole mass ($\xi_{RM}$) is weak. We also find a tight correlation with a similar slope between the soft X-ray luminosity and broad emission line luminosity for both radio-loud and radio-quiet AGNs, which implies that their soft X-ray emission is unbeamed and probably related to the accretion process. With the current larger sample of AGNs, we are able to study the redshift evolution of the black hole fundamental plane relation for both radio-loud and radio-quiet subsamples. We find that there is no clear evidence of evolution for radio-quiet AGNs, while for radio-loud ones there is a weak trend where $\xi_{RM}$ decreases as the redshift increases. This may be understood in part as due to the observed evolution of the radio spectral index as a function of redshift. Finally, we discuss the relativistic beaming effect and some other uncertainties related to the black hole fundamental plane (Abridged).

Some extensions of the Standard Model provide Dark Matter candidate particles with sub-GeV mass. These Light Dark Matter particles have been considered for example in Warm Dark Matter scenarios (e.g. the keV scale sterile neutrino, axino or gravitino). Moreover MeV scale DM candidates have been proposed in supersymmetric models and as source of the 511 keV line from the Galactic center. In this paper the possibility of direct detection of a Light Dark Matter candidate is investigated considering the inelastic scattering processes on the electron or on the nucleus targets. Some theoretical arguments are developed and related phenomenological aspects are discussed. Allowed volumes and regions for the characteristic phenomenological parameters of the considered scenarios are derived from the DAMA/NaI annual modulation data.

We present observations of four z>= SDSS quasars at 350 micron with the SHARC-II bolometer camera on the Caltech Submillimeter Observatory. These are among the deepest observations that have been made by SHARC-II at 350 micron, and three quasars are detected at >=3 sigma significance, greatly increasing the sample of 350 micron (corresponds to rest frame wavelengths of <60 micron at z>=5), detected high-redshift quasars. The derived rest frame far-infrared (FIR) emission in the three detected sources is about five to ten times stronger than that expected from the average SED of the local quasars given the same 1450A luminosity. Combining the previous submillimeter and millimeter observations at longer wavelengths, the temperatures of the FIR-emitting warm dust from the three quasar detections are estimated to be in the range of 39 to 52 K. Additionally, the FIR-to-radio SEDs of the three 350 micron detections are consistent with the emission from typical star forming galaxies. The FIR luminosities are ~10^{13} L_solar and the dust masses are >= 10^{8}M_solar. These results confirm that huge amounts of warm dust can exist in the host galaxies of optically bright quasars as early as z~6. The universe is so young at these epochs (~1 Gyr) that a rapid dust formation mechanism is required. We estimate the size of the FIR dust emission region to be about a few kpc, and further provide a comparison of the SEDs among different kinds of dust emitting sources to investigate the dominant dust heating mechanism.

Some extensions of the Standard Model provide Dark Matter candidate particles which can have a dominant coupling with the lepton sector of the ordinary matter. Thus, such Dark Matter candidate particles ($\chi^{0}$) can be directly detected only through their interaction with electrons in the detectors of a suitable experiment, while they are lost by experiments based on the rejection of the electromagnetic component of the experimental counting rate. These candidates can also offer a possible source of the 511 keV photons observed from the galactic bulge. In this paper this scenario is investigated. Some theoretical arguments are developed and related phenomenological aspects are discussed. Allowed intervals and regions for the characteristic phenomenological parameters of the considered model and of the possible mediator of the interaction are also derived considering the DAMA/NaI data.

The channeling effect of low energy ions along the crystallographic axes and planes of NaI(Tl) crystals is discussed in the framework of corollary investigations on WIMP Dark Matter candidates. In fact, the modeling of this existing effect implies a more complex evaluation of the luminosity yield for low energy recoiling Na and I ions. In the present paper related phenomenological arguments are developed and possible implications are discussed at some extent.

We have detected emission by the CO 5-4 and 6-5 rotational transitions at $z = 5.7722\pm 0.0006$ from the host galaxy of the SDSS quasar J0927+2001 using the Plateau de Bure interferometer. The peak line flux density for the CO 5-4 line is $0.72 \pm 0.09$ mJy, with a line FWHM = $610 \pm 110$ km s$^{-1}$. The implied molecular gas mass is $(1.6 \pm 0.3) \times 10^{10}$ M$_\odot$. We also detect the 90 GHz continuum at $0.12 \pm 0.03$ mJy, consistent with a 47K dust spectrum extrapolated from higher frequencies. J0927+2001 is the second example of a huge molecular gas reservoir within the host galaxy of a quasar within 1 Gyr of the big bang. Observations of J0927+2001 are consistent with a massive starburst coeval with a bright quasar phase in the galaxy, suggesting the rapid formation of both a super-massive black hole through accretion, and the stellar host spheroid, at a time close to the end of cosmic reionization.

We present millimeter and radio observations of 13 SDSS quasars at reshifts z~6. We observed eleven of them with the Max-Planck Millimeter Bolometer Array (MAMBO-2) at the IRAM 30m-telescope at 250 GHz and all of them with the Very Large Array (VLA) at 1.4 GHz. Four sources are detected by MAMBO-2 and six are detected by the VLA at >=3 sigma level. These sources, together with another 6 published in previous papers,yield a submillimeter/millimeter and radio observed SDSS quasar sample at z~6. We use this sample to investigate the far-infrared (FIR) andradio properties of optically bright quasars in the early universe. We compare this sample to lower redshift samples of quasars observed inthe submillimeter and millimeter wavelengths ((sub)mm), and find that the distribution of the FIR to B band optical luminosity ratio (L_FIR/L_B) is similar from z~2 to 6. We find a weak correlation between the FIR luminosity (L_FIR) and B band optical luminosity (L_B) byincluding the (sub)mm observed samples at all redshifts. Some strong (sub)mm detections in the z~6 sample have radio-to-FIR ratios within the range defined by star forming galaxies, which suggests possible co-eval star forming activity with the powerful AGN in these sources. We calculate the rest frame radio to optical ratios (R*_1.4=L_{v, 1.4GHz}/L_{v, 4400A}) for all of the VLA observed sources in the z~6 quasar sample. Only one radio detection in this sample, J083643.85+005453.3, has R*_1.4~40 and can be considered radio loud. There are no strong radio sources (R*_1.4>=100) among these SDSS quasars at z~6. These data are consistent with, although do not set strong constraints on, a decreasing radio-loud quasar fraction with increasing redshift.

We discuss observations of the first galaxies, within cosmic reionization, at centimeter and millimeter wavelengths. We present a summary of current observations of the host galaxies of the most distant QSOs ($z \sim 6$). These observations reveal the gas, dust, and star formation in the host galaxies on kpc-scales. These data imply an enriched ISM in the QSO host galaxies within 1 Gyr of the big bang, and are consistent with models of coeval supermassive black hole and spheroidal galaxy formation in major mergers at high redshift. Current instruments are limited to studying truly pathologic objects at these redshifts, meaning hyper-luminous infrared galaxies ($L_{FIR} \sim 10^{13}$ L$_\odot$). ALMA will provide the one to two orders of magnitude improvement in millimeter astronomy required to study normal star forming galaxies (ie. Ly-$\alpha$ emitters) at $z \sim 6$. ALMA will reveal, at sub-kpc spatial resolution, the thermal gas and dust -- the fundamental fuel for star formation -- in galaxies into cosmic reionization.

We have searched for HI 21cm absorption toward the two brightest radio AGN at high redshift, J0924--2201 at $z = 5.20$, and J0913+5919 at $z=5.11$, using the Giant Meter Wave Radio Telescope (GMRT). These data set a 3$\sigma$ upper limit to absorption of $< 30%$ at 40 km s$^{-1}$ resolution for the 30 mJy source J0913+5919, and $< 3%$ for the 0.55 Jy source J0924--2201 at 20 km s$^{-1}$ resolution. For J0924--2201, limits to broader lines at the few percent level are set by residual spectral baseline structure. For J0924--2201 the column density limit per 20 km s $^{-1}$ channel is: N(HI) $< 2.2\times 10^{18} \rm T_s$ cm$^{-2}$ over a velocity range of -700 km s$^{-1}$ to $+1180$ km s$^{-1}$ centered on the galaxy redshift determined through CO emission, assuming a covering factor of one. For J0913+5919 the column density limit per 40 km s$^{-1}$ channel is: N(HI) $< 2.2\times 10^{19} \rm T_s$ cm$^{-2}$ within $\pm 2400$ km s$^{-1}$ of the optical redshift. These data rule out any cool, high column density HI clouds within roughly $\pm 1000$ km s$^{-1}$ of the galaxies, as are often seen in Compact Steep Spectrum radio AGN, or clouds that might correspond to residual gas left over from cosmic reionization.

Because of the contamination of jets, using the size -- continuum luminosity relation can overestimate the broad line region (BLR) size and black hole mass for radio-loud AGNs. We propose a new relation between the BLR size and $H_{\beta}$ emission line luminosity and present evidences for using it to get more accurate black hole masses of radio-loud AGNs. For extremely radio-loud AGNs such as blazars with weak/absent emission lines, we suggest to use the fundamental plane relation of their elliptical host galaxies to estimate the central velocity dispersions and black hole masses, if their host galaxies can be mapped.

We present observations at 1.4 and 250 GHz of the $z\sim 5.7$ Ly$\alpha$ emitters (LAE) in the COSMOS field found by Murayama et al.. At 1.4 GHz there are 99 LAEs in the lower noise regions of the radio field. We do not detect any individual source down to 3$\sigma$ limits of $\sim 30\mu$Jy beam$^{-1}$ at 1.4 GHz, nor do we detect a source in a stacking analysis, to a 2$\sigma$ limit of $2.5\mu$Jy beam$^{-1}$. At 250 GHz we do not detect any of the 10 LAEs that are located within the central regions of the COSMOS field covered by MAMBO ($20' \times 20'$) to a typical 2$\sigma$ limit of $S_{250} < 2$mJy. The radio data imply that there are no low luminosity radio AGN with $L_{1.4} > 6\times 10^{24}$ W Hz$^{-1}$ in the LAE sample. The radio and millimeter observations also rule out any highly obscured, extreme starbursts in the sample, ie. any galaxies with massive star formation rates $> 1500$ M$_\odot$ year$^{-1}$ in the full sample (based on the radio data), or 500 M$_\odot$ year$^{-1}$ for the 10% of the LAE sample that fall in the central MAMBO field. The stacking analysis implies an upper limit to the mean massive star formation rate of $\sim 100$ M$_\odot$ year$^{-1}$.

The relationship between the emission-line equivalent width and the continuum luminosity, so called the Baldwin effect, plays an important role in studying the physics of the broad line region of AGNs. Using the archived ultraviolet spectra obtained by IUE, HST, and HUT in 1978-2002, we investigated the intrinsic CIV Baldwin effect of a well-studied Seyfert 1 galaxy NGC 4151. Both its continuum flux and CIV emission-line flux varied about two orders of magnitude in more than two decades, making it one of the best targets for studying the slope variation of the Baldwin effect. We fitted the CIV line profile of the 490 archived UV spectra of NGC 4151 with the same model consisting of a few Gaussian components, and derived the slope in the log-log plot for the total flux of CIV emission line against the UV continuum flux in different observation epochs.We found that the slope is not constant for NGC 4151, as it varies from 0.58 in the highest flux epoch to 0.83 in the lowest flux epoch. The slope evidently decreases as the continuum flux increases, which reinforces the previous findings of the non-constant slope in the H_beta Baldwin effect of NGC 5548 and the CIV Baldwin effect of Fairall 9. Our result suggests that such a non-constant slope may not be unusual for AGNs. Its physical origin is probably related to the different non-linear responses of the emission line to the variable ionizing continuum caused by the different accretion modes at different luminosity levels. We briefly discuss the effects of various absorption components in the CIV line profile of NGC 4151 and argue that the slope variation is not driven mainly by the absorption effect. Intensive, long-term and high resolution spectral observations of strongly variable AGNs are needed in the future to confirm our result.

The determination of the central black hole mass is crucial to the understanding of AGN physics. In this paper we briefly review some methods that are currently used to estimate the black hole mass of AGNs. Particularly we demonstrate the importance of two correlations: one between the black hole mass and the stellar velocity dispersion and another one between the broad line region (BLR) size and the optical continuum luminosity. Besides applying these relations in deriving black hole masses of various types of AGNs, we also employed the fundamental plane of elliptical galaxies to estimate the central velocity dispersions of AGN host galaxies, and then the black hole masses of AGNs including BL Lac objects. In addition, we derived another empirical relation between the BLR size and H_beta emission line luminosity from AGNs with the BLR size measured by reverberation mapping studies, and argued that more accurate black hole masses of extremely radio-loud AGNs could be obtained with it than using the usual R-L_{5100\AA} relation because of the jet contribution to the optical continuum. Finally we pointed out that black hole mass estimation is very much helpful to determine the accretion rate and understand the accretion process in AGNs.

We derived the black hole fundamental plane relationship among the 1.4GHz radio luminosity (L_r), 0.1-2.4keV X-ray luminosity (L_X), and black hole mass (M) from a uniform broad line SDSS AGN sample including both radio loud and radio quiet X-ray emitting sources. We found in our sample that the fundamental plane relation has a very weak dependence on the black hole mass, and a tight correlation also exists between the Eddington luminosity scaled X-ray and radio luminosities for the radio quiet subsample. Additionally, we noticed that the radio quiet and radio loud AGNs have different power-law slopes in the radio--X-ray non-linear relationship. The radio loud sample displays a slope of 1.39, which seems consistent with the jet dominated X-ray model. However, it may also be partly due to the relativistic beaming effect. For radio quiet sample the slope of the radio--X-ray relationship is about 0.85, which is possibly consistent with the theoretical prediction from the accretion flow dominated X-ray model. We briefly discuss the reason why our derived relationship is different from some previous works and expect the future spectral studies in radio and X-ray bands on individual sources in our sample to confirm our result.

Based on the measured sizes of broad line region of the reverberation-mapping AGN sample, two new empirical relations are introduced to estimate the central black hole masses of radio-loud high-redshift ($z > 0.5$) AGNs. First, using the archival $IUE/HST$ spectroscopy data at UV band for the reverberation-mapping objects, we obtained two new empirical relations between the BLR size and \Mg/\C emission line luminosity. Secondly, using the newly determined black hole masses of the reverberation-mapping sample for calibration, two new relationships for determination of black hole mass with the full width of half maximum and the luminosity of \Mg/\C line are also found. We then apply the relations to estimate the black hole masses of AGNs in Large Bright Quasar Surveyq and a sample of radio-loud quasars. For the objects with small radio-loudness, the black hole mass estimated using the $R_{\rm BLR} - L_{\eMg/\eC}$ relation is consistent with that from the $R_{BLR} - L_{3000\AA/1350 \AA}$ relation. But for radio-loud AGNs, the mass estimated from the $R_{BLR} - L_{\eMg/\eC}$ relation is systematically lower than that from the continuum luminosity $L_{3000\AA/1350\AA}$. Because jets could have significant contributions to the UV/optical continuum luminosity of radio-loud AGNs, we emphasized again that for radio-loud AGNs, the emission line luminosity may be a better tracer of the ionizing luminosity than the continuum luminosity, so that the relations between the BLR size and UV emission line luminosity should be used to estimate the black hole masses of high redshift radio-loud AGNs.

An empirical relation between the broad line region (BLR) size and optical continuum luminosity is often adopted to estimate the BLR size and then the black hole mass of AGNs. However, optical luminosity may not be a good indicator of photoionizing luminosity for extremely radio-loud AGNs because the jets usually contribute significantly to the optical continuum. Therefore, the black hole masses derived for blazar-type AGNs with this method are probably overestimated. Here we first derived a tight empirical relation between the BLR size and the H$_\beta$ emission line luminosity, $R(\rm{light-days})= 24.05(L_{H_\beta}/10^{42} ergs s^{-1})^{0.68}$, from a sample of 34 AGNs with the BLR size estimated with the reverberation mapping technique. Then we applied this relation to estimate the black hole masses of some AGNs and found that for many extremely radio-loud AGNs the black hole masses obtained with the $R-L_{H_\beta}$ relation are systematically lower than those derived previously with the $R-L_{5100\AA}$ relation, while for radio-quiet and slightly radio-loud AGNs the results obtained with these two methods are almost the same. The difference of black hole masses estimated with these two relations increases with the radio-loudness for extremely radio-loud AGNs, which is consistent with the fact that their equivalent widths of H$_{\beta}$ emission line become smaller at higher radio-loudness. If the small H$_{\beta}$ equivalent widths of extremely radio-loud AGNs are indeed caused by the beaming effect, we argue that the optical emission line luminosity may be a better tracer of ionizing luminosity for blazar-type AGNs and the black hole mass derived with the $R-L_{H_\beta}$ relation are probably more accurate.