We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

SN 2018hti is a Type I superluminous supernova (SLSN~I) with an absolute $g$-band magnitude of $-22.2$ at maximum brightness, discovered in a metal-poor galaxy at a redshift of 0.0612. We present extensive photometric and spectroscopic observations of this supernova, covering the phases from $\sim -35$ days to more than +340 days from the $r$-band maximum. Combining our $BVgri$-band photometry with {\it Swift} UVOT optical/ultraviolet photometry, we calculated the peak luminosity as $\sim 3.5\times10^{44}$ erg s$^{-1}$. Modeling the observed light curve reveals that the luminosity evolution of SN 2018hti can be produced by an ejecta mass of 5.8 $M_\odot$ and a magnetar with a magnetic field of $B=1.8\times10^{13}$~G having an initial spin period of $P_0=1.8$ ms. Based on such a magnetar-powered scenario and a larger sample, a correlation between the spin of the magnetar and the kinetic energy of the ejecta can be inferred for most SLSNe~I, suggesting a self-consistent scenario. Like for other SLSNe~I, the host galaxy of SN 2018hti is found to be relatively faint ($M_{g} = -17.75$ mag) and of low metallicity ($Z=0.3~Z_\odot$), with a star-formation rate of 0.3 $M_\odot$ yr$^{-1}$. According to simulation results of single-star evolution, SN 2018hti could originate from a massive, metal-poor star with a zero-age main sequence (ZAMS) mass of 25--40 $M_\odot$, or from a less massive rotating star with $M_\mathrm{ZAMS} \approx 16$--25 $M_\odot$. For the case of a binary system, its progenitor could also be a star with $M_\mathrm{ZAMS} \gtrsim 25$ $M_\odot$.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.

The detection of binary black hole coalescences by LIGO/Virgo has aroused the interest in primordial black holes (PBHs), because they could be both the progenitors of these black holes and a compelling candidate of dark matter (DM). PBHs are formed soon after the enhanced scalar perturbations re-enter horizon during radiation dominated era, which would inevitably induce gravitational waves as well. Searching for such scalar induced gravitational waves (SIGWs) provides an elegant way to probe PBHs. We perform the first direct search for the signals of SIGWs accompanying the formation of PBHs in North American Nanohertz Observatory for Gravitational waves (NANOGrav) 11-year data set. No statistically significant detection has been made, and hence we place a stringent upper limit on the abundance of PBHs at $95\%$ confidence level. In particular, less than one part in a million of the total DM mass could come from PBHs in the mass range of $[2 \times 10^{-3}, 7\times 10^{-1}] \Msun$.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

We present a study on the radial profiles of D4000,luminosity-weighted stellar ages $\tau_L$,and luminosity-weighted stellar metallicities $[Z/H]_L$ of 3654 nearby galaxies($0.01<z<0.15$)using the IFU spectroscopic data from the MaNGA survey available in the SDSS DR15,in an effort to explore the connection between median stellar population radial gradients($\nabla$D4000,$\nabla\tau_L,\nabla[Z/H]_L$)out to~$1.5R_e$ and various galaxy properties,including stellar mass($M_\star$),specific star formation rate(sSFR),morphologies,and local environment. We find that $M_\star$ is the single most predictive physical property for$\nabla$D4000 and$\nabla[Z/H]_L$. The most predictive properties for $\nabla\tau_L$ are sSFR,and to a lesser degree,$M_\star$. The environmental parameters,including local galaxy overdensities and central-satellite division,have virtually no correlation with stellar population radial profiles for the whole sample,but the $\nabla$D4000 of star-forming satellite galaxies with$M_\star\lesssim 10^{10}M_\odot$exhibit a significant positive correlation with galaxy overdensities. Galaxies with lower sSFR have on average steeper negative stellar population gradients,and this sSFR dependence is stronger for more massive star-forming galaxies. The negative correlation between the median stellar population gradients and$M_\star$ are best described largely as segmented relationships, whereby median gradients of galaxies with$\log M_\star\lesssim 10$(with the exact value depending on sSFR)have much weaker mass dependence than galaxies with higher$M_\star$. While the dependence of the radial gradients of ages and metallicities on T-Types and central stellar mass surface densities are generally not significant,galaxies with later T-Types or lower central mass densities tend to have significantly lower D4000,younger$\tau_L$ and lower$[Z/H]_L$ across the radial ranges probed in this study.

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50" and a grid spacing of 30". The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

In this paper, we study the motion of test particles around two exact charged black-hole solutions in Einstein-{\AE}ther theory. Specifically, we first consider the quasi-periodic oscillations (QPOs) and their resonances generated by the particle moving in the Einstein-{\AE}ther black hole and then turn to study the periodic orbits of the massive particles. For QPOs, we drop the usually adopted assumptions $\nu_U=\nu_\theta$, $\nu_L=\nu_r$, and $\nu_U/\nu_L=3/2$ with $\nu_U$ ($\nu_L$) and $\nu_r$ ($\nu_\theta$) being the upper (lower) frequency of QPOs and radial (vertical) epicyclic frequency of the orbiting particles, respectively. Instead, we put-forward a new working ansatz for which the Keplerian radius is much closer to that of the innermost stable circular orbit and explore in details the effects of the {\ae}ther field on the frequencies of QPOs. We then realize good curves for the frequencies of QPOs, which fit to data of three microquasars very well by ignoring any effects of rotation and magnetic fields. For the study of the periodic orbits of massive particle moving in the two Einstein-{\AE}ther black hole solutions, we first consider the geodesic equations in the two black holes and solve the geodesic equation for one of black hole (the black hole with \ae{}ther parameter $c_{14}=0$ but $c_{123}\neq 0$) analytically. The innermost stable circular orbits are also analyzed and we find the isco radius increases with increasing $c_{13}$ for the first type black hole while decreases with increasing $c_{14}$ for the second one. We also obtain several periodic orbits and find that they share similar taxonomy schemes as the periodic equatorial orbits in the Schwarzschild/Kerr metrics. The equations for null geodesics are also briefly considered, where we study circular photon orbits and bending angles for gravitational lensing.

Gravitational accretion accumulates the original mass, and this process is crucial for us to understand the initial phases of star formation. Using the specific infall profiles in optically thick and thin lines, we searched the clumps with infall motion from the Milky Way Imaging Scroll Painting (MWISP) CO data in previous work. In this study, we selected 133 sources of them as a sub-sample for further research and identification. The excitation temperatures of these sources are between 7.0 and 38.5 K, while the H_2 column densities are between 10^21 and 10^23 cm^-2. We have observed optically thick lines HCO+ (1-0) and HCN (1-0) using the DLH 13.7-m telescope, and found 56 sources of them with blue profile and no red profile in these two lines, which are likely to have infall motions, with the detection rate of 42\%. It suggests that using CO data to restrict sample can effectively improve the infall detection rate. Among these confirmed infall sources, there are 43 associated with Class 0/I young stellar objects (YSOs), and 13 are not. These 13 sources are probably associated with the sources in earlier evolutionary stage. By comparison, the confirmed sources which are associated with Class 0/I YSOs have higher excitation temperatures and column densities, while the other sources are colder and have lower column densities. Most infall velocities of the sources we confirmed are between 10^-1 to 10^0 km s^-1, which is consistent with previous studies.

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

We present a pipeline for geomorphological analysis that uses structure from motion (SfM) and deep learning on close-range aerial imagery to estimate spatial distributions of rock traits (diameter, size, and orientation) along a tectonic fault scarp. Unpiloted aircraft systems (UAS) have enabled acquisition of high-resolution imagery at close range, revolutionizing domains such as infrastructure inspection, precision agriculture, and disaster response. Our pipeline leverages UAS-based imagery to help scientists gain a better understanding of tectonic surface processes. We start by using SfM on aerial imagery to produce georeferenced orthomosaics and digital elevation models (DEM), then a human expert annotates rocks on a set of image tiles sampled from the orthomosaics. These annotations are used to train a deep neural network to detect and segment individual rocks in the whole site. This pipeline automatically extracts semantic information (rock boundaries) on large volumes of unlabeled, high-resolution aerial imagery, which allows subsequent structural analysis and shape descriptors to result in estimates of rock diameter, size and orientation. We present results of two experiments conducted along a fault scarp in the Volcanic Tablelands near Bishop, California. We conducted the first experiment with a hexrotor and a multispectral camera to produce a DEM and five spectral orthomosaics in red, green, blue, red edge (RE), and near infrared (NIR). We then trained deep neural networks with different input channel combinations to study the most effective learning method for inference. In the second experiment, we deployed a DJI Phantom 4 Pro equipped with an RGB camera, and focused on the spatial difference of rock-trait histograms in a larger area. Although presented in the context of geology, our pipeline can be extended to a variety of geomorphological analysis tasks in other domains.

The origin of chemically peculiar stars and non-zero eccentricity in evolved close binaries have been long-standing problems in binary stellar evolution. Answers to these questions may trace back to an intense mass transfer during AGB binary phase. In this work, we use AstroBEAR to solve the 3D radiation-hydrodynamic equations and calculate the mass transfer rate in asymptotic-giant-branch (AGB) binaries that undergo the wind-Roche-lobe-overflow or Bondi-Hoyle-Lyttleton (BHL) accretion. MESA produces the density and temperature of the boundary condition of the AGB star. To improve the resolution of the dynamics of a circumbinary disk, we implement an azimuthal angle-dependent 3D radiation transfer. We consider optically thin cooling and obtain the number density of the coolants by solving Saha equations. One of the goals of this work is to illustrate the transition from the wind-Roche-lobe-overflow to BHL accretion. Both circumbinary disks and spiral structure outflows can appear in the simulations. Circumbinary disks may form when the optical thickness in the equatorial region increases. The increase of the optical thickness is due to the deflected wind. The resulting mass transfer efficiency in our models is up to a factor of 8 times higher than what the standard BHL accretion scenario predicts, and the outflow gains up to 91% of its initial angular momentum when it reaches 1.3 binary separations. Consequently, some AGB binaries may undergo orbit shrinkage, and some will expand. The high mass transfer efficiency is closely related to the presence of the circumbinary disks.

We use data from galaxy-galaxy weak lensing to perform a novel test on General Relativity (GR). In particular, we impose strong constraints using the torsional (teleparallel) formulation of gravity in which the deviation from GR is quantified by a single parameter $\alpha$, an approximation which is always valid at low-redshift Universe and weak gravitational fields. We calculate the difference in the deflection angle and eventually derive the modified Excess Surface Density profile, which is mainly affected at small scales. Hence, confronting the predictions with weak lensing data from the Sloan Digital Sky Survey Data Release 7 we obtain the upper bound on the deviation parameter, which, expressed via the dimensionless percentage in the universe energy content, reads as ${\rm{log_{10}}}\Omega_\alpha\le -18.52_{-0.42}^{+0.80}$({\it stat.})$_{-0.39}^{+1.50}$({\it sys.}) with systematics mainly arising from the modelling. To our knowledge, this is the first time that GR is verified at such an accuracy at the corresponding scales.

We analytically calculate the scalar induced gravitational waves (SIGWs) and find a log-dependent slope of SIGW in the infrared regions $(f<f_c)$, namely $n_{\mathrm{GW}}(f)=3-2/\ln(f_c/f)$, and $n_{\mathrm{GW}}(f)=2-2/\ln(f_c/f)$ near the peak if the power spectrum of scalar curvature perturbation is quite narrow, where $f_c$ is roughly the frequency at the peak of SIGW. Such a log-dependent slope can be taken as a new template for distinguishing SIGW from other sources.

We apply the inverse-Gaussianization method proposed in \citealt{arXiv:1607.05007} to fast produce weak lensing convergence maps and investigate the peak statistics, including the peak height counts and peak steepness counts, in these mocks. We find that the distribution of peak height and steepness is in good agreement with the simulation. The difference is $\lesssim 20\%$ for these peak statistics in the maps at source redshift $z_s=1$. Besides, the loss of off-diagonal elements in peak covariance motivates us to consider the super sample variance in weak lensing peak statistics. We propose correction methods to effectively recover the (anti-)correlation among different bins by adding scatters in the mean value of these mocks. Finally, as an example of the application, we adopt the improved inverse-Gaussianization method with the above improvement to fast generate 40,000 mocks to calculate precision matrices between the power spectrum and peak statistics.

In this letter we calculate the scalar induced gravitational waves (SIGWs) accompanying the formation of primordial black hole during the radiation dominated era in three different gauges, i.e. synchronous gauge, Newton gauge and uniform curvature gauge, and we find that the energy density spectra of SIGWs, $\ogw(k)$, are identical in these three different gauges.

Precise measurements of the energy spectra of cosmic rays (CRs) show various kinds of features deviating from single power-laws, which give very interesting and important implications on their origin and propagation. Previous measurements from a few balloon and space experiments indicate the existence of spectral softenings around 10 TV for protons (and probably also for Helium nuclei). Very recently, the DArk Matter Particle Explorer (DAMPE) measurement about the proton spectrum clearly reveals such a softening with a high significance. Here we study the implications of these new measurements, as well as the groundbased indirect measurements, on the origin of CRs. We find that a single component of CRs fails to fit the spectral softening and the air shower experiment data simultaneously. In the framework of multiple components, we discuss two possible scenarios, the multiple source population scenario and the background plus nearby source scenario. Both scenarios give reasonable fits to the wide-band data from TeV to 100 PeV energies. Considering the anisotropy observations, the nearby source model is favored.

We use the multiconfiguration Dirac-Hartree-Fock (MCDHF) method combined with the relativistic configuration interaction (RCI) approach (GRASP2K) to provide a consistent set of transition energies and radiative transition data for the lower $n =3$ states in all Cl-like ions of astrophysical importance, from \ion{Cr}{8} to \ion{Zn}{14}. We also provide excitation energies calculated for \mbox{Fe X} using the many-body perturbation theory (MBPT, implemented within FAC). The comparison of the present MCDHF results with MBPT and with the available experimental energies indicates that the theoretical excitation energies are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Detailed comparisons for Fe~X and Ni~XII highlight discrepancies in the experimental energies found in the literature. Several new identifications are proposed.

We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

We study the evolution of the hydrogen ionization field during the epoch of reionization (EoR) using semi-numerical simulations. By calculating the Minkowski functionals (MFs) of the 21cm brightness temperature field, which provides a quantitative description of topology of the neutral and ionized regions, we find that the reionization process can be divided into five stages, each with different topological structures, corresponding to isolated ionized regions ("bubbles"); connected ionized region ("ionized fibers") which percolate through the cosmic volume; a sponge-like configuration with intertwined neutral and ionized regions after the overlap of bubbles; connected neutral regions ("neutral fibers") when the ionized region begin to surround the remaining neutral region; and isolated neutral regions ("islands"), before the whole space is ionized except for the rare dense clumps in galaxies and minihalos and the reionization process is completed. We use the MFs and the size statistics on the ionized or neutral regions to distinguish these different stages, and find the neutral fractions at which the transition of the different stages occur. At the later stages of reionization the neutral regions are more isolated than the ionized regions, this is the motivation for the island model description of reionization. We compare the late stage topological evolution in the island model and the bubble model, showing that in the island model the neutral fibers are more easily broken into multiple pieces due to the ionizing background.

We report new infrared measurements of the supermassive black hole at the Galactic Center, Sgr A*, over a decade that was previously inaccessible at these wavelengths. This enables a variability study that addresses variability timescales that are ten times longer than earlier published studies. Sgr A* was initially detected in the near-infrared with adaptive optics observations in 2002. While earlier data exists in form of speckle imaging (1995 - 2005), Sgr A* was not detected in the initial analysis. Here, we improved our speckle holography analysis techniques. This has improved the sensitivity of the resulting speckle images by up to a factor of three. Sgr A* is now detectable in the majority of epochs covering 7 years. The brightness of Sgr A* in the speckle data has an average observed K magnitude of 16.0, which corresponds to a dereddened flux density of $3.4$ mJy. Furthermore, the flat power spectral density (PSD) of Sgr A* between $\sim$80 days and 7 years shows its uncorrelation in time beyond the proposed single power-law break of $\sim$245 minutes. We report that the brightness and its variability is consistent over 22 years. This analysis is based on simulations using Witzel et al. (2018) model to characterize infrared variability from 2006 to 2016. Finally, we note that the 2001 periapse of the extended, dusty object G1 had no apparent effect on the near-infrared emission from accretion flow onto Sgr A*. The result is consistent with G1 being a self-gravitating object rather than a disrupting gas cloud.

The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy, twice the level of the previously observed peak flux from $>13,000$ measurements over 130 nights with the VLT and Keck Telescopes. We also observe a factor of 75 change in flux over a 2-hour time span with no obvious color changes between 1.6 $\mu$m and 2.1 $\mu$m. The distribution of flux variations observed this year is also significantly different than the historical distribution. Using the most comprehensive statistical model published, the probability of a single night exhibiting peak flux levels observed this year, given historical Keck observations, is less than $0.3\%$. The probability to observe the flux levels similar to all 4 nights of data in 2019 is less than $0.05\%$. This increase in brightness and variability may indicate a period of heightened activity from Sgr A* or a change in its accretion state. It may also indicate that the current model is not sufficient to model Sgr A* at high flux levels and should be updated. Potential physical origins of Sgr A*'s unprecedented brightness may be from changes in the accretion-flow as a result of the star S0-2's closest passage to the black hole in 2018 or from a delayed reaction to the approach of the dusty object G2 in 2014. Additional multi-wavelength observations will be necessary to both monitor Sgr A* for potential state changes and to constrain the physical processes responsible for its current variability.

Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc.Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.

General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $\Upsilon$. Our result, $\Upsilon=0.88 \pm 0.17$, is consistent with General Relativity ($\Upsilon=1$) and excludes a Newtonian model ($\Upsilon=0$ ) with a statistical significance of 5 $\sigma$.

Primordial gravitational waves generated during inflation lead to the B-mode polarization in the cosmic microwave background and a stochastic gravitational wave background in the Universe. We will explore the current constraint on the tilt of primordial gravitational-wave spectrum, and forecast how the future observations can improve the current constraint.

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 possibility that primordial black holes (PBHs) represent all of the dark matter (DM) in the Universe and explain the coalescences of binary black holes detected by LIGO/Virgo has attracted a lot of attention. PBHs are generated by the enhancement of scalar perturbations which inevitably produce the induced gravitational waves (GWs). We calculate the induced GWs up to the third-order correction which not only enhances the amplitude of induced GWs, but also extends the cutoff frequency from $2k_*$ to $3k_*$. Such effects of the third-order correction lead to an around $10\%$ increase of the signal-to-noise ratio (SNR) for both LISA and pulsar timing array (PTA) observations, and significantly widen the mass range of PBHs in the stellar mass window accompanying detectable induced GWs for PTA observations including IPTA, FAST and SKA. On the other hand, the null detections of the induced GWs by LISA and PTA experiments will exclude the possibility that all of the DM is comprised of PBHs and the GW events detected by LIGO/Virgo are generated by PBHs.

We present extensive ground-based and $Hubble~Space~Telescope$ ($HST$) photometry of the highly reddened, very nearby type Ia supernova (SN Ia) 2014J in M82, covering the phases from 9 days before to about 900 days after the $B$-band maximum. SN 2014J is similar to other normal SNe Ia near the maximum light, but it shows flux excess in the $B$ band in the early nebular phase. This excess flux emission can be due to light scattering by some structures of circumstellar materials located at a few 10$^{17}$ cm, consistent with a single degenerate progenitor system or a double degenerate progenitor system with mass outflows in the final evolution or magnetically driven winds around the binary system. At t$\sim$+300 to $\sim$+500 days past the $B$-band maximum, the light curve of SN 2014J shows a faster decline relative to the $^{56}$Ni decay. Such a feature can be attributed to the significant weakening of the emission features around [Fe III] $\lambda$4700 and [Fe II] $\lambda$5200 rather than the positron escape as previously suggested. Analysis of the $HST$ images taken at t$>$600 days confirms that the luminosity of SN 2014J maintains a flat evolution at the very late phase. Fitting the late-time pseudo-bolometric light curve with radioactive decay of $^{56}$Ni, $^{57}$Ni and $^{55}$Fe isotopes, we obtain the mass ratio $^{57}$Ni/$^{56}$Ni as $0.035 \pm 0.011$, which is consistent with the corresponding value predicted from the 2D and 3D delayed-detonation models. Combined with early-time analysis, we propose that delayed-detonation through single degenerate scenario is most likely favored for SN 2014J.

Hemispherical power asymmetry has emerged as a new challenge to cosmology in early universe. While the cosmic microwave background (CMB) measurements indicated the asymmetry amplitude $A \simeq 0.07$ at the CMB scale $k_{\rm CMB}\simeq 0.0045\,{\rm Mpc}^{-1}$, the high-redshift quasar observations found no significant deviation from statistical isotropy. This conflict can be reconciled in some scale-dependent asymmetry models. We put forward a new parameterization of scale-dependent asymmetric power spectrum, inspired by a multi-speed inflation model. The 21-cm power spectrum from the epoch of reionization can be used to constrain the scale-dependent hemispherical asymmetry. We demonstrate that an optimum, multi-frequency observation by the Square Kilometre Array (SKA) Phase 2 can impose a constraint on the amplitude of the power asymmetry anomaly at the level of $\Delta A \simeq 0.2$ at $0.056 \lesssim k_{\rm 21cm} \lesssim 0.15 \,{\rm Mpc}^{-1}$. This limit may be further improved by an order of magnitude as $\Delta A \simeq 0.01$ with a cosmic variance limited experiment such as the Omniscope.

We investigate how the next generation gravitational-wave (GW) detectors, such as Einstein Telescope (ET) and Cosmic Explorer (CE), can be used to distinguish primordial black holes (PBHs) from astrophysical black holes (ABHs). Since a direct detection of sub-solar mass black holes can be taken as the smoking gun for PBHs, we figure out the detectable limits of the abundance of sub-solar mass PBHs in cold dark matter by the targeted search for sub-solar mass PBH binaries and binaries containing a sub-solar mass PBH and a super-solar mass PBH, respectively. On the other hand, according to the different redshift evolutions of merger rate for PBH binaries and ABH binaries, we forecast the detectable event rate distributions for the PBH binaries and ABH binaries by ET and CE respectively, which can serve as a method to distinguish super-solar mass PBHs from ABHs.

The advent of gravitational wave (GW) and multi-messenger astronomy has stimulated the research on the formation mechanisms of binary black holes (BBHs) observed by LIGO/Virgo. In literature, the progenitors of these BBHs could be stellar-origin black holes (sBHs) or primordial black holes (PBHs). In this paper we calculate the Stochastic Gravitational-Wave Background (SGWB) from BBHs, covering the astrophysical and primordial scenarios separately, together with the one from binary neutron stars (BNSs). Our results indicate that PBHs contribute a stronger SGWB than that from sBHs, and the total SGWB from both BBHs and BNSs has a high possibility to be detected by the future observing runs of LIGO/Virgo and LISA. On the other hand, the SGWB from BBHs and BNSs also contributes an additional source of confusion noise to LISA's total noise curve, and then weakens LISA's detection abilities. For instance, the detection of massive black hole binary (MBHB) coalescences is one of the key missions of LISA, and the largest detectable redshift of MBHB mergers can be significantly reduced.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

We present the results based on the optical $R$-band polarization observations of 280 stars distributed towards the dark globule LDN\,1225. {\it Gaia} data release 2 parallaxes along with the polarization data of $\sim$200 stars have been used to (a) constrain the distance of LDN\,1225 as 830$\pm$83~pc, (b) determine the contribution of interstellar polarization (ISP), and (c) characterize the dust properties and delineate the magnetic field (B-field) morphology of LDN\,1225. We find that B-fields are more organized and exhibit a small dispersion of 12$\degr$. Using the $^{12}$CO molecular line data from the Purple Mountain Observatory (PMO), along with the column density, dispersion in B-fields, we estimate B-field strength to be $\sim$56\,$\pm$\,10\,$\mu$G, magnetic to turbulence pressure to be $\sim$3\,$\pm$\,2, and the mass-to-magnetic flux ratio (in units of critical value) to be~$<$\,1. These results indicate the dominant role of B-fields in comparison to turbulence and gravity in rendering the cloud support. B-fields are aligned parallel to the low-density parts (traced by $^{12}$CO map) of the cloud, in contrast they are neither parallel nor perpendicular to the high-density core structures (traced by $^{13}$CO and C$^{18}$O maps). LDN\,1225 hosts two 70\,$\mu$m sources which seem to be of low-mass Class 0 sources. The total-to-selective extinction derived using optical and near-infrared photometric data is found to be anomalous ($R_{V}$~$=$~3.4), suggesting dust grain growth in LDN\,1225. Polarization efficiency of dust grains follows a power-law index of $-$0.7 inferring that optical polarimetry traces B-fields in the outer parts of the cloud.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

The Milky Way Imaging Scroll Painting (MWISP) project is an unbiased Galactic plane CO survey for mapping regions of l=-10deg to +250deg and |b|<5.2deg with the 13.7 m telescope of the Purple Mountain Observatory. The legacy survey aims to observe the 12CO, 13CO, and C18O (J=1-0) lines simultaneously with full sampling using the nine-beam Superconducting SpectroScopic Array Receiver (SSAR) system with an instantaneous bandwidth of 1 GHz. In this paper, the completed 250 deg^2 data from l=+25.8deg to +49.7deg are presented with a grid spacing of 30" and a typical rms noise level of ~0.5 K for 12CO at the channel width of 0.16 km/s and ~0.3 K for 13CO and C18O at 0.17 km/s. The high quality data with moderate resolution (50"), uniform sensitivity, and high spatial dynamic range, allow us to investigate the details of molecular clouds (MCs) traced by the three CO isotope lines. Three interesting examples are briefly investigated, including distant Galactic spiral arms traced by CO emission with VLSR<0 km/s, the bubble-like dense gas structure near the HII region W40, and the MCs distribution perpendicular to the Galactic plane.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

At intermediate and high redshifts, measurements of galaxy star-formation rates are usually based on rest-frame ultraviolet (UV) data. A correction for dust attenuation, A_UV, is needed for these measurements. This correction is typically inferred from UV spectral slopes (beta) using an equation known as "Meurer's Relation." In this paper, we study this relation at a redshift of 1.5 using images and photometric measurements in the rest-frame UV (HST) through mid-infrared (Spitzer). It is shown that massive star-forming galaxies (above 10^10 Msun) have dust corrections that are dependent on their inclination to the line-of-sight. Edge-on galaxies have higher A_UV and infrared excess (IRX=L(IR)/L(UV)) than face-on galaxies at a given beta. Interestingly, dust corrections for low-mass star-forming galaxies do not depend on inclination. This is likely because more massive galaxies have more disk-like shapes/kinematics, while low-mass galaxies are more prolate and have more disturbed kinematics. To account for an inclination-dependent dust correction, a modified Meurer's Relation is derived: A_UV=4.43+1.99 beta - 1.73 (b/a-0.67), where b/a is the galaxy axis ratio. This inclination-dependence of A_UV can be explained by a two-component model of the dust distribution inside galaxies. In such a model, the dust attenuation of edge-on galaxies has a higher contribution from a "mixture" component (dust uniformly mixed with stars in the diffuse interstellar medium), and a lower contribution from a "birth cloud" component (near-spherical dust shells surrounding young stars in H II regions) than that of face-on galaxies. The difference is caused by the larger path-lengths through disks at higher inclinations.

Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.

Up to now several gravitational-wave events from the coalescences of black hole binaries have been reported by LIGO/VIRGO, and imply that black holes should have an extended mass function. We work out the merger rate distribution of primordial-black-hole binaries with a general mass function by taking into account the torques by all primordial black holes and linear density perturbations. In the future, many more coalescences of black hole binaries are expected to be detected, and the one-dimensional and two-dimensional merger rate distributions will be crucial for reconstructing the mass function of primordial black holes.

The morphologies of planetary nebula have long been believed to be due to wind shaping processes in which a fast wind from the central star impacts a previously ejected envelope. Asymmetries assumed to exist in the slow wind envelope lead to inertial confinement shaping the resulting interacting winds flow. We present new results demonstrating the effectiveness of Common Envelope Evolution (CEE) at producing aspherical envelopes which, when impinged upon by a spherical fast stellar wind, produce highly bipolar, jet-like outflows. We have run two simple cases using the output of a single PHANTOM SPH CEE simulation. Our work uses the Adaptive Mesh Refinement code AstroBEAR to track the interaction of the fast wind and CEE ejecta allowing us to follow the morphological evolution of the outflow lobes at high resolution in 3-D. Our two models bracket low and high momentum output fast winds. We find the interaction leads to highly collimated bipolar outflows. In addition, the bipolar morphology depends on the fast wind momentum injection rate. With this dependence comes the initiation of significant symmetry breaking between the top and bottom bipolar lobes. Our simulations, though simplified, confirm the long-standing belief that CEE can plan a major role in PPN and PN shaping.

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

We assemble 207 variable quasars from the Sloan Digital Sky Survey, all with at least 3 observations, to analyze C IV narrow absorption doublets, and obtain 328 C IV narrow absorption line systems. We find that 19 out of 328 C IV narrow absorption line systems were changed by $\mid\Delta W_r^{\lambda1548}\mid\ge3\sigma_{\Delta W_r^{\lambda1548}}$ on timescales from 15.9 to 1477 days at rest-frame. Among the 19 obviously variable C IV systems, we find that (1) 14 systems have relative velocities $\upsilon_r>0.01c$ and 4 systems have $\upsilon_r>0.1c$, where $c$ is the speed of light; (2) 13 systems are accompanied by other variable C IV systems; (3) 9 systems were changed continuously during multiple observations; and (4) 1 system with $\upsilon_{\rm r}=16~862$ km/s was enhanced by $\Delta W_r^{\lambda1548}=2.7\sigma_{\Delta W_r^{\lambda1548}}$ in 0.67 day at rest-frame. The variations of absorption lines are inversely correlated with the changes in the ionizing continuum. We also find that large variations of C IV narrow absorption lines are difficult to form on a short timescale.

We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of prolate, oblate and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low m and high redshifts, and oblate at high M and low redshifts, qualitatively consistent with van der Wel et al. (2014), implying that galaxies tend to evolve from prolate to oblate. These results are consistent with the predictions from simulations (Ceverino et al. 2015, Tomassetti et al. 2016) that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy center baryon dominated. We give probabilities of a galaxy's being prolate, oblate or spheroidal as a function of its M, redshift, projected b/a and a, which can facilitate target selections of galaxies with specific shapes at hight redshifts. We also give predicted optical depths of galaxies, which are qualitatively consistent with the expected correlation that AV should be higher for edge-on disk galaxies in each log a slice at low redshift and high mass bins.

Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

Shanghai Astronomical Observatory has upgraded its DiFX cluster to 420 CPU cores and a 432-TB storage system at the end of 2014. An international network connection for the raw data transfer has also been established. The routine operations for IVS sessions including CRF, AOV, and APSG series began in early 2015. In addition to the IVS observations, the correlator is dedicated to astrophysical and astrometric programs with the Chinese VLBI Network and international joint VLBI observations. It also worked with the new-built Tianma 65-m radio telescope and successfully found fringes as high as at X/Ka and Q bands in late 2015. A more powerful platform is planned for the high data rate and massive data correlation tasks in the future.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

Because of vast volume of data being produced by today's scientific simulations and experiments, lossy data compressor allowing user-controlled loss of accuracy during the compression is a relevant solution for significantly reducing the data size. However, lossy compressor developers and users are missing a tool to explore the features of scientific datasets and understand the data alteration after compression in a systematic and reliable way. To address this gap, we have designed and implemented a generic framework called Z-checker. On the one hand, Z-checker combines a battery of data analysis components for data compression. On the other hand, Z-checker is implemented as an open-source community tool to which users and developers can contribute and add new analysis components based on their additional analysis demands. In this paper, we present a survey of existing lossy compressors. Then we describe the design framework of Z-checker, in which we integrated evaluation metrics proposed in prior work as well as other analysis tools. Specifically, for lossy compressor developers, Z-checker can be used to characterize critical properties of any dataset to improve compression strategies. For lossy compression users, Z-checker can detect the compression quality, provide various global distortion analysis comparing the original data with the decompressed data and statistical analysis of the compression error. Z-checker can perform the analysis with either coarse granularity or fine granularity, such that the users and developers can select the best-fit, adaptive compressors for different parts of the dataset. Z-checker features a visualization interface displaying all analysis results in addition to some basic views of the datasets such as time series. To the best of our knowledge, Z-checker is the first tool designed to assess lossy compression comprehensively for scientific datasets.

Using 3D radiation-hydrodynamic simulations and analytic theory, we study the orbital evolution of asymptotic-giant-branch (AGB) binary systems for various initial orbital separations and mass ratios, and thus different initial accretion modes. The time evolution of binary separations and orbital periods are calculated directly from the averaged mass loss rate, accretion rate and angular momentum loss rate. We separately consider spin-orbit synchronized and zero spin AGB cases. We find that the the angular momentum carried away by the mass loss together with the mass transfer can effectively shrink the orbit when accretion occurs via wind-Roche-lobe overflow. In contrast, the larger fraction of mass lost in Bondi-Hoyle-Lyttleton accreting systems acts to enlarge the orbit. Synchronized binaries tend to experience stronger orbital period decay in close binaries. We also find that orbital period decay is faster when we account for the nonlinear evolution of the accretion mode as the binary starts to tighten. This can increase the fraction of binaries that result in common envelope, luminous red novae, Type Ia supernovae and planetary nebulae with tight central binaries. The results also imply that planets in the the habitable zone around white dwarfs are unlikely to be found.

Because of the vast volume of data being produced by today's scientific simulations, lossy compression allowing user-controlled information loss can significantly reduce the data size and the I/O burden. However, for large-scale cosmology simulation, such as the Hardware/Hybrid Accelerated Cosmology Code (HACC), where memory overhead constraints restrict compression to only one snapshot at a time, the lossy compression ratio is extremely limited because of the fairly low spatial coherence and high irregularity of the data. In this work, we propose a pattern-matching (similarity searching) technique to optimize the prediction accuracy and compression ratio of SZ lossy compressor on the HACC data sets. We evaluate our proposed method with different configurations and compare it with state-of-the-art lossy compressors. Experiments show that our proposed optimization approach can improve the prediction accuracy and reduce the compressed size of quantization codes compared with SZ. We present several lessons useful for future research involving pattern-matching techniques for lossy compression.

We present here the results of the first part of the VLBI Ecliptic Plane Survey (VEPS) program. The goal of the program is to find all compact sources within $7.5^\circ$ of the ecliptic plane which are suitable as calibrators for anticipated phase referencing observations of spacecraft and determine their positions with accuracy at the 1.5~nrad level. We run the program in two modes: the search mode and the refining mode. In the search mode, a complete sample of all sources brighter than 50 mJy at 5 GHz listed in the Parkes-MIT-NRAO (PMN) and Green Bank 6~cm (GB6) catalogs, except those previously detected with VLBI, is observed. In the refining mode, the positions of all ecliptic plane sources, including those found in the search mode, are improved. By October 2016, thirteen 24-hr sessions that targeted all sources brighter than 100~mJy have been observed and analyzed. Among 3320 observed target sources, 555 objects have been detected. We also conducted a number of follow-up VLBI experiments in the refining mode and improved the positions of 249 ecliptic plane sources.

Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth force obtained by using 19 years of observations of the two best measured short-period stars S0-2 and S0-38 ; (ii) an upper limit on the secular advance of the argument of the periastron for the star S0-2 ; (iii) a sensitivity analysis showing that the relativistic redshift of S0-2 will be measured after its closest approach to the black hole in 2018.

We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy connection exists beyond the local universe. To this purpose, we only need a statistically representative, rather than a complete, sample of low-mass galaxies, which enables our study to $z\gtrsim1.5$. For each low-mass galaxy, we measure the projected distance ($d_{proj}$) to its nearest massive neighbor (M*$>10^{10.5} M_\odot$) within a redshift range. At a given redshift and M*, the environmental quenching effect is considered to be observed if the $d_{proj}$ distribution of QGs ($d_{proj}^Q$) is significantly skewed toward lower values than that of star-forming galaxies ($d_{proj}^{SF}$). For galaxies with $10^{8} M_\odot < M* < 10^{10} M_\odot$, such a difference between $d_{proj}^Q$ and $d_{proj}^{SF}$ is detected up to $z\sim1$. Also, about 10\% of the quenched galaxies in our sample are located between two and four virial radii ($R_{Vir}$) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, $d_{proj}^Q / R_{Vir}$, decreases with satellite M* at $M* \lesssim 10^{9.5} M_\odot$, but increases with satellite M* at $M* \gtrsim 10^{9.5} M_\odot$. This trend suggests a smooth, if any, transition of the quenching timescale around $M* \sim 10^{9.5} M_\odot$ at $0.5<z<1.0$.

This paper uses radial colour profiles to infer the distributions of dust, gas and star formation in z=0.4-1.4 star-forming main sequence galaxies. We start with the standard UVJ-based method to estimate dust extinction and specific star formation rate (sSFR). By replacing J with I band, a new calibration method suitable for use with ACS+WFC3 data is created (i.e. UVI diagram). Using a multi-wavelength multi-aperture photometry catalogue based on CANDELS, UVI colour profiles of 1328 galaxies are stacked in stellar mass and redshift bins. The resulting colour gradients, covering a radial range of 0.2--2.0 effective radii, increase strongly with galaxy mass and with global $A_V$. Colour gradient directions are nearly parallel to the Calzetti extinction vector, indicating that dust plays a more important role than stellar population variations. With our calibration, the resulting $A_V$ profiles fall much more slowly than stellar mass profiles over the measured radial range. sSFR gradients are nearly flat without central quenching signatures, except for $M_*>10^{10.5} M_{\odot}$, where central declines of 20--25 per cent are observed. Both sets of profiles agree well with previous radial sSFR and (continuum) $A_V$ measurements. They are also consistent with the sSFR profiles and, if assuming a radially constant gas-to-dust ratio, gas profiles in recent hydrodynamic models. We finally discuss the striking findings that SFR scales with stellar mass density in the inner parts of galaxies, and that dust content is high in the outer parts despite low stellar-mass surface densities there.

A wide bandwidth, dual polarized, modified four-square antenna is presented as a feed antenna for radio astronomical measurements. A linear array of these antennas is used as a line-feed for cylindrical reflectors for Tianlai, a radio interferometer designed for 21~cm intensity mapping. Simulations of the feed antenna beam patterns and scattering parameters are compared to experimental results at multiple frequencies across the 650 - 1420 MHz range. Simulations of the beam patterns of the combined feed array/reflector are presented as well.

We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star incurs a short period of outburst. In all cases, the secondary accretes part of the ejected mass and also influences the mass-loss rate of the AGB star. The ejected mass may remain gravitationally bound to the binary system and form a circumbinary disk, or contribute to an accretion disk around the secondary. In other cases, the ejecta will escape the binary system. The accretion rate on to the secondary changes non-linearly with binary separation. In our closest binary simulations, our models exemplify the wind Roche lobe overflow while in our wide binary cases, the mass transfer exhibits Bondi-Hoyle accretion. The morphologies of the outflows in the binary systems are varied. The variety may provide clues to how the late AGB phase influences planetary nebulae shaping. We employ the adaptive-mesh-refinement code ASTROBEAR for our simulations and include ray-tracing, radiation transfer, cooling and dust formation. To attain the highest computational efficiency and the most stable results, all simulations are run in the corotating frame.

We report the detection of a curved magnetic field in the ring-like shell of the bubble N4, derived from near-infrared polarization of reddened diskless stars located behind this bubble. The magnetic field in the shell is curved and parallel to the ring-like shell, and its strength is estimated to be $\sim120\,\mu$G in the plane of the sky. The magnetic field strength in the shell is significantly enhanced compared to the local field strength. We calculate the mass-to-flux ratio for the submillimeter clumps in the shell and find that they are all magnetically subcritical. Our results demonstrate that the magnetic field strengthens as the interstellar medium is compressed into a shell, and suggest that the magnetic field has the potential to hinder star formation triggered by \ion{H}{2} region expansion.

Binary stars can interact via mass transfer when one member (the primary) ascends onto a giant branch. The amount of gas ejected by the binary and the amount of gas accreted by the secondary over the lifetime of the primary influence the subsequent binary phenomenology. Some of the gas ejected by the binary will remain gravitationally bound and its distribution will be closely related to the formation of planetary nebulae. We investigate the nature of mass transfer in binary systems containing an AGB star by adding radiative transfer to the AstroBEAR AMR Hydro/MHD code.

Recent observations of the L2 Puppis system suggest that this Mira-like variable may be in the early stages of forming a bipolar planetary nebula (PN). As one of nearest and brightest AGB stars, thought be a binary, L2 Puppis serves as a benchmark object for studying the late-stages of stellar evolution. We perform global, three-dimensional, adaptive-mesh-refinement hydrodynamic simulations of the L2 Puppis system with AstroBEAR. We use the radiative transfer code RADMC-3D to construct the broad-band spectral-energy-distribution (SED) and synthetic observational images from our simulations. Given the reported binary parameters, we are able to reproduce the current observational data if a short pulse of dense material is released from the AGB star with a velocity sufficient to escape the primary but not the binary. Such a situation could result from a thermal pulse, be induced by a periastron passage of the secondary, or could be launched if the primary ingests a planet.

Protostars are mostly found in star-forming regions, where the natal molecular gas still remains. In about 5' west of the molecular bubble N4, N4W is identified as a star-forming clump hosting three Class II (IRS\,1\,--\,3), and one Class I (IRS\,4) young stellar objects (YSOs), as well as a submillimeter source SMM1. The near-IR polarization imaging data of N4W reveal two infrared reflection nebulae close to each other, which are in favor of the outflows of IRS\,1 and IRS\,2. The bipolar mid-IR emission centered on IRS\,4 and the arc-like molecular gas shell are lying on the same axis, indicating a bipolar molecular outflow from IRS\,4. There are two dust temperature distributions in N4W. The warmer one is widely distributed and has a temperature $T_\mathrm{d}\gtrsim28\,\mathrm{K}$, with the colder one from the embedded compact submillimeter source SMM1. N4W's mass is estimated to be $\sim2.5\times10^3\,M_\odot$, and the mass of SMM1 is $\sim5.5\times10^2\,M_\odot$ at $T_\mathrm{d}=15\,\mathrm{K}$, calculated from the CO\,$1-0$ emission and $870\,\mu$m dust continuum emission, respectively. Based on the estimates of bolometric luminosity of IRS\,1\,--\,4, these four sources are intermediate-mass YSOs at least. SMM1 is gravitationally bound, and is capable of forming intermediate-mass stars or even possibly massive stars. The co-existence of the IR bright YSOs and the submillimeter source represents potential sequential star formation processes separated by $\sim0.5$\,Myr in N4W. This small age spread implies that the intermediate-mass star formation processes happening in N4W are almost coeval.

The possibility that mass ejected during Asymptotic Giant Branch (AGB) stellar evolution phases falls back towards the star has been suggested in applications ranging from the formation of accretion disks to the powering of late-thermal pulses. In this paper, we seek to explicate the properties of fallback flow trajectories from mass-loss events. We focus on a transient phase of mass ejection with sub-escape speeds, followed by a phase of a typical AGB wind. We solve the problem using both hydrodynamic simulations and a simplified one-dimensional analytic model that matches the simulations. For a given set of initial wind characteristics, we find a critical shell velocity that distinguishes between "shell fallback" and "shell escape". We discuss the relevance of our results for both single and binary AGB stars. In particular, we discuss how our results help to frame further studies of fallback as a mechanism for forming the substantial population of observed post-AGB stars with dusty disks.

The bright soft X-ray transient Nova Muscae 1991 was intensively observed during its entire 8-month outburst using the Large Area Counter (LAC) onboard the Ginga satellite. Recently, we obtained accurate estimates of the mass of the black hole primary, the orbital inclination angle of the system, and the distance. Using these crucial input data and Ginga X-ray spectra, we have measured the spin of the black hole using the continuum-fitting method. For four X-ray spectra of extraordinary quality we have determined the dimensionless spin parameter of the black hole to be a/M = 0.63 (-0.19, +0.16) (1 sigma confidence level), a result that we confirm using eleven additional spectra of lower quality. Our spin estimate challenges two published results: It is somewhat higher than the value predicted by a proposed relationship between jet power and spin; and we find that the spin of the black hole is decidedly prograde, not retrograde as has been claimed.

The optical counterpart of the black-hole soft X-ray transient Nova Muscae 1991 has brightened by $\Delta{V}\approx0.8$ mag since its return to quiescence 23 years ago. We present the first clear evidence that the brightening of soft X-ray transients in quiescence occurs at a nearly linear rate. This discovery, and our precise determination of the disk component of emission obtained using our $simultaneous$ photometric and spectroscopic data, have allowed us to identify and accurately model archival ellipsoidal light curves of the highest quality. The simultaneity, and the strong constraint it provides on the component of disk emission, is a key element of our work. Based on our analysis of the light curves, and our earlier measurements of the mass function and mass ratio, we have obtained for Nova Muscae 1991 the first accurate estimates of its systemic inclination $i=43.2^{+2.1}_{-2.7}$ deg, and black hole mass $M=11.0^{+2.1}_{-1.4}\ M_\odot$. Based on our determination of the radius of the secondary, we estimate the distance to be $D=4.95^{+0.69}_{-0.65}$ kpc. We discuss the implications of our work for future dynamical studies of black-hole soft X-ray transients.

We carried out observations toward the giant molecular cloud W 37 with the $J = 1 - 0$ transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O using the 13.7 m single-dish telescope at the Delingha station of Purple Mountain Observatory. Based on the three CO lines, we calculated the column densities, cloud masses for the molecular clouds with radial velocities at around $+20 \mathrm{km s}^{-1}$. The gas mass of W 37, calculated from $^{13}$CO emission, is $1.7\times10^5 M_\odot$, above the criteria of giant molecular cloud. The dense ridge of W 37 is a dense filament, which is supercritical in linear mass ratio. Dense clumps found by C$^{18}$O emission are aligned along the dense ridge with a regular interval about 2.8 pc, similar to the clump separation caused by large-scale `sausage instability'. We confirm the identification of the giant molecular filament (GMF) G 18.0-16.8 by \cite{2014A&A...568A..73R} and find a new giant filament, G16.5-15.8, located in the west 0.8 degree of G 18.0-16.8. Both GMFs are not gravitationally bound, as indicated by their low linear mass ratio ($\sim80 M_\odot \mathrm{pc}^{-1}$). We compared the gas temperature map with the dust temperature map from \emph{Herschel} images, and find similar structures. The spatial distributions of class I objects and the dense clumps is reminiscent of triggered star formation occurring in the northwestern part of W 37, which is close to NGC 6611.

We assemble 3524 quasars from Sloan Digital Sky Survey (SDSS) with repeated observations to search for variations of narrow C IV1548,1551 and Mg II2796,2803 absorption doublets in spectral regions shortward of 7000 Ang at the observed frame, which corresponds to time-scales of about 150 ~ 2643 days at quasar rest frame. In these quasar spectra, we detect 3580 C IV absorption systems with z_{abs} = 1.5188 ~ 3.5212, and 1809 Mg II absorption systems with z_{abs} = 0.3948 ~ 1.7167. In term of the absorber velocity (beta) distribution at quasar rest frame, we find a substantial number of C IV absorbers with beta<0.06, which might be connected to the absorptions of quasar outflows. The outflow absorptions peak at v~2000 km/s and drop rapidly below the peak value. Among 3580 C IV absorption systems, 52 systems (~ 1.5%) show obvious variations in equivalent widths at the absorber rest frame (Wr): 16 enhanced, 16 emerged, 12 weaken, and 8 disappeared systems, respectively. We find that changes in Wr548 are neither related to time-scales of the two SDSS observations, nor to absorber velocities at the quasar rest frame. Variable absorptions in low-ionization species are important to constraint the physical conditions of absorbing gas. There are two variable Mg II absorption systems measured from SDSS spectra detected by Hacker et al. However, in our Mg II$ absorption sample, we find that neither shows variable absorption with confident levels of >4sigma for lambda2796 lines and >3sigma for lambda2803 lines.

Inspired by the $f(R)$ non-linear massive gravity, we propose a new kind of modified gravity model, namely $f(T)$ non-linear massive gravity, by adding the dRGT mass term reformulated in the vierbein formalism, to the $f(T)$ theory. We then investigate the cosmological evolution of $f(T)$ massive gravity, and constrain it by using the latest observational data. We find that it slightly favors a crossing of the phantom divide line from the quintessence-like phase ($w_{de} > -1$) to the phantom-like one ($w_{de} < -1$) as redshift decreases.

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.

We use the parameterized post-Newtonian (PPN) formalism to explore the weak field approximation of teleparallel gravity non-minimally coupling to a scalar field $\phi$, with arbitrary coupling function $\omega(\phi)$ and potential $V(\phi)$. We find that all the PPN parameters are identical to general relativity (GR), which makes this class of theories compatible with the Solar System experiments. This feature also makes the theories quite different from the scalar-tensor theories, which might be subject to stringent constraints on the parameter space, or need some screening mechanisms to pass the Solar System experimental constraints.

We investigate the surroundings of the hypercompact HII region M17 UC1 to probe the physical properties of the associated young stellar objects and the environment of massive star formation. Five of the seven point sources in this region show $L$-band excess emission. Geometric match is found between the H_2 emission and near-IR polarized light in the vicinity of IRS5A, and between the diffuse mid-IR emission and near-IR polarization north of UC1. The H_2 emission is typical for dense PDRs, which are FUV pumped initially and repopulated by collisional de-excitation. The spectral types of IRS5A and B273A are B3-B7 V/III and G4-G5 III, respectively. The observed infrared luminosity L_IR in the range 1-20 micron is derived for three objects; we obtain 2.0x10^3 L_\sun for IRS5A, 13 L_\sun for IRS5C, and 10 L_\sun for B273A. IRS5 might be a young quadruple system. Its primary star IRS5A is confirmed to be a high-mass protostellar object (~ 9 M_\sun, ~1x10^5 yrs); it might have terminated accretion due to the feedback from the stellar activities (radiation pressure, outflow) and the expanding HII region of M17. UC1 might also have terminated accretion because of the expanding hypercompact HII region ionized by itself. The disk clearing process of the low-mass YSOs in this region might be accelerated by the expanding HII region. The outflows driven by UC1 are running in south-north with its northeastern side suppressed by the expanding ionization front of M17; the blue-shifted outflow lobe of IRS5A is seen in two types of tracers along the same line of sight in the form of H_2 emission filament and mid-emission. The H_2 line ratios probe the properties of M17 SW PDR, which is confirmed to have a clumpy structure with two temperature distributions: warm, dense molecular clumps with n_H>10^5 cm^-3 and T~575 K and cooler atomic gas with n_H~3.7x10^3-1.5x10^4 cm-3 and T~50-200 K.

As the second work in a series of papers aiming to detect absorption systems in the quasar spectra of the Baryon Oscillation Spectroscopic Survey, we continue the analysis of Paper I by expanding the quasar sample to those quasars with zem>2.4. This yields a sample of 21,963 appropriate quasars to search for narrow C IV absorptions with Wr>=0.2A for both lines. There are 9708 quasars with at least one appropriate absorption system imprinted on their spectra. From these spectra, we detect 13,919 narrow Civabsorption systems whose absorption redshifts cover a range of zabs=1.8784 - 4.3704. In this paper and Paper 1, we have selected 37,241 appropriate quasars with median SNR>= 4 and 1.54<= zem <=5.16 to visually analyze narrow C IV absorption doublets one by one. A total of 15,999 quasars are found to have at least one appropriate absorption system imprinted on their spectra. From these 15,999 quasar spectra, we have detected 23,336 appropriate C IV absorption systems with Wr>=0.2A whose absorption redshifts cover a range of zabs=1.4544 - 4.3704. The largest values of Wr are 3.19 A for the \lambda1548 absorption line and 2.93 A for the lambda1551 absorption line, respectively. We find that only a few absorbers show large values of Wr. About 1.1% absorbers of the total absorbers have Wr>=2.0 A.

We have assembled absorption systems by visually identifying $C~IV\lambda\lambda1548,1551$ absorption doublets in the quasar spectra of the Baryon Oscillation Spectroscopic Survey (BOSS) one by one. This paper is the first of the series work. In this paper, we concern quasars with relatively low redshifts and high signal-to-noise ratios for their spectra, and hence we limit our analysis on quasars with $z_{em}\le2.4$ and on the doublets with $W_r\lambda1548\ge0.2$ \AA. Out of the more than 87,000 quasars in the Data Release 9, we limit our search to 10,121 quasars that have the appropriate redshifts and spectra with high enough signal-to-noise ratios to identify narrow C IV absorption lines. Among them, 5,442 quasars are detected to have at least one $C~IV\lambda\lambda1548,1551$ absorption doublet. We obtain a catalog containing 8,368 $C~IV\lambda\lambda1548,1551$ absorption systems, whose redshifts are within $z_{abs}=1.4544$ - $2.2805$. In this catalog, about $33.7\%$ absorbers have $0.2$ \AA$\le W_r\lambda1548<0.5$ \AA, about $45.9\%$ absorbers have $0.5$ \AA$\le W_r\lambda1548<1.0$ \AA, about $19.2\%$ absorbers have $1.0$ \AA$\le W_r\lambda1548<2.0$ \AA, and about $1.2\%$ absorbers have $W_r\lambda1548\ge2.0$ \AA.

The $L_p - E_p$ relation is a well-known relation in gamma-ray bursts. Its implication remains unclear. We propose to investigate the underlying mechanisms of this relation by considering the corresponding kinetic and dynamic mechanisms separately. In this way, one can tell how much the kinetic or dynamic mechanism contributes to the index of the relationship. Our analysis gives rise to several conclusions. (1) The index of the kinetic effect in the $L_p - E_p$ relation can simply be derived from the theory of special relativity, which is generally larger than 2, depending on the situation concerned. (2) The index of the dynamic effect in the relation can be deduced from observation once a model of jets is adopted. According to current GRB data, we find: the dynamic effect alone tends to make an anti-correlation between $L_p$ and $E_p$; in terms of statistics, the dynamic effect is obviously smaller than the kinetic effect; in the situation of jets with moving discrete radio clouds which moves directly towards the observer, the index of the dynamic effect is currently constrained within $(-1.6, -1)$, while in other situations of jets, the constrains are different; both internal and external shocks can account for the current data.

Gamma-ray bursts (GRBs) are believed to belong to two classes and they are conventionally divided according to their durations. This classification scheme is not satisfied due to the fact that duration distributions of the two classes are heavily overlapped. We collect a new sample (153 sources) of GRBs and investigate the distribution of the logarithmic deviation of the $E_p$ value from the Amati relation. The distribution possesses an obvious bimodality structure and it can be accounted for by the combination of two Gaussian curves. Based on this analysis, we propose to statistically classify GRBs in the well-known $E_{p}$ vs. $E_{iso}$ plane with the logarithmic deviation of the $E_p$ value. This classification scheme divides GRBs into two groups: Amati type bursts and non-Amati type bursts. While Amati type bursts well follow the Amati relation, non-Amati type bursts do not. It shows that most Amati type bursts are long duration bursts and the majority of non-Amati type bursts are short duration bursts. In addition, it reveals that Amati type bursts are generally more energetic than non-Amati type bursts. An advantage of the new classification is that the two kinds of burst are well distinguishable and therefore their members can be identified in certainty.

In this letter, we present for the first time, the discovery of the disappearance of a narrow Mg II $\lambda\lambda2796,2803$ absorption system from the spectra of quasar SDSS J165501.31+260517.4 ($z_{\rm e}=1.8671$). This absorber is located at $z_{\rm abs} =1.7877$, and has a velocity offset of $8,423\rm ~km~s^{-1}$ with respect to the quasar. According to the velocity offset and the line variability, this narrow Mg II $\lambda\lambda2796,2803$ absorption system is likely intrinsic to the quasar. Since the corresponding UV continuum emission and the absorption lines of another narrow Mg II $\lambda\lambda2796,2803$ absorption system at $z_{\rm abs}=1.8656$ are very stable, we think that the disappearance of the absorption system is unlikely to be caused by the change in ionization of absorption gas. Instead, it likely arises from the motion of the absorption gas across the line of sight.

In this paper, we report on two C IV narrow absorption systems emerging from the spectrum of the quasar SDSS J095254.10 +021932.8, which are located at $z_{\rm abs}=2.0053$ and $z_{\rm abs}=1.8011$, respectively. These features have velocity offsets of $\rm \sim 13,800$ and $\rm \sim 34,700~km~s^{-1}$ with respect to the quasar, and show obviously partial coverage to the background emission sources. The two C IV absorption systems are imprinted on the spectrum obtained by the Sloan Digital Sky Survey (SDSS) III on 2011 March 25, which cannot be observed from the spectrum obtained by SDSS-I/II on 2000 December 30. The time interval of the two SDSS observations is 1186.4 d at the quasar rest frame. Because the continuum radiation, broad emission lines and mass accretion rate are stable for the two observations, together with the time interval of observations, we believe that the observed spectral variations are less likely to be caused by the changes in the ionization states of the absorbing gas or by a new outflow, arising from the situation where there was previously no outflow at all. Instead, they probably arise from the motions of multiple streaming gas across the sightline of the quasar.

Quasar J152645.61+193006.7 ($z_{\rm e}=2.5771$) was observed by the Sloan Digital Sky Survey (SDSS) on 2006 May 31, and again on 2011 April 9. The time interval of the SDSS two observations is 497.4 days at the quasar rest frame. From the spectra of this quasar, we detect a phenomenon of disappearance of a narrow $\rm Mg II \lambda\lambda2796,2803$ absorption system with a velocity up to $\rm 166,129 km s^{-1}$ with respect to the quasar. This disappearance event can be caused by changes in the ionization state of absorbing gas or by the bulk motion of the absorbing gas across the quasar sightline. The coverage fraction analysis shows that this absorber partially covers the background emission sources with an effective coverage fraction of $C_{\rm f}=0.40\pm0.06$. The time variation analysis and the coverage fraction analysis imply that this absorber might be intrinsic to the quasar. However, the scenario of a cosmologically separated foreground object located at $z=0.9170$ accounting for the phenomenon cannot be ruled out according to current available data.

The current accelerated expansion of our universe could be due to an unknown energy component with negative pressure (dark energy) or a modification to general relativity (modified gravity). On the other hand, recently warm dark matter (WDM) remarkably rose as an alternative of cold dark matter (CDM). Obviously, it is of interest to distinguish these different types of models. In fact, many attempts have been made in the literature. However, in the present work, we show that WDM, modified gravity and coupled CDM form a trinity, namely, they are indistinguishable by using the cosmological observations of both cosmic expansion history and growth history. Therefore, to break this degeneracy, the other complementary probes beyond the ones of cosmic expansion history and growth history are required.

We investigate the proposed high-mass young stellar object (YSO) candidate CEN34, thought to be associated with the star forming region M17. Its optical to near-infrared (550-2500 nm) spectrum reveals several photospheric absorption features, such as H{\alpha}, Ca triplet and CO bandheads but lacks any emission lines. The spectral features in the range 8375-8770{\AA} are used to constrain an effective temperature of 5250\pm250 (early-/mid-G) and a surface gravity of 2.0\pm0.3 (supergiant). The spectral energy distribution of CEN34 resembles the SED of a high-mass YSO or an evolved star. Moreover, the observed temperature and surface gravity are identical for high-mass YSOs and evolved stars. The radial velocity relative to LSR (V_LSR) of CEN34 as obtained from various photospheric lines is of the order of -60 km/s and thus distinct from the +25 km/s found for several OB stars in the cluster and for the associated molecular cloud. The SED modeling yields ~ 10^{-4} M_sun of circumstellar material which contributes only a tiny fraction to the total visual extinction (11 mag). In the case of a YSO, a dynamical ejection process is proposed to explain the V_LSR difference between CEN34 and M17. Additionally, to match the temperature and luminosity, we speculate that CEN34 had accumulated the bulk of its mass with accretion rate > 4x10^{-3} M_sun/yr in a very short time span (~ 10^3 yrs), and currently undergoes a phase of gravitational contraction without any further mass gain. However, all the aforementioned characteristics of CEN34 are compatible with an evolved star of 5-7 M_sun and an age of 50-100 Myrs, most likely a background post-AGB star with a distance between 2.0 kpc and 4.5 kpc. We consider the latter classification as the more likely interpretation. Further discrimination between the two possible scenarios should come from the more strict confinement of CEN34's distance.

Although $\Lambda$CDM model is very successful in many aspects, it has been seriously challenged. Recently, warm dark matter (WDM) remarkably rose as an alternative of cold dark matter (CDM). In the literature, many attempts have been made to determine the equation-of-state parameter (EoS) of WDM. However, in most of the previous works, it is usually assumed that the EoS of dark matter (DM) is constant (and usually the EoS of dark energy is also constant). Obviously, this assumption is fairly restrictive. It is more natural to assume a variable EoS for WDM (and dark energy). In the present work, we try to constrain the EoS of variable WDM with the current cosmological observations. We find that the best fits indicate WDM, while CDM is still consistent with the current observational data. However, $\Lambda$CDM is still better than WDM models from the viewpoint of goodness-of-fit. So, in order to distinguish WDM and CDM, the further observations on the small/galactic scale are required. On the other hand, in this work we also consider WDM whose EoS is constant, while the role of dark energy is played by various models. We find that the cosmological constraint on the constant EoS of WDM is fairly robust.

After carefully cross-identifying a previously discovered GALEX-selected Lyman Break Galaxy (LBG) candidates one-to-one with their optical counterparts in the field of the CDF-S, we re-estimate their photometric redshifts using multi-wavelength data from UV, optical to NIR. We refine a new updated sample of 383 LBGs at $0.7\la z \la 1.4$. Most LBGs are classified as starburst and irregular types. Ages spread from several Myr to 1.5Gyr. Their dust-corrected star formation rates (SFRs) and stellar masses ($M_*$) are from $4\my$ to $220\my$ and from $2.3\times 10^8 \msun$ to $4 \times 10^{11} \msun$. The rest-frame FUV luminosity function of LBGs are presented. LBGs of irregular types mainly distribute along the "main sequence" of star forming galaxies while most LBGs of starburst types locate in the starburst region. A "downsizing" effect is clearly found and LBGs distribute in the "blue" cloud. HST images in F606W ($V$ band) and F850LP ($z$ band) are taken from the GEMS and GOODS-S surveys. SExtractor and GALFIT are applied to get their morphological parameters. A morphological sample of 142 LBGs with reliable results of $\sersic$ and sizes in both bands is defined. We find that LBGs at $z\sim 1 $ are dominated by disk-like galaxies. Correlations between photometric and morphological properties of LBGs are investigated. Strong correlations between their half-light radii and $M_*$, i.e., size-stellar mass relations, are found in both bands. Physical connections between correlations and the "downsizing" effect are discussed.

We present the analysis of Spitzer/IRAC and NIR imaging observation of AFGL 5157, an active star forming region. In the IRAC images, this region shows strong PAH emissions in channel 4 and H$_2$ emissions in channel 2. Many of the H$_2$ features are aligned to form jet-like structures. Three bipolar jets in the NH$_3$ core region and a couple of jets northwest of the core have been identified. We identify the possible driving agents of the bipolar jets and show them to be very young. An embedded cluster has been detected in the NH$_3$ core; many members in the cluster show their SEDs increasing from JHK bands toward longer wavelength, indicative of their early evolutionary stages. Millimeter and sub-millimeter continuum emissions are found to coincide spatially with these presumable Class 0/I sources, in the NH$_3$ core and the NW subregion. The existence of H$_2$ bipolar jets and very young stellar objects suggests that star formation is still going on at present epoch in these subregions. Combining the information from previous studies, we propose a sequential star formation scenario in the whole AFGL 5157 region.

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.

Based on the joint-observations of the radio broadband spectral emissions of solar eclipse on August 1, 2008 at Jiuquan (total eclipse) and Huairou (partial eclipse) at the frequencies of 2.00 -- 5.60 GHz (Jiuquan), 2.60 -- 3.80 GHZ (Chinese solar broadband radiospectrometer, SBRS/Huairou), and 5.20 -- 7.60 GHz (SBRS/Huairou), the authors assemble a successive series of broadband spectrum with a frequency of 2.60 -- 7.60 GHz to observe the solar eclipse synchronously. This is the first attempt to analyze the solar eclipse radio emission under the two telescopes located at different places with broadband frequencies in the periods of total and partial eclipse. With these analyses, the authors made a new semiempirical model of the coronal plasma density of the quiet Sun and made a comparison with the classic models.

We have collected all available spectra and photometric data from SDSS catalogue for bright AGNs complied from the first three months of the Fermi large area telescope all-sky survey. Based on the 106 high-confidence and 11 low-confidence association bright AGN list, the photometry data are collected from SDSS DR7 for 28 sources (12 BL Lacs and 16 FSRQs), two of which are low-confidence association bright AGNs. Among these 28 SDSS photometric sources, SDSS spectra are available for 20 sources (6 BL Lacs and 14 FSRQs). The black hole mass M_BH and the broad line region (BLR) luminosity were obtained for 14 FSRQs by measuring the line-width and strength of broad emission lines from SDSS spectra. The broad emission lines measurements of five FSRQs are presented for the first time in this work. The optical continuum emission of these 14 FSRQs is found to be likely dominated by the nonthermal jet emission through comparing the relationship between the broad Mg II line and continuum luminosity to that of radio quiet AGNs. The black hole mass of 14 FSRQs ranges from 10^8.2 to 10^9.9 solar mass, with most of sources larger than 10^9 solar mass. The Eddington ratio L_bol/L_Edd ranges from 10^-1.5 to ~ 1. This implies that the optically thin, geometrically thick accretion disk may exist in these FSRQs.

The relationship between the black hole mass and velocity dispersion indicated with [O III] line width is investigated for a sample of 87 flat-spectrum radio quasars (FSRQs) selected from SDSS DR3 quasar catalogue. We found the M_bh - sigma_[O III] relation is deviated from Tremaine et al. relation for nearby inactive galaxies, with a larger black hole mass at given velocity dispersion. There is no strong evidence of cosmology evolution in M_bh - sigma_[O III] relation up to z~0.8. A significant correlation between the [O III] luminosity and Broad Line Region (BLR) luminosity is found. When transferring the [O III] luminosity to Narrow Line Region (NLR) luminosity, the BLR luminosity is, on average, larger than NLR one by about one order of magnitude. We found a strong correlation between the synchrotron peak luminosity and NLR luminosity, which implies a tight relation between the jet physics and accretion process.

We constructed a sample of 185 Flat Spectrum Radio Quasars (FSRQs) by cross-correlating the Shen et al.'s SDSS DR3 X-ray quasar sample with FIRST and GB6 radio catalogues. From the spectrum energy distribution (SED) constructed using multi-band (radio, UV, optical, Infrared and X-ray) data, we derived the synchrotron peak frequency and peak luminosity. The black hole mass and the broad line region (BLR) luminosity (then the bolometric luminosity) were obtained by measuring the line-width and strength of broad emission lines from SDSS spectra. We define a subsample of 118 FSRQs, of which the nonthermal jet emission is thought to be dominated over the thermal emission from accretion disk and host galaxy. For this subsample, we found 25 FSRQs having synchrotron peak frequency > 10^{15} Hz, which is higher than the typical value for FSRQs. While only a weak anti-correlation is found between the synchrotron peak frequency and peak luminosity, it becomes significant when combining with the Wu et al.'s sample of 170 BL Lac objects. At similar peak frequency, the peak luminosity of FSRQs with $\nupeak > 10^{15}$ Hz is systematically higher than that of BL Lac objects, with some FSRQs out of the range covered by BL Lac objects. Although high $\nupeak$ are found in some FSRQs, they do not reach the extreme value of BL Lacs. For the subsample of 118 FSRQs, we found significant correlations between the peak luminosity and black hole mass, the Eddington ratio, and the BLR luminosity, indicating that the jet physics may be tightly related with the accretion process.

The supernova yields of r-process elements are obtained as a function of the mass of their progenitor stars from the abundance patterns of extremely metal-poor stars on the left-side [Ba/Mg]-[Mg/H] boundary with a procedure proposed by Tsujimoto and Shigeyama. The ejected masses of r-process elements associated with stars of progenitor mass $M_{ms}\leq18M_{\odot}$ are infertile sources and the SNe II with 20$M_{\odot}\leq M_{ms}\leq 40M_{\odot}$are the dominant source of r-process nucleosynthesis in the Galaxy. The ratio of these stars 20$M_{\odot}\leq M_{ms}\leq40M_{\odot}$ with compared to the all massive stars is about $\sim$18%. In this paper, we present a simple model that describes a star's [r/Fe] in terms of the nucleosynthesis yields of r-process elements and the number of SN II explosions. Combined the r-process yields obtained by our procedure with the scatter model of the Galactic halo, the observed abundance patterns of the metal-poor stars can be well reproduced