We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

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.

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.

High-redshift blazars are valuable tools to study the early Universe. So far only a handful of $\gamma$-ray blazars have been found at redshifts above 3. $\gamma$-ray signals are detected in the direction of PMN J2219-2719 ($z=3.63$) and PMN J2321-0827 ($z=3.16$) by analyzing the 10-year $Fermi$-LAT Pass 8 data. PMN J2219-2719 is not distinguished from the background in the global analysis. During the 5-month epoch, the TS value is 47.8 and the flux is more than 10 times of the 10-year averaged flux. In addition, the angular distance between the $\gamma$-ray position and the radio position of PMN J2219-2719 is only ${0.04}^{\circ}$. Moreover, the $\gamma$-ray and infrared light curves of long time scale are very similar, which support the association between the $\gamma$-ray source and PMN J2219-2719. The global analysis of PMN J2321-0827 suggest a new $\gamma$-ray source, during the flare phase, the TS value is 61.4 and the $\gamma$-ray flux increased significantly. The association probability suggests that PMN J2321-0827 may be the counterpart of the new $\gamma$-ray source. In the future, the number of high-redshift $\gamma$-ray sources will increase by combining $Fermi$-LAT and the upcoming Large Synoptic Survey Telescope.

To understand the formation and quenching processes of local massive red spiral galaxies with $M_{\ast} > 10^{10.5}M_{\odot}$, we perform a statistical analysis of their spectroscopic and structural properties, and compare them with elliptical and blue spiral galaxies of similar mass. The sample was selected from the stellar mass catalog of galaxies in SDSS DR7, according to their locations on the u-r color-stellar mass diagram. We find that red spirals harbor compact cores with high stellar mass surface densities measured by $\Sigma_1$ and they are bulge-dominated. Particularly, the red spirals, especially their bulges follow the $\Sigma_1$-$M_{\ast}$ ridgeline for quenched galaxies. Furthermore, the red spirals show similarly large central D$_n(4000)$, high [Mg/Fe] and dark matter halo mass to ellipticals. These results suggest that the bulges of red spirals formed within a short timescale before redshift ~ 1-2 and were quenched via a fast mode, similar to ellipticals. Careful examinations of the optical morphologies reveal that ~70% of red spirals show strong bars, rings/shells and even merging features, which suggests that interactions or mergers might have played an important role in the formation of red spirals. In contrast, most of the massive blue spirals have completely different spectral and structural properties from red spirals. However, the blue spirals with high $\Sigma_1$ ($\Sigma_1 > 10^{9.5} M_\odot \, {\rm kpc}^{-2}$) show similar structural and morphological properties, as well as similar halo mass and HI mass to red spirals. We discuss rejuvenation from red to blue as a possible explanation for these high $\Sigma_1$ blue spirals.

The Chern-Simons term, through which the cosmic Axion-like field couples to the electromagnetic field, has the effect to rotate CMB polarization directions and to break the CPT symmetry. This rotation will change the CMB power spectra, no matter isotropic or anisotropic the rotation angle is. In this paper we revisit this issue by further considering the correlations between the (anisotropic) rotation angle $\alpha$ and the CMB temperature and (unrotated) $E$ polarization fields. These correlations could be generated in the Axion-like models with nonzero potential under the adiabatic initial condition. We first investigate how these correlations contribute further modifications to the CMB power spectra, then calculate the CMB bispectra for the temperature and rotated polarization fields. These bispectra would vanish if the $T\alpha$ and $E\alpha$ correlations are absent. So, they are useful in searching for CPT violation and the $T\alpha$ and $E\alpha$ correlations arisen in the Axion-like models.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

In this work, we present a new method to estimate cosmological parameters accurately based on Artificial Neural Network (ANN), and a code called ECoPANN (Estimating Cosmological Parameters with ANN) is developed to achieve parameter inference. We test the ANN method by estimating the basic parameters of the concordance cosmological model using the simulated temperature power spectrum of the cosmic microwave background (CMB). The results show that the ANN performs excellently on best-fit values and errors of parameters, as well as correlations between parameters when compared with that of the Markov chain Monte Carlo (MCMC) method. Besides, for a well trained ANN model, it is capable of estimating parameters for multiple experiments that have different precisions, which can greatly reduce the consumption of time and computing resources for parameter inference. Moreover, the well trained ANN is capable of discovering new potential physics in both the current and future higher precision observations. In addition, we extend the ANN to a multi-branch network to achieve joint constraint on parameters. We test the multi-branch network using the simulated temperature and polarization power spectra of CMB, type Ia supernovae, and baryon acoustic oscillation, and almost obtain the same results as the MCMC method. Therefore, we propose that the ANN can provide an alternative way to accurate and fast estimate cosmological parameters, and ECoPANN can be applied to the research of cosmology and even other broader scientific fields.

GB 1508+5714 is a high-redshift blazar (z = 4.3) and a spectrally soft gamma-ray source has been found in its direction. By analyzing 11.4-yr Fermi-LAT data, significant long-term variability of the gamma-ray source is confirmed. More importantly, a gamma-ray emission flare appeared in an epoch of several tens of days in year 2018, when the flux is about four times of the value from the global fit. Meanwhile, optical flares displayed in both i- and r-bands from the Zwicky Transient Facility light curves. Detections of the simultaneous gamma-ray and optical brightening provide a decisive proof to pin down the association between the gamma-ray source and GB 1508+5714, which makes it the first identified gamma-ray blazar beyond redshift 4. Future multiwavelength campaigns are needed to further investigate its jet properties, and, upcoming wide-deep-fast optical telescopes together with the gamma-ray sky surveyors will shed lights on the role that the AGN jets play in the early cosmic time.

In this work, we propose a new nonparametric approach for reconstructing a function from observational data using an Artificial Neural Network (ANN), which has no assumptions about the data and is a completely data-driven approach. We test the ANN method by reconstructing functions of the Hubble parameter measurements $H(z)$ and the distance-redshift relation $D_L(z)$ of Type Ia supernovae. We find that both $H(z)$ and $D_L(z)$ can be reconstructed with high accuracy. Furthermore, we estimate cosmological parameters using the reconstructed functions of $H(z)$ and $D_L(z)$ and find the results are consistent with those obtained using the observational data directly. Therefore, we propose that the function reconstructed by ANN can represent the actual distribution of observational data and can be used for parameter estimation in further cosmological research. In addition, we present a new strategy for training and evaluating the neural network, and a code for reconstructing functions using ANN has been developed and is available for download.

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

We present observations of HCN 3-2 emissions towards 37 local galaxies using 10-m Submillimeter Telescope (SMT). HCN 3-2 emission is detected in 23 galaxies. The correlation of infrared luminosity (LIR) and the luminosity of HCN 3-2 line emission measured in our sample is fitted with a slope of 1.11 and correlation coefficient of 0.91, which follows the linear correlation found in other dense gas tracers in the literatures. Although molecular gas above a certain volume density threshold (i.e., $n_{\rm H_2}\geq$ 10$^4$ cm$^{-3}$) statistically gave similar relation with infrared luminosity, the large scatter of HCN 3-2/HCN1-0 ratios for galaxies with different LIR indicates that dense gas masses estimated from the line luminosities of only one transition of dense gas tracers should be treated with caution for individual galaxies.

The Search for Extraterrestrial Intelligence (SETI) attempts to address the possibility of the presence of technological civilizations beyond the Earth. Benefiting from high sensitivity, large sky coverage, an innovative feed cabin for China's Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed the SETI first observations with FAST's newly commisioned 19-beam receiver; we report preliminary results in this paper. Using the data stream produced by the SERENDIP VI realtime multibeam SETI spectrometer installed at FAST, as well as its off-line data processing pipelines, we identify and remove four kinds of radio frequency interference(RFI): zone, broadband, multi-beam, and drifting, utilizing the Nebula SETI software pipeline combined with machine learning algorithms. After RFI mitigation, the Nebula pipeline identifies and ranks interesting narrow band candidate ET signals, scoring candidates by the number of times candidate signals have been seen at roughly the same sky position and same frequency, signal strength, proximity to a nearby star or object of interest, along with several other scoring criteria. We show four example candidates groups that demonstrate these RFI mitigation and candidate selection. This preliminary testing on FAST data helps to validate our SETI instrumentation techniques as well as our data processing pipeline.

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored.

We propose that the periodic fast radio bursts of FRB 180916.J0158+65 are sourced by axion emission (mass $m_{a} \sim 10^{-14}$ eV) from cosmic superstrings. Some of the emitted axions are converted to photons by magnetic fields as they travel along the line of sight to Earth. An impulsive burst of axion emission generates a photon signal typically lasting for milliseconds and varying with frequency in the observed manner. We find a range of parameters in our cosmic string network model consistent with the properties of FRB 180916.J0158+65. We suggest followup gravitational wave observations to test our model.

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $\sim 1.0\ {\rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($\sim 4.0$ mm) being a factor of $\sim10$ larger than other vortex regions. Both large ($\sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($\sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

The WorldWide Telescope(WWT) is a scientific visualization platform which can browse deep space images, star catalogs, and planetary remote sensing data from different observation facilities in a three-dimensional virtual scene. First launched and then open-sourced by Microsoft Research, the WWT is now managed by the American Astronomical Society (AAS). Hierarchical Progressive Survey (HiPS) is an astronomical data release scheme proposed by Centre de Donn\'ees astronomiques de Strasbourg (CDS) and has been accepted as a recommendation by International Virtual Observatory Alliance (IVOA). The HiPS solution has been adopted widely by many astronomical institutions for data release. Since WWT selected Hierarchical Triangular Mesh (HTM) as the standard for data visualization in the early stage of development, data released by HiPS cannot be visualized in WWT, which significantly limits the application of WWT. This paper introduces the implementation method for HiPS dataset visualization in WWT, and introduces HiPS data projection, mesh rendering, and data index implementation in WWT. Taking Chang'E-2 lunar probe data as an example, this paper introduces how to convert planetary remote sensing data into a HiPS dataset and integrate it into WWT. This paper also compares the efficiency and memory consumption of WWT loading its native data and HiPS data, and illustrates the application of HiPS in scientific data visualization and science education in WWT.

As part of our "Ultra-Fast Astronomy" (UFA) program, we present initial measurements of silicon photomultiplier detector (SiPM) performance for application to measurement of variability of optical astronomical targets on millisecond and shorter time scales. We mounted two different SiPM detectors at the focal plane of the 0.7 meter Nazarbayev University Transient Telescope at the Assy-Turgen Astrophysical Observatory (NUTTelA-TAO), with no filter in front of the detector, and measured their performance. The $3mm\times3mm$ SiPM detectors had a field of view of $2.2716'\times2.2716'$ for each SiPM channel. During the nights of 28-29/10/2019, we measured sky background, bright stars, and an artificial source with a 100Hz flashing frequency. We compared detected SiPM counts with Gaia satellite G band flux values to show that our SiPM detectors have linear response. With our two SiPMs (models S14520-3050VS and S14160-3050HS), we measured a dark current of $\sim$130 and $\sim$85 kilo count per second (kcps), and a sky background of $\sim$191 and $\sim$203 kcps, respectively. We measured an intrinsic crosstalk of 10.34$\%$ and 10.52$\%$ for the two detectors, respectively. We derive a 5$\sigma$ sensitivity of 14.7 and 15.2 Gaia G-band magnitude for 200ms exposures, for the two detectors, respectively. We derive, for transient events in 10$\mu s$ time window, a sensitivity of 20 detected photons or 6.1 Gaia G-band magnitudes. For shorter timescales , e.g. 100ns, our detection limit is limited by crosstalk to 12 detected photons, which correspond to fluence of $\sim$133 photons per square meter within 100ns.

Liu and colleagues recently proposed an elliptical accretion disk model for tidal disruption events (TDEs). They showed that the accretion disks of optical/UV TDEs are large and highly eccentric and suggested that the broad optical emission lines with complex and diverse profiles originate in the cool eccentric accretion disk of random inclination and orientation. In this paper, we calculate the radiation efficiency of the elliptical accretion disk and investigate the implications for the observations of TDEs. We compile the observational data of the peak bolometric luminosity and the total radiation energy after peak brightness of 18 TDE sources and compare the observations to the expectations from the elliptical accretion disk model. Our results show that the observations are well consistent with the expectations and that the majority of the orbital energy of the stellar debris is advected onto the black hole (BH) without being converted into radiation. By comparing the observations and the expectations, we derive the masses of the disrupted stars and BHs of the TDEs. The BH masses obtained in this paper are well consistent with those calculated with the M_BH-sigma_* relation. Our results provide an effective method to measure the masses of BHs from thousands or more of the TDEs to be discovered in the ongoing and next-generation sky surveys, no matter whether the BHs are located at the centers of the galactic nuclei or wandering in disks and halos.

We report the detection of an emission feature at the 12 sigma level with FWHM line width of about 450 km/s toward the nearest quasi-stellar object, QSO Mrk 231. Based on observations with the IRAM 30 m telescope and the NOEMA Interferometer, the 11-10 transition of molecular oxygen is the likely origin of line with rest frequency close to 118.75 GHz. The velocity of the O2 emission in Mrk 231 coincides with the red wing seen in CO emission, suggesting that it is associated with the outflowing molecular gas, located mainly at about ten kpc away from the central AGN. This first detection of extragalactic molecular oxygen provides an ideal tool to study AGN-driven molecular outflows on dynamic time scales of tens of Myr. O2 may be a significant coolant for molecular gas in such regions affected by AGN-driven outflows. New astrochemical models are needed to explain the implied high molecular oxygen abundance in such regions several kpc away from the center of galaxies.

High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations.

The accretion-powered X-ray pulsar GX 301-2 was observed with the balloon-borne X-Calibur hard X-ray polarimeter during late December 2018, with contiguous observations by the NICER X-ray telescope, the Swift X-ray Telescope and Burst Alert Telescope, and the Fermi Gamma-ray Burst Monitor spanning several months. The observations detected the pulsar in a rare apastron flaring state coinciding with a significant spin-up of the pulsar discovered with the Fermi GBM. The X-Calibur, NICER, and Swift observations reveal a pulse profile strongly dominated by one main peak, and the NICER and Swift data show strong variation of the profile from pulse to pulse. The X-Calibur observations constrain for the first time the linear polarization of the 15-35 keV emission from a highly magnetized accreting neutron star, indicating a polarization degree of (27+38-27)% (90% confidence limit) averaged over all pulse phases. We discuss the spin-up and the X-ray spectral and polarimetric results in the context of theoretical predictions. We conclude with a discussion of the scientific potential of future observations of highly magnetized neutron stars with the more sensitive follow-up mission XL-Calibur.

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$\beta$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$\beta$, He II, He I, and Fe II lines with respect to the varying 5100\AA continuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

Galaxy mergers and interactions are expected to play a significant role leading to offsets between gas and stellar motions in galaxies. Herein we crossmatch galaxies in MaNGA MPL-8 with the Dark Energy Spectroscopic Instrument (DESI) Legacy Surveys and identify 538 galaxies with merging/interacting features to investigate their position angle offsets ($\Delta$PA) between gas and stellar rotation. We find that there is a much higher merging/interacting fraction in misaligned galaxies (30$^\circ$ $\leqslant$ $\Delta$PA $<$150$^\circ$) than that in co-rotators ($\Delta$PA $<$ 30$^\circ$). This result corroborates that merging/interacting is one of the sources to produce misaligned galaxies and recent merging events can contribute to maximum 40\% of such misalignment. Furthermore, the marginal merging/interacting fraction in star-forming (SF) counter-rotators ($\Delta$PA $\geqslant$ 150$^\circ$) indicates the minor role of merging in the origin of SF counter-rotators. In addition, the slightly smaller merging/interacting fraction in non star-forming (non-SF) counter-rotators (0.14$\pm$0.039) than that in non-SF misaligned galaxies (0.22$\pm$0.031) agrees with that counter-rotating is a stable state, where tidal features disappear. Finally, the ratio of co-rotators to counter-rotators in non-SF merging/interacting galaxies is about 8:1, much larger than the prediction from the isotropic merging (1:1), which supports the speculation that gas and stars prefer to be aligned during merging, as the orbital angular momentum transfers to gas and stellar spin.

We present program objectives and specifications for the first generation Ultra-Fast Astronomy (UFA) observatory which will explore a new astrophysical phase space by characterizing the variability of the optical (320 nm - 650 nm) sky in the millisecond to nanosecond timescales. One of the first objectives of the UFA observatory will be to search for optical counterparts to fast radio bursts (FRB) that can be used to identify the origins of FRB and probe the epoch of reionization and baryonic matter in the interstellar and intergalactic mediums. The UFA camera will consist of two single-photon resolution fast-response detector 16x16 arrays operated in coincidence mounted on the 0.7 meter Nazarbayev University Transient Telescope at the Assy-Turgen Astrophysical Observatory (NUTTelA-TAO) located near Almaty, Kazakhstan. We are currently developing two readout systems that can measure down to the microsecond and nanosecond timescales and characterizing two silicon photomultipliers (SiPM) and one photomultiplier tube (PMT) to compare the detectors for the UFA observatory and astrophysical observations in general.

Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks which include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of sub-mm sized grains often appears in the dust emission at the outer edge of the dust disk.

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector ${\pi}_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this ${\pi}_{\rm E}$ determination with the angular Einstein radius $\theta_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $\mu_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

We performed high-sensitivity mapping observations of several complex organic molecules around Sagittarius B2 with ARO 12m telescope at 3-mm wavelength. Based on their spatial distribution, molecules can be classified as either "extended" that detected not only in Sgr B2(N) and Sgr B2(M), or "compact" that only detected toward or near to Sgr B2(N) and Sgr B2(M). The "extended" molecules including glycolaldehyde (CH2OHCHO), methyl formate (CH3OCHO), formic acid (t-HCOOH), ethanol (C2H5OH) and methyl amine (CH3NH2), while the "compact" molecules including dimethyl ether (CH3OCH3), ethyl cyanide (C2H5CN), and amino acetonitrile (H2NCH2CN). These "compact" molecules are likely produced under strong UV radiation, while "extended" molecules are likely formed under low-temperature, via gas-phase or grain surface reactions. The spatial distribution of "warm" CH2OHCHO at 89 GHz differ from the spatial distribution of "cold" CH2OHCHO observed at 13 GHz. We found evidence for an overabundance of CH2OHCHO compared to that expected from the gas-phase model, which indicates that grain-surface reactions are necessary to explain the origin of CH2OHCHO in Sagittarius B2. Grain-surface reactions are also needed to explain the correlation between the abundances of "cold" CH2OHCHO and C2H5OH. These results demonstrate the importance of grain-surface chemistry in the production of complex organic molecules.

This report summarizes the present status of neutrino non-standard interactions (NSI). After a brief overview, several aspects of NSIs are discussed, including connection to neutrino mass models, model-building and phenomenology of large NSI with both light and heavy mediators, NSI phenomenology in both short- and long-baseline neutrino oscillation experiments, neutrino cross-sections, complementarity of NSI with other low- and high-energy experiments, fits with neutrino oscillation and scattering data, DUNE sensitivity to NSI, effective field theory of NSI, as well as the relevance of NSI to dark matter and cosmology. We also discuss the open questions and interesting future directions that can be pursued by the community at large. This report is based on talks and discussions during the Neutrino Theory Network NSI workshop held at Washington University in St. Louis from May 29-31, 2019 (https://indico.cern.ch/event/812851/)

Given observations of the standard candles and the cosmic chronometers, we apply Pad\'{e} parameterization to the comoving distance and the Hubble paramter to find how stringent the constraint is set to the curvature parameter by the data. A weak informative prior is introduced in the modeling process to keep the inference away from the singularities. Bayesian evidence for different order of Pad\'{e} parameterizations is evaluated during the inference to select the most suitable parameterization in light of the data. The data we used prefer a parameterization form of comoving distance as $D_{01}(z)=\frac{a_0 z}{1+b_1 z}$ as well as a competitive form $D_{02}(z)=\frac{a_0 z}{1+b_1 z + b_2 z^2}$. Similar constraints on the spatial curvature parameter are established by those models and given the Hubble constant as a byproduct: $\Omega_k = 0.25^{+0.14}_{-0.13}$ (68\% confidence level [C.L.]), $H_0 = 67.7 \pm 2.0$ km/s/Mpc (68\% C.L.) for $D_{01}$, and $\Omega_k = -0.01 \pm 0.13$ (68\% C.L.), $H_0 = 68.8 \pm 2.0$ km/s/Mpc (68\% C.L.) for $D_{02}$. The evidence of different models demonstrates the qualitative analysis of the Pad\'{e} parameterizations for the comoving distance.

This work uses a combination of a variational auto-encoder and generative adversarial network to compare different dark energy models in light of observations, e.g., the distance modulus from type Ia supernovae. The network finds an analytical variational approximation to the true posterior of the latent parameters in the models, yielding consistent model comparison results with those derived by the standard Bayesian method, which suffers from a computationally expensive integral over the parameters in the product of the likelihood and the prior. The parallel computational nature of the network together with the stochastic gradient descent optimization technique leads to an efficient way to compare the physical models given a set of observations. The converged network also provides interpolation for a dataset, which is useful for data reconstruction.

In this work, a nearby luminous AGN sample is selected from HST, where only sources with both X-ray emission observed by \textit{Chandra/XMM-Newton} and radio flux detected by VLA/VLBA/VLBI/MERLIN are adopted to keep high precision. We get a sample of 30 luminous AGNs finally, which consist of 11 RLAGN and 19 RQAGN. It is found that the relationship between $R_{\rm {UV}}$ and $\alpha_{\rm {ox}}$, which was firstly reported by Li & Xie (2017) in LLAGN, and other relationships are all absent in RLAGN, probably due to the complex physical process therein. Our results indicate that the X-ray emission from jet should play an important role in RLAGN and further support the transition of accretion mode between LLAGN and RLAGN. On the other hand, the traditional relationships in RQAGN, such as $\alpha_{\rm {ox}}$ and $\lambda$, $\Gamma$ and $\lambda$, are found to be well consistent with previous works.

We investigate the impact of a highly eccentric 10 $M_{\rm \oplus}$ (where $M_{\rm \oplus}$ is the Earth mass) planet embedded in a dusty protoplanetary disk on the dust dynamics and its observational implications. By carrying out high-resolution 2D gas and dust two-fluid hydrodynamical simulations, we find that the planet's orbit can be circularized at large radii. After the planet's orbit is circularized, partial gap opening and dust ring formation happen close to the planet's circularization radius, which can explain the observed gaps/rings at the outer region of disks. When the disk mass and viscosity become low, we find that an eccentric planet can even open gaps and produce dust rings close to the pericenter and apocenter radii before its circularization. This offers alternative scenarios for explaining the observed dust rings and gaps in protoplanetary disks. A lower disk viscosity is favored to produce brighter rings in observations. An eccentric planet can also potentially slow down the dust radial drift in the outer region of the disk when the disk viscosity is low ($\alpha \lesssim2\times10^{-4}$) and the circularization is faster than the dust radial drift.

Cosmological CPT violation will rotate the polarized direction of CMB photons, convert partial CMB E mode into B mode and vice versa. It will generate non-zero EB, TB spectra and change the EE, BB, TE spectra. This phenomenon gives us a way to detect the CPT-violation signature from CMB observations, and also provides a new mechanism to produce B mode polarization. In this paper, we perform a global analysis on tensor-to-scalar ratio $r$ and polarization rotation angles based on current CMB datasets with both low $\ell$ (Planck, BICEP2/Keck Array) and high $\ell$ (POLARBEAR, SPTpol, ACTPol). Benefited from the high precision of CMB data, we obtain the isotropic rotation angle $\bar{\alpha} = -0.01^\circ \pm 0.37^\circ $ at 68% C.L., the variance of the anisotropic rotation angles $C^{\alpha}(0)<0.0032\,\mathrm{rad}^2$, the scale invariant power spectrum $D^{\alpha\alpha}_{\ell \in [2, 350]}<4.71\times 10^{-5} \,\mathrm{rad}^2$ and $r<0.057$ at 95% C.L.. Our result shows that with the polarization rotation effect, the 95% upper limit on $r$ gets tightened by 17%.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We show that the Deep Underground Neutrino Experiment (DUNE), with significant but feasible new efforts, has the potential to deliver world-leading results in solar neutrinos. With a 100~kton-year exposure, DUNE could detect $\gtrsim 10^5$ signal events above 5~MeV electron energy. Separate precision measurements of neutrino-mixing parameters and the $^8$B flux could be made using two detection channels ($\nu_e + \, ^{40}$Ar and $\nu_{e,\mu,\tau} + e^-$) and the day-night effect ($> 10 \sigma$). New particle physics may be revealed through the comparison of solar neutrinos (with matter effects) and reactor neutrinos (without), which is discrepant by $\sim 2 \sigma$ (and could become $5.6 \sigma$). New astrophysics may be revealed through the most precise measurement of the $^8$B flux (to 2.5\%) and the first detection of the {\it hep} flux (to 11\%). {\it DUNE is required:} No other experiment, even proposed, has been shown capable of fully realizing these discovery opportunities.

Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2-D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is obtained by including multiple dust fluids in a single grid-based hydrodynamic simulation and solving the Smoluchowski equation for dust coagulation on top of solving for the hydrodynamic evolution. We find that fragmentation of dust aggregates trapped in a pressure bump outside of the planetary gap leads to an enhancement in density of small grains. We compare the results obtained from the full coagulation treatment to the commonly used, fixed dust size approach and to previously applied, less computationally intensive methods for including dust coagulation. We find that the full coagulation results cannot be reproduced using the fixed-size treatment, but some can be mimicked using a relatively simple method for estimating the characteristic dust size in every grid cell.

In this work, we revisit the size-luminosity relation of the extended narrow line regions (ENLRs) using a large sample of nearby active galactic nuclei (AGN) from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The ENLRs ionized by the AGN are identified through the spatially resolved BPT diagram, which results in a sample of 152 AGN. By combining our AGN with the literature high-luminosity quasars, we found a tight log-linear relation between the size of the ENLR and the AGN [O III]{\lambda}5007{\AA} luminosity over four orders of magnitude of the [O III] luminosity. The slope of this relation is 0.42 $\pm$ 0.02 which can be explained in terms of a distribution of clouds photoionized by the AGN. This relation also indicates the AGN have the potential to ionize and heat the gas clouds at a large distance from the nuclei without the aids of outflows and jets for the low-luminosity Seyferts.

We present Submillimeter Array (SMA) observations of seven massive molecular clumps which are dark in the far-infrared for wavelengths up to 70 $\mu$m. Our 1.3 mm continuum images reveal 44 dense cores, with gas masses ranging from 1.4 to 77.1 M$_{\odot}$. Twenty-nine dense cores have masses greater than 8 M$_{\odot}$ and the other fifteen dense cores have masses between 1.4 and 7.5 M$_{\odot}$. Assuming the core density follows a power-law in radius $\rho \propto r^{-b}$, the index $b$ is found to be between 0.6 and 2.1 with a mean value of 1.3. The virial analysis reveals that the dense cores are not in virial equilibrium. CO outflow emission was detected toward 6 out of 7 molecular clumps and associated with 17 dense cores. For five of these cores, CO emissions appear to have line-wings at velocities of greater than 30 km s$^{-1}$ with respect to the source systemic velocity, which indicates that most of the clumps harbor protostars and thus are not quiescent in star formation. The estimated outflow timescale increase with core mass, which likely indicates that massive cores have longer accretion timescale than that of the less massive ones. The fragmentation analysis shows that the mass of low-mass and massive cores are roughly consistent with thermal and turbulent Jeans masses, respectively.

The hierarchical galaxy formation model predicts supermassive black hole binaries (SMBHBs) in galactic nuclei. Due to the gas poor environment and the limited spatial resolution in observations they may hide in the center of many a galaxy. However, a close encounter of a star with one of the supermassive black holes (SMBHs) may tidally disrupt it to produce a tidal disruption event (TDE) and temporarily light up the SMBH. In a previous work, we investigated with direct N-BODY simulations the evolution of TDE rates of SMBHB systems in galaxy mergers of equal mass. In this work we extend to unequal mass mergers. Our results show that, when two SMBHs are far away from each other, the TDE rate of each host galaxy is similar as in an isolated galaxy. As the two galaxies and their SMBHs separation shrinks, the TDE rate is increasing gradually until it reaches a maximum shortly after the two SMBHs become bound. In this stage, the averaged TDE rate can be enhanced by several times to an order of magnitude relative to isolated single galaxies. Our simulations show that the dependence of the TDE accretion rate on the mass ratio in this stage can be well fitted by power law relations for both SMBHs. After the bound SMBHB forms, the TDE rate decreases with its further evolution. We also find that in minor mergers TDEs of the secondary SMBH during and after the bound binary formation are mainly contributed by stars from the other galaxy.

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario.

We report two-dimensional spectroscopic analysis of massive red spiral galaxies ($M_{*}$ $>$ 10$^{10.5}$ $M_{\odot}$) and compare them to blue spiral and red elliptical galaxies above the same mass limit based on the public SDSS DR15 MaNGA observations. We find that the stellar population properties of red spiral galaxies are more similar to those of elliptical galaxies than to blue spiral galaxies. Red spiral galaxies show a shallow mass-weighted age profile, and they have higher stellar metallicity and Mgb/${\rm \langle Fe \rangle}$ across the whole 1.5$R_{\rm e}$ as compared to blue spirals, but all these properties are close to those of elliptical galaxies. One scenario to explain this is that red spirals form as remnants of very gas-rich major mergers that happened above $z$$\sim$1.

Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHa 330, are unlikely to be induced by the candidate vortices in them.

We characterized the S13360-3050CS Multi-Pixel Photon Counter (MPPC), a silicon photomultiplier (SiPM) manufactured by Hamamatsu Photonics K.K.. Measurements were obtained inside a light tight dark box using 365 nm, 400 nm, 525 nm, 660 nm, 720 nm, 810 nm, and 900 nm light-emitting diodes (LED) and the Citiroc 1A front-end evaluation system manufactured by Weeroc. At a 2.95V over voltage, we measured a dark count rate of 5.07x$10^{5}$ counts per second at 26$^{\circ}$C, crosstalk probability of 8.7$\%$, photon detection efficiency of 36$\%$ at 400 nm, linear range of 1.8x$10^{7}$ photons per second, and saturation at 5x$10^8$ photons per second. The S13360-3050CS MPPC is a candidate detector for the Ultra-Fast Astronomy (UFA) telescope which will characterize the optical sky in the millisecond to nanosecond timescales using two SiPM arrays operated in coincidence mounted on the 0.7 meter Nazarbayev University Transient Telescope at the Assy-Turgen Astrophysical Observatory (NUTTelA-TAO) located near Almaty, Kazakhstan. One objective of the UFA telescope will be to search for optical counterparts to fast radio bursts (FRB) that can be used to identify the origins of FRB and probe the epoch of reionization and baryonic matter in the interstellar and intergalactic mediums.

We present the photometric observations of blazars S5 0716+714 and 3C 273 with high temporal resolution (30--60s) in the $I$ or $R$ bands. The observations were performed with a 1.02 m optical telescope from 2007 March 07 to 2012 May 16. The $F$-test, one-way analysis of variance (ANOVA) test, and $z$-transformed discrete correlation function (ZDCF) cross-correlation zero lag test are used to search for intra-day variability (IDV). Four and five reliable IDVs survive three tests for S5 0716+714 and 3C 273, respectively. IDVs are found for S5 0716+714 and 3C 273. A flare on 2008 May 08 has $\Delta I \approx$ 0.06$\pm$0.01 mag in a duration of 0.54 hr for S5 0716+714. A flare on 2011 May 10 shows $\Delta R \approx$ 0.05$\pm$0.01 mag in a duration of 0.40 hr for 3C 273. Sharp dips appear on 2011 May 9 for 3C 273, and show $\Delta R \approx$ 0.05$\pm$0.01 mag. Under the assumptions that the IDV is tightly connected to black hole mass, $M_{\bullet}$, and that the flare durations are representative of the minimum characteristic timescales, we can estimate upper bounds to $M_{\bullet}$. In the case of the Kerr black holes, $M_{\bullet} \la 10^{8.91} M_{\odot}$ and $M_{\bullet} \la 10^{9.02} M_{\odot}$ are given for S5 0716+714 and 3C 273, respectively. These mass measurements are consistent with those measurements reported in the literatures. Also, we discuss the origins of optical variations found in this work.

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

With the accumulation of big data of CME observations by coronagraphs, automatic detection and tracking of CMEs has proven to be crucial. The excellent performance of convolutional neural network in image classification, object detection and other computer vision tasks motivates us to apply it to CME detection and tracking as well. We have developed a new tool for CME Automatic detection and tracking with MachinE Learning (CAMEL) techniques. The system is a three-module pipeline. It is first a supervised image classification problem. We solve it by training a neural network LeNet with training labels obtained from an existing CME catalog. Those images containing CME structures are flagged as CME images. Next, to identify the CME region in each CME-flagged image, we use deep descriptor transforming to localize the common object in an image set. A following step is to apply the graph cut technique to finely tune the detected CME region. To track the CME in an image sequence, the binary images with detected CME pixels are converted from cartesian to polar coordinate. A CME event is labeled if it can move in at least two frames and reach the edge of coronagraph field of view. For each event, a few fundamental parameters are derived. The results of four representative CMEs with various characteristics are presented and compared with those from four existing automatic and manual catalogs. We find that CAMEL can detect more complete and weaker structures, and has better performance to catch a CME as early as possible.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk surface density profile follows a power-law of the form $\rho_r \propto r^{-4}$. For the disk density profiles that can lead to a small reduced $\chi^2$ of goodness-of-fit, we find that there is a clear linear correlation between the disk mass and the power-law index $\gamma$ in the equation of disk density profile. This suggests that the $^{12}$CO disk of LkCa 15 is optically thick and we can fit its $^{12}$CO radial intensity profile using either a lower disk mass with a smaller $\gamma$ or a higher disk mass with a bigger $\gamma$. By comparing the $^{12}$CO channel maps of the best-fit model with disk models with higher or lower masses, we find that a disk mass of $\sim$ 0.1 $M_{\odot}$ can best reproduce the observed morphology of the $^{12}$CO channel maps. The dust continuum map at 0.87 mm of the LkCa 15 disk shows an inner cavity of the similar size of the best-fit gas model, but its out edge is at $\sim$ 200 AU, much smaller than the fitted gas disk. Such a discrepancy between the outer edges of the gas and dust disks is consistent with dust drifting and trapping models.

We present detector characterization of a state-of-the-art near-infrared (950nm - 1650 nm) Discrete Avalanche Photodiode detector (NIRDAPD) 5x5 array. We designed an experimental setup to characterize the NIRDAPD dark count rate, photon detection efficiency (PDE), and non-linearity. The NIRDAPD array was illuminated using a 1050 nm light-emitting diode (LED) as well as 980 nm, 1310 nm, and 1550 nm laser diodes. We find a dark count rate of 3.3x10$^6$ cps, saturation at 1.2x10$^8$ photons per second, a photon detection efficiency of 14.8% at 1050 nm, and pulse detection at 1 GHz. We characterized this NIRDAPD array for a future astrophysical program that will search for technosignatures and other fast (>1 Ghz) astrophysical transients as part of the Pulsed All-sky Near-infrared Optical Search for Extraterrestrial Intelligence (PANOSETI) project. The PANOSETI program will consist of an all-sky optical (350 - 800 nm) observatory capable of observing the entire northern hemisphere instantaneously and a wide-field NIR (950 - 1650 nm) component capable of drift scanning the entire sky in 230 clear nights. PANOSETI aims to be the first wide-field fast-time response near-infrared transient search.

In order to search for the line-like signals in the Fermi-LAT data, we have analyzed totally 49152 regions of interest (ROIs) that cover the whole sky. No ROI displays a line signal with test statistic (TS) value above 25, while for 50 ROIs weak line-like excesses with ${\rm TS}>16$ are presented. The intrinsic significances of these potential signals are further reduced by the large trial factor introduced in such kind of analysis. For the largest TS value of 24.3 derived in our analysis, the corresponding global significance is only $0.54\sigma$. We thus do not find any significant line-like signal and set up constraints on the cross section of dark matter annihilating to gamma-ray lines, $\left<\sigma{v}\right>_{\gamma\gamma}$.

We performed a survey in the SiO $J=5\rightarrow4$ line toward a sample of 199 Galactic massive star-forming regions at different evolutionary stages with the SMT 10 m and CSO 10.4 m telescopes. The sample consists of 44 infrared dark clouds (IRDCs), 86 protostellar candidates, and 69 young \HII\ regions. We detected SiO $J=5\rightarrow4$ line emission in 102 sources, with a detection rate of 57\%, 37\%, and 65\% for IRDCs, protostellar candidates, and young \HII\ regions, respectively. We find both broad line with Full Widths at Zero Power (FWZP) $>$ 20 \kms and narrow line emissons of SiO in objects at various evolutionary stages, likely associated with high-velocity shocks and low-velocity shocks, respectively. The SiO luminosities do not show apparent differences among various evolutionary stages in our sample. We find no correlation between the SiO abundance and the luminosity-to-mass ratio, indicating that the SiO abundance does not vary significantly in regions at different evolutionary stages of star formation.

The Deep Underground Neutrino Experiment (DUNE) could be revolutionary for MeV neutrino astrophysics, because of its huge detector volume, unique event reconstruction capabilities, and excellent sensitivity to the $\nu_e$ flavor. However, its backgrounds are not yet known. A major background is expected due to muon spallation of argon, which produces unstable isotopes that later beta decay. We present the first comprehensive study of MeV spallation backgrounds in argon, detailing isotope production mechanisms and decay properties, analyzing beta energy and time distributions, and proposing experimental cuts. We show that above a nominal detection threshold of 5-MeV electron energy, the most important backgrounds are --- surprisingly --- due to low-A isotopes, such as Li, Be, and B, even though high-A isotopes near argon are abundantly produced. We show that spallation backgrounds can be powerfully rejected by simple cuts, with clear paths for improvements. We compare these background rates to rates of possible MeV astrophysical neutrino signals in DUNE, including solar neutrinos (detailed in a companion paper [Capozzi et al. arXiv:1808.08232 [hep-ph]]), supernova burst neutrinos, and the diffuse supernova neutrino background. Further, to aid trigger strategies, in the Appendixes we quantify the rates of single and multiple MeV events due to spallation, radiogenic neutron capture, and other backgrounds, including through pileup. Our overall conclusion is that DUNE has high potential for MeV neutrino astrophysics, but reaching this potential requires new experimental initiatives.

The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.

We analyze the intrinsic velocity dispersion properties of 648 star-forming galaxies observed by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, to explore the relation of intrinsic gas velocity dispersions with star formation rates (SFRs), SFR surface densities ($\rm{\Sigma_{SFR}}$), stellar masses and stellar mass surface densities ($\rm{\Sigma_{*}}$). By combining with high z galaxies, we found that there is a good correlation between the velocity dispersion and the SFR as well as $\rm{\Sigma_{SFR}}$. But the correlation between the velocity dispersion and the stellar mass as well as $\rm{\Sigma_{*}}$ is moderate. By comparing our results with predictions of theoretical models, we found that the energy feedback from star formation processes alone and the gravitational instability alone can not fully explain simultaneously the observed velocity-dispersion/SFR and velocity-dispersion/$\rm{\Sigma_{SFR}}$ relationships.

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

The flavor composition of high-energy astrophysical neutrinos is a rich observable. However, present analyses cannot effectively distinguish particle showers induced by $\nu_e$ versus $\nu_\tau$. We show that this can be accomplished by measuring the intensities of the delayed, collective light emission from muon decays and neutron captures, which are, on average, greater for $\nu_\tau$ than for $\nu_e$. This new technique would significantly improve tests of the nature of astrophysical sources and of neutrino properties. We discuss the promising prospects for implementing it in IceCube and other detectors.

The formation of neutron stars (NSs), both from collapses of massive stars and mergers of compact objects, can be usually indicated by bright transients emitted from explosively-ejected material. In particular, if the newborn NSs can rotate at a millisecond period and have a sufficiently high magnetic field, then the spin-down of the NSs would provide a remarkable amount of energy to the emitting material. As a result, super-luminous supernovae could be produced in the massive stellar collapse cases, while some unusual fast evolving and luminous optical transients could arise from the cases of NS mergers and accretion-induced collapses of white dwarfs. In all cases, if the dipolar magnetic fields of the newborn NSs can be amplified to be as high as $10^{15}$ G, a relativistic jet could be launched and then a gamma-ray burst can be produced as the jet successfully breaks out from the surrounding nearly-isotropic ejected material.

Light, weakly-coupled dark sectors may be naturally decoupled in the early universe and enter equilibrium with the Standard Model bath during the epoch of primordial nucleosynthesis. The equilibration and eventual decoupling of dark sector states modifies the expansion rate of the universe, which alters the predicted abundances of the light elements. This effect can be encompassed in a time-varying contribution to $N_{\mathrm{eff}}$, the effective number of neutrino species, such that $N_{\mathrm{eff}}$ during nucleosynthesis differs from its measured value at the time of recombination. We investigate the impact of such variations on the light element abundances with model-independent templates for the time-dependence of $N_{\mathrm{eff}}$ as well as in specific models where a dark sector equilibrates with neutrinos or photons. We find that significant modifications of the expansion rate are consistent with the measured abundances of light nuclei, provided that they occur during specific periods of nucleosynthesis. In constraining concrete models, the relative importance of the cosmic microwave background and primordial nucleosynthesis is highly model-dependent.

We present the study of a sample of seven "bare" active galactic nuclei (AGN) observed with Suzaku. We interpret the spectrum of these sources with a relativistic reflection component and we employ our model RELXILL_NK to test the Kerr nature of their supermassive black holes. We constrain the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$, in which the Kerr metric is recovered when $\alpha_{13} = \alpha_{22} = 0$. All our measurements are consistent with the hypothesis that the spacetime geometry around these supermassive objects is described by the Kerr solution. For some sources, we obtain quite strong constraints on $\alpha_{13}$ and $\alpha_{22}$ when compared to those found in our previous studies. We discuss the systematic uncertainties in our tests and the implications of our results.

Retrieving the mass of a gravitational-wave (GW) source is a fundamental but difficult problem because the mass is degenerate with redshift. In astronomy, three types of redshift exist, namely cosmological, Doppler, and gravitational redshift, but the latter two are normally too weak to affect the observation. In this Letter, we show that the current astrophysical models allow binary black holes (BBHs) to merge within $10$ Schwarzschild radii of a supermassive black hole (SMBH). We find that in this case both the Doppler and gravitational redshift are significant, and in the most extreme condition they could increase the "apparent" black-hole mass and distance by a factor of $1.9-3.4$. We show that such a factor is consistent with the distribution in the distance-mass diagram of the ten BBHs detected so far by LIGO/Virgo. We also discuss the difficulties of this redshift scenario caused by the low event rate predicted by the current models, as well the potential solutions.

We use the spinor helicity formalism to study KLT-like relations for the inflationary graviton four-point correlation function. New features are observed in this correlation function compared to the graviton scattering amplitude in flat spacetime. After obtaining the general momentum dependence, collinear, squeezed and collapsed limits are considered to further study the features of the correlation function, and the relation to the corresponding flat space scattering amplitude.

We present results of millimeter line observations toward four local gas-rich galaxies (NGC 3079, NGC 4258, NGC 6240 and VII Zw 31) with the IRAM 30 meter millimeter telescope. More than 33 lines in these four sources were detected, including normal dense gas tracers (HCN 1-0, HCO$^+$ 1-0, and C$_2$H 1-0, etc) and their isotopic species. H$^{13}$CN (1-0) and H$^{13}$CO$^+$ (1-0) are detected for the first time in NGC 4258. Optical depths of HCN 1-0 and HCO$^{+}$ 1-0 were estimated with detected isotopic lines in NGC 4258, which were 4.1 and 2.6, respectively. HC$_3$N $J=29-28$, which requires high volume density and high temperature to excite, was detected in NGC 6240. High ratios of HCO$^+$/HCN in NGC 4258 and NGC 6240 imply that this ratio might not be a perfect diagnostic tool between AGN and starburst environments, due to contamination/combination of both processes. The low HC$_3$N/HCN line ratios with less than 0.15 in NGC 4258, NGC 6240 and the non-detection of HC$_3$N line in NGC 3079 and VII Zw 31 indicates that these four galaxies are HC$_3$N-poor galaxies. The variation of fractional abundance of CN in different types of galaxies is large.

In solar flares and other astrophysical systems, a major challenge for solving particle acceleration problem associated with magnetic reconnection is the enormous scale separation between kinetic scales and observed reconnection scale. Because of this, it has been difficult to draw any definite conclusions by just using kinetic simulations. Particle acceleration model that solves energetic particle transport equation can capture the main acceleration physics found in kinetic simulations, and thus provide a practical way to make observable predictions and directly compare model results with observations. Here we study compression particle acceleration in magnetic reconnection by solving Parker (diffusion-advection) transport equation using velocity and magnetic fields from two-dimensional high-Lundquist-number magnetohydrodynamics (MHD) simulations of a low-$\beta$ reconnection layer. We show that the compressible reconnection layer can give significant particle acceleration, leading to the formation of power-law particle energy distributions. We analyze the acceleration rate and find that the acceleration in the reconnection layer is a mixture of first order and second order Fermi processes. When including a guide field, we find the spectrum becomes steeper and both the power-law cutoff energy and maximum particle energy decrease as plasma becomes less compressible. This model produces 2D particle distribution that one can use to generate radiation map and directly compare with solar flare observations. This provides a framework to explain particle acceleration at large-scale astrophysical reconnection sites, such as solar flares.

Gravitational waves emitted from compact binary coalescence can be subject to wave diffraction if they are gravitationally lensed by an intervening mass clump whose Schwarzschild timescale matches the wave period. Waves in the ground-based frequency band $f\sim 10$--$10^3\,$Hz are sensitive to clumps with masses $M_E \sim 10^2$--$10^3\,M_\odot$ enclosed within the impact parameter. These can be the central parts of low mass $M_L \sim 10^3$--$10^6\,M_\odot$ dark matter halos, which are predicted in Cold Dark Matter scenarios but are challenging to observe. Neglecting finely-tuned impact parameters, we focus on lenses aligned generally on the Einstein scale for which multiple lensed images may not form in the case of an extended lens. In this case, diffraction induces amplitude and phase modulations whose sizes $\sim 10\%$--$20\%$ are small enough so that standard matched filtering with unlensed waveforms do not degrade, but are still detectable for events with high signal-to-noise ratio. We develop and test an agnostic detection method based on dynamic programming, which does not require a detailed model of the lensed waveforms. For pseudo-Jaffe lenses aligned up to the Einstein radius, we demonstrate that a pair of fully upgraded aLIGO/Virgo detectors can extract diffraction imprints from binary black hole mergers out to $z_s \sim 0.2$--$0.3$. The prospect will improve dramatically for a third-generation detector for which binary black hole mergers out to $z_s \sim 2$--$4$ will all become valuable sources.

Several nearby protoplanetary disks have been observed to display large scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas over-densities and are expected to concentrate dust particles with Stokes number around unity. This process might catalyze the formation of planetesimals. Whereas recent observations showed that dust crescent are indeed regions where millimeter-size particles have abnormally high concentration relative to the gas and smaller grains, no observations have yet shown that the gas within the crescent region counter-rotates with respect to the protoplanetary disk. Here we investigate the detectability of anticyclonic features through measurement of the line-of-sight component of the gas velocity obtained with ALMA. We carry out 2D hydrodynamic simulations and 3D radiative transfer calculation of a protoplanetary disk characterized by a vortex created by the tidal interaction with a massive planet. As a case study, the disk parameters are chosen to mimic the IRS 48 system, which has the most prominent crescent observed to date. We generate synthetic ALMA observations of both the dust continuum and 12CO emission around the frequency of 345 GHz. We find that the anticyclonic features of vortex are weak but can be detected if both the source and the observational setup are properly chosen. We provide a recipe for maximizing the probability to detect such vortex features and present an analysis procedure to infer their kinematic properties.

We report the detection of a transient $\gamma$-ray source in the direction of B3 1428+422 ($z=4.72$) by analyzing the 110-month {\it Fermi}-LAT Pass 8 data. The new transient $\gamma-$ray source is far away from the Galactic plane and has a rather soft spectrum, in agreement with being a high redshift blazar. We suggest that the newly discovered transient is the $\gamma$-ray counterpart of B3 1428+422, which could be the {\it most distant} GeV source detected so far. The detection of a group of such distant $\gamma-$ray blazars will be helpful to reconstruct the evolution of the luminosity function and to study the extragalactic background light at such high redshifts.

ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($\alpha\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened by one planet depend on input parameters, and apply our results to real systems. Gap locations (equivalently, spacings) are the easiest metric to use when making comparisons between theory and observations, as positions can be robustly measured. We fit the locations of gaps empirically as functions of planet mass and disk aspect ratio. We find that the locations of the double gaps in HL Tau and TW Hya, and of all three gaps in HD 163296, are consistent with being opened by a sub-Saturn mass planet. This scenario predicts the locations of other gaps in HL Tau and TW Hya, some of which appear consistent with current observations. We also show how the Rossby wave instability may develop at the edges of several gaps and result in multiple dusty vortices, all caused by one planet. A planet as low in mass as Mars may produce multiple dust gaps in the terrestrial planet forming region.

The radio-X-ray slope in the fundamental plane of radio-loud active galactic nuclei (AGNs) is found to be steeper compared with that of radio-quiet AGNs in previous works. In this work, we reinvestigate the fundamental plane in radio-loud AGNs by compiling a sample of 13 low-excitation radio galaxies (LERG) from the 3CR radio galaxies, for the reason that the accretion mode in LERG is believed to be a radiatively inefficient accretion flow. All the sources in our sample possess the data available at both the 5 GHz core radio luminosity detected by VLA/VLBI/VLBA and the core X-ray luminosity detected by Chandra/XMM-Newton. Surprisingly, we find the slope in the fundamental plane ($\log L_{\rm R}=0.52 \log L_{\rm X}+ 0.84 \log M_{\rm BH} + 10.84$) of LERG is well consistent with that reported by \citet{m2003}. However, the normalization is found to be shifted by about 0.7 dex, which can be due to the difference on magnetic field strength in different objects. A shallower slope of $L_{\rm R}-L_{\rm X}$ relation ($L_{\rm R}\sim L_{\rm X}^{0.63}$) is also given by our sample, which demonstrates that the X-ray emission in LERG may come from accretion disc instead of jets as suggested by previous works.

Early-type galaxies (ETGs) are composed of two distinct populations: high-mass and low-mass, which are likely to be built via gas-poor merging and gas-rich merging/accretion, respectively. However, it is difficult to directly associate low-mass ETGs with gas-rich processes, because currently they are gas poor with no signs of ongoing star formation. We report a discovery of an ETG (SDSS J142055.01+400715.7) with Mstellar=10^10 Msun that offers direct evidence for gas-rich merging as the origin of low-mass ETGs. The integrated properties of the galaxy are consistent with a typical low-mass ETG, but the outer and inner regions show distinct dispersion- and rotation-dominated kinematics, respectively. There are some tidal features surrounding the galaxy. These two facts suggest very recent galaxy merging. Furthermore, the inner disk harbors on-going star formation, indicating the merging to be gas rich. This type of galaxy is rare but it may be a demonstration of the role the transient phase of gas-rich merging plays in making a low-mass ETG.

The groundbreaking discovery of the optical transient AT2017gfo associated with GW170817 opens a unique opportunity to study the physics of double neutron star (NS) mergers. We argue that the standard interpretation of AT2017gfo as being powered by radioactive decays of r-process elements faces the challenge of simultaneously accounting for the peak luminosity and peak time of the event, as it is not easy to achieve the required high mass, and especially the low opacity of the ejecta required to fit the data. A plausible solution would be to invoke an additional energy source, which is probably provided by the merger product. We consider energy injection from two types of the merger products: (1) a post-merger black hole powered by fallback accretion; and (2) a long-lived NS remnant. The former case can only account for the early emission of AT2017gfo, with the late emission still powered by radioactive decay. In the latter case, both early- and late-emission components can be well interpreted as due to energy injection from a spinning-down NS, with the required mass and opacity of the ejecta components well consistent with known numerical simulation results. We suggest that there is a strong indication that the merger product of GW170817 is a long-lived (supramassive or even permanently stable), low magnetic field NS. The result provides a stringent constraint on the equations of state of NSs.

On 2017 September 22 the IceCube neutrino observatory detected a track-like, very-high-energy event (IceCube-170922A) that is spatially associated with TXS 0506+056, a quasar at a redshift of $z=0.3365$. This source is characterized by the increased acitivies in a very wide energy range (from radio to TeV) during these days. To investigate the possible connection of the PeV neutrino emission with the GeV activity of blazars, in this work we select 116 bright sources and analyze their lightcurves and spectra. We focus on the sources displaying GeV activities. Among these blazars, TXS 0506+056 seems to be typical in many aspects but is distinguished by the very strong GeV activties. We suggest to search for neutrino outburst in the historical data of IceCube, as recently done for TXS 0506+056, from the directions of these more energetic and harder blazars with strong GeV activities.

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

In this work, we search for $\gamma$-ray emission from the directions of some nearby Milky Way dwarf spheroidal galaxies (dSphs) and candidates with the publicly-available Pass 8 data of Fermi-LAT. Our sample includes 12 sources with the distances $<50$ kpc. Very weak $\gamma$-ray excesses ($\sim 2\sigma$) are found in some dSphs/candidates, consistent with that reported in the previous literature. Intriguingly, the peak test statistic (TS) value of the weak emission from the direction of Reticulum II rises continually. If interpreted as dark matter (DM) annihilation, the peak TS value is 13.5 for the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ and the DM mass of $m_\chi \sim 16$ GeV. The combination of all these nearby sources yields a more significant (with local significance $> 4\sigma$) $\gamma$-ray signal.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

In this paper, we study the strong gravitational lensing of gravitational waves (GWs) from a statistical perspective, with particular focus on the high frequency GWs from stellar binary black hole coalescences. These are most promising targets for ground-based detectors such as Advanced Laser Interferometer Gravitational Wave Observatory (aLIGO) and the proposed Einstein Telescope (ET) and can be safely treated under the geometrical optics limit for GW propagation. We perform a thorough calculation of the lensing rate, by taking account of effects caused by the ellipticity of lensing galaxies, lens environments, and magnification bias. We find that in certain GW source rate scenarios, we should be able to observe strongly lensed GW events once per year ($\sim1~\text{yr}^{-1}$) in the aLIGO survey at its design sensitivity; for the proposed ET survey, the rate could be as high as $\sim80~\text{yr}^{-1}$. These results depend on the estimate of GW source abundance, and hence can be correspondingly modified with an improvement in our understanding of the merger rate of stellar binary black holes. We also compute the fraction of four-image lens systems in each survey, predicting it to be $\sim30$ per cent for the aLIGO survey and $\sim6$ per cent for the ET survey. Finally, we evaluate the possibility of missing some images due to the finite survey duration, by presenting the probability distribution of lensing time delays. We predict that this selection bias will be insignificant in future GW surveys, as most of the lens systems ($\sim90$ per cent) will have time delays less than $\sim1$ month, which will be far shorter than survey durations.

We report simultaneous observations at different energy bands in radio, optical, UV, X-rays and $\gamma$ rays of the flat-spectrum radio-quasar 3C 345. We built the light curve of the source at different frequencies from 2008, the beginning of the \textit{Fermi} all-sky survey, to 2016, using new data and public archives. In particular we obtained several optical spectra, to study the behavior of emission lines and the continuum in different activity states and to derive the black hole mass. 3C 345 showed two flaring episodes in 2009, which occurred simultaneously in $\gamma$ ray, optical/UV and X-rays, and were later followed in radio. The source shows an inverse Compton dominated spectral energy distribution, which moved from higher to lower frequencies from the high to the low state. The reverberation of emission lines during one outburst event allowed us to constrain the location of production of $\gamma$ rays very close to the broad-line region, and possibly in the jet-base. We report the observation of an increased accretion after the outburst, possibly induced by the decrease of magnetic field intensity with respect to the low state.

We investigate two sympathetic filament eruptions observed by the New Vacuum Solar Telescope (NVST) on 2015 October 15. The full picture of the eruptions is obtained from the corresponding SDO/AIA observations. The two filaments start from active region NOAA 12434 in the north and end in one large quiescent filament channel in the south. The left filament erupts firstly, followed by the right filament eruption about 10 minutes later. Clear twist structure and rotating motion are observed in both filaments during the eruption. Both eruptions are failed, since the filaments firstly rise up, then flow towards the south and merge into the southern large quiescent filament. We also observe repeating activations of mini filaments below the right filament after its eruption. Using magnetic field models constructed based on SDO/HMI magnetograms by flux rope insertion method, we find that the left filament eruption is likely to be triggered by kink instability, while weakening of overlying magnetic fields due to magnetic reconnection at an X-point between the two filament systems might play an important role in the onset of the right filament eruption.

We study five sequential solar flares (SOL2015-08-07) occurring in Active Region 12396 observed with the Goode Solar Telescope (GST) at the BBSO, complemented by IRIS and SDO observations. The main flaring region is an arch filament system (AFS) consisting of multiple bundles of dark filament threads enclosed by semi-circular flare ribbons. We study the magnetic configuration and evolution of the active region by constructing coronal magnetic field models based on SDO/HMI magnetograms using two independent methods, i.e., the nonlinear force-free field (NLFFF) extrapolation and the flux rope insertion method. The models consist of multiple flux ropes with mixed signs of helicity, i.e., positive (negative) in the northern (southern) region, which is consistent with the GST observations of multiple filament bundles. The footprints of quasi-separatrix layers (QSLs) derived from the extrapolated NLFFF compare favorably with the observed flare ribbons. An interesting double-ribbon fine structure located at the east border of the AFS is consistent with the fine structure of the QSL's footprint. Moreover, magnetic field lines traced along the semi-circular footprint of a dome-like QSL surrounding the AFS are connected to the regions of significant helicity and Poynting flux injection. The maps of magnetic twist show that positive twist became dominant as time progressed, which is consistent with the injection of positive helicity before the flares. We hence conclude that these circular shaped flares are caused by 3D magnetic reconnection at the QSLs associated with the AFS possessing mixed signs of helicity.

Recent Atacama Large Millimeter and Submillimeter Array (ALMA) observations of the protoplanetary disk around the Herbig Ae star HD 163296 revealed three depleted dust gaps at 60, 100 and 160 au in the 1.3 mm continuum as well as CO depletion in the middle and outer dust gaps. However, no CO depletion was found in the inner dust gap. To examine the planet--disk interaction model, we present results of two-dimensional two fluid (gas + dust) hydrodynamic simulations coupled with three-dimensional radiative transfer simulations. In order to fit the high gas-to-dust ratio of the first gap, we find the Shakura--Sunyaev viscosity parameter $\alpha$ must be very small ($\lesssim 10^{-4}$) in the inner disk. On the other hand, a relatively large $\alpha$ ($\sim 7.5\times 10^{-3}$) is required to reproduce the dust surface density in the outer disk. We interpret the variation of $\alpha$ as an indicator of the transition from an inner dead zone to the outer magnetorotational instability (MRI) active zone. Within $\sim 100$ au, the HD 163296 disk's ionization level is low, and non-ideal magnetohydrodynamic (MHD) effects could suppress the MRI, so the disk can be largely laminar. The disk's ionization level gradually increases toward larger radii, and the outermost disk ($r > 300$ au) becomes turbulent due to MRI. Under this condition, we find that the observed dust continuum and CO gas line emissions can be reasonably fit by three half-Jovian-mass planets (0.46, 0.46 and 0.58 $M_\textrm{J}$) at 59, 105 and 160 au, respectively.

We focus on a series of $f(R)$ gravity theories in Palatini formalism to investigate the probabilities of producing the late-time acceleration for the flat Friedmann-Robertson-Walker (FRW) universe. We apply statefinder diagnostic to these cosmological models for chosen series of parameters to see if they distinguish from one another. The diagnostic involves the statefinder pair $\{r,s\}$, where $r$ is derived from the scale factor $a$ and its higher derivatives with respect to the cosmic time $t$, and $s$ is expressed by $r$ and the deceleration parameter $q$. In conclusion, we find that although two types of $f(R)$ theories: (i) $f(R) = R + \alpha R^m - \beta R^{-n}$ and (ii) $f(R) = R + \alpha \ln R - \beta$ can lead to late-time acceleration, their evolutionary trajectories in the $r-s$ and $r-q$ planes reveal different evolutionary properties, which certainly justify the merits of statefinder diagnostic. Additionally, we utilize the observational Hubble parameter data (OHD) to constrain these models of $f(R)$ gravity. As a result, except for $m=n=1/2$ of (i) case, $\alpha=0$ of (i) case and (ii) case allow $\Lambda$CDM model to exist in 1$\sigma$ confidence region. After adopting statefinder diagnostic to the best-fit models, we find that all the best-fit models are capable of going through deceleration/acceleration transition stage with late-time acceleration epoch, and all these models turn to de-Sitter point ($\{r,s\}=\{1,0\}$) in the future. Also, the evolutionary differences between these models are distinct, especially in $r-s$ plane, which makes the statefinder diagnostic more reliable in discriminating cosmological models.

As an essential part of China’s Gravitational Waves Program, the Ali CMB Polarization Telescope (AliCPT) is a ground-based experiment aiming at the Primordial Gravitational Waves (PGWs) by measuring B-mode polarization of Cosmic Microwave Background (CMB). First proposed in 2014 and currently in fast construction phase, AliCPT is China’s first CMB project that plans for commissioning in 2019. Led by the Institute of High Energy Physics (IHEP) under the Chinese Academy of Sciences (CAS), the project is a worldwide collaboration of more than fifteen universities and research institutes. Ali CMB Project is briefly introduced.

The axion-like particle (ALP)-photon mixing in the magnetic field around $\gamma$-ray sources or along the line-of-sight could induce oscillation between photons and ALPs, which then causes irregularities in the $\gamma$-ray spectra. In this work we try to search for such spectral irregularities in the spectrum of PKS $2155-304$ using 8.6 years of the Fermi Large Area Telescope (Fermi-LAT) data. No significant evidence for the presence of ALP-photon oscillation is obtained, and the parameter space of ALPs is constrained. The exclusion region sensitively depends on the poorly known magnetic field of host galaxy cluster of PKS $2155-304$. If the magnetic field is as high as $\sim 10~{\rm \mu G}$, the "hole"-like parameter region allowed in Ref.~\cite{Fermi} can be ruled out.

We study the rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital types, we can classify tidal disruption events (TDEs) into three main categories: eccentric, parabolic, and hyperbolic TDEs. Respective TDEs are characterized by two critical values of the orbital eccentricity: the lower critical eccentricity is the one below which the stars on eccentric orbits cause the finite, intense accretion, and the higher critical eccentricity above which the stars on hyperbolic orbits cause no accretion. Moreover, we find that the parabolic TDEs are divided into three subclasses: precisely parabolic, marginally eccentric, and marginally hyperbolic TDEs. We analytically derive that the mass fallback rate of the marginally eccentric TDEs can be flatter and slightly higher than the standard fallback rate proportional to $t^{-5/3}$, whereas it can be flatter and lower for the marginally hyperbolic TDEs. We confirm by N-body experiments that only few eccentric, precisely parabolic, and hyperbolic TDEs can occur in a spherical stellar system with a single intermediate-mass to supermassive black hole. A substantial fraction of the stars approaching to the black hole would cause the marginally eccentric or marginally hyperbolic TDEs.

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

With many large science equipment constructing and putting into use, astronomy has stepped into the big data era. The new method and infrastructure of big data processing has become a new requirement of many astronomers. Cloud computing, Map/Reduce, Hadoop, Spark, etc. many new technology has sprung up in recent years. Comparing to the high performance computing(HPC), Data is the center of these new technology. So, a new computing architecture infrastructure is necessary, which can be shared by both HPC and big data processing. Based on Astronomy Cloud project of Chinese Virtual Observatory (China-VO), we have made much efforts to optimize the designation of the hybrid computing platform. which include the hardware architecture, cluster management, Job and Resource scheduling.

Public observatory project is playing more and more important role in science popularization education and scientific research, and many amateur astronomers also have began to build their own observatories in remote areas. As a result of the limitation of technical condition and construction funds for amateur astronomers, their system often breaks down, and then a stable remote unattended operation system becomes very critical. Hardware connection and control is the basic and core part in observatory design. Here we propose a conception of engineering hardware design for public observatory operation as a bridge between observatory equipment and observation software. It can not only satisfy multiple observation mode requirement, but also save cost.

IAU Inter-Commission B2-C1-C2 WG Data-driven Astronomy Education and Public Outreach (DAEPO) was launched officially in April 2017. With the development of many mega-science astronomical projects, for example CTA, DESI, EUCLID, FAST, GAIA, JWST, LAMOST, LSST, SDSS, SKA, and large scale simulations, astronomy has become a Big Data science. Astronomical data is not only necessary resource for scientific research, but also very valuable resource for education and public outreach (EPO), especially in the era of Internet and Cloud Computing. IAU WG Data-driven Astronomy Education and Public Outreach is hosted at the IAU Division B (Facilities, Technologies and Data Science) Commission B2 (Data and Documentation), and organized jointly with Commission C1 (Astronomy Education and Development), Commission C2 (Communicating Astronomy with the Public), Office of Astronomy for Development (OAD), Office for Astronomy Outreach (OAO) and several other non IAU communities, including IVOA Education Interest Group, American Astronomical Society Worldwide Telescope Advisory Board, Zooniverse project and International Planetarium Society. The working group has the major objectives to: Act as a forum to discuss the value of astronomy data in EPO, the advantages and benefits of data driven EPO, and the challenges facing to data driven EPO; Provide guidelines, curriculum, data resources, tools, and e-infrastructure for data driven EPO; Provide best practices of data driven EPO. In the paper, backgrounds, current status and working plans in the future are introduced. More information about the WG is available at: http://daepo.china-vo.org/

Recently, gamma-ray emission in the direction of Coma, with a TS value of $\sim 40$, has been reported. In this work we will discuss the possibility of such a residual emission coming from dark matter annihilation. Our results show that the gamma-ray emission within the Coma region is also spatially correlated to the mass distribution derived from weak gravitational lensing measurements very well. However the dark matter models are not supported by the spectral analysis results and constraints by observations of other targets. Thus we derive the upper limits of the dark matter annihilation cross section according to the observation of the Coma region.

High-redshift blazars are one of the most powerful monsters in the universe and $\gamma$-ray variability carries crucial information of their relativistic jets. In this work we present results of the first systematical temporal analysis of {\it Fermi}-LAT data of all known seven $\gamma$-ray blazars beyond redshift 3. Significant long-term $\gamma$-ray variability is found from five sources in monthly $\gamma$-ray light curves, in which three of them are reported for the first time. Furthermore, intraday $\gamma$-ray variations are detected from NVSS J053954$-$283956 and NVSS J080518$+$614423. Doubling variability timescale of the former source is limited as short as $\lesssim$ 1 hour (at the source frame). Together with variability amplitude over one order of magnitude, NVSS J053954$-$283956 is the most distant $\gamma$-ray flaring blazar so far. Meanwhile, intraday optical variability of NVSS J163547$+$362930 is found based on archival PTF/iPTF light curve. Benefited from multiwavelength activity of these sources, constraints of their Doppler factors as well as locations of $\gamma$-ray radiation region and indications for the SDSS high redshift jetted active galactic nuclei deficit are discussed.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

SPIDER is a balloon-borne instrument designed to map the polarization of the millimeter-wave sky at large angular scales. SPIDER targets the B-mode signature of primordial gravitational waves in the cosmic microwave background (CMB), with a focus on mapping a large sky area with high fidelity at multiple frequencies. SPIDER's first longduration balloon (LDB) flight in January 2015 deployed a total of 2400 antenna-coupled Transition Edge Sensors (TESs) at 90 GHz and 150 GHz. In this work we review the design and in-flight performance of the SPIDER instrument, with a particular focus on the measured performance of the detectors and instrument in a space-like loading and radiation environment. SPIDER's second flight in December 2018 will incorporate payload upgrades and new receivers to map the sky at 285 GHz, providing valuable information for cleaning polarized dust emission from CMB maps.

The Cosmic Microwave Background (CMB) Polarization plays an important role in current cosmological studies. CMB B-mode polarization is the most effective probe to primordial gravitational waves (PGWs) and a test of the inflation as well as other theories of the early universe such as bouncing and cyclic universe. So far, major ground-based CMB polarization experiments are located in the southern hemisphere.Recently, China has launched the Ali CMB Polarization Telescope (AliCPT) in Tibetan Plateau to measure CMB B mode polarization and detect the PGWs in northern hemisphere. AliCPT include two stages, the first one is to build a telescope at the 5250m site (AliCPT-1) and the second one is to have a more sensitive telescope at a higher altitude of about 6000m (AliCPT-2). In this paper, we report the atmospherical conditions, sky coverage and the current infrastructure associated with AliCPT. We analyzed the reanalysis data from MERRA-2 together with radiosonde data from the Ali Meteorological Service and found that the amount of water vapor has a heavy seasonal variation and October to March is the suitable observation time. We also found 95/150 GHz to be feasible for AliCPT-1 and higher frequencies to be possible for AliCPT-2. Then we analyzed the observable sky and the target fields, and showed that Ali provides us a unique opportunity to observe CMB with less foreground contamination in the northern hemisphere and is complementary to the existed southern CMB experiments. Together with the developed infrastructure, we point out that Ali opens a new window for CMB observation and will be one of the major sites in the world along with Antarctic and Atacama.

We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne SPIDER instrument. These FPUs are vital to SPIDER's primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a 16 x 16 grid of corrugated silicon feedhorns coupled to an array of aluminum-manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1,530 polarization sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be ~3 pW at 300 mK with a less than 6% variation across each array at one standard deviation. These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the SPIDER-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We report new searches for the solar axions and galactic axion-like dark matter particles, using the first low-background data from PandaX-II experiment at China Jinping Underground Laboratory, corresponding to a total exposure of about $2.7\times 10^4$ kg$\cdot$day. No solar axion or galactic axion-like dark matter particle candidate has been identified. The upper limit on the axion-electron coupling ($g_{Ae}$) from the solar flux is found to be about $4.35 \times 10^{-12}$ in mass range from $10^{-5}$ to 1 keV/$c^2$ with 90\% confidence level, similar to the recent LUX result. We also report a new best limit from the $^{57}$Fe de-excitation. On the other hand, the upper limit from the galactic axions is on the order of $10^{-13}$ in the mass range from 1 keV/$c^2$ to 10 keV/$c^2$ with 90\% confidence level, slightly improved compared with the LUX.

We report the detection of widespread CH$_2$OHCHO and HOCH$_2$CH$_2$OH emission in Galactic center giant molecular cloud Sagittarius B2 using the Shanghai Tianma 65m Radio Telescope. Our observations show for the first time that the spatial distribution of these two important prebiotic molecules extends over 15 arc-minutes, corresponding to a linear size of approximately 36 pc. These two molecules are not just distributed in or near the hot cores. The abundance of these two molecules seems to decrease from the cold outer region to the central region associated with star-formation activity. Results present here suggest that these two molecules are likely to form through a low temperature process. Recent theoretical and experimental studies demonstrated that prebiotic molecules can be efficiently formed in icy grain mantles through several pathways. However, these complex ice features cannot be directly observed, and most constraints on the ice compositions come from millimeter observations of desorbed ice chemistry products. These results, combined with laboratory studies, strongly support the existence of abundant prebiotic molecules in ices.

We present our observations of the optical intra-day variability (IDV) in $\gamma$-ray BL Lac object Mrk 501. The observations were run with the 1.02 m and 2.4 m optical telescopes at Yunnan Observatories from 2005 April to 2012 May. The light curve at the $R$ band on 2010 May 15 passes both variability tests (the $F$ test and the ANOVA test). A flare within the light curve on 2010 May 15 has a magnitude change $\Delta m = 0.03 \pm 0.005_{\rm{stat}} \pm 0.007_{\rm{sys}}$ mag, \textbf{a darkening timescale of $\tau_{\rm{d}}=$ 26.7 minutes}, and an amplitude of IDV $Amp=2.9\% \pm0.7\%$. A decline \textbf{described by 11 consecutive flux measurements} within the flare can be fitted linearly with a Pearson's correlation coefficient $r = 0.945$ at the confidence level of $> 99.99\%$. Under the assumptions that the IDV is tightly connected to the mass of the black hole, \textbf{and that the flare duration, being two times $\tau_{\rm{d}}$, is representative of the minimum characteristic timescale, we can derive upper bounds to the mass of the black hole}. In the case of the Kerr black hole, the timescale of $\Delta t_{\rm{min}}^{\rm{ob}}=$ 0.89 hours gives $M_{\bullet}\la 10^{9.20} M_{\odot}$, which is consistent with measurements reported in the literature. This agreement indicates that the hypothesis about $M_{\bullet}$ and $\Delta t_{\rm{min}}^{\rm{ob}}$ is consistent with the measurements/data.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33<\ell<307$. No other limits exist over this full range of angular scales, and SPIDER improves upon the previous limit by several orders of magnitude, providing 95% C.L. constraints on $\ell (\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K ^2$ to 255 $\mu K ^2$ at 150 GHz for a thermal CMB spectrum. As linear CMB polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to obtain even stronger constraints on circular polarization.

By considering the realization of the emergent universe scenario in Eddington-inspired Born-Infeld (EiBI) theory, we study the stability of the Einstein static universe filled with perfect fluid in EiBI theory against both the homogeneous and inhomogeneous scalar perturbations in this work. We find that in both the spatially flat and closed cases, the emergent universe scenario is no longer viable, since the Einstein static universe cannot be stable against both the homogeneous and inhomogeneous scalar perturbations simultaneously. However, the emergent universe scenario survives in the spatially open case, while the Einstein static universe can be stable under some conditions.

It has been argued for years that the accretion mode changes from bright active galactic nuclei (AGNs) to low-luminosity AGNs (LLAGNs) at a rough dividing point of bolometric Eddington ratio $\lambda \sim 10^{-2}$. In this work, we strengthen this scenario through investigation of the relationship between the radio loudness $R_{\rm UV}$ and the optical-to-X-ray spectral index $\alpha_{\rm ox}$ in LLAGNs with $10^{-6} \lesssim \lambda \lesssim 10^{-3}$. We compile from literature a sample of 32 LLAGNs, consisting 18 LINERs and 14 low Eddington ratio Seyfert galaxies, and observe a strong negative $R_{\rm UV}$--$\alpha_{\rm ox}$ relationship, with large scatter in both $R_{\rm UV}$ and $\alpha_{\rm ox}$. We further demonstrate that this negative correlation, and the additional two negative relationships reported in literature ($R_{\rm UV}$--$\lambda$ and $\alpha_{\rm ox}$--$\lambda$ correlations), can be understood consistently and comprehensively under the truncated accretion--jet model, the model that has been applied successfully applied to LLAGNs. We argue that the scatter in the observations are (mainly) due to the spread in the viscosity parameter $\alpha$ of a hot accretion flow, a parameter that potentially can serve as a diagnose of the strength and/or configuration of magnetic fields in accretion flows.

We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal continuum for direct comparison with observations. In low viscosity gas ($\alpha\lesssim10^{-4}$), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on system parameters, additional rings may manifest for a single planet. A double gap located at tens of AUs from a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ~0".03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems like HL Tau (D5 and D6) and TW Hya (at 37 and 43 AU); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.

Significant undulations appear in the light curve of a recently discovered super-luminous supernova (SLSN) SN 2015bn after the first peak, while the underlying profile of the light curve can be well explained by a continuous energy supply from a central engine, possibly the spin-down of a millisecond magnetar. We propose that these undulations are caused by an intermittent pulsed energy supply, indicating an energetic flare activity of the central engine of the SLSN. Many post-burst flares were discovered during X-ray afterglow observations of Gamma-Ray Bursts (GRBs). We find that the SLSN flares described here approximately obey the empirical correlation between the luminosity and time scale of GRB flares, extrapolated to the relevant longer time scales of SLSN flares. This confirms the possible connection between these two different phenomena as previously suggested.

Milky Way-like galaxies are predicted to host a very large number of dark matter subhalos. Some massive and nearby subhalos could generate detectable gamma-rays, appearing as unidentified, spatially-extended and stable gamma-ray sources. We search for such sources in the third Fermi Large Area Telescope source List (3FGL) and report the identification of a new candidate, 3FGL J1924.8-1034. With the Fermi-LAT Pass 8 data, we find that 3FGL J1924.8-1034 is spatially-extended at a high confidence level of $5.4\sigma$, with a best-fit extension radius of $\sim0.15^{\circ}$. No significant variability has been found and its gamma-ray spectrum is well fitted by the dark matter annihilation into $b\bar{b}$ with a mass of $\sim 43$ GeV. All these facts make 3FGL J1924.8-1034 a possible dark matter subhalo candidate. However, due to the limited angular resolution, the possibility of that the spatial extension of 3FGL J1924.8-1034 is caused by the contamination from the other un-resolved point source can not be ruled out.

The ELUCID project aims to build a series of realistic cosmological simulations that reproduce the spatial and mass distribution of the galaxies as observed in the Sloan Digital Sky Survey (SDSS). This requires powerful reconstruction techniques to create constrained initial conditions. We test the reconstruction method by applying it to several $N$-body simulations. We use 2 medium resolution simulations from each of which three additional constrained $N$-body simulations were produced. We compare the resulting friend of friend catalogs by using the particle indexes as tracers, and quantify the quality of the reconstruction by varying the main smoothing parameter. The cross identification method we use proves to be efficient, and the results suggest that the most massive reconstructed halos are effectively traced from the same Lagrangian regions in the initial conditions. Preliminary time dependence analysis indicates that high mass end halos converge only at a redshift close to the reconstruction redshift. This suggests that, for earlier snapshots, only collections of progenitors may be effectively cross-identified.

We identify a largely model-independent signature of dark matter interactions with nucleons and electrons. Dark matter in the local galactic halo, gravitationally accelerated to over half the speed of light, scatters against and deposits kinetic energy into neutron stars, heating them to infrared blackbody temperatures. The resulting radiation could potentially be detected by the James Webb Space Telescope, the Thirty Meter Telescope, or the European Extremely Large Telescope. This mechanism also produces optical emission from neutron stars in the galactic bulge, and X-ray emission near the galactic center, because dark matter is denser in these regions. For GeV - PeV mass dark matter, dark kinetic heating would initially unmask any spin-independent or spin-dependent dark matter-nucleon cross-sections exceeding $2 \times 10^{-45}$ cm$^2$, with improved sensitivity after more telescope exposure. For lighter-than-GeV dark matter, cross-section sensitivity scales inversely with dark matter mass because of Pauli blocking; for heavier-than-PeV dark matter, it scales linearly with mass as a result of needing multiple scatters for capture. Future observations of dark sector-warmed neutron stars could determine whether dark matter annihilates in or only kinetically heats neutron stars. Because inelastic inter-state transitions of up to a few GeV would occur in relativistic scattering against nucleons, elusive inelastic dark matter like pure Higgsinos can also be discovered.

By fitting the bolometric light curves of 31 super-luminous supernovae (SLSNe) with the magnetar engine model, we derive the ejecta masses and magnetar parameters for these SLSNe. The lower boundary of magnetic field strengths of SLSN magnetars can be set just around the critical field strength $B_{\rm c}$ of electron Landau quantization. In more details, SLSN magnetars can further be divided into two subclasses of magnetic fields of $\sim(1-5)B_{\rm c}$ and $\sim(5-10)B_{\rm c}$, respectively. It is revealed that these two subclasses of magnetars are just associated with the slow-evolving and fast-evolving bolometric light curves of SLSNe. In comparison, the magnetars harbored in gamma-ray bursts (GRBs) and associated hypernovae are usually inferred to have much higher magnetic fields with a lower boundary about $\sim10B_{\rm c}$. This robustly suggests that it is the magnetic fields that play the crucial role in distinguishing SLSNe from GRBs/hypernovae. The rotational energy of SLSN magnetars are found to be correlated with the masses of supernova ejecta, which provides a clue to explore the nature of their progenitors. Moreover, the distribution of ejecta masses of SLSNe is basically intermediate between those of normal core-collapse supernovae and hypernovae. This could indicate an intrinsic connection among these different stellar explosions.

The distinctive morphology of head-tail radio galaxies reveals strong interactions between the radio jets and their intra-cluster environment, the general consensus on the morphology origin of head-tail sources is that radio jets are bent by violent intra-cluster weather. We demonstrate in this paper that such strong interactions provide a great opportunity to study the jet properties and also the dynamics of intra-cluster medium (ICM). By three-dimensional magnetohydrodynamical simulations, we analyse the detailed bending process of a magnetically dominated jet, based on the magnetic tower jet model. We use stratified atmospheres modulated by wind/shock to mimic the violent intra-cluster weather. Core sloshing is found to be inevitable during the wind-cluster core interaction, which induces significant shear motion and could finally drive ICM turbulence around the jet, making it difficult for jet to survive. We perform detailed comparison between the behaviour of pure hydrodynamical jets and magnetic tower jet, and find that the jet-lobe morphology could not survive against the violent disruption in all of our pure hydrodynamical jet models. On the other hand, the head-tail morphology is well reproduced by using a magnetic tower jet model bent by wind, in which hydrodynamical instabilities are naturally suppressed and the jet could always keep its integrity under the protection of its internal magnetic fields. Finally, we also check the possibility for jet bending by shock only. We find that shock could not bend jet significantly, so could not be expected to explain the observed long tails in head-tail radio galaxies.

Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.

New constraints are presented on the spin-dependent WIMP-nucleon interaction from the PandaX-II experiment, using a data set corresponding to a total exposure of 3.3$\times10^4$ kg-days. Assuming a standard axial-vector spin-dependent WIMP interaction with $^{129}$Xe and $^{131}$Xe nuclei, the most stringent upper limits on WIMP-neutron cross sections for WIMPs with masses above 10 GeV/c$^{2}$ are set in all dark matter direct detection experiments. The minimum upper limit of $4.1\times 10^{-41}$ cm$^2$ at 90\% confidence level is obtained for a WIMP mass of 40 GeV/c$^{2}$. This represents more than a factor of two improvement on the best available limits at this and higher masses. These improved cross-section limits provide more stringent constraints on the effective WIMP-proton and WIMP-neutron couplings.

We discuss an in-situ evaluation of the $^{85}$Kr, $^{222}$Rn, and $^{220}$Rn background in PandaX-I, a 120-kg liquid xenon dark matter direct detection experiment. Combining with a simulation, their contributions to the low energy electron-recoil background in the dark matter search region are obtained.

Virtual Observatory (VO) is a data-intensively online astronomical research and education environment, which takes advantages of advanced information technologies to achieve seamless and global access to astronomical information. AstroCloud is a cyber-infrastructure for astronomy research initiated by Chinese Virtual Observatory (China-VO) project, and also a kind of physical distributed platform which integrates lots of tasks such as telescope access proposal management, data archiving, data quality control, data release and open access, cloud based data processing and analysis. It consists of five application channels, i.e. observation, data, tools, cloud and public and is acting as a full lifecycle management system and gateway for astronomical data and telescopes. Physically, the platform is hosted in six cities currently, i.e. Beijing, Nanjing, Shanghai, Kunming, Lijiang and Urumqi, and serving more than 17 thousand users. Achievements from international Virtual Observatories and Cloud Computing are adopted heavily. In the paper, backgrounds of the project, architecture, Cloud Computing environment, key features of the system, current status and future plans are introduced.

X-ray polarimetry has seen a growing interest in recent years. Improvements in detector technology and focusing X-ray optics now enable sensitive astrophysical X-ray polarization measurements. These measurements will provide new insights into the processes at work in accreting black holes, the emission of X-rays from neutron stars and magnetars, and the structure of AGN jets. X-Calibur is a balloon-borne hard X-ray scattering polarimeter. An X-ray mirror with a focal length of 8 m focuses X-rays onto the detector, which consists of a plastic scattering element surrounded by Cadmium-Zinc-Telluride detectors, which absorb and record the scattered X-rays. Since X-rays preferentially scatter perpendicular to their polarization direction, the polarization properties of an X-ray beam can be inferred from the azimuthal distribution of scattered X-rays. A close alignment of the X-ray focal spot with the center of the detector is required in order to reduce systematic uncertainties and to maintain a high photon detection efficiency. This places stringent requirements on the mechanical and thermal stability of the telescope structure. During the flight on a stratospheric balloon, X-Calibur makes use of the Wallops Arc-Second Pointer (WASP) to point the telescope at astrophysical sources. In this paper, we describe the design, construction, and test of the telescope structure, as well as its performance during a 25-hour flight from Ft. Sumner, New Mexico. The carbon fiber-aluminum composite structure met the requirements set by X-Calibur and its design can easily be adapted for other types of experiments, such as X-ray imaging or spectroscopic telescopes.

Supermassive black hole binaries (SMBHBs) are productions of the hierarchical galaxy formation model. There are many close connections between central SMBH and its host galaxy because the former plays very important roles on the formation and evolution of a galaxy. For this reason, the evolution of SMBHBs in merging galaxies is an essential problem. Since there are many discussions about SMBHB evolution in gas rich environment, we focus on the quiescent galaxy, using tidal disruption as a diagnostic tool. Our study is based on a series of numerical large particle number direct N-body simulations for dry major mergers. According to the simulation results, the evolution can be divided into three phases. In phase I, the tidal disruption rate for two well separated SMBHs in merging system has similar level to single SMBH in isolate galaxy. After two SMBHs getting close enough to form a bound binary in phase II, the disruption rate can be enhanced for ~ 2 order of magnitudes within a short time. This "boosted" disruption stage finishes after the SMBHB evolving to compact binary system in phase III, corresponding to a drop back of disruption rate to a level of a few times higher than that in Phase I. How to correctly extrapolate our N-body simulation results to reality, and implications of our results to observations, are discussed too.

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after $10^4$ orbits, or $2.5$ Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

We present spectral line mapping observations toward four massive star-forming regions (Cepheus A, DR21S, S76E and G34.26+0.15), with the IRAM 30 meter telescope at 2 mm and 3 mm bands. Totally 396 spectral lines from 51 molecules, one helium recombination line, ten hydrogen recombination lines, and 16 unidentified lines were detected in these four sources. An emission line of nitrosyl cyanide (ONCN, 14$_{0,14}$-13$_{0,13}$) was detected in G34.26+0.15, as first detection in massive star-forming regions. We found that the $c$-C$_{3}$H$_{2}$ and NH$_{2}$D show enhancement in shocked regions as suggested by evidences of SiO and/or SO emission. Column density and rotational temperature of CH$_{3}$CN were estimated with the rotational diagram method for all four sources. Isotope abundance ratios of $^{12}$C/$^{13}$C were derived using HC$_{3}$N and its $^{13}$C isotopologue, which were around 40 in all four massive star-forming regions and slightly lower than the local interstellar value ($\sim$65). $^{14}$N/$^{15}$N and $^{16}$O/$^{18}$O abundance ratios in these sources were also derived using double isotopic method, which were slightly lower than that in local interstellar medium. Except for Cep A, $^{33}$S/$^{34}$S ratio in the other three targets were derived, which were similar to that in the local interstellar medium. The column density ratios of N(DCN)/N(HCN) and N(DCO$^{+}$)/N(HCO$^{+}$) in these sources were more than two orders of magnitude higher than the elemental [D]/[H] ratio, which is 1.5$\times$10$ ^{-5}$. Our results show the later stage sources, G34.26+0.15 in particular, present more molecular species than earlier stage ones. Evidence of shock activity is seen in all stages studied.

Observations of solar wind turbulence indicate the existence of multi-scale pressure-balanced structures (PBSs) in the solar wind. In this work, we conduct a numerical simulation to investigate multi-scale PBSs and in particular their formation in compressive MHD turbulence. By the use of a higher order Godunov code Athena,a driven compressible turbulence with an imposed uniform guide field is simulated. The simulation results show that both the magnetic pressure and the thermal pressure exhibit a turbulent spectrum with a Kolmogorov-like power law, and that in many regions of the simulation domain they are anti-correlated. The computed wavelet cross-coherence spectrum of the magnetic pressure and the thermal pressure, as well as their space series, indicate the existence of multi-scale PBSs, with the small PBSs being embedded in the large ones. These multi-scale PBSs are likely to be related with the highly oblique-propagating slow-mode waves, as the traced multi-scale PBS is found to be traveling in a certain direction at a speed consistent with that predicted theoretically for a slow-mode wave propagating in the same direction.

N-body simulations predict that galaxies at the Milky Way scale host a large number of dark matter (DM) subhalos. Some of these subhalos, if they are massive enough or close enough to the Earth, might be detectable in $\gamma$ rays due to the DM annihilation. 3FGL J2212.5+0703, an unidentified gamma-ray source, has been suggested to be the counterpart candidate of a DM subhalo by Bertoni et al. (2015, 2016). In this work we analyze the Fermi-LAT Pass 8 data of 3FGL J2212.5+0703 to independently test the DM subhalo hypothesis of this source. In order to suppress the possible contamination from two nearby very-bright blazars, we just take into account the front-converting gamma-rays which have better angular resolutions than that of the back-converting photons. In addition to the spatial distribution analysis, we have extended the spectrum analysis down to the energies of $\sim 100$ MeV, and thoroughly examined the variability of the emission during the past 8 years. We confirm that 3FGL J2212.5+0703 is a steady and spatially-extended gamma-ray emitter at a high confidence level. The spectrum is well consistent with that expected from DM annihilation into $b\bar{b}$. The introduction of a phenomenological LogParabola spectrum just improves the fit slightly. All these results suggest that 3FGL J2212.5+0703 could be indicative of a DM subhalo.

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

To explain the accelerated expansion of our universe, many dark energy models and modified gravity theories have been proposed so far. It is argued in the literature that they are difficult to be distinguished on the cosmological scales. Therefore, it is well motivated to consider the relevant astrophysical phenomena on (or below) the galactic scales. In this work, we study the stability of self-gravitating differentially rotating galactic disks in $f(T)$ theory, and obtain the local stability criteria in $f(T)$ theory, which are valid for all $f(T)$ theories satisfying $f(T=0)=0$ and $f_T (T=0)\not=0$, if the adiabatic approximation and the weak field limit are considered. The information of the function $f(T)$ is mainly encoded in the parameter $\alpha\equiv 1/f_T(T=0)$. We find that the local stability criteria in $f(T)$ theory are quite different from the ones in Newtonian gravity, general relativity, and other modified gravity theories such as $f(R)$ theory. We consider that this might be a possible hint to distinguish $f(T)$ theory from general relativity and other modified gravity theories on (or below) the galactic scales.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

III Zw 2 is the prototype of radio-intermediate quasars. Although there is the evidence of possessing strong jet, significant $\gamma$-ray emission has not been reported before. In this work, we carry out a detailed analysis of the latest {\it Fermi}-LAT {\it Pass} 8 data. No significant $\gamma$-ray signal has been detected in the time-averaged 7-year {\it Fermi}-LAT data of III Zw 2, we however have identified two distinct $\gamma$-ray flares with isotropic luminosities of $\sim 10^{45}$ erg $\rm s^{-1}$. Multiwavelength data analysis (including also the optical photometric observations from Yunnan Observatories) are presented and the main finding is the simultaneous optical and $\gamma$-ray flares of III Zw 2 appearing in Nov. 2009. Violent $\gamma$-ray variability with doubling timescale of 2.5 hours was detected in another $\gamma$-ray flare in May 2010, for which the 3-hour $\gamma$-ray peak flux is $\sim 250$ times of the average flux in 7 years. Rather similar behaviors are observed in blazars and the blazar model can reasonably reproduce the spectral energy distribution of III Zw 2 in a wide energy range, strongly suggesting that its central engine resembles that of blazars. In view of its core which shares radio similarities with young radio sources together with the weak extended radio lobe emission, we suggest that III Zw 2 harbors a recurrent activity core and serves as a valuable target for investigating the fuelling and triggering of the activity in radio loud AGNs.

We investigate two dark energy cosmological models (i.e., the $\Lambda$CDM and $\phi$CDM models) with massive neutrinos assuming two different neutrino mass hierarchies in both the spatially flat and non-flat scenarios, where in the $\phi$CDM model the scalar field possesses an inverse power-law potential, $V(\phi)\propto {\phi}^{-\alpha}$ ($\alpha>0$). Cosmic microwave background data from Planck 2015, baryon acoustic oscillations data from 6dFGS, SDSS-MGS, BOSS-LOWZ and BOSS CMASS-DR11, the JLA compilation of Type Ia supernova apparent magnitude observations, and the Hubble Space Telescope $H_0$ prior, are jointly employed to constrain the model parameters. We first determine constraints assuming three species of degenerate massive neutrinos. In the spatially flat (non-flat) $\Lambda$CDM model, the sum of neutrino masses is bounded as $\Sigma m_{\nu} < 0.165 (0.299)$ eV at 95% confidence level (CL). Correspondingly, in the flat (non-flat) $\phi$CDM model, we find $\Sigma m_{\nu} < 0.164 (0.301)$ eV at 95% CL. The inclusion of spatial curvature as a free parameter results in a significant broadening of confidence regions for $\Sigma m_{\nu}$ and other parameters. In the scenario where the total neutrino mass is dominated by the heaviest neutrino mass eigenstate, we can obtain the similar conclusions as those obtained in the degenerate neutrino mass scenario. In addition, the results show that the bounds on $\Sigma m_{\nu}$ based on two different neutrino mass hierarchies have insignificant differences in the spatially flat case for both the $\Lambda$CDM and $\phi$CDM models, however, the corresponding differences are larger in the non-flat case.

The optical/ultraviolet (UV) variability of quasars has been discovered to be correlated with other quasar properties, such as luminosity, black hole mass and rest-frame wavelength. However,the origin of variability has been a puzzle so far. In this work, we upgrade the accretion disc model (Li & Cao 2008), which assumed the variability is caused by the change of global mass accretion rate, by constraining the disc size to match the viscous timescale of accretion disc to the variability timescale observed and by including the irradiation/X-ray reprocessing to make the emitted spectrum become steeper. We find this hybrid model can reproduce the observed bluer-when-brighter trend quite well, which is used to validate the theoretical model by several works recently. The traditional correlation between the variability amplitude and rest-frame wavelength can also be well fitted by our model. In addition, a weak positive correlation between variability amplitude and black hole mass is present, qualitatively consistent with recent observations.

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances.

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

Very recently the Dark Energy Survey (DES) Collaboration has released their second group of Dwarf spheroidal (dSph) galaxy candidates. With the publicly-available Pass 8 data of Fermi-LAT we search for $\gamma-$ray emissions from the directions of these eight newly discovered dSph galaxy candidates. No statistically significant $\gamma-$ray signal has been found in the combined analysis of these sources. With the empirically estimated J-factors of these sources, the constraint on the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ is comparable to the bound set by the joint analysis of fifteen previously known dSphs with kinematically constrained J-factors for the dark matter mass $m_\chi>250$ GeV. In the direction of Tucana III (DES J2356-5935), one of the nearest dSph galaxy candidates that is $\sim 25$ kpc away, there is a weak $\gamma-$ray signal and its peak test statistic (TS) value for the dark matter annihilation channel $\chi\chi\rightarrow \tau^{+}\tau^{-1}$ is $\approx 6.7$ at $m_\chi \sim 15$ GeV. The significance of the possible signal likely increases with time. More data is highly needed to pin down the physical origin of such a GeV excess.

A supra-massive neutron star (NS) spinning extremely rapidly could survive from a merger of NS-NS binary. The spin-down of this remnant NS that is highly magnetized would power the isotropic merger ejecta to produce a bright mergernova emission in ultraviolet/optical bands. Before the mergernova, the early interaction between the NS wind and the ejecta can drive a forward shock propagating outwards into the ejecta. As a result, a remarkable amount of heat can be accumulated behind the shock front and the final escaping of this heat can produce a shock breakout emission. We describe the dynamics and thermal emission of this shock with a semi-analytical model. It is found that sharp and luminous breakout emission appears mainly in soft X-rays with a luminosity of $\sim10^{45}~\rm erg~s^{-1}$ at a few hours after the merger, by leading the mergernova emission as a precursor. Therefore, detection of such an X-ray precursor would provide a smoking-gun evidence for identifying NS-powered mergernovae and distinguishing them from the radioactive-powered ones (i.e., kilonovae or macronovae). The discovery of NS-powered mergernovae would finally help to confirm the gravitational wave signals due to the mergers and the existence of supra-massive NSs.

As the cyber-infrastructure for Astronomical research from Chinese Virtual Observatory (China-VO) project, AstroCloud has been archived solid progresses during the last one year. Proposal management system and data access system are re-designed. Several new sub-systems are developed, including China-VO PaperData, AstroCloud Statics and Public channel. More data sets and application environments are integrated into the platform. LAMOST DR1, the largest astronomical spectrum archive was released to the public using the platform. The latest progresses will be introduced.

The Large sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) is the largest optical telescope in China. In last four years, the LAMOST telescope has published four editions data (pilot data release, data release 1, data release 2 and data release 3). To archive and release these data (raw data, catalog, spectrum etc), we have set up a data cycle management system, including the transfer of data, archiving, backup. And through the evolution of four software versions, mature established data release system.

We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices of disk parameters, assuming an already settled dust distribution and no planet migration. We can obtain quite adequate fits to the observed dust emission using three planets with masses 0.35, 0.17, and 0.26 $M_{Jup}$ at 13.1, 33.0, and 68.6 AU, respectively. Implications for the planet formation as well as the limitations of this scenario are discussed.

The study of reionization history plays an important role in understanding the evolution of our universe. It is commonly believed that the intergalactic medium (IGM) in our universe are fully ionized today, however the reionizing process remains to be mysterious. A simple instantaneous reionization process is usually adopted in modern cosmology without direct observational evidence. However, the history of ionization fraction, $x_e(z)$ will influence cosmic microwave background (CMB) observables and constraints on optical depth $\tau$. With the mocked future data sets based on featured reionization model, we find the bias on $\tau$ introduced by instantaneous model can not be neglected. In this paper, we study the cosmic reionization history in a model independent way, the so called principle component analysis (PCA) method, and reconstruct $x_e (z)$ at different redshift $z$ with the data sets of Planck, WMAP 9 years temperature and polarization power spectra, combining with the baryon acoustic oscillation (BAO) from galaxy survey and type Ia supernovae (SN) Union 2.1 sample respectively. The results show that reconstructed $x_e(z)$ is consistent with instantaneous behavior, however, there exists slight deviation from this behavior at some epoch. With PCA method, after abandoning the noisy modes, we get stronger constraints, and the hints for featured $x_e(z)$ evolution could become a little more obvious.

In this paper, we used standard rulers and standard candles (separately and jointly) to explore five popular dark energy models under assumption of spatial flatness of the Universe. As standard rulers, we used a data set comprising 118 galactic-scale strong lensing systems (individual standard rulers if properly calibrated for the mass density profile) combined with BAO diagnostics (statistical standard ruler). Supernovae Ia served asstandard candles. Unlike in the most of previous statistical studies involving strong lensing systems, we relaxed the assumption of singular isothermal sphere (SIS) in favor of its generalization: the power-law mass density profile. Therefore, along with cosmological model parameters we fitted the power law index and its first derivative with respect to the redshift (thus allowing for mass density profile evolution). It turned out that the best fitted $\gamma$ parameters are in agreement with each other irrespective of the cosmological model considered. This demonstrates that galactic strong lensing systems may provide a complementary probe to test the properties of dark energy. Fits for cosmological model parameters which we obtained are in agreement with alternative studies performed by the others. Because standard rulers and standard candles have different parameter degeneracies, combination of standard rulers and standard candles gives much more restrictive results for cosmological parameters. At last, we attempted at model selection based on information theoretic criteria (AIC and BIC). Our results support the claim, that cosmological constant model is still the best one and there is no (at least statistical) reason to prefer any other more complex model.

Cosmic-ray muons and especially their secondaries break apart nuclei ("spallation") and produce fast neutrons and beta-decay isotopes, which are backgrounds for low-energy experiments. In Super-Kamiokande, these beta decays are the dominant background in 6--18 MeV, relevant for solar neutrinos and the diffuse supernova neutrino background. In a previous paper, we showed that these spallation isotopes are produced primarily in showers, instead of in isolation. This explains an empirical spatial correlation between a peak in the muon Cherenkov light profile and the spallation decay, which Super-Kamiokande used to develop a new spallation cut. However, the muon light profiles that Super-Kamiokande measured are grossly inconsistent with shower physics. We show how to resolve this discrepancy and how to reconstruct accurate profiles of muons and their showers from their Cherenkov light. We propose a new spallation cut based on these improved profiles and quantify its effects. Our results can significantly benefit low-energy studies in Super-Kamiokande, and will be especially important for detectors at shallower depths, like the proposed Hyper-Kamiokande.

We have constrained the charge-mass ($\varepsilon-m$) phase space of millicharged particles through the simulation of the rotational evolution of neutron stars, where an extra slow-down effect due to the accretions of millicharged dark matter particles is considered. For a canonical neutron star of $M=1.4~M_{\odot}$ and $R=10~{\rm km}$ with typical magnetic field strength $B_{0}=10^{12}$ G, we have shown an upper limit of millicharged particles, which is compatible with recently experimental and observational bounds. Meanwhile, we have also explored the influences on the $\varepsilon-m$ phase space of millicharged particles for different magnetic fields $B_{0}$ and dark matter density $\rho_{\rm{DM}}$ in the vicinity of the neutron star.

Previous works suggested that the state transitions in an X-ray binary can be triggered by accreting inverse magnetic field from its companion star. A key point of this mechanism is the accretion and magnification of large-scale magnetic fields from outer boundary of a thin disk. However, how such a process can be realized is still an open question. In this work, we check this issue in a realistic X-ray binary system. According to our calculations, a quite strong initial magnetic field $B\sim 10^2-10^3$ G is required in order to assure that the large-scale magnetic field can be effectively dragged inward and magnified with the accretion of gas. Thus, such a picture probably can be present in high-mass X-ray binaries possessing strong stellar magnetic field, e.g., Cyg X-1.

We report the results of a weakly-interacting massive particle (WIMP) dark matter search using the full 80.1\;live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting low-mass WIMPs, achieving a photon detection efficiency of 9.6\%. With a fiducial liquid xenon target mass of 54.0\,kg, no significant excess event were found above the expected background. A profile likelihood analysis confirms our earlier finding that the PandaX-I data disfavor all positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent bound on the low mass WIMP is set at WIMP mass below 10\,GeV/c$^2$, demonstrating that liquid xenon detectors can be competitive for low-mass WIMP searches.

Aiming at exploring the nature of dark energy, we use thirty-six observational Hubble parameter data (OHD) in the redshift range $0 \leqslant z \leqslant 2.36$ to make a cosmological model-independent test of the two-point $Omh^2(z_{2};z_{1})$ diagnostic. In $\Lambda$CDM, we have $Omh^2 \equiv \Omega_{m}h^2$, where $\Omega_{m}$ is the matter density parameter at present. We bin all the OHD into four data points to mitigate the observational contaminations. By comparing with the value of $\Omega_{m}h^2$ which is constrained tightly by the Planck observations, our results show that in all six testing pairs of $Omh^2$ there are two testing pairs are consistent with $\Lambda$CDM at $1\sigma$ confidence level (CL), whereas for another two of them $\Lambda$CDM can only be accommodated at $2\sigma$ CL. Particularly, for remaining two pairs, $\Lambda$CDM is not compatible even at $2\sigma$ CL. Therefore it is reasonable that although deviations from $\Lambda$CDM exist for some pairs, cautiously, we cannot rule out the validity of $\Lambda$CDM. We further apply two methods to derive the value of Hubble constant $H_0$ utilizing the two-point $Omh^2(z_{2};z_{1})$ diagnostic. We obtain $H_0 = 71.23\pm1.54$ ${\mathrm{km \ s^{-1} \ Mpc^{-1}}}$ from inverse variance weighted $Omh^2$ value (method (I)) and $H_0 = 69.37\pm1.59$ ${\mathrm{km \ s^{-1} \ Mpc^{-1}}}$ that the $Omh^2$ value originates from Planck measurement (method (II)), both at $1\sigma$ CL. Finally, we explore how the error in OHD propagate into $w(z)$ at certain redshift during the reconstruction of $w(z)$. We argue that the current precision on OHD is not sufficient small to ensure the reconstruction of $w(z)$ in an acceptable error range, especially at the low redshift

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $\gamma$-ray emission from 3C 275.1. The $\gamma$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $\gamma$-ray emission from other LDQs. 3C 275.1 is the known $\gamma$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $\gamma$-ray data. The parsec jet model, however, is favored by the potential $\gamma$-ray variability at the timescale of months. We suggest that some dimmer $\gamma$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $\gamma$-ray background.

Searching for the signal of primordial gravitational waves in the B-modes (BB) power spectrum is one of the key scientific aims of the cosmic microwave background (CMB) polarization experiments. However, this could be easily contaminated by several foreground issues, such as the thermal dust emission. In this paper we study another mechanism, the cosmic birefringence, which can be introduced by a CPT-violating interaction between CMB photons and an external scalar field. Such kind of interaction could give rise to the rotation of the linear polarization state of CMB photons, and consequently induce the CMB BB power spectrum, which could mimic the signal of primordial gravitational waves at large scales. With the recent polarization data of BICEP2 and the joint analysis data of BICEP2/Keck Array and Planck, we perform a global fitting analysis on constraining the tensor-to-scalar ratio $r$ by considering the polarization rotation angle which can be separated into a background isotropic part and a small anisotropic part. Since the data of BICEP2 and Keck Array experiments have already been corrected by using the "self-calibration" method, here we mainly focus on the effects from the anisotropies of CMB polarization rotation angle. We find that including the anisotropies in the analysis could slightly weaken the constraints on $r$, when using current CMB polarization measurements. We also simulate the mock CMB data with the BICEP3-like sensitivity. Very interestingly, we find that if the effects of the anisotropic polarization rotation angle can not be taken into account properly in the analysis, the constraints on $r$ will be dramatically biased. This implies that we need to break the degeneracy between the anisotropies of the CMB polarization rotation angle and the CMB primordial tensor perturbations, in order to measure the signal of primordial gravitational waves accurately.

We provide a general analysis on the properties of emitting material of some rapidly evolving and luminous transients discovered recently with the Pan-STARRS1 Medium Deep Survey. It is found that these transients are probably produced by a low-mass non-relativistic outflow that is continuously powered by a newly born, rapidly spinning, and highly magnetized neutron star. Such a system could originate from an accretion-induced collapse of a white dwarf or a merger of a neutron star-neutron star binary. Therefore, observations to these transients would be helpful for constraining white dwarf and neutron star physics and/or for searching and identifying gravitational wave signals from the mergers.

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

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

Chinese Spectral Radio Heliograph (CSRH) is an advanced aperture synthesis solar radio heliograph, developed by National Astronomical Observatories, Chinese Academy of Sciences independently. It consists of 100 reflector antennas, which are grouped into two antenna arrays (CSRH-I and CSRH-II) for low and high frequency bands respectively. The frequency band of CSRH-I is 0.4-2GHz and for CSRH-II, the frequency band is 2-15GHz. In the antenna and feed system, CSRH uses an Eleven feed to receive signals coming from the Sun, the radiation pattern with lower side lobe and back lobe of the feed is well radiated. The characteristics of gain G and antenna noise temperature T effect the quality of solar radio imaging. For CSRH, measured G is larger than 60 dBi and $ T $ is less than 120K, after CSRH-I was established, we have successfully captured a solar radio burst between 1.2-1.6GHz on November 12, 2010 through this instrument and this event was confirmed through the observation of Solar Broadband Radio Spectrometer (SBRS) at 2.84GHz and Geostationary Operational Environmental Satellite (GOES). In addition, an image obtained from CSRH-I clearly reveals the profile of the solar radio burst. The other observational work is the imaging of Fengyun-2E geosynchronous satellite which is assumed to be a point source. This data processing method indicates that, the method of deleting errors about dirty image could be used for processing other surface sources.

Crucial questions about solar and supernova neutrinos remain unanswered. Super-Kamiokande has the exposure needed for progress, but detector backgrounds are a limiting factor. A leading component is the beta decays of isotopes produced by cosmic-ray muons and their secondaries, which initiate nuclear spallation reactions. Cuts of events after and surrounding muon tracks reduce this spallation decay background by $\simeq 90\%$ (at a cost of $\simeq 20\%$ deadtime), but its rate at 6--18 MeV is still dominant. A better way to cut this background was suggested in a Super-Kamiokande paper [Bays {\it et al.}, Phys.~Rev.~D {\bf 85}, 052007 (2012)] on a search for the diffuse supernova neutrino background. They found that spallation decays above 16 MeV were preceded near the same location by a peak in the apparent Cherenkov light profile from the muon; a more aggressive cut was applied to a limited section of the muon track, leading to decreased background without increased deadtime. We put their empirical discovery on a firm theoretical foundation. We show that almost all spallation decay isotopes are produced by muon-induced showers and that these showers are rare enough and energetic enough to be identifiable. This is the first such demonstration for any detector. We detail how the physics of showers explains the peak in the muon Cherenkov light profile and other Super-K observations. Our results provide a physical basis for practical improvements in background rejection that will benefit multiple studies. For solar neutrinos, in particular, it should be possible to dramatically reduce backgrounds at energies as low as 6 MeV.

We restudy the issue of cross-correlations between broad-line and jet emission variations, and aim to locate the position of radio (and gamma-ray) emitting region in jet of active galactic nuclei (AGNs). Considering the radial profiles of the radius and number density of clouds in a spherical broad-line region (BLR), we derive new formulae connecting jet emitting position $R_{\rm{jet}}$ to time lag $\tau_{\rm{ob}}$ between broad-line and jet emission variations, and BLR radius. Also, formulae are derived for a disk-like BLR and a spherical shell BLR. The model-independent FR/RSS method is used to estimate $\tau_{\rm{ob}}$. For 3C 120, positive lags of about 0.3 yr are found between the 15 GHz emission and the H$\beta$, H$\gamma$ and He II $\lambda 4686$ lines, including broad-line data in a newly published paper, indicating the line variations lead the 15 GHz ones. Each of the broad-line light curves corresponds to a radio outburst. $R_{\rm{jet}}=$1.1--1.5 parsec (pc) are obtained for 3C 120. For 3C 273, a common feature of negative time lags is found in the cross-correlation functions between light curves of radio emission and the Balmer lines, and as well Ly$\alpha \/\ \lambda 1216$ and C~{\sc iv} $\lambda 1549$ lines. $R_{\rm{jet}}=$ 1.0--2.6 pc are obtained for 3C 273. The estimated $R_{\rm{jet}}$ are comparable for 3C 120 and 3C 273, and the gamma-ray emitting positions will be within $\sim$ 1--3 pc from the central engines. Comparisons show that the cloud number density and radius radial distributions and the BLR structures only have negligible effects on $R_{\rm{jet}}$.

We perform 3D relativistic ideal MHD simulations to study the collisions between high-$\sigma$ (Poynting-flux-dominated) blobs which contain both poloidal and toroidal magnetic field components. This is meant to mimic the interactions inside a highly variable Poynting-flux-dominated jet. We discover a significant electromagnetic field (EMF) energy dissipation with an Alfv\'enic rate with the efficiency around 35\%. Detailed analyses show that this dissipation is mostly facilitated by the collision-induced magnetic reconnection. Additional resolution and parameter studies show a robust result that the relative EMF energy dissipation efficiency is nearly independent of the numerical resolution or most physical parameters in the relevant parameter range. The reconnection outflows in our simulation can potentially form the multi-orientation relativistic mini-jets as needed for several analytical models. We also find a linear relationship between the $\sigma$ values before and after the major EMF energy dissipation process. Our results give support to the proposed astrophysical models that invoke significant magnetic energy dissipation in Poynting-flux-dominated jets, such as the internal collision-induced magnetic reconnection and turbulence (ICMART) model for GRBs, and reconnection triggered mini-jets model for AGNs.

Cadmium Zinc Telluride and Cadmium Telluride are the detector materials of choice for the detection of X-rays in the X-ray energy band E >= 5keV with excellent spatial and spectral resolution and without cryogenic cooling. Owing to recent breakthroughs in grazing incidence mirror technology, next-generation hard X-ray telescopes will achieve angular resolution between 5 and 10 arc seconds - about an order of magnitude better than that of the NuSTAR hard X-ray telescope. As a consequence, the next generation of X-ray telescopes will require pixelated X-ray detectors with pixels on a grid with a lattice constant of <= 250um. Additional detector requirements include a low energy threshold of less than 5keV and an energy resolution of less than one keV. The science drivers for a high angular-resolution X-ray mission include studies and measurements of black hole spins, the cosmic evolution of super-massive black holes, active galactic nuclei feedback, and the behaviour of matter at very high densities. In this contribution, we report on our R&D studies with the goal to optimise small-pixel Cadmium Zinc Telluride and Cadmium Telluride detectors.

In this paper we use the current and future cosmic microwave background (CMB) experiments to test the Charge-Parity-Time Reversal (CPT) symmetry. We consider a CPT-violating interaction in the photon sector $\mathcal{L}_{\rm cs}\sim p_\mu A_\nu \tilde{F}^{\mu\nu}$ which gives rise to a rotation of the polarization vectors of the propagating CMB photons. By combining the nine-year WMAP, BOOMERanG 2003 and BICEP1 observations, we obtain the current constraint on the isotropic rotation angle $\bar{\alpha} = -2.12 \pm 1.14$ ($1\sigma$), indicating an about $2\sigma$ significance of the CPT violation. Here, we particularly take the systematic errors of CMB measurements into account. Then, we study the effects of the anisotropies of the rotation angle [$\Delta{\alpha}({\bf \hat{n}})$] on the CMB polarization power spectra in detail. Due to the small effects, the current CMB polarization data can not constrain the related parameters very well. We obtain the 95\% C.L. upper limit of the variance of the anisotropies of the rotation angle $C^\alpha(0) < 0.035$ from all the CMB datasets. More interestingly, including the anisotropies of rotation angle could lower the best fit value of $r$ and relax the tension on the constraints of $r$ between BICEP2 and Planck. Finally, we investigate the capabilities of future Planck polarization measurements on $\bar{\alpha}$ and $\Delta{\alpha}({\bf \hat{n}})$. Benefited from the high precision of Planck data, the constraints of the rotation angle can be significantly improved.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

Data access and interoperability module connects the observation proposals, data, virtual machines and software. According to the unique identifier of PI (principal investigator), an email address or an internal ID, data can be collected by PI's proposals, or by the search interfaces, e.g. conesearch. Files associated with the searched results could be easily transported to cloud storages, including the storage with virtual machines, or several commercial platforms like Dropbox. Benefitted from the standards of IVOA (International Observatories Alliance), VOTable formatted searching result could be sent to kinds of VO software. Latter endeavor will try to integrate more data and connect archives and some other astronomical resources.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences). The ultimate goal of this project is to provide a comprehensive end-to-end astronomy research environment where several independent systems seamlessly collaborate to support the full lifecycle of the modern observational astronomy based on big data, from proposal submission, to data archiving, data release, and to in-situ data analysis and processing. In this paper, the architecture and key designs of the AstroCloud platform are introduced, including data access middleware, access control and security framework, extendible proposal workflow, and system integration mechanism.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences){\url{http://astrocloud.china-vo.org}}\citep{O8-5_Cui_adassxxiv}. To archive the astronomical data in China, we present the implementation of the astronomical data archiving system (ADAS). Data archiving and quality control are the infrastructure for the AstroCloud. Throughout the data of the entire life cycle, data archiving system standardized data, transferring data, logging observational data, archiving ambient data, And storing these data and metadata in database. Quality control covers the whole process and all aspects of data archiving.

Using high-resolution, multiple-passband Hubble Space Telescope images spanning the entire optical/near-infrared wavelength range, we obtained a statistically complete sample, $U$-band selected sample of 846 extended star clusters across the disk of the nearby starburst galaxy M82. Based on careful analysis of their spectral energy distributions, we determined their galaxy-wide age and mass distributions. The M82 clusters exhibit three clear peaks in their age distribution, thus defining a relatively young, log(t/yr) < 7.5, an intermediate-age, log(t/yr) $\in$ [7.5, 8.5], and an old sample, log(t/yr) > 8.5. Comparison of the completeness-corrected mass distributions offers a firm handle on the galaxy's star cluster disruption history. The most massive star clusters in the young and old samples are (almost) all concentrated in the most densely populated central region, while the intermediate-age sample's most massive clusters are more spatially dispersed, which may reflect the distribution of the highest-density gas throughout the galaxy's evolutionary history, combined with the solid-body nature of the galaxy's central region.

We consider the origin of compact, short-period, Jupiter-mass planets. We propose that their diverse structure is caused by giant impacts of embryos and super-Earths or mergers with other gas giants during the formation and evolution of these hot Jupiters. Through a series of numerical simulations, we show that typical head-on collisions generally lead to total coalescence of impinging gas giants. Although extremely energetic collisions can disintegrate the envelope of gas giants, these events seldom occur. During oblique and moderately energetic collisions, the merger products retain higher fraction of the colliders' cores than their envelopes. They can also deposit considerable amount of spin angular momentum to the gas giants and desynchronize their spins from their orbital mean motion. We find that the oblateness of gas giants can be used to infer the impact history. Subsequent dissipation of stellar tide inside the planets' envelope can lead to runaway inflation and potentially a substantial loss of gas through Roche-lobe overflow. The impact of super-Earths on parabolic orbits can also enlarge gas giant planets' envelope and elevates their tidal dissipation rate over $\sim $ 100 Myr time scale. Since giant impacts occur stochastically with a range of impactor sizes and energies, their diverse outcomes may account for the dispersion in the mass-radius relationship of hot Jupiters.

We carried out two-dimensional high-resolution simulations to study the effect of dust feedback on the evolution of vortices induced by massive planets in protoplanetary disks. Various initial dust to gas disk surface density ratios ($0.001$ -- $0.01$) and dust particle sizes (Stokes number $4\times10^{-4}$ -- $0.16$) are considered. We found that while dust particles migrate inwards, vortices are very effective in collecting them. When dust density becomes comparable to gas density within the vortex, a dynamical instability is excited and it alters the coherent vorticity pattern and destroys the vortex. This dust feedback effect is stronger with higher initial dust/gas density ratio and larger dust grain. Consequently, we found that the disk vortex lifetime can be reduced up to a factor of 10. We discuss the implications of our findings on the survivability of vortices in protoplanetary disks and planet formation.

We report on the first dark-matter (DM) search results from PandaX-I, a low threshold dual-phase xenon experiment operating at the China Jinping Underground Laboratory. In the 37-kg liquid xenon target with 17.4 live-days of exposure, no DM particle candidate event was found. This result sets a stringent limit for low-mass DM particles and disfavors the interpretation of previously-reported positive experimental results. The minimum upper limit, $3.7\times10^{-44}$\,cm$^2$, for the spin-independent isoscalar DM-particle-nucleon scattering cross section is obtained at a DM-particle mass of 49\,GeV/c$^2$ at 90\% confidence level.

In this paper, we set the new limits on the equation of state parameter (EoS) of dark energy with the observations of cosmic microwave background radiation (CMB) from Planck satellite, the type Ia supernovae from Pan-STARRS and the baryon acoustic oscillation (BAO). We consider two parametrization forms of EoS: a constant $w$ and time evolving $w(a)=w_0+w_a(1-a)$. The results show that with a constant EoS, $w=-1.141\pm{0.075}$ ($68\%~C.L.$), which is consistent with $\Lambda$CDM at about $2\sigma$ confidence level. For a time evolving $w(a)$ model, we get $w_0=-1.09^{+0.16}_{-0.18}$ ($1\sigma~C.L.$), $w_a=-0.34^{+0.87}_{-0.51}$ ($1\sigma~C.L.$), and in this case $\Lambda$CDM can be comparable with our observational data at $1\sigma$ confidence level. In order to do the parametrization independent analysis, additionally we adopt the so called principal component analysis (PCA) method, in which we divide redshift range into several bins and assume $w$ as a constant in each redshift bin (bin-w). In such bin-w scenario, we find that for most of the bins cosmological constant can be comparable with the data, however, there exists few bins which give $w$ deviating from $\Lambda$CDM at more than $2\sigma$ confidence level, which shows a weak hint for the time evolving behavior of dark energy. To further confirm this hint, we need more data with higher precision.

Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulationsof the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal component that is continuously being injected into the domain. This flux injection mimics the poloidal currents driven by the anode-cathode voltage drop in the experiment. The injected toroidal field stretches the poloidal fields to large distances, while forming a collimated jet along with several other key features. Detailed comparisons between 3D MHD simulations and experimental measurements provide a comprehensive description of the interplay among magnetic force, pressure and flow effects. In particular, we delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms. With suitably chosen parameters that are derived from experiments, the jet in the simulation agrees quantitatively with the experimental jet in terms of magnetic/kinetic/inertial energy, total poloidal current, voltage, jet radius, and jet propagation velocity. Specifically, the jet velocity in the simulation is proportional to the poloidal current divided by the square root of the jet density, in agreement with both the experiment and analytical theory. This work provides a new and quantitative method for relating experiments, numerical simulations and astrophysical observation, and demonstrates the possibility of using terrestrial laboratory experiments to study astrophysical jets.

Among $\sim 160$ of the multiple exoplanetary systems confirmed, about $30\%$ of them have neighboring pairs with a period ratio $\leq 2$. A significant fraction of these pairs are around mean motion resonance (MMR), more interestingly, peak around 2:1 and 3:2, with a clear absence of more closely packed MMRs with period ratios less than 4:3, regardless of planet masses. Here we report numerical simulations demonstrating that such MMR behavior places important constraints on the disk evolution stage out of which the observed planets formed. Multiple massive planets (with mass $\geq 0.8$ $M_{\rm Jup}$) tend to end up with a 2:1 MMR mostly independent of the disk masses but low-mass planets (with mass $\leq 30$ $M_{\oplus}$) can have MMRs larger than 4:3 only when the disk mass is quite small, suggesting that the observed dynamical architecture of most low-mass-planet pairs was established late in the disk evolution stage, just before it was dispersed completely.

We present a model-independent determination of the curvature parameter $\Omega_k$ by using the Hubble parameter $H(z)$ and angular diameter distance $D_A(z)$ from the recent baryon acoustic oscillation (BAO) measurements. Each $H(z)$ and $D_A(z)$ pair from a BAO measurement can constrain a curvature parameter. The accuracy of the curvature measurement improves with increased redshift of $H(z)$ and $D_A(z)$ data. By using the $H(z)$ and $D_A(z)$ pair derived from BAO Lyman $\alpha$ forest measurement at $z=2.36$, the $\Omega_k$ is confined to be -0.05$\pm$0.06, which is consistent with the curvature $-0.037^{+0.044}_{-0.042}$ constrained by the nine-year WMAP data only. Considering future BAO meausurements, at least one order of magnitude improvement of this curvature measurement could be expected.

Recent observations of large-scale asymmetric features in protoplanetary disks suggest that large-scale vortices exist in such disks. Massive planets are known to be able to produce deep gaps in protoplanetary disks. The gap edges could become hydrodynamically unstable to the Rossby wave/vortex instability and form large-scale vortices. In this study we examine the long term evolution of these vortices by carrying out high-resolution two dimensional hydrodynamic simulations that last more than $10^4$ orbits (measured at the planet's orbit). We find that the disk viscosity has a strong influence on both the emergence and lifetime of vortices. In the outer disk region where asymmetric features are observed, our simulation results suggest that the disk viscous $\alpha$ needs to be low $\sim 10^{-5 }$ - $10^{-4}$ to sustain vortices to thousands and up to $10^{4}$ orbits in certain cases. The chance of finding a vortex feature in a disk then decreases with smaller planet orbital radius. For $\alpha \sim 10^{-3}$ or larger, even planets with masses of 5 Jupiter-masses will have difficulty either producing or sustaining vortices. We have also studied the effects of different disk temperatures and planet masses. We discuss the implications of our findings on current and future protoplanetary disk observations.

The recently released BICEP2 data detected the primordial B-mode polarization in the Cosmic Microwave Background (CMB) map which strongly supports for a large tensor-to-scalar ratio, and thus, is found to be in tension with the Planck experiment with no evidence of primordial gravitational waves. Such an observational tension, if confirmed by forthcoming measurements, would bring a theoretical challenge for the very early universe models. To address this issue, we in the present paper revisit a single field inflation model, which includes a modulated potential. We show that this inflation model can give rise to a sizable negative running behavior for the spectral index of primordial curvature perturbation and a large tensor-to-scalar ratio. Applying these properties, our model can nicely explain the combined Planck and BICEP2 observations. To examine the validity of analytic calculations, we numerically confront the predicted temperature and B-mode power spectra with the latest CMB observations and explicitly show that our model is consistent with the current data.

Galaxy mergers play a key role in the evolution of galaxies and the growth of their central supermassive black holes (SMBHs). A search for (active) SMBH binaries (SMBHBs) at the centers of the merger remnants is currently ongoing. Perhaps the greatest challenge is to identify the inactive SMBHBs, which might be the most abundant, but are also the most difficult to identify. Liu et al. predicted characteristic drops in the light curves of tidal disruption events (TDEs), caused by the presence of a secondary SMBH. Here, we apply that model to the light curve of the optically inactive galaxy SDSS J120136.02+300305.5, which was identified as a candidate TDE with XMM-Newton. We show that the deep dips in its evolving X-ray light curve can be well explained by the presence of a SMBHB at its core. A SMBHB model with a mass of the primary of M_BH = 10^7 M_sun, a mass ratio q=0.08, and a semimajor axis a_b = 0.6 mpc is in good agreement with the observations. Given that primary mass, introducing an orbital eccentricity is needed, with e_b = 0.3. Alternatively, a lower mass primary of M_BH = 10^6 M_sun in a circular orbit fits the light curve well. Tight binaries like this one, which have already overcome the "final parsec problem," are prime sources of gravitational wave radiation once the two SMBHs coalesce. Future transient surveys, which will detect TDEs in large numbers, will place tight constraints on the SMBHB fraction in otherwise non-active galaxies.

A maximum jet efficiency line $R\sim 25$ ($R=L_{\rm jet}/L_{\rm bol}$), found in FRII radio galaxies by Fernandes et al., was extended to cover the full range of jet power by Punsly. Recent general relativistic magnetohydrodynamic (GRMHD) simulations on jet formation mainly focused on the enhancement of jet power. In this work, we suggest that the jet efficiency could be very high even for conventional jet power if the radiative efficiency of disk were much smaller. We adopt the model of a thin disk with magnetically driven winds to investigate the observational high efficiency jets in FRII radio galaxies. It is found that the structure of a thin disk can be significantly altered by the feedback of winds. The temperature of disk gradually decreases with increasing magnetic field; the disk density, surface density and pressure also change enormously. The lower temperature and higher surface density in inner disk result in the rapid decrease of radiative efficiency. Thus, the jet efficiency is greatly improved even the jet power is conventional. Our results can explain the observations quite well. A theoretical maximum jet efficiency $R \sim 1000$ suggested by our calculations is large enough to explain all the high jet efficiency in observations even considering the episodic activity of jets.

We present high spatial resolution X-ray spectroscopy of supernova remnant G292.0+1.8 with the {\sl Chandra} observations. The X-ray emitting region of this remnant was divided into 25 $\times$ 25 pixels with a scale of 20$\arcsec$ $\times$ 20$\arcsec$ each. Spectra of 324 pixels were created and fitted with an absorbed one component non-equilibrium ionization model. With the spectral analysis results we obtained maps of absorbing column density, temperature, ionization age, and the abundances for O, Ne, Mg, Si, S, and Fe. The abundances of O, Ne and Mg show tight correlations between each other in the range of about two orders of magnitude, suggesting them all from explosive C/Ne burning. Meanwhile, the abundances of Si and S are also well correlated, indicating them to be the ashes of explosive explosive O-burning or incomplete Si-burbing. The Fe emission lines are not prominent among the whole remnant, and its abundance are significantly deduced, indicating that the reverse shock may have not propagated to the Fe-rich ejecta. Based on relative abundances of O, Ne, Mg, Si and Fe to Si, we suggest a progenitor mass of $25-30 M_{\odot}$ for this remnant.

Star formation can be triggered by compression from wind or supernova driven shock waves that sweep over molecular clouds. Because these shocks will likely contain processed elements, triggered star formation has been proposed as an explanation for short lived radioactive isotopes (SLRI) in the Solar System. Previous studies have tracked the triggering event to the earliest phases of collapse and have focused on the shock properties required for both successful star formation and mixing of SLRI's. In this paper, we use Adaptive Mesh Refinement (AMR) simulation methods, including sink particles, to simulate the full collapse and subsequent evolution of a stable Bonnor-Ebert sphere subjected to a shock and post-shock wind. We track the flow of the cloud material after a star (a sink particle) has formed. For non-rotating clouds we find robust triggered collapse and little bound circumstellar material remaining around the post-shock collapsed core. When we add initial cloud rotation we observe the formation of disks around the collapsed core which then interact with the post-shock flow. Our results indicate that these circumstellar disks are massive enough to form planets and are long-lived, in spite of the ablation driven by post-shock flow ram pressure. As a function of the initial conditions, we also track the time evolution of the accretion rates and particle mixing between between the ambient wind and cloud material. The latter is maximized for cases of highest mach number.

When muons travel through matter, their energy losses lead to nuclear breakup ("spallation") processes. The delayed decays of unstable daughter nuclei produced by cosmic-ray muons are important backgrounds for low-energy astrophysical neutrino experiments, e.g., those seeking to detect solar neutrino or Diffuse Supernova Neutrino Background (DSNB) signals. Even though Super-Kamiokande has strong general cuts to reduce these spallation-induced backgrounds, the remaining rate before additional cuts for specific signals is much larger than the signal rates for kinetic energies of about 6 -- 18 MeV. Surprisingly, there is no published calculation of the production and properties of these backgrounds in water, though there are such studies for scintillator. Using the simulation code FLUKA and theoretical insights, we detail how muons lose energy in water, produce secondary particles, how and where these secondaries produce isotopes, and the properties of the backgrounds from their decays. We reproduce Super-Kamiokande measurements of the total background to within a factor of 2, which is good given that the isotope yields vary by orders of magnitude and that some details of the experiment are unknown to us at this level. Our results break aggregate data into component isotopes, reveal their separate production mechanisms, and preserve correlations between them. We outline how to implement more effective background rejection techniques using this information. Reducing backgrounds in solar and DSNB studies by even a factor of a few could help lead to important new discoveries.

In this paper, we study the cosmological constraints from the measurements of Hubble parameters---$H(z)$ data. Here, we consider two kinds of $H(z)$ data: the direct $H_0$ probe from the Hubble Space Telescope (HST) observations of Cepheid variables with $H_0=73.8\pm2.4$ ${\rm km\,s^{-1}\,Mpc^{-1}}$ and several measurements on the Hubble parameter at high redshifts $H(z)$. Employing Markov Chain Monte Carlo method, we also combine the WMAP nine-year data (WMAP9), the baryon acoustic oscillations (BAO) and type Ia supernovae (SNIa) "Union2.1" compilation to determine the cosmological parameters, such as the equation of state (EoS) of dark energy $w$, the curvature of the universe $\Omega_k$, the total neutrino mass $\sum{m_\nu}$, the effective number of neutrinos $N_{\rm eff}$, and the parameters associated with the power spectrum of primordial fluctuations. These $H(z)$ data provide extra information on the accelerate rate of our Universe at high redshifts. Therefore, adding these $H(z)$ data significantly improves the constraints on cosmological parameters, such as the number of relativistic species. Moreover, we find that direct prior on $H_0$ from HST can also give good constraints on some parameters, due to the degeneracies between these parameters and $H_0$.

The absence of thermal instability in the high/soft state of black hole X-ray binaries, in disagreement with the standard thin disk theory, is a long-standing riddle for theoretical astronomers. We have tried to resolve this question by studying the thermal stability of a thin disk with magnetically driven winds in the $\dot{M}- \Sigma$ plane. It is found that disk winds can greatly decrease the disk temperature and thus help the disk become more stable at a given accretion rate. The critical accretion rate $\dot{M}_{\rm crit}$ corresponding to the thermal instability threshold is increased significantly in the presence of disk winds. For $\alpha=0.01$ and $B_{\rm \phi}=10B_{\rm _p}$, the disk is quite stable even for a very weak initial poloidal magnetic field [$\beta_{\rm p,0}\sim 2000, \beta_{\rm p}=(P_{\rm {gas}}+P_{\rm rad})/(B_{\rm p}^2/8\pi)$]. But when $B_{\rm \phi}=B_{\rm _p}$ or $B_{\rm \phi}=0.1B_{\rm _p}$, a somewhat stronger (but still weak) field ($\beta_{\rm p,0}\sim 200$ or $\beta_{\rm p,0}\sim 20$) is required to make the disk stable. Nevertheless, despite the great increase of $\dot{M}_{\rm crit}$, the luminosity threshold corresponding to instability remains almost constant or even decreases slowly with increasing $\dot{M}_{\rm crit}$ due to the decrease of gas temperature. The advection and diffusion timescales of the large-scale magnetic field threading the disk are also investigated in this work. We find that the advection timescale can be smaller than the diffusion timescale in a disk with winds, because the disk winds take away most of the gravitational energy released in the disk, resulting in the decrease of the magnetic diffusivity $\eta$ and the increase of the diffusion timescale.

We present a systematic method for evaluation of perturbation observables in non-canonical single-field inflation models within the slow-roll approximation, which allied with field redefinitions enables predictions to be established for a wide range of models. We use this to investigate various non-canonical inflation models, including Tachyon inflation and DBI inflation. The Lambert $W$ function will be used extensively in our method for the evaluation of observables. In the Tachyon case, in the slow-roll approximation the model can be approximated by a canonical field with a redefined potential, which yields predictions in better agreement with observations than the canonical equivalents. For DBI inflation models we consider contributions from both the scalar potential and the warp geometry. In the case of a quartic potential, we find a formula for the observables under both non-relativistic and relativistic behaviour of the scalar DBI inflaton. For a quadratic potential we find two branches in the non-relativistic case, determined by the competition of model parameters, while for the relativistic case we find consistency with results already in the literature. We present a comparison to the latest Planck satellite observations. Most of the non-canonical models we investigate, including the Tachyon, are better fits to data than canonical models with the same potential, but we find that DBI models in the slow-roll regime have difficulty in matching the data.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

We present direct astrophysical N-body simulations with up to a few million bodies using our parallel MPI/CUDA code on large GPU clusters in China, Ukraine and Germany, with different kinds of GPU hardware. These clusters are directly linked under the Chinese Academy of Sciences special GPU cluster program in the cooperation of ICCS (International Center for Computational Science). We reach about the half the peak Kepler K20 GPU performance for our phi-GPU code [2], in a real application scenario with individual hierarchically block time-steps with the high (4th, 6th and 8th) order Hermite integration schemes and a real core-halo density structure of the modeled stellar systems. The code and hardware are mainly used to simulate star clusters [23, 24] and galactic nuclei with supermassive black holes [20], in which correlations between distant particles cannot be neglected.

Relativistic, magnetized jets are observed to propagate to very large distances in many Active Galactic Nuclei (AGN). We use 3D relativistic MHD (RMHD) simulations to study the propagation of Poynting flux-driven jets in AGN. These jets are assumed already being launched from the vicinity ($\sim 10^3$ gravitational radii) of supermassive black holes. Jet injections are characterized by a model described in Li et al. (2006) and we follow the propagation of these jets to ~ parsec scales. We find that these current-carrying jets are always collimated and mildly relativistic. When $\alpha$, the ratio of toroidal-to-poloidal magnetic flux injection, is large the jet is subject to non-axisymmetric current-driven instabilities (CDI) which lead to substantial dissipation and reduced jet speed. However, even with the presence of instabilities, the jet is not disrupted and will continue to propagate to large distances. We suggest that the relatively weak impact by the instability is due to the nature of the instability being convective and the fact that the jet magnetic fields are rapidly evolving on Alfv\'enic timescale. We present the detailed jet properties and show that far from the jet launching region, a substantial amount of magnetic energy has been transformed into kinetic energy and thermal energy, producing a jet magnetization number $\sigma < 1$. In addition, we have also studied the effects of a gas pressure supported "disk" surrounding the injection region and qualitatively similar global jet behaviors were observed. We stress that jet collimation, CDIs, and the subsequent energy transitions are intrinsic features of current-carrying jets.

We discuss a new mechanism which can be responsible for the origin of the primordial perturbation in inflationary models, the inhomogeneous DBI reheating scenario. Light DBI fields fluctuate during inflation, and finally create the density perturbations through modulation of the inflation decay rate. In this note, we investigate the curvature perturbation and its non-Gaussianity from this new mechanism. Presenting generalized expressions for them, we show that the curvature perturbation not only depends on the particular process of decay but is also dependent on the sound speed $c_s$ from the DBI action. More interestingly we find that the non-Gaussianity parameter $f_{NL}$ is independent of $c_s$. As an application we exemplify some decay processes which give a viable and detectable non-Gaussianity. Finally we find a possible connection between our model and the DBI-Curvaton mechanism.

In the paper, we investigate correlation between broad-line and radio variations for broad-line radio galaxy 3C 120. By the z-transformed discrete correlation function method and the model-independent flux randomization/random subset selection (FR/RSS) Monte Carlo method, we find that the broad Hbeta line variations lead the 15 GHz variations. The FR/RSS method shows that the Hbeta line variations lead the radio variations by a factor of tau_{ob}=0.34 +/- 0.01 yr. This time lag can be used to locate the position of emitting region of radio outbursts in jet, on the order of ~ 5 light-years, from the central engine. This distance is much larger than the size of broad-line region. The large separation of the radio outburst emitting region from the broad-line region will observably influence the gamma-ray emission in 3C 120.

In this paper we present a catalog which includes 141 bright candidates ($\leq10.27$ mag, V band) showing the infrared (IR) excess at 22 $\mu$m. Of which, 38 stars are known IR excess stars or disk, 23 stars are double or multiple stars and 4 are Be stars. While the remaining more than 70 stars are identified as the 22 $\mu$m excess candidates in our work. The criterion of selecting candidates is $K_s-[22]_{\mu m}$. All these candidates are selected from \emph{WISE} All-sky data cross-correlated with \emph{Hipparcos} Main Catalog and the likelihood-ratio technique is employed. Considering the effect of background, we introduce the \emph{IRAS} 100 $\mu$m level to exclude the high background. We also estimated the coincidence probability of these sources. In addition, we presented the optical to mid-infrared SEDs and optical images of all the candidates, and gave the observed optical spectra of 6 stars with NAOC's 2.16-m telescope. To measure for the dust amount around each star, the fractional luminosity is also provided. We also test whether our method of selecting IR excess stars can be used to search for extra-solar planets, we cross-matched our catalog with known IR-excess stars having planets but none is matched. Finally, we give the fraction of stars showing IR-excess for different spectral type of main-sequence stars.

We study the interaction of strong shock waves with magnetized clumps. Previous numerical work focused on the simplified scenario in which shocked clumps are immersed in a globally uniform magnetic field that extends through both the clump and the ambient medium. Here we consider the complementary circumstance in which the field is completely self-contained within the clumps. This could arise naturally during clump formation via dynamical or thermal instabilities for example as magnetic field pinches off from the ambient medium. Using our AMR MHD code AstroBEAR, we carry out a series of simulations with magnetized clumps that have different self-contained magnetic field configurations. We find that the clump and magnetic evolution are sensitive to the fraction of magnetic field aligned with versus perpendicular to the shock normal. The relative strength of magnetic pressure and tension in the different field configurations allows us to analytically understand the different cases of post-shock evolution. We also show how turbulence and the mixing it implies depends of the initial field configuration and suggest ways in which observed shock-clump morphology may be used as a proxy for identifying internal field topologies a posteriori.

Context. The steep-spectrum radio quasars (SSRQs) are powerful radio sources, with thermal emission from accretion disk and jet nonthermal emission likely both contributing in the Ultraviolet (UV)/optical luminosity, however the former may play a dominant role. While the UV/optical variability of SSRQs has been poorly studied, little is known on the mechanism of their variability. Aims. We investigate the mechanism of the UV/optical variability of SSRQs. Methods. A sample of eighteen SSRQs has been established in SDSS Stripe 82 region in our previous works, in which the flux and spectral variability have been studied. In this work, we construct the flux-flux diagram using SDSS u and i multi-epoch data for these eighteen SSRQs. The standard accretion disk model is used to fit the flux-flux variations, in order to explore the variability mechanism. Results. The model fit to flux-flux diagram are tuned with fixed black hole mass and varying accretion rate. We found that the flux-flux diagram of all our SSRQs can be qualitatively described by the standard accretion disk model with change of the accretion rate. Although non-thermal jet power-law emission can also qualitatively reproduce the variability, the reasonable accretion rates and black hole masses required to fit the flux-flux variations suggest that a disk emission with variable accretion rate is a plausible description of the data.

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

The ephemeris of Phoebe, the ninth satellite of Saturn, is not very accurate. Previous dynamical models were usually too simplified, the astrometry is heterogeneous and, the Saturn's ephemeris itself is an additionnal source of error. The aim is to improve Phoebe's ephemeris by using a large set of observations, correcting some systematic errors and updating the dynamical model. The dynamical model makes use of the most recent ephemeris of planets and Saturnian satellites. The astrometry of Phoebe is improved by using a compilation of ground-based and space-based observations and by correcting the bias in stellar catalogues used for the reduction. We present an accurate ephemeris of Phoebe with residuals of 0.45 arcsec and with an estimated accuracy of Phoebe's position of less that 100 km on 1990-2020 period.

A measurement of the absolute fluorescence yield of the 337 nm nitrogen band, relevant to ultra-high energy cosmic ray (UHECR) detectors, is reported. Two independent calibrations of the fluorescence emission induced by a 120 GeV proton beam were employed: Cherenkov light from the beam particle and calibrated light from a nitrogen laser. The fluorescence yield in air at a pressure of 1013 hPa and temperature of 293 K was found to be $Y_{337} = 5.61\pm 0.06_{stat} \pm 0.21_{syst}$ photons/MeV. When compared to the fluorescence yield currently used by UHECR experiments, this measurement improves the uncertainty by a factor of three, and has a significant impact on the determination of the energy scale of the cosmic ray spectrum.

We extend the ModeCode software of Mortonson, Peiris and Easther to enable numerical computation of perturbations in K-inflation models, where the scalar field no longer has a canonical kinetic term. Focussing on models where the kinetic and potential terms can be separated into a sum, we compute slow-roll predictions for various models and use these to verify the numerical code. A Markov chain Monte Carlo analysis is then used to impose constraints from WMAP7 data on the addition of a term quadratic in the kinetic energy to the Lagrangian of simple chaotic inflation models. For a quadratic potential, the data do not discriminate against addition of such a term, while for a quartic (\lambda \phi^4) potential inclusion of such a term is actually favoured. Overall, constraints on such a term from present data are found to be extremely weak.

In this paper we investigate the possible direct, non-gravitational interaction between holographic dark energy (HDE) and dark matter. Firstly, we start with two simple models with the interaction terms $Q \propto \rho_{dm}$ and $Q \propto \rho_{de}$, and then we move on to the general form $Q \propto \rho_m^\alpha\rho_{de}^\beta$. The cosmological constraints of the models are obtained from the joint analysis of the present Union2.1+BAO+CMB+$H_0$ data. We find that the data slightly favor an energy flow from dark matter to dark energy, although the original HDE model still lies in the 95.4% confidence level (CL) region. For all models we find $c<1$ at the 95.4% CL. We show that compared with the cosmic expansion, the effect of interaction on the evolution of $\rho_{dm}$ and $\rho_{de}$ is smaller, and the relative increment (decrement) amount of the energy in the dark matter component is constrained to be less than 9% (15%) at the 95.4% CL. By introducing the interaction, we find that even when $c<1$ the big rip still can be avoided due to the existence of a de Sitter solution at $z\rightarrow-1$. We show that this solution can not be accomplished in the two simple models, while for the general model such a solution can be achieved with a large $\beta$, and the big rip may be avoided at the 95.4% CL.

We explore the cosmological implications of five modified gravity (MG) models by using the recent cosmological observational data, including the recently released SNLS3 type Ia supernovae sample, the cosmic microwave background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7-yr observations, the baryon acoustic oscillation results from the Sloan Digital Sky Survey data release 7, and the latest Hubble constant measurement utilizing the Wide Field Camera 3 on the Hubble Space Telescope. The MG models considered include the Dvali-Gabadadze-Porrati(DGP) model, two $f(R)$ models, and two $f(T)$ models. We find that compared with the $\Lambda$CDM model, MG models can not lead to a appreciable reduction of the $\chi^2_{min}$. The analysis of AIC and BIC shows that the simplest cosmological constant model($\Lambda$CDM) is still most preferred by the current data, and the DGP model is strongly disfavored. In addition, from the observational constraints, we also reconstruct the evolutions of the growth factor in these models. We find that the current available growth factor data are not enough to distinguish these MG models from the $\Lambda$CDM model.

Giant X-ray cavities lie in some active galactic nuclei (AGNs) locating in central galaxies of clusters, most of these cavities are thought to be inflated by jets of AGNs. The jets can be either powered by rotating black holes or the accretion disks surrounding black holes, or both. In this work, we choose the most energetic cavity, MS 0735+7421, with stored energy ~ 10^62 erg, to constrain the jet formation mechanisms and the evolution of the central massive black hole in this source. The bolometric luminosity of the AGN in this cavity is ~ 10^(-5) L_Edd, however, the mean power of the jet required to inflate the cavity is estimated as ~ 0.02 L_Edd, which implies that the source has experienced strong outbursts previously. During outbursts, the jet power and the mass accretion rate should be significantly higher than its present values. We construct an accretion disk model, in which the angular momentum and energy carried away by jets is properly included, to calculate the spin and mass evolution of the massive black hole. In our calculations, different jet formation mechanisms are employed, and we find that the jets generated with the Blandford-Znajek (BZ) mechanism are unable to produce the giant cavity with ~ 10^62 erg in this source. Only the jets accelerated with the combination of the Blandford-Payne (BP) and BZ mechanisms can successfully inflate such a giant cavity, if the magnetic pressure is close to equipartition with the total (radiation+gas) pressure of the accretion disk. For dynamo generated magnetic field in the disk, such an energetic giant cavity can be inflated by the magnetically driven jets only if the initial black hole spin parameter a_0 > 0.95. Our calculations show that the final spin parameter a of the black hole is always ~ 0:9 - 0.998 for all the computational examples which can provide sufficient energy for the cavity of MS 0735+7421.

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. The coalescence of the SMBHBs is a distinct source of gravitational wave (GW) radiation. The detections of the strong GW radiation and their possible electromagnetic counterparts are essential. Numerical relativity suggests that the post-merger supermassive black hole (SMBH) gets a kick velocity up to 4000 km/s due to the anisotropic GW radiations. Here we investigate the dynamical co-evolution and interaction of the recoiling SMBHs and their galactic stellar environments with one million direct N-body simulations including the stellar tidal disruption by the recoiling SMBHs. Our results show that the accretion of disrupted stars does not significantly affect the SMBH dynamical evolution. We investigate the stellar tidal disruption rates as a function of the dynamical evolution of oscillating SMBHs in the galactic nuclei. Our simulations show that most of stellar tidal disruptions are contributed by the unbound stars and occur when the oscillating SMBHs pass through the galactic center. The averaged disruption rate is ~10^{-6} M_\odot yr^{-1}, which is about an order of magnitude lower than that by a stationary SMBH at similar galactic nuclei. Our results also show that a bound star cluster is around the oscillating SMBH of about ~ 0.7% the black hole mass. In addition, we discover a massive cloud of unbound stars following the oscillating SMBH. We also investigate the dependence of the results on the SMBH masses and density slopes of the galactic nuclei.

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the amplification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar plasma. The blast wave is driven by injecting high pressure in the simulation domain. The interstellar magnetic field can be amplified by two different processes, occurring in different regions. One is facilitated by the fluid vorticity generated by the ``rippled" shock front interacting with the background turbulence. The resulting turbulent flow keeps amplifying the magnetic field, consistent with earlier work \citep{Giacalone2007}. The other process is facilitated by the growth of the Rayleigh-Taylor instability at the contact discontinuity between the ejecta and the shocked medium. This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We investigate the dependence of the amplification on numerical parameters such as grid-cell size and on various physical parameters. We show the magnetic field has a characteristic radial profile that the downstream magnetic field gets progressively stronger away from the shock. This is because the downstream magnetic field needs a finite time to reach the efficient amplification, and will get further amplified in the Rayleigh-Taylor region. In our simulation we do not observe a systematic strong magnetic field within a small distance to the shock. This indicates that if the magnetic-field amplification in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional instabilities, plasma kinetics and/or cosmic ray effect may need to be considered to explain the strong magnetic field in supernova remnants.

Context. Multi-wavelength observations of gamma-ray burst (GRB) afterglows provide important information about the activity of their central engines and their environments. In particular, the short timescale variability, such as bumps and/or rebrightening features visible in the multi-wavelength light curves, is still poorly understood. Aims. We analyze the multi-wavelength observations of the GRB100219A afterglow at redshift 4.7. In particular, we attempt to identify the physical origin of the late achromatic flares/bumps detected in the X-ray and optical bands. Methods. We present ground-based optical photometric data and Swift X-ray observations on GRB100219A. We analyzed the temporal behavior of the X-ray and optical light curves, as well as the X-ray spectra. Results. The early flares in the X-ray and optical light curves peak simultaneously at about 1000 s after the burst trigger, while late achromatic bumps in the X-ray and optical bands appear at about 20000 s after the burst trigger. These are uncommon features in the afterglow phenomenology. Considering the temporal and spectral properties, we argue that both optical and X-ray emissions come from the same mechanism. The late flares/bumps may be produced by late internal shocks from long-lasting activity of the central engine. An off-axis origin for a structured jet model is also discussed to interpret the bump shapes. The early optical bump can be interpreted as the afterglow onset, while the early X-ray flare could be caused by the internal activity. GRB 100219A exploded in a dense environment as revealed by the strong attenuation of X-ray emission and the optical-to-X-ray spectral energy distribution.

Interfaces between hot and cold magnetized plasmas exist in various astrophysical contexts, for example where hot outflows impinge on an ambient interstellar medium (ISM). It is of interest to understand how the structure of the magnetic field spanning the interface affects the temporal evolution of the temperature gradient. Here we explore the relation between the magnetic field topology and the heat transfer rate by adding various fractions of tangled vs. ordered field across a hot-cold interface allow the system to evolve to a steady state. We find a simple mathematical relation for the rate of heat conduction as a function of the initial ratio of ordered to tangled field across the interface. We discuss potential implications for the astrophysical context of magnetized wind blown bubbles (WBB) around evolved stars.

The physical nature of compressible turbulence is of fundamental importance in a variety of astrophysical settings. We present the first direct evidence that mean kinetic energy cascades conservatively beyond a transitional "conversion" scale-range despite not being an invariant of the compressible flow dynamics. We use high-resolution three-dimensional simulations of compressible hydrodynamic turbulence on $512^3$ and $1024^3$ grids. We probe regimes of forced steady-state isothermal flows and of unforced decaying ideal gas flows. The key quantity we measure is pressure dilatation cospectrum, $E^{PD}(k)$, where we provide the first numerical evidence that it decays at a rate faster than $k^{-1}$ as a function of wavenumber. This is sufficient to imply that mean pressure dilatation acts primarily at large-scales and that kinetic and internal energy budgets statistically decouple beyond a transitional scale-range. Our results suggest that an extension of Kolmogorov's inertial-range theory to compressible turbulence is possible.

In this paper, a new method is proposed to estimate the broad-line region sizes of UV lines $R^{\rm{uv}}_{\rm{BLR}}$. It is applied to 3C 273. First, we derive the time lags of radio emission relative to broad emission lines Ly$\alpha$ and C IV by the ZDCF method. The broad lines lag the 5, 8, 15, 22 and 37 GHz emission. The measured lags $\tau^{\rm{uv}}_{\rm{ob}}$ are of the order of years. For a given line, $\tau^{\rm{uv}}_{\rm{ob}}$ decreases as the radio frequency increases. This trend results from the radiative cooling of relativistic electrons. Both UV lines have a lag of $\tau^{\rm{uv}}_{\rm{ob}}=-2.74^{+0.06}_{-0.25}$ yr relative to the 37 GHz emission. These results are consistent with those derived from the Balmer lines in Paper I. Second, we derive the time lags of the lines Ly$\alpha$, CIV, H$\gamma$, H$\beta$ and H$\alpha$ relative to the 37 GHz emission by the FR/RSS Monte Carlo method. The measured lags are $\tau_{\rm{ob}}=-3.40^{+0.31}_{-0.05}$, $-3.40^{+0.41}_{-0.14}$, $-2.06^{+0.36}_{-0.92}$, $-3.40^{+1.15}_{-0.20}$ and $-3.56^{+0.35}_{-0.18}$ yr for the lines Ly$\alpha$, CIV, H$\gamma$, H$\beta$ and H$\alpha$, respectively. These estimated lags are consistent with those derived by the ZDCF method within the uncertainties. Based on the new method, we derive $R^{\rm{uv}}_{\rm{BLR}}=2.54^{+0.71}_{-0.35}$--$4.01^{+0.90}_{-1.16}$ and $2.54^{+0.80}_{-0.43}$--$4.01^{+0.98}_{-1.24}$ light-years for the Ly$\alpha$ and CIV lines, respectively. Considering the uncertainties, these estimated sizes are consistent with those obtained in the classical reverberation mapping for the UV lines and the Balmer lines. This indicates that their emitting regions are not separated so large as in the classical mapping of the UV and optical lines. These results seem to depart from the stratified ionization structures obtained in the classical mapping.

Theory and simulations suggest that magnetic fields from radio jets and lobes powered by their central super massive black holes can be an important source of magnetic fields in the galaxy clusters. This is paper II in a series of studies where we present self-consistent high-resolution adaptive mesh refinement cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial masses ranging from 1 $\times$ 10$^{14}$ to 2 $\times$ 10$^{15}$ M$_{\odot}$. In this work we examine the effects of the mass and merger history on the final magnetic properties. We find that the evolution of magnetic fields is qualitatively similar to those of previous studies. In most clusters, the injected magnetic fields can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process with hierarchical mergers, while the amplification history and the magnetic field distribution depend on the cluster formation and magnetism history. This can be very different for different clusters. The total magnetic energies in these clusters are between 4 $\times$ 10$^{57}$ and $10^{61}$ erg, which is mainly decided by the cluster mass, scaling approximately with the square of the total mass. Dynamically older relaxed clusters usually have more magnetic fields in their ICM. The dynamically very young clusters may be magnetized weakly since there is not enough time for magnetic fields to be amplified.

We probe the cosmic acceleration by using the recently released SNLS3 sample of 472 type Ia supernovae. Combining this type Ia supernovae dataset with the cosmic microwave background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7-yr observations, the baryon acoustic oscillation results from the Sloan Digital Sky Survey data release 7, and the Hubble constant measurement from the Wide Field Camera 3 on the Hubble Space Telescope, we measure the dark energy equation of state $w$ and the deceleration parameter $q$ as functions of redshift by using the Chevallier-Polarski-Linder parametrization. Our result is consistent with a cosmological constant at 1$\sigma$ confidence level, without evidence for the recent slowing down of the cosmic acceleration. Furthermore, we consider three binned parametrizations ($w$ is piecewise constant in redshift $z$) based on different binning methods. The similar results are obtained, i.e., the $\Lambda$CDM model is still nicely compatible with current observations.

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

Interior to the gaseous envelopes of Saturn, Uranus, and Neptune, there are high-density cores with masses larger than 10 Earth masses. According to the conventional sequential accretion hypothesis, such massive cores are needed for the onset of efficient accretion of their gaseous envelopes. However, Jupiter's gaseous envelope is more massive and core may be less massive than those of Saturn. In order to account for this structural diversity and the super-solar metallicity in the envelope of Jupiter and Saturn, we investigate the possibility that they may have either merged with other gas giants or consumed several Earth-mass proto-planetary embryos during or after the rapid accretion of their envelope. In general, impinging sub-Earth-mass planetesimals disintegrate in gas giants' envelopes deposit heavy elements well outside the cores and locally suppress the convection. Consequently, their fragments sediment to promote the growth of cores. Through a series of numerical simulations, we show that it is possible for colliding super-Earth-mass embryos to reach the cores of gas giants. Direct parabolic collisions also lead to the coalescence of gas giants and merging of their cores. In these cases, the energy released from the impact leads to vigorous convective motion throughout the envelope and the erosion of the cores. This dichotomy contributes to the observed dispersion in the internal structure and atmospheric composition between Jupiter and Saturn and other gas giant planets and elsewhere.

The cosmological model of dark energy interacting with cold dark matter without coupling to the baryonic matter, is studied in the background of both classical Einstein and loop quantum cosmology. We consider two types of interacting models. In the former model, the interaction is a linear combination of the densities of two dark sectors, while in the latter model, the interaction with a constant transfer rate depends only on the density of cold dark matter. It is shown that the dynamical results in loop quantum cosmology are different from those in classical Einstein cosmology for both two kinds of interacting models. Moreover, the form of the interaction affects significantly the dynamical results in both kinds of cosmology.

The metric $f(R)$ theories of gravity are generalized to five-dimensional spacetimes. By assuming a hypersurface-orthogonal Killing vector field representing the compact fifth dimension, the five-dimensional theories are reduced to their four-dimensional formalism. Then we study the cosmology of a special class of $f(R)=\alpha R^m$ models in a spatially flat FRW spacetime. It is shown that the parameter $m$ can be constrained to a certain range by the current observed deceleration parameter, and its lower bound corresponds to the Kaluza-Klein theory. It turns out that both expansion and contraction of the extra dimension may prescribe the smooth transition from the deceleration era to the acceleration era in the recent past as well as an accelerated scenario for the present universe. Hence five-dimensional $f(R)$ gravity can naturally account for the present accelerated expansion of the universe. Moreover, the models predict a transition from acceleration to deceleration in the future, followed by a cosmic recollapse within finite time. This differs from the prediction of the five-dimensional Brans-Dicke theory but is in consistent with a recent prediction based on loop quantum cosmology.

There is a significant nonlinear correlation between the Eddington ratio (L_bol=L_Edd) and the Eddington-scaled kinetic power (L_kin=L_Edd) of jets in low luminosity active galactic nuclei (AGNs) (Merloni & Heinz). It is believed that these low luminosity AGNs contain advection dominated accretion flows (ADAFs). We adopt the ADAF model developed by Li & Cao, in which the global dynamics of ADAFs with magnetically driven outflows is derived numerically, to explore the relation between bolometric luminosity and kinetic power of jets. We find that the observed relation, L_kin/L_Edd ~ (L_bol=L_Edd)^0.49, can be well reproduced by the model calculations with reasonable parameters for ADAFs with magnetically driven outflows. Our model calculations is always consistent with the slope of the correlation independent of the values of the parameters adopted. Compared with the observations, our results show that over 60% of the accreted gas at the outer radius escapes from the accretion disc in a wind before the gas falls into the black holes. The observed correlation between Eddington-scaled kinetic power and Bondi power can also be qualitatively reproduced by our model calculations. Our results show that the mechanical efficiency varies from 10^-2 ~ 10^-3, which is roughly consistent with that required in AGN feedback simulations.

The class of exotic Jupiter-mass planets that orbit very close to their parent stars were not explicitly expected before their discovery. The recently found transiting planet WASP-12b has a mass Mp = 1.4(+/-0.1) Jupiter masses (MJ), a mean orbital distance of only 3.1 stellar radii (meaning it is subject to intense tidal forces), and a period of 1.1 days. Its radius 1.79(+/- 0.09) RJ is unexpectedly large and its orbital eccentricity 0.049(+/-0:015) is even more surprising as such close orbits are in general quickly circularized. Here we report an analysis of its properties, which reveals that the planet is losing mass to its host star at a rate ~ 10^-7 MJ yr^-1. The planets surface is distorted by the stars gravity and the light curve produced by its prolate shape will differ by about ten per cent from that of a spherical planet. We conclude that dissipation of the stars tidal perturbation in the planets convective envelope provides the energy source for its large volume. We predict up to 10mJy CO band-head (2.292 micron) emission from a tenuous disk around the host star, made up of tidally stripped planetary gas. It may also contain a detectable resonant super-Earth, as a hypothetical perturber that continually stirs up WASP-12b's eccentricity.

The Palatini $f(R)$ gravity, is able to probably explain the late time cosmic acceleration without the need for dark energy, is studied. In this paper, we investigate a number of $f(R)$ gravity theories in Palatini formalism by means of statefinder diagnosis. We consider two types of $f(R)$ theories: (i) $f(R)=R+\alpha R^{m}-\beta R^{-n}$ and (ii) $f(R)=R+\alpha ln R+\beta$. We find that the evolutionary trajectories in the $s-r$ and $q-r$ planes for various types of the Palatini $f(R)$ theories reveal different evolutionary properties of the universe. Additionally, we use the observational $H(z)$ data to constrain models of $f(R)$ gravity.

We carry out 2-D high resolution numerical simulations of type I planet migration with different disk viscosities. We find that the planet migration is strongly dependent on disk viscosities. Two kinds of density wave damping mechanisms are discussed. Accordingly, the angular momentum transport can be either viscosity dominated or shock dominated, depending on the disk viscosities. The long term migration behavior is different as well. Influences of the Rossby vortex instability on planet migration are also discussed. In addition, we investigate very weak shock generation in inviscid disks by small mass planets and compare the results with prior analytic results.

Supermassive black hole binaries (SMBHBs) are products of galaxy mergers, and are important in testing Lambda cold dark matter cosmology and locating gravitational-wave-radiation sources. A unique electromagnetic signature of SMBHBs in galactic nuclei is essential in identifying the binaries in observations from the IR band through optical to X-ray. Recently, the flares in optical, UV, and X-ray caused by supermassive black holes (SMBHs) tidally disrupting nearby stars have been successfully used to observationally probe single SMBHs in normal galaxies. In this Letter, we investigate the accretion of the gaseous debris of a tidally disrupted star by a SMBHB. Using both stability analysis of three-body systems and numerical scattering experiments, we show that the accretion of stellar debris gas, which initially decays with time $\propto t^{-5/3}$, would stop at a time $T_{\rm tr} \simeq \eta T_{\rm b}$. Here, $\eta \sim0.25$ and $T_{\rm b}$ is the orbital period of the SMBHB. After a period of interruption, the accretion recurs discretely at time $T_{\rm r} \simeq \xi T_b$, where $\xi \sim 1$. Both $\eta$ and $\xi$ sensitively depend on the orbital parameters of the tidally disrupted star at the tidal radius and the orbit eccentricity of SMBHB. The interrupted accretion of the stellar debris gas gives rise to an interrupted tidal flare, which could be used to identify SMBHBs in non-active galaxies in the upcoming transient surveys.

We study the global dynamics of advection-dominated accretion flows (ADAFs) with magnetically driven outflows. A fraction of gases in the accretion flow is accelerated into the outflows, which leads to decreasing of the mass accretion rate in the accretion flow towards the black hole. We find that the r-dependent mass accretion rate is close to a power-law one, m_dot r^s, as assumed in the advection-dominated inflow-outflow solution (ADIOS), in the outer region of the ADAF, while it deviates significantly from the power-law r-dependent accretion rate in the inner region of the ADAF. It is found that the structure of the ADAF is significantly changed in the presence of the outflows. The temperatures of the ions and electrons in the ADAF decreases in the presence of outflows, as a fraction of gravitational power released in the ADAF is tapped to accelerate the outflows.

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($\mu$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

We describe 2D hydrodynamic simulations of the migration of low-mass planets ($\leq 30 M_{\oplus}$) in nearly laminar disks (viscosity parameter $\alpha < 10^{-3}$) over timescales of several thousand orbit periods. We consider disk masses of 1, 2, and 5 times the minimum mass solar nebula, disk thickness parameters of $H/r = 0.035$ and 0.05, and a variety of $\alpha$ values and planet masses. Disk self-gravity is fully included. Previous analytic work has suggested that Type I planet migration can be halted in disks of sufficiently low turbulent viscosity, for $\alpha \sim 10^{-4}$. The halting is due to a feedback effect of breaking density waves that results in a slight mass redistribution and consequently an increased outward torque contribution. The simulations confirm the existence of a critical mass ($M_{cr} \sim 10 M_{\oplus}$) beyond which migration halts in nearly laminar disks. For $\alpha \ga 10^{-3}$, density feedback effects are washed out and Type I migration persists. The critical masses are in good agreement with the analytic model of Rafikov (2002). In addition, for $\alpha \la 10^{-4}$ steep density gradients produce a vortex instability, resulting in a small time-varying eccentricity in the planet's orbit and a slight outward migration. Migration in nearly laminar disks may be sufficiently slow to reconcile the timescales of migration theory with those of giant planet formation in the core accretion model.

Disk self-gravity could play an important role in the dynamic evolution of interaction between disks and embedded protoplanets. We have developed a fast and accurate solver to calculate the disk potential and disk self-gravity forces for disk systems on a uniform polar grid. Our method follows closely the method given by Chan et al. (2006), in which an FFT in the azimuthal direction is performed and a direct integral approach in the frequency domain in the radial direction is implemented on a uniform polar grid. This method can be very effective for disks with vertical structures that depend only on the disk radius, achieving the same computational efficiency as for zero-thickness disks. We describe how to parallelize the solver efficiently on distributed parallel computers. We propose a mode-cutoff procedure to reduce the parallel communication cost and achieve nearly linear scalability for a large number of processors. For comparison, we have also developed a particle-based fast tree-code to calculate the self-gravity of the disk system with vertical structure. The numerical results show that our direct integral method is at least two order of magnitudes faster than our optimized tree-code approach.

We have performed nonlinear ideal magnetohydrodynamic simulations of the long term evolution of a magnetized low-density "bubble" plasma formed by a radio galaxy in a stratified cluster medium. It is found that about 3.5% of the initial magnetic energy remains in the bubble after $\sim 8 \times 10^{9}$ years, and the initial magnetic bubble expansion is adiabatic. The bubble can survive for at least $8 \times 10^9$ years due to the stabilizing effect of the bubble magnetic field on Rayleigh-Taylor and Kelvin-Holmholtz instabilities, possibly accounting for "ghost cavities" as observed in Perseus-A\@. A filament structure spanning about 500 kpc is formed along the path of bubble motion. The mean value of the magnetic field inside this structure is $\sim 0.88$ $\mu$G at $\sim8\times10^9$ years. Finally, the initial bubble momentum and rotation have limited influence on the long term evolution of the bubble.

We present cosmological magnetohydrodynamic simulations of the formation of a galaxy cluster with magnetic energy feedback from an active galactic nuclei (AGN). We demonstrate that X-ray cavities can be produced by the magnetically dominated jet-lobe system that is supported by a central axial current. The cavities are magnetically dominated and their morphology is determined jointedly by the magnetic fields and the background cluster pressure profile. The expansion and motion of the cavities are driven initially by the Lorentz force of the magnetic fields, and the cavities only become buoyant at late stages ($> 500$ Myr). We find that up to $80%-90%$ of the injected magnetic energy goes into doing work against the hot cluster medium, heating it, and lifting it in the cluster potential.

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magnetohydrodynamical simulations of Population III star formation including the Biermann Battery effect. We find that the Population III stars formed including this effect are both qualitatively and quantitatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe peak magnetic fields of ~10^-9 G in the center of our star-forming halo at z ~ 17.55. The magnetic fields created by the Biermann Battery effect are predominantly formed early in the evolution of the primordial halo at low density and large spatial scales, and then grow through compression and by shear flows. The fields seen in this calculation are never large enough to be dynamically important (with beta >= 10^{15} at all times), and should be considered the minimum possible fields in existence during Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational instability at higher densities and smaller spatial scales.

We apply magnetohydrodynamic (MHD) modeling to the radio galaxy Hercules A for investigating the jet-driven shock, jet/lobe transition, wiggling, and magnetic field distribution associated with this source. The model consists of magnetic tower jets in a galaxy cluster environment, which has been discussed in a series of our papers. The profile of underlying ambient gas plays an important role in jet-lobe morphology. The balance between the magnetic pressure generated by axial current and the ambient gas pressure can determine the lobe radius. The jet body is confined jointly by the external pressure and gravity inside the cluster core radius R_c, while outside R_c it expands radially to form fat lobes in a steeply decreasing ambient thermal pressure gradient. The current-carrying jets are responsible for generating a strong, tightly wound helical magnetic field. This magnetic configuration will be unstable against the current-driven kink mode and it visibly grows beyond R_c where a separation between the jet forward and return currents occurs. The reversed pinch profile of global magnetic field associated with the jet and lobes produces projected B-vector distributions aligned with the jet flow and the lobe edge. AGN-driven shock powered by the expanding magnetic tower jet surrounds the jet/lobe structure and heats the ambient ICM. The lobes expand subsonically; no obvious hot spots are produced at the heads of lobes. Several key features in our MHD modeling may be qualitatively supported by the observations of Hercules A.

In this paper, we study a curvaton model described by a Dirac-Born-Infeld-like action. We calculate the spectrum of curvature perturbation induced by DBI-curvaton and estimate its non-Gaussianity. We find that in the limit of low sound speed the amplitude of curvature perturbation is enhanced more than that in DBI inflation. This result also means that the inflationary scale with DBI-curvaton may be lower than that in usual curvaton model. In addition, we also find that the local non-Gaussianity level in DBI-curvaton is generally about 10 times larger than that in usual curvaton model, which is interesting for current observations. This work in some sense explores a new application of brane-cosmology in inflation.

Some different correlations between optical-UV variability and other quasar properties, such as luminosity, black hole mass and rest-frame wavelength, were discovered. The positive correlation between optical-UV variability amplitude and black hole mass was first found by Wold et al., and this was confirmed by Wilhite et al. We suggest that the accretion disk model can explain these correlations, provided the optical-UV variability is triggered by the change of accretion rate. The disk temperature of accretion discs decreases with increasing black hole mass, which leads to systematical spectral shape difference with black hole mass even if the black hole is accreting at the same rate m_dot (m_dot = M_dot / M_dotEdd). The observed positive correlation between optical-UV variability and black hole mass can be well reproduced by our model calculations, if the mean accretion rate m_dot0 ~ 0.1 with variation of m_delta ~ 0.4 - 0.5 m_dot0. We also found that the observed correlations of optical-UV variability amplitude with luminosity or rest-frame wavelength can be qualitatively explained by this accretion disc model.

Nonlinear ideal magnetohydrodynamic (MHD) simulations of the propagation and expansion of a magnetic "bubble" plasma into a lower density, weakly-magnetized background plasma are presented. These simulations mimic the geometry and parameters of the Plasma Bubble Expansion Experiment (PBEX) [A. G. Lynn, Y. Zhang, S. C. Hsu, H. Li, W. Liu, M. Gilmore, and C. Watts, Bull. Amer. Phys. Soc. {\bf 52}, 53 (2007)], which is studying magnetic bubble expansion as a model for extra-galactic radio lobes. The simulations predict several key features of the bubble evolution. First, the direction of bubble expansion depends on the ratio of the bubble toroidal to poloidal magnetic field, with a higher ratio leading to expansion predominantly in the direction of propagation and a lower ratio leading to expansion predominantly normal to the direction of propagation. Second, an MHD shock and a trailing slow-mode compressible MHD wavefront are formed ahead of the bubble as it propagates into the background plasma. Third, the bubble expansion and propagation develop asymmetries about its propagation axis due to reconnection facilitated by numerical resistivity and to inhomogeneous angular momentum transport mainly due to the background magnetic field. These results will help guide the initial experiments and diagnostic measurements on PBEX.

Magnetic fields play an important role in almost all astrophysical phenomena including star formation. But due to the difficulty in analytic modeling and observation, magnetic fields are still poorly studied and numerical simulation has become a major tool. We have implemented a cosmological magnetohydrodynamics package for Enzo which is an AMR hydrodynamics code designed to simulate structure formation. We use the TVD solver developed by S. Li as the base solver. In addition, we employ the constrained transport(CT) algorithm as described by D. Balsara. For interpolation magnetic fields to fine grids we used a divergence free quadratic reconstruction, also described by Balsara. We present results from several test problems including MHD caustics, MHD pancake and galaxy cluster formation with magnetic fields. We also discuss possible applications of our AMR MHD code to first star research.

In the present paper, we study the non-Gaussianity of multiple-field inflation model using the method of the field equation. We start from reviewing the background and the perturbation theory of multiple-field inflation, and then derive the Klein-Gorden equation for the perturbations at second order. Afterward, we calculate the tree-level bispectrum of the fields' perturbations and finally give the corresponding parameter $f_{NL}$ for the curvature perturbation $\zeta$ in virtue of the $\delta N$ formalism. We also compare our result with the one already obtained from the Lagrangian formalism, and find they are consistent. This work may help us understand perturbation theory of inflation more deeply.

(abridged) Search for planets around main-sequence (MS) stars more massive than the Sun is hindered by their hot and rapidly spinning atmospheres. This obstacle has been sidestepped by radial-velocity surveys of those stars on their post-MS evolutionary track (G sub-giant and giant stars). Preliminary observational findings suggest a deficiency of short-period hot Jupiters around the observed post MS stars, although the total fraction of them with known planets appears to increase with their mass. Here we consider the possibility that some very close- in gas giants or a population of rocky planets may have either undergone orbital decay or been engulfed by the expanding envelope of their intermediate-mass host stars. If such events occur during or shortly after those stars' main sequence evolution when their convection zone remains relatively shallow, their surface metallicity can be significantly enhanced by the consumption of one or more gas giants. We show that stars with enriched veneer and lower-metallicity interior follow slightly modified evolution tracks as those with the same high surface and interior metallicity. As an example, we consider HD149026, a marginal post MS 1.3 Msun star. We suggest that its observed high (nearly twice solar) metallicity may be confined to the surface layer as a consequence of pollution by the accretion of either a planet similar to its known 2.7-day-period Saturn-mass planet, which has a 70 Mearth compact core, or a population of smaller mass planets with a comparable total amount of heavy elements. It is shown that an enhancement in surface metallicity leads to a reduction in effective temperature, in increase in radius and a net decrease in luminosity. The effects of such an enhancement are not negligible in the determinations of the planet's radius based on the transit light curves.

We argue the behavior of Poynting flux-dominated outflows from AGN in the galactic cluster systems by performing three-dimensional MHD simulations within the framework of the "magnetic tower" model. Of particular interests are the structure of MHD waves, the cylindrical radial force balance, the (de)collimation, and the stability properties of magnetic tower jets. Transition between the jet/lobe and the formation of wiggling jet by growing current-driven instability are discussed.

We present a numerical method for spatially 1.5-dimensional and time-dependent studies of accretion disks around black holes, that is originated from a combination of the standard pseudo-spectral method and the adaptive domain decomposition method existing in the literature, but with a number of improvements in both the numerical and physical senses. In particular, we introduce a new treatment for the connection at the interfaces of decomposed subdomains, construct an adaptive function for the mapping between the Chebyshev-Gauss-Lobatto collocation points and the physical collocation points in each subdomain, and modify the over-simplified 1-dimensional basic equations of accretion flows to account for the effects of viscous stresses in both the azimuthal and radial directions. Our method is verified by reproducing the best results obtained previously by Szuszkiewicz & Miller on the limit-cycle behavior of thermally unstable accretion disks with moderate viscosity. A new finding is that, according to our computations, the Bernoulli function of the matter in such disks is always and everywhere negative, so that outflows are unlikely to originate from these disks. We are encouraged to study the more difficult case of thermally unstable accretion disks with strong viscosity, and wish to report our results in a subsequent paper.

Extrasolar planets found with radial velocity surveys have masses ranging from several Earth to several Jupiter masses. While mass accretion onto protoplanetary cores in weak-line T-Tauri disks may eventually be quenched by a global depletion of gas, such a mechanism is unlikely to have stalled the growth of some known planetary systems which contain relatively low-mass and close-in planets along with more massive and longer period companions. Here, we suggest a potential solution for this conundrum. In general, supersonic infall of surrounding gas onto a protoplanet is only possible interior to both of its Bondi and Roche radii. At a critical mass, a protoplanet's Bondi and Roche radii are equal to the disk thickness. Above this mass, the protoplanets' tidal perturbation induces the formation of a gap. Although the disk gas may continue to diffuse into the gap, the azimuthal flux across the protoplanets' Roche lobe is quenched. Using two different schemes, we present the results of numerical simulations and analysis to show that the accretion rate increases rapidly with the ratio of the protoplanet's Roche to Bondi radii or equivalently to the disk thickness. In regions with low geometric aspect ratios, gas accretion is quenched with relatively low protoplanetary masses. This effect is important for determining the gas-giant planets' mass function, the distribution of their masses within multiple planet systems around solar type stars, and for suppressing the emergence of gas-giants around low mass stars.

In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy flow that affects the baryon dynamics. We have developed an Eulerian-based cosmological magnetohydrodynamics code (CosmoMHD) with modern shock capturing schemes to study the formation and evolution of cosmic structures in the presence of magnetic fields. The code solves the ideal MHD equations as well as the non-equilibrium rate equations for multiple species, the Vlasov equation for dynamics of collisionless particles, the Poisson's equation for the gravitational potential field and the equation for the evolution of the intergalactic ionizing radiation field. In addition, a detailed star formation prescription and feedback processes are implemented. Several methods for solving the MHD by high-resolution schemes with finite-volume and finite-difference methods are implemented. The divergence-free condition of the magnetic fields is preserved at a level of computer roundoff error via the constraint transport method. We have also implemented a high-resolution method via dual-equation formulations to track the thermal energy accurately in very high Mach number or high Alfven-Mach number regions. Several numerical tests have demonstrated the efficacy of the proposed schemes.

Stability properties of ``magnetic tower'' jets propagating in the gravitationally stratified background have been examined by performing three-dimensional magnetohydrodynamic simulations. The current-carrying, Poynting flux-dominated magnetic tower jet, which possesses a highly wound helical magnetic field, is subject to the current-driven instability (CDI). We find that, under general physical conditions including small perturbations in the initial background profiles, the propagating magnetic tower jets develop the non-axisymmetric, $m=1$ kink mode of the CDI. The kink mode grows on the local Alfv\'en crossing time scale. In addition, two types of kink modes appear in the system. At the central region where external thermal pressure confinement is strong, only the internal kink mode is excited and will grow. A large distance away from the central region where the external thermal pressure becomes low, the external kink mode is observed. As a result, the exterior of magnetic tower jets will be deformed into a large-scale wiggled structure. We also discuss extensively the different physical processes that contribute to the overall stability properties of the magnetic tower jets. Specifically, when the jet propagates in an initially unperturbed background, we find that they can survive the kink mode beyond the point predicted by the well-known Kruskal-Shafranov (K-S) criterion. The stabilization in this case comes mainly from the dynamical relaxation of magnetic twists during the propagation of magnetic towers; the magnetic pitch is reduced and the corresponding K-S critical wavelength becomes longer as the tower jet proceeds. Furthermore, we show that the pressure-driven and Kelvin-Helmholtz instabilities do not occur in the magnetic tower jets.

Based on a new approach on modeling the magnetically dominated outflows from AGNs (Li et al. 2006), we study the propagation of magnetic tower jets in gravitationally stratified atmospheres (such as a galaxy cluster environment) in large scales ($>$ tens of kpc) by performing three-dimensional magnetohydrodynamic (MHD) simulations. We present the detailed analysis of the MHD waves, the cylindrical radial force balance, and the collimation of magnetic tower jets. As magnetic energy is injected into a small central volume over a finite amount of time, the magnetic fields expand down the background density gradient, forming a collimated jet and an expanded ``lobe'' due to the gradually decreasing background density and pressure. Both the jet and lobes are magnetically dominated. In addition, the injection and expansion produce a hydrodynamic shock wave that is moving ahead of and enclosing the magnetic tower jet. This shock can eventually break the hydrostatic equilibrium in the ambient medium and cause a global gravitational contraction. This contraction produces a strong compression at the head of the magnetic tower front and helps to collimate radially to produce a slender-shaped jet. At the outer edge of the jet, the magnetic pressure is balanced by the background (modified) gas pressure, without any significant contribution from the hoop stress. On the other hand, along the central axis of the jet, hoop stress is the dominant force in shaping the central collimation of the poloidal current. The system, which possesses a highly wound helical magnetic configuration, never quite reaches a force-free equilibrium state though the evolution becomes much slower at late stages. The simulations were performed without any initial perturbations so the overall structures of the jet remain mostly axisymmetric.

We report a detailed spectroscopic abundance analysis for a sample of 18 F-K dwarfs of the young open cluster IC 4665. Stellar parameters and element abundances of Li, O, Mg, Si, Ca, Ti, Cr, Fe and Ni have been derived using the spectroscopic synthesis tool SME (Spectroscopy Made Easy). Within the measurement uncertainties the iron abundance is uniform with a standard deviation of 0.04 dex. No correlation is found between the iron abundance and the mass of the stellar convective zone, and between the Li abundance and the Fe abundance. In other words, our results do not reveal any signature of accretion and therefore do not support the scenario that stars with planets (SWPs) acquire their on the average higher metallicity compared to field stars via accretion of metal-rich planetary material. Instead the higher metallicity of SWPs may simply reflect the fact that planet formation is more efficient in high metallicity environs. However, since that many details of the planet system formation processes remain poorly understood, further studies are needed for a final settlement of the problem of the high metallicity of SWPs. The standard deviation of [Fe/H] deduced from our observations, taken as an upper limit on the metallicity dispersion amongst the IC 4665 member stars, has been used to constrain proto-planetary disk evolution, terrestrial and giant planets formation and evolution processes. Our results do not support the possibility that the migration of gas giants and the circularization of terrestrial planets' orbits are regulated by their interaction with a residual population of planetesimals and dust particles.

We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimately lead to the excitation of a secondary instability. We also demonstrate that very high numerical resolution is required to both follow the potential vorticity changes and identify the location where the secondary instability is first excited. Low-resolution results are shown to give the wrong location. We establish the robustness of a secondary instability and its impact on the torque onto the planet. After the saturation of the instability, the disk shows large-scale non-axisymmetry, causing the torque on the planet to oscillate with large amplitude. The impact of the oscillating torque on the protoplanet's migration remains to be investigated.

We obtain global solutions of radiatively inefficiently accretion flows around black holes. Whether and where convection develops in a flow are self-consistently determined with the mixing-length theory. The solutions can be divided into three types according to the strength of normal viscosity. Type I solution corresponds to large viscosity parameter $\alpha \ga 0.1$, which is purely advection-dominated and with no convection, and has been extensively studied in the literature. Type II solution is for moderate $\alpha \sim 0.01$, which has a three-zone structure. The inner zone is advection-dominated, the middle zone is convection-dominated and ranges from a few tens to a few thousands of gravitational radii, and the outer zone is convectively stable and matches outward a Keplerian disc. The net energy flux throughout the flow is inward as in type I solution. Type III solution which is for small $\alpha \la 0.001$ consists of two zones as Abramowicz et al. suggested previously: an inner advection-dominated zone and an outer convection-dominated zone, separated at a radius of a few tens of gravitational radii. This type of solution has an outward net energy flux. In both type II and III solutions the radial density profile is between the 1/2 law of self-similar convection-dominated accretion flow model and the 3/2 law of self-similar advection-dominated accretion flow model, and the efficiency of energy release is found to be extremely low. Our results are in good agreement with those of recent numerical simulations.