Recently solar flares exhibiting a double J-shaped ribbons in the lower solar atmosphere have been paid increasing attention in the context of extending the two-dimensional standard flare model to three dimensions, as motivated by the spatial correlation between photospheric current channels and flare ribbons. Here we study the electric currents through the photospheric area swept by flare ribbons (termed synthesized ribbon area or SRA), with a sample of 71 two-ribbon flares, of which 36 are J-shaped. Electric currents flowing through one ribbon are highly correlated with those through the other, therefore belonging to the same current system. The non-neutrality factor of this current system is independent of the flare magnitude, implying that both direct and return currents participate in flares. J-shaped flares are distinct from non-J-shaped flares in the following aspects: 1) Electric current densities within J-shaped SRA are significantly smaller than those within non-J-shaped SRA, but J-shaped SRA and its associated magnetic flux are also significantly larger. 2) Electric currents through SRA are positively correlated with the flare magnitude, but J-shaped flares show stronger correlation than non-J-shaped flares. 3) The majority (75\%) of J-shaped flares are eruptive, while the majority (86\%) of non-J-shaped flares are confined; accordingly, hosting active regions of J-shaped flares are more likely to be sigmoidal than non-J-shaped flares. Thus, J-shaped flares constitute a distinct subset of two-ribbon flares, probably the representative of eruptive ones. Further, we found that combining SRA and its associated magnetic flux has the potential to differentiate eruptive from confined flares.

We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses.

The VRI passbands light curves of V344 Lac were presented and analyzed by using the latest version of the W-D code. The observed spectrum reveals that V344 Lac is not an A3 type but would be a later F type star according to the yielded temperature. The results of solution show that V344 Lac is an A-subtype contact binary, with a mediate photometric mass ratio of 0.387(0.003) and a mediate contact factor of 44.6(3.0)%. Based on the parallax given by Gaia, the parameters of the components are estimated as: M1 = 1.16 Ms, M2 = 0.45 Ms, R1 = 1.31 Rs, R2 = 0.88 Rs, L1 = 2.512 Ls, L2 = 1.057 Ls. The period investigation indicates that V344 Lac may have an eccentric orbital oscillation, with P3 = 12.4(0.5) yr, A3 = 0.0020(0.0002) d, and e = 0.38(0.16). Analysis shows such oscillation would be caused by a magnetic activity which can be explained by the Applegate mechanism. Meanwhile, according to the value of l3 and the estimated physical parameters of V344 Lac, the mass of the third companion may be 0.79 Ms. This third body could be a wide company.

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$.

We investigate a new mechanism to create large curvature perturbations on small scales due to parameter resonance in a single-field inflationary model with a small periodic structure upon the potential. After reentering the horizon, the amplified curvature perturbations can lead to observable primordial black holes as well as stochastic gravitational waves. The mass of primordial black holes and frequency of the induced gravitational waves depend on the model parameters. The resulted primordial black hole could constitute all dark matter or a fraction of dark matter in the universe, and corresponding stochastic gravitational waves fall in the frequency band measurable for the pulsar timing array and the space-based gravitational wave detectors.

We investigate two-field inflationary models in which scalar cosmological pertubations are generated via a spectator field nonminimally coupled to gravity, with the particular emphasis on curvaton scenarios. The principal advantage of these models is in the possibility to tune the spectator spectral index via the nonminimal coupling. Our models naturally yield red spectrum of the adiabatic perturbation demanded by observations. We study how the nonminimal coupling affects the spectrum of the curvature perturbation generated in the curvaton scenarios. In particular we find that for small, negative nonminimal couplings the spectral index gets a contribution that is negative and linear in the nonminimal coupling. Since in this way the curvature spectrum becomes redder, some of curvaton scenarios can be saved, which would otherwise be ruled out. In the power law inflation we find that a large nonminimal coupling is excluded since it gives the principal slow roll parameter that is of the order of unity. Finally, we point out that nonminimal coupling can affect the postinflationary growth of the spectator perturbation, and in this way the effectiveness of the curvaton mechanism.

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

In this paper, we present a numerical method to study the predicted lightcurves as a function of viewing angle. We extrapolate the fitting formulae for the mass and velocity of tidal dynamical ejecta across a wide range of mass ratio validated with 66 simulations and use them in the calculations of the kilonova lightcurves. The calculated peak luminosity of a BH-NS merger kilonova is typically about a few times $10^{41}\ {\rm erg\ s^{-1}}$, which is always $\lesssim4.5\times10^{41}\ {\rm erg\ s^{-1}}$. This corresponds to the AB absolute magnitudes fainter than $\sim -15\ {\rm mag}$ in optical and $\sim -16\ {\rm mag}$ in infrared. Since the projected photosphere area of the dynamical ejecta is much larger than that of the disk wind outflows, the dynamical ejecta usually contribute to the majority of the kilonova emission from BH-NS mergers. The fitted blackbody temperature and the shape of the observed multi-band lightcurves are insensitive to the line of sight. The peak time of the observed multi-band lightcurves, affected by the light propagation effect, is related to the relative motion direction between the dynamical ejecta and the observer. The observed luminosity varies with the projected photosphere area determined by the viewing angles. However, the predicted peak luminosity only varies by a factor of $\sim (2 - 3)$ (or by $\sim1\ {\rm mag}$) for different viewing angles. When the short-duration gamma-ray burst afterglow is taken into account, for an on-axis geometry, the kilonova emission is usually outshone by the afterglow emission and can be only observed in the redder bands, especially in the $K$-band at late times. Compared with GW170817/AT2017gfo, the BH-NS merger kilonovae are optically dim but possibly infrared bright. At the same epoch after the merger, the blackbody fitting temperature of the BH-NS merger kilonovae is lower than that of GW170817/AT2017gfo.

LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility and remanent magnetic moment from the test masses and disturb them from their geodesic movement. LISA Pathfinder carried on-board a dedicated magnetic measurement subsystem with noise levels of 10 $ \rm nT \ Hz^{-1/2}$ from 1 Hz down to 1 mHz. In this paper we report on the magnetic measurements throughout LISA Pathfinder operations. We characterise the magnetic environment within the spacecraft, study the time evolution of the magnetic field and its stability down to 20 $\mu$Hz, where we measure values around 200 $ \rm nT \ Hz^{-1/2}$ and identify two different frequency regimes, one related to the interplanetary magnetic field and the other to the magnetic field originating inside the spacecraft. Finally, we characterise the non-stationary component of the fluctuations of the magnetic field below the mHz and relate them to the dynamics of the solar wind.

Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels \emph{Swift} data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult to be interpreted within the framework of a black hole (BH) central engine, but are consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supra-massive magnetar into a black hole when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature would show up, if the fall-back accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in \emph{Swift} sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow light curves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be well interpreted with our proposed model. Future observations are likely to discover more similar events, which could offer more information of the properties of the magnetar as well as the newborn BH.

Inspired by \cite{Jiang:2018uce}, we propose a similar curvaton mechanism whose realization occurs in preheating process, in which the effective mass is running (its potential consists of coupling part and exponential part whose contribution is subdominant comparing to the coupling part). The production of curvaton contains the cases of narrow resonance and broad resonances whose criteria comes via the spectral index of curvaton. Since the inflationary potential is chaotic inflation (quadratic potential), it could smoothly transit into the preheating process. Once the entropy perturbation transferred into curvature perturbation, we will use $\delta N$ formalism to investigate its validity. By neglecting the contribution of exponential potential of curvaton, we calculate power spectrum $P_\zeta$ and non linear Non-Gaussian parameter $f_{NL}$. Our calculation analytically shows that these two observables are independent of potential of inflaton. Finally, as the curvaton almost decay (inflaton field vanishes), the exponential potential will be approaching a constant of order of cosmological constant, which may play a role of dark energy.

We test the Emergent Gravity(EG) theory using the galaxy-galaxy lensing technique based on SDSS DR7 data. In the EG scenario, we do not expect color dependence of the galaxy sample in the 'apparent dark matter' predicted by EG, which is exerted only by the baryonic mass. If the baryonic mass is similar, then the predicted lensing profiles from the baryonic mass should be similar according to the EG, regardless of the color of the galaxy sample. We use the stellar mass of the galaxy as a proxy of its baryonic mass. We divide our galaxy sample into 5 stellar mass bins, and further classify them as red and blue subsamples in each stellar mass bin. If we set halo mass and concentration as free parameters, $\Lambda$CDM is favored by our data in terms of the reduced $\chi^2$ while EG fails to explain the color dependence of ESDs from the galaxy-galaxy lensing measurement.

The orientation, chirality, and dynamics of solar eruptive filaments is a key to understanding the magnetic field of coronal mass ejections (CMEs) and therefore to predicting the geoeffectiveness of CMEs arriving at Earth. However, confusion and contention remain over the relationship between the filament chirality, magnetic helicity, and sense of rotation during eruption. To resolve the ambiguity in observations, in this paper, we used stereoscopic observations to determine the rotation direction of filament apex and the method proposed by Chen et al. (2014) to determine the filament chirality. Our sample of 12 eruptive active-region filaments establishes a strong one-to-one relationship, i.e., during the eruption, sinistral/dextral filaments (located in the southern/northern hemisphere) rotate clockwise/counterclockwise when viewed from above, and corroborates a weak hemispheric preference, i.e., a filament and related sigmoid both exhibit a forward (reverse) S shape in the southern (northern) hemisphere, which suggests that the sigmoidal filament is associated with a low-lying magnetic flux rope with its axis dipped in the middle. As a result of rotation, the projected S shape of a filament is anticipated to be reversed during eruption.

We consider gravitational radiation and electromagnetic radiation from point mass binary with electric charges in a Keplerian orbit, and calculate the merger rate distribution of primordial black hole binaries with charges and a general mass function by taking into account gravitational torque and electromagnetic torque by the nearest primordial black hole. We apply the formalism to the extremal charged case and find that $\alpha=-(m_i+m_j)^2\partial^2 \ln {\cal R}(m_i,m_j)/\partial m_i \partial m_j=12/11$, which is independent of the mass function.

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

In this paper, we present the study for the energy reservoir powering the light curves (LCs) of PS1-12cil and SN 2012aa which are superluminous and luminous supernovae (SNe), respectively. The multi-band and bolometric LCs of these two SNe show unusual secondary bumps after the main peaks. The two-peaked LCs cannot be explained by any simple energy-source models (e.g., the $^{56}$Ni cascade decay model, the magnetar spin-down model, and the ejecta-circumstellar medium interaction model). Therefore, we employ the $^{56}$Ni plus ejecta-circumstellar medium (CSM) interaction (CSI) model, the magnetar plus CSI model, and the double CSI model to fit their bolometric LCs, and find that both these two SNe can be explained by the double CSI model and the magnetar plus CSI model. Based on the modeling, we calculate the the time when the shells were expelled by the progenitors: provided that they were powered by double ejecta-shell CSI, the inner and outer shells might be expelled $\sim 0.2-3.6$ and $\sim 2-25$ years before the explosions of the SNe, respectively; the shells were expelled $\sim 2-20$ years before the explosions of the SNe if they were powered by magnetars plus CSI.

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$\sigma$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

The joint detection of GW 170817 and GRB 170817A indicated that at least a fraction of short gamma ray bursts (SGRBs) originate from binary neutron star (BNS) mergers. One possible remnant of a BNS merger is a rapidly rotating, strongly magnetized neutron star, which has been discussed as one possible central engine for GRBs. For a rapidly rotating magnetar central engine, the deposition of the rotation energy into the ejecta launched from the merger could lead to bright radio emission. The brightness of radio emission years after an SGRB would provide an estimation of the kinetic energy of ejecta and hence, a possible constraint on the BNS merger product. We perform a more detailed calculation on the brightness of radio emission from the interaction between the merger ejecta and circumburst medium in the magnetar scenario, invoking several important physical processes such as generic hydrodynamics, relativistic effects, and the deep Newtonian phase. We use the model to constrain the allowed parameter space for 15 SGRBs that have late radio observations. Our results show that an injected energy $E_{\rm inj} \sim 10^{52}$ erg is allowed for all the cases, which suggests that the possibility of a supra-massive or hyper-massive neutron star remnant is not disfavored by the available radio data.

We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 \mu m}$) and 22 $\mu$m ($L_{\rm 22 \mu m}$), as well as the dependence of the residuals on a variety of galaxy properties. The correlation with $L_{\rm 12 \mu m}$ is tighter and more linear, but galaxies with relatively low stellar masses and blue colors fall significantly below the mean $L_{\rm CO(1-0)}-L_{\rm 12\mu m}$ relation. We propose a new estimator of the CO(1-0) luminosity (and thus the total molecular gas mass) that is a linear combination of three parameters: $L_{\rm 12 \mu m}$, $M_\ast$ and $g-r$. We show that, with a scatter of only 0.18 dex in log $(L_{\rm CO(1-0)})$, this estimator provides unbiased estimates for galaxies of different properties and types. An immediate application of this estimator to a compiled sample of galaxies with only CO(J=2-1) observations yields a distribution of the CO(J=2-1) to CO(J=1-0) luminosity ratios ($R21$) that agrees well with the distribution of real observations, in terms of both the median and the shape. Application of our estimator to the current MaNGA sample reveals a gas-poor population of galaxies that are predominantly early-type and show no correlation between molecular gas-to-stellar mass ratio and star formation rate, in contrast to gas-rich galaxies. We also provide alternative estimators with similar scatters, based on $r$ and/or $z$ band luminosities instead of $M_\ast$. These estimators serve as cheap and convenient $M_{\rm mol}$ proxies to be potentially applied to large samples of galaxies, thus allowing statistical studies of gas-related processes of galaxies.

We present the Star cluster Hunting Pipeline (SHiP) which can identify star clusters in GAIA DR2 data, and establish a star cluster catalog for the Galactic disk. A Friend of Friend based cluster finder method is used to identify star clusters using 5-dimensional stellar parameters, $l, b, \varpi, \mu_\alpha\cos\delta$, and $\mu_\delta$. Our new catalog contains 2443 star cluster candidates identified from disk stars located within $|b|=25^\circ$ and with $G<18$ mag. An automatic isochrone fitting scheme is applied to all cluster candidates. With a combination of parameters obtained from isochrone fitting, we classify cluster candidates into three classes (Class 1, 2 and 3). Class 1 clusters are the most probable star cluster candidates with the most stringent criteria. Most of these clusters are nearby (within 4 kpc). Our catalog is cross-matched with three Galactic star cluster catalogs, Kharchenko et al. (2013), Cantat-Gaudin et al. (2018, 2019), and Bica et al. (2019). The proper motion and parallax of matched star clusters are in good agreement with these earlier catalogs. We discover 76 new star cluster candidates that are not listed in these 3 catalogs. The majority of these are clusters older than log(age/yr)=8.0, and are located in the inner disk with $|b|<5^\circ$. The recent discovery of new star clusters suggests that current Galactic star cluster catalogs are still incomplete. Among the Class 1 cluster candidates, we find 56 candidates for star cluster groups. The pipeline, the catalog and the member list containing all candidates star clusters and star cluster groups have been made publicly available.

The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,\mu\textrm{Ns}$ to 230$\,\mu\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.

It is well-known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 days, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar+56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar+56Ni or cooling+magnetar+56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar+56Ni is very high.

As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM), we present a measurement of the mass of the supermassive black hole (SMBH) in the nearby early-type galaxy NGC 0383 (radio source 3C 031). This measurement is based on Atacama Large Millimeter/sub-milimeter Array (ALMA) cycle 4 and 5 observations of the 12CO(2-1) emission line with a spatial resolution of 58x32pc2 (0."18x0."1). This resolution, combined with a channel width of 10 km/s, allows us to well resolve the radius of the black hole sphere of influence (measured as R_SOI = 316pc = 0."98), where we detect a clear Keplerian increase of the rotation velocities. NGC 0383 has a kinematically-relaxed, smooth nuclear molecular gas disc with weak ring/spiral features. We forward-model the ALMA data cube with the Kinematic Molecular Simulation (KinMS) tool and a Bayesian Markov Chain Monte Carlo method to measure a SMBH mass of (4.2+/-0.7)x10^9 Msun, a F160W-band stellar mass-to-light ratio that varies from 2.8+/-0.6 Msun/Lsun in the centre to 2.4+/-0.3 Msun/Lsun at the outer edge of the disc and a molecular gas velocity dispersion of 8.3+/-2.1 km/s (all 3-sigma uncertainties). We also detect unresolved continuum emission across the full bandwidth, consistent with synchrotron emission from an active galactic nucleus. This work demonstrates that low-J CO emission can resolve gas very close to the SMBH (~140,000 Schwarzschild radii) and hence that the molecular gas method is highly complimentary to megamaser observations as it can probe the same emitting material.

We develop a formalism to calculate the merger rate density of primordial black hole binaries with a general mass function, by taking into account the merger history of primordial black holes. We apply the formalism to three specific mass functions, monochromatic, power-law and log-normal cases. In the former case, the merger rate density is dominated by the single-merger events, while in the latter two cases, the contribution of the multiple-merger events on the merger rate density can not be ignored. The effects of the merger history on the merger rate density depend on the mass function.

The quasi-thermal components found in many Fermi gamma-ray bursts (GRBs) imply that the photosphere emission indeed contributes to the prompt emission of many GRBs. But whether the observed spectra empirically fitted by the Band function or cutoff power law, especially the spectral and peak energy ($E_{p}$) evolutions can be explained by the photosphere emission model alone needs further discussion. In this work, we investigate in detail the time-resolved spectra and $E_{p}$ evolutions of photospheric emission from a structured jet, with an inner-constant and outer-decreasing angular Lorentz factor profile. Also, a continuous wind with a time-dependent wind luminosity has been considered. We show that the photosphere spectrum near the peak luminosity is similar to the cutoff power-law spectrum. The spectrum can have the observed average low-energy spectral index $\alpha $ $ \sim -1$, and the distribution of the low-energy spectral index in our photosphere model is similar to that observed ($-2\lesssim $ $\alpha \lesssim 0$). Furthermore, the two kinds of spectral evolutions during the decay phase, separated by the width of the core ($\theta _{c}$), are consistent with the time-resolved spectral analysis results of several Fermi multi-pulse GRBs and single-pulse GRBs, respectively. Also, for this photosphere model we can reproduce the two kinds of observed $E_{p}$ evolution patterns rather well. Thus, by considering the photospheric emission from a structured jet, we reproduce the observations well for the GRBs best fitted by the cutoff power-law model for the peak-flux spectrum or the time-integrated spectrum.

We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in which this balancing force is applied continuously, with the advantage that the acceleration noise during free fall is measured in the absence of the actuation force, thus eliminating associated noise and force calibration errors. The differential acceleration noise measurement presented here with the free-fall mode agrees with noise measured with the continuous actuation scheme, representing an important and independent confirmation of the LPF result. An additional measurement with larger actuation forces also shows that the technique can be used to eliminate actuation noise when this is a dominant factor.

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements.

In this work, we investigate the formation of a magnetic flux rope (MFR) above the central polarity inversion line (PIL) of NOAA Active Region 12673 during its early emergence phase. Through analyzing the photospheric vector magnetic field, extreme ultraviolet (EUV) and ultraviolet (UV) images, extrapolated three-dimensional (3D) non-linear force-free fields (NLFFFs), as well as the photospheric motions, we find that with the successive emergence of different bipoles in the central region, the conjugate polarities separate, resulting in collision between the non-conjugated opposite polarities. Nearly-potential loops appear above the PIL at first, then get sheared and merge at the collision locations as evidenced by the appearance of a continuous EUV sigmoid on 2017 September 4, which also indicates the formation of an MFR. The 3D NLFFFs further reveal the gradual buildup of the MFR, accompanied by the appearance of two elongated bald patches (BPs) at the collision locations and a very low-lying hyperbolic flux tube configuration between the BPs. The final MFR has relatively steady axial flux and average twist number of around $2.1\times 10^{20}$~Mx and -1.5, respective. Shearing motions are found developing near the BPs when the collision occurs, with flux cancellation and UV brightenings being observed simultaneously, indicating the development of a process named as "collisional shearing" (firstly identified by Chintzoglou et al. 2019). The results clearly show that the MFR is formed by "collisional shearing", i.e., through shearing and flux cancellation driven by the collision between non-conjugated opposite polarities during their emergence.

With the aim of understanding the physical mechanisms of confined flares, we selected 18 confined flares during 2011-2017, and classified the confined flares into two types based on their different dynamic properties and magnetic configurations. "Type I" of confined flares are characterized by slipping reconnection, strong shear, and stable filament. "Type II" flares have nearly no slipping reconnection, and have a configuration in potential state after the flare. Filament erupts but is confined by strong strapping field. "Type II" flares could be explained by 2D MHD models while "type I" flares need 3D MHD models. 7 flares of 18 ($\sim$39 \%) belong to "type I" and 11 ($\sim$61 \%) are "type II" confined flares. The post-flare loops (PFLs) of "type I" flares have a stronger non-potentiality, however, the PFLs in "type II" flares are weakly sheared. All the "type I" flares exhibit the ribbon elongations parallel to the polarity inversion line (PIL) at speeds of several tens of km s$^{-1}$. For "type II" flares, only a small proportion shows the ribbon elongations along the PIL. We suggest that different magnetic topologies and reconnection scenarios dictate the distinct properties for the two types of flares. Slipping agnetic reconnections between multiple magnetic systems result in "type I" flares. For "type II" flares, magnetic reconnections occur in anti-parallel magnetic fields underlying the erupting filament. Our study shows that "type I" flares account for more than one third of the overall large confined flares, which should not be neglected in further studies.

Space very long baseline interferometry (VLBI) has unique applications in high-resolution imaging of fine structure of astronomical objects and high-precision astrometry due to the key long space-Earth or space-space baselines beyond the Earth's diameter. China has been actively involved in the development of space VLBI in recent years. This review briefly summarizes China's research progress in space VLBI and the future development plan.

LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall accuracy in many ways. A dedicated temperature measurement subsystem, with noise levels in 10$\,\mu$K$\,$Hz$^{-1/2}$ down to $1\,$mHz was part of the diagnostics unit on board LPF. In this paper we report on the temperature measurements throughout mission operations, characterize the thermal environment, estimate transfer functions between different locations and report temperature stability (and its time evolution) at frequencies as low as 10$\,\mu$Hz, where typically values around $1\,$K$\,$Hz$^{-1/2}$ were measured.

A second emission enhancement in warm coronal extreme-ultraviolet (EUV) lines (about 2-7 MK) during some solar flares is known as the EUV late phase. Imaging observations confirm that the late phase emission originates from a set of longer or higher loops than the main flare loops. Nevertheless, some questions remain controversial: What is the relationship between these two loop systems? What is the heating source of late phase emission, a heating accompany the main phase heating or occuring quite later? In this paper, we present clear evidence for heating source in a late-phase solar flare: magnetic reconnection of overlying field in a quadrupolar magnetic configuration. The event is triggered by an erupted core structure that eventually leads to a coronal mass ejection (CME). Cusp feature and its shrinkage motion high in the late-phase emission region are the manifestation of the later phase reconnection following the main flare reconnection. Using the enthalpy-based thermal evolution of loops (EBTEL) model, we reasonably reproduce the late-phase emissions in some EUV lines. We suggest that a continuous additional heating is responsible for the appearance of the elongated EUV late phase.

To investigate the factors that control the success and/or failure of solar eruptions, we study the magnetic field and 3-Dimensional (3D) configuration of 16 filament eruptions during 2010 July - 2013 February. All these events, i.e., erupted but failed to be ejected to become a coronal mass ejection (CME), are failed eruptions with the filament maximum height exceeding $100 Mm$. The magnetic field of filament source regions is approximated by a potential field extrapolation method. The filament 3D configuration is reconstructed from three vantage points by the observations of STEREO Ahead/Behind and SDO spacecraft. We calculate the decay index at the apex of these failed filaments and find that in 7 cases, their apex decay indexes exceed the theoretical threshold ($n_{crit} = 1.5$) of the torus instability. We further determine the orientation change or rotation angle of each filament top during the eruption. Finally, the distribution of these events in the parameter space of rotation angle versus decay index is established. Four distinct regimes in the parameter space are empirically identified. We find that, all the torus-unstable cases (decay index $n > 1.5$), have a large rotation angles ranging from $50^\circ - 130^\circ$. The possible mechanisms leading to the rotation and failed eruption are discussed. These results imply that, besides the torus instability, the rotation motion during the eruption may also play a significant role in solar eruptions.

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

SN~2007D is a nearby (redshift $z = 0.023146$), luminous Type Ic supernova (SN) having a narrow light curve (LC) and high peak luminosity. Previous research based on the assumption that it was powered by the $^{56}$Ni cascade decay suggested that the inferred $^{56}$Ni mass and the ejecta mass are $\sim 1.5$M$_{\odot}$ and $\sim 3.5$M$_{\odot}$, respectively. In this paper, we employ some multiband LC models to model the $R$-band LC and the color ($V-R$) evolution of SN~2007D to investigate the possible energy sources powering them. We find that the pure $^{56}$Ni model is disfavored; the multiband LCs of SN~2007D can be reproduced by a magnetar whose initial rotational period $P_{0}$ and magnetic field strength $B_p$ are $7.28_{-0.21}^{+0.21}$ (or $9.00_{-0.42}^{+0.32}$) ms and $3.10_{-0.35}^{+0.36}\times 10^{14}$ (or $2.81_{-0.44}^{+0.43}\times 10^{14}$) G, respectively. By comparing the spectrum of SN~2007D with that of some superluminous SNe (SLSNe), we find that it might be a luminous SN like several luminous ``gap-filler" optical transients that bridge ordinary and SLSNe, rather than a genuine SLSN.

Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions.

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

We develop an analytic formalism for computing the merger rate of primordial black hole binaries with a general mass function by taking into account the torques by the surrounding primordial black holes and linear density perturbations. We find that $\alpha=-(m_{i}+m_{j})^2\partial^{2} \ln {\mathcal R}(m_{i},m_{j})/\partial m_{i} \partial m_{j}=36/37$ is independent of the mass function. Moreover, the ratio of the merger rate density of primordial black hole binaries by taking into account the torques by the surrounding primordial black holes to by the nearest primordial black hole is independent of the masses of binaries.

Since the 2017 Nobel Prize in Physics was awarded for the observation of gravitational waves, it is fair to say that the epoch of gravitational wave astronomy (GWs) has begun. However, a number of interesting sources of GWs can only be observed from space. To demonstrate the feasibility of the Laser Interferometer Space Antenna (LISA), a future gravitational wave observatory in space, the LISA Pathfinder satellite was launched on December, 3rd 2015. Measurements of the spurious forces accelerating an otherwise free-falling test mass, and detailed investigations of the individual subsystems needed to achieve the free-fall, have been conducted throughout the mission. This overview article starts with the purpose and aim of the mission, explains satellite hardware and mission operations and ends with a summary of selected important results and an outlook towards LISA. From the LISA Pathfinder experience, we can conclude that the proposed LISA mission is feasible.

We present high angular resolution (0.3" or 37 pc) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the CO(2-1) line emission from a central disc in the early-type galaxy NGC 524. This disc is shown to be dynamically relaxed, exhibiting ordered rotation about a compact 1.3mm continuum source, which we identify as emission from an active supermassive black hole (SMBH). There is a hole at the centre of the disc slightly larger than the SMBH sphere of influence. An azimuthal distortion of the observed velocity field is found to be due to either a position angle warp or radial gas flow over the inner 2.5". By forward-modelling the observations, we obtain an estimate of the SMBH mass of $4.0^{+3.5}_{-2.0}\times10^8\,\mathrm{M_\odot}$, where the uncertainties are at the $3\sigma$ level. The uncertainties are dominated by the poorly constrained inclination and the stellar mass-to-light ratio of this galaxy, and our measurement is consistent with the established correlation between SMBH mass and stellar velocity dispersion. Our result is roughly half that of the previous stellar dynamical measurement, but is consistent within the uncertainties of both. We also present and apply a new tool for modelling complex molecular gas distributions.

In our previous work, we propose a cross spectrum based method to extract single pulse signals from RFI contaminated data, which is originated from geodetic VLBI postprocessing. This method fully utilizes fringe phase information of the cross spectrum and hence maximizes signal power, however the localization was not discussed in that work yet. As the continuation of that work, in this paper, we further study how to localize single pulses using astrometric solving method. Assuming that the burst is a point source, we derive the burst position by solving a set of linear equations given the relation between residual delay and offset to a priori position. We find that the single pulse localization results given by both astrometric solving and radio imaging are consistent within 3 sigma level. Therefore we claim that it is possible to derive the position of a single pulse with reasonable precision based on only 3 or even 2 baselines with 4 milliseconds integration. The combination of cross spectrum based detection and the localization proposed in this work then provide a thorough solution for searching single pulse in VLBI observation. According to our calculation, our pipeline gives comparable accuracy as radio imaging pipeline. Moreover, the computational cost of our pipeline is much smaller, which makes it more practical for FRB search in regular VLBI observation. The pipeline is now publicly available and we name it as "VOLKS", which is the acronym of "VLBI Observation for frb Localization Keen Searcher".

The science operations of the LISA Pathfinder mission has demonstrated the feasibility of sub-femto-g free-fall of macroscopic test masses necessary to build a LISA-like gravitational wave observatory in space. While the main focus of interest, i.e. the optical axis or the $x$-axis, has been extensively studied, it is also of interest to evaluate the stability of the spacecraft with respect to all the other degrees of freedom. The current paper is dedicated to such a study, with a focus set on an exhaustive and quantitative evaluation of the imperfections and dynamical effects that impact the stability with respect to its local geodesic. A model of the complete closed-loop system provides a comprehensive understanding of each part of the in-loop coordinates spectra. As will be presented, this model gives very good agreements with LISA Pathfinder flight data. It allows one to identify the physical noise source at the origin and the physical phenomena underlying the couplings. From this, the performances of the stability of the spacecraft, with respect to its geodesic, are extracted as a function of frequency. Close to $1 mHz$, the stability of the spacecraft on the $X_{SC}$, $Y_{SC}$ and $Z_{SC}$ degrees of freedom is shown to be of the order of $5.0\ 10^{-15} m\ s^{-2}/\sqrt{Hz}$ for X and $4.0 \ 10^{-14} m\ s^{-2}/\sqrt{Hz}$ for Y and Z. For the angular degrees of freedom, the values are of the order $3\ 10^{-12} rad\ s^{-2}/\sqrt{Hz}$ for $\Theta_{SC}$ and $3\ 10^{-13} rad\ s^{-2}/\sqrt{Hz}$ for $H_{SC}$ and $\Phi_{SC}$.

Recent wide-field surveys discovered new types of peculiar optical transients that showed diverse behaviors of the evolution of photospheric properties. We develop a general theory of homologous explosions with constant opacity, paying special attention on the evolution of the photospheric radius $R_{\rm ph}$. We find that regardless of the density distribution profile, $R_{\rm ph}$ always increases early on and decreases at late times. This result does not depend on the radiation and cooling processes inside the ejecta.The general rising/falling behavior of $R_{\rm ph}$ can be used to quickly diagnose whether the source originates from a supernova-like explosion. The shape of the $R_{\rm ph}$ evolution curve depends on the density profile, so the observations may be used to directly diagnose the density profile as well as the temperature profile of the ejecta. All the well-monitored supernovae show such a $R_{\rm ph}$ rising/falling behavior, which is consistent with our theory. The recently discovered peculiar transient AT2018cow showed a continuous decay of $R_{\rm ph}$, which is disfavored to be of a supernova-like explosion origin. Our result therefore supports the interpretation of this transient as a tidal disruption event.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

Magnetic flux ropes (MFRs) are believed to be the core structure in solar eruptions, nevertheless, their formation remains intensely debated. Here we report a rapid buildup process of an MFR-system during a confined X2.2 class flare occurred on 2017 September 6 in NOAA AR 12673, three hours after which the structure erupted to a major coronal mass ejection (CME) accompanied by an X9.3 class flare. For the X2.2 flare, we do not find EUV dimmings, separation of its flare ribbons, or clear CME signatures, suggesting a confined flare. For the X9.3 flare, large-scale dimmings, separation of its flare ribbons, and a CME show it to be eruptive. By performing a time sequence of nonlinear force-free fields (NLFFFs) extrapolations we find that: until the eruptive flare, an MFR-system was located in the AR. During the confined flare, the axial flux and the lower bound of the magnetic helicity for the MFR-system were dramatically enhanced by about 86% and 260%, respectively, although the mean twist number was almost unchanged. During the eruptive flare, the three parameters were all significantly reduced. The results evidence the buildup and release of the MFR-system during the confined and the eruptive flare, respectively. The former may be achieved by flare reconnection. We also calculate the pre-flare distributions of the decay index above the main polarity inversion line (PIL) and find no significant difference. It indicates that the buildup of the magnetic flux and helicity of the MFR-system may play a role in facilitating its final eruption.

The Intermediate Palomar Transient Factory reported the discovery of an unusual type II-P supernova iPTF14hls. Instead of a ~100-day plateau as observed for ordinary type II-P supernovae, the light curve of iPTF14hls has at least five distinct peaks, followed by a steep decline at ~1000 days since discovery. Until 500 days since discovery, the effective temperature of iPTF14hls is roughly constant at 5000-6000K . In this paper we propose that iPTF14hls is likely powered by intermittent fallback accretion. It is found that the light curve of iPTF14hls can be well fit by the usual t^{-5/3} accretion law until ~1000 days post discovery when the light curve transitions to a steep decline. To account for this steep decline, we suggest a power-law density profile for the late accreted material, rather than the constant profile as appropriated for the t^{-5/3} accretion law. Detailed modeling indicates that the total fallback mass is ~0.2M_{sun}, with an ejecta mass M_{ej}~21M_{sun}. We find the third peak of the light curve cannot be well fit by the fallback model, indicating that there could be some extra rapid energy injection. We suggest that this extra energy injection may be a result of a magnetic outburst if the central object is a neutron star. These results indicate that the progenitor of iPTF14hls could be a massive red supergiant.

The gravitational wave event GW170817 associated with the short gamma-ray burst (GRB) 170817A confirms that binary neutron star (BNS) mergers are one of the origins of short GRBs. The associated kilonova emission, radioactively powered by nucleosynthesized heavy elements, was also detected. Black hole-neutron star (BH-NS) mergers have been argued to be another promising origin candidate of short GRBs and kilonovae. Numerical simulations show that the ejecta in BH-NS mergers is geometrically much more anisotropic than the BNS merger case. In this paper, we investigate observational signatures of kilonova emission from the anisotropic ejecta in BH-NS mergers. We find that a bump appears on the bolometric luminosity light curve due to the inhomogeneous mass distribution in the latitudinal direction. The decay slope of the single-band light curve becomes flatter and the spectrum also deviates from a single-temperature blackbody radiation spectrum due to the gradient in the velocity distribution of the ejecta. Future detection or non-detection of such signatures would be useful to test the mass ejection geometry in BH-NS mergers.

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

We study the $R^p$ inflationary model of [Muller:1989rp] for $p>2$ using the result of Ref. [Motohashi:2014tra]. After reproducing the observable quantities: the power spectral index $n_s$, its corresponding running $\alpha=\frac{dn_s}{d\ln(k)}$ and the tensor to scalar ration $r$ in terms of e-folding number $N$ and $p$, we show that $R^p$ inflation model is still alive as $p$ is from $2$ to $2.02$. In this range, our calculation confirms that $n_s$ and $r$ agree with observations and $\alpha$ is of order $10^{-4}$ which needs more precise observational constraints. We find that, as the value of $p$ increases, all $n_s$, $r$ and $|\alpha|$ decrease. However, the precise interdependence between these observables is such that this class of models can in principle be tested by the next generation of dedicated satellite CMB probes.

LISA Pathfinder (LPF) was a European Space Agency mission with the aim to test key technologies for future space-borne gravitational-wave observatories like LISA. The main scientific goal of LPF was to demonstrate measurements of differential acceleration between free-falling test masses at the sub-femto-g level, and to understand the residual acceleration in terms of a physical model of stray forces, and displacement readout noise. A key step toward reaching the LPF goals was the correct calibration of the dynamics of LPF, which was a three-body system composed by two test-masses enclosed in a single spacecraft, and subject to control laws for system stability. In this work, we report on the calibration procedures adopted to calculate the residual differential stray force per unit mass acting on the two test-masses in their nominal positions. The physical parameters of the adopted dynamical model are presented, together with their role on LPF performance. The analysis and results of these experiments show that the dynamics of the system was accurately modeled and the dynamical parameters were stationary throughout the mission. Finally, the impact and importance of calibrating system dynamics for future space-based gravitational wave observatories is discussed.

The first gravitational-wave event from the merger of a binary neutron star system (GW170817) was detected recently. The associated short gamma-ray burst (GRB 170817A) has a low isotropic luminosity ($\sim 10^{47}$ erg s$ ^{-1}$) and a peak energy $E_{p}\sim 145$ keV during the initial main emission between -0.3 and 0.4 s. The origin of this short GRB is still under debate, but a plausible interpretation is that it is due to the off-axis emission from a structured jet. We consider two possibilities. First, since the best-fit spectral model for the main pulse of GRB 170817A is a cutoff power law with a hard low-energy photon index ($\alpha =-0.62_{-0.54}^{+0.49} $), we consider an off-axis photosphere model. We develop a theory of photosphere emission in a structured jet and find that such a model can reproduce a low-energy photon index that is softer than a blackbody through enhancing high-latitude emission. The model can naturally account for the observed spectrum. The best-fit Lorentz factor along the line of sight is $\sim 20$, which demands that there is a significant delay between the merger and jet launching. Alternatively, we consider that the emission is produced via synchrotron radiation in an optically thin region in an expanding jet with decreasing magnetic fields. This model does not require a delay of jet launching but demands a larger bulk Lorentz factor along the line of sight. We perform Markov Chain Monte Carlo fitting to the data within the framework of both models and obtain good fitting results in both cases.

Multi-band observations of the fast radio burst (FRB) 131104 show that this burst may be associated with a gamma-ray transient entitled Swift J0644.5-5111. Follow-up observations for potential X-ray and radio counterparts of FRB 131104/Swift J0644.5-5111 got null results and provided the upper limits of the emission flux at 5.5 GHz, 7.5 GHz, $U$-band, and X-ray band. By assuming this association and using these upper limits, environmental properties (the fraction of energyin a magnetic field $\varepsilon_{\rm B}$ and the number density $n$) of the progenitor system of FRB 131104/Swift J0644.5-5111 were constrained in the context of the standard afterglow model that neglects the non-relativistic effect and jet effect by several groups. In this paper, we adopt a unified afterglow model that takes into account the non-relativistic effect and jet effect and use the upper limits of four bands (5.5 GHz, 7.5 GHz, $U$-band, and X-ray) to obtain more stringent constraints on the parameter space spanned by $\varepsilon_{\rm B}$ and $n$. We thus suggest that FRB 131104/Swift J0644.5-5111 might originate from a black hole-neutron star merger event. Moreover, we calculate multi-band emissions from a kilonova powered by the radioactivity of $r$-process elements synthesized in the ejected neutron-rich material and find that the $U$-band emission from the putative kilonova is significantly lower than the upper limit of the observations.

The successful joint observation of the gravitational wave event GW170817 and its multi-wavelength electromagnetic counterparts first enables human to witness a definite merger event of two neutron stars (NSs). This historical event confirms the origin of short-duration gamma-ray bursts (GRBs), and in particular, identifies the theoretically-predicted kilonova phenomenon that is powered by radioactive decays of $r$-process heavy elements. However, whether a long-lived remnant NS could be formed during this merger event remains unknown, although such a central engine has been suggested by afterglow observations of some short-duration GRBs. By invoking this long-lived remnant NS, we here propose a model of hybrid energy sources for the kilonova AT2017gfo associated with GW 170817. While the early emission of AT2017gfo is still powered radioactively as usually suggested, its late emission is primarily caused by delayed energy injection from the remnant NS. In our model, only one single opacity is required and an intermediate value of $\kappa\simeq0.97\,\rm cm^2g^{-1}$ is revealed, which could be naturally provided by lanthanide-rich ejecta that is deeply ionized by the emission from a wind of the NS. These self-consistent results indicate that a long-lived remnant NS, which must own a very stiff equation of state, had been formed during the merger event of GW170817. This provides a very stringent constraint on the strong interaction in nuclear-quark matter. It is further implied that such GW events could provide a probe of the early spin and magnetic evolutions of NSs, e.g., the burying of surface magnetic fields.

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

We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies.

Galactic cosmic-ray (GCR) energy spectra observed in the inner heliosphere are modulated by the solar activity, the solar polarity and structures of solar and interplanetary origin. A high counting rate particle detector (PD) aboard LISA Pathfinder (LPF), meant for subsystems diagnostics, was devoted to the measurement of galactic cosmic-ray and solar energetic particle integral fluxes above 70 MeV n$^{-1}$ up to 6500 counts s$^{-1}$. PD data were gathered with a sampling time of 15 s. Characteristics and energy-dependence of GCR flux recurrent depressions and of a Forbush decrease dated August 2, 2016 are reported here. The capability of interplanetary missions, carrying PDs for instrument performance purposes, in monitoring the passage of interplanetary coronal mass ejections (ICMEs) is also discussed.

In this paper, we analyzed a failed and a successful eruption that initiated from the same polarity inversion line within NOAA AR 11387 on December 25, 2011. They both started from a reconnection between sheared arcades, having distinct pre-eruption conditions and eruption details: before the failed one, the magnetic fields of the core region had a weaker non-potentiality; the external fields had a similar critical height for torus instability, a similar local torus-stable region, but a larger magnetic flux ratio (of low corona and near-surface region) as compared to the successful one. During the failed eruption, a smaller Lorentz force impulse was exerted on the outward ejecta; the ejecta had a much slower rising speed. Factors that might lead to the initiation of the failed eruption are identified: 1) a weaker non-potentiality of the core region, and a smaller Lorentz force impulse gave the ejecta a small momentum; 2) the large flux ratio, and the local torus-stable region in the corona provided strong confinements that made the erupting structure regain an equilibrium state.

We have measured the radial profiles of isophotal ellipticity ($\varepsilon$) and disky/boxy parameter A$_4$ out to radii of about three times the semi-major axes for $\sim4,600$ star-forming galaxies (SFGs) at intermediate redshifts $0.5<z<1.8$ in the CANDELS/GOODS-S and UDS fields. Based on the average size versus stellar-mass relation in each redshift bin, we divide our galaxies into Small SFGs (SSFGs), i.e., smaller than average for its mass, and Large SFGs (LSFGs), i.e., larger than average. We find that, at low masses ($M_{\ast} < 10^{10}M_{\odot}$), the SSFGs generally have nearly flat $\varepsilon$ and A$_4$ profiles for both edge-on and face-on views, especially at redshifts $z>1$. Moreover, the median A$_4$ values at all radii are almost zero. In contrast, the highly-inclined, low-mass LSFGs in the same mass-redshift bins generally have monotonically increasing $\varepsilon$ with radius and are dominated by disky values at intermediate radii. These findings at intermediate redshifts imply that low-mass SSFGs are not disk-like, while low-mass LSFGs appear to harbour disk-like components flattened by significant rotation. At high masses ($M_{\ast} > 10^{10}M_{\odot}$), highly-inclined SSFGs and LSFGs both exhibit a general, distinct trend for both $\varepsilon$ and A$_4$ profiles: increasing values with radius at lower radii, reaching maxima at intermediate radii, and then decreasing values at larger radii. Such a trend is more prevalent for more massive ($M_{\ast} > 10^{10.5}M_{\odot}$) galaxies or those at lower redshifts ($z<1.4$). The distinct trend in $\varepsilon$ and A$_4$ can be simply explained if galaxies possess all three components: central bulges, disks in the intermediate regions, and halo-like stellar components in the outskirts.

In this paper, we investigate two hydrogen-poor superluminous supernovae (SLSNe) iPTF15esb and iPTF13dcc whose light curves (LCs) show significant deviation from the smooth rise and fall. The LC of iPTF15esb exhibits two peaks and a post-peak plateau, and furthermore the late-time spectrum of iPTF15esb shows a strong, broad H$\alpha$ emission line. The early-time LC of iPTF13dcc shows a long duration bump followed by the second peak. Here we propose an ejecta-circumstellar medium (CSM) interaction model involving multiple shells/winds and use it to explain the LCs of iPTF15esb and iPTF13dcc. We find that the theoretical LCs reproduced by this model can well match the observations of iPTF15esb and iPTF13dcc. Based on our results, we infer that the progenitors have undergone multiple violent mass-loss processes before the SN explosion. In addition, we find that the variation trend of our inferred densities of the shells is consistent with that predicted by the stellar mass-loss history before an SN explosion. Further investigations for other bumpy SLSNe/SNe would shed light on their nature and provide a probe for the mass-loss history of their progenitors.

It is discussed that the formation of the spiral galaxies is driven by the cosmic background rotation, not a result of an isolated evolution proposed by the density wave theory. To analyze the motions of the galaxies, a simple double particle galaxy model is considered and the Coriolis force formed by the rotational background is introduced. The numerical analysis shows that not only the trajectory of the particle is the spiral shape, but also the relationship between the velocity and the radius reveals both the existence of spiral arm and the change of the arm number. In addition, the results of the three-dimensional simulation also give the warped structure of the spiral galaxies, and shows that the disc surface of the warped galaxy, like a spinning coins on the table, exists a whole overturning movement. Through the analysis, it can be concluded that the background environment of the spiral galaxies have a large-scale rotation, and both the formation and evolution of hurricane-like spiral galaxies are driven by this background rotation.

The search of fast radio burst (FRB) is a hot topic in current radio astronomy study. In this work, we carry out single pulse search for a VLBI pulsar observation data set using both auto spectrum and cross spectrum search method. The cross spectrum method is first proposed in Liu et al. (2018), which maximizes the signal power by fully utilizing the fringe phase information of the baseline cross spectrum. The auto spectrum search method is based on the popular pulsar software package PRESTO, which extracts single pulses from the auto spectrum of each station. According to our comparison, the cross spectrum method is able to enhance the signal power and therefore extract single pules from highly RFI contaminated data, which makes it possible to carry out FRB search in regular VLBI observations with the presence of RFIs.

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

In this paper, we present the detailed analysis of recurrent homologous jets originating from an emerging negative magnetic flux at the edge of an Active Region. The observed jets show multi-thermal features. Their evolution shows high consistence with the characteristic parameters of the emerging flux, suggesting that with more free magnetic energy, the eruptions tend to be more violent, frequent and blowout-like. The average temperature, average electron number density and axial speed are found to be similar for different jets, indicating that they should have been formed by plasmas from similar origins. Statistical analysis of the jets and their footpoint region conditions reveals a strong positive relationship between the footpoint-region total 131 {\AA} intensity enhancement and jets' length/width. Stronger linearly positive relationships also exist between the total intensity enhancement/thermal energy of the footpoint regions and jets' mass/kinetic/thermal energy, with higher cross-correlation coefficients. All the above results, together, confirm the direct relationship between the magnetic reconnection and the jets, and validate the important role of magnetic reconnection in transporting large amount of free magnetic energy into jets. It is also suggested that there should be more free energy released during the magnetic reconnection of blowout than of standard jet events.

Test mass charging caused by cosmic rays will be a significant source of acceleration noise for space-based gravitational wave detectors like LISA. Operating between December 2015 and July 2017, the technology demonstration mission LISA Pathfinder included a bespoke monitor to help characterise the relationship between test mass charging and the local radiation environment. The radiation monitor made in situ measurements of the cosmic ray flux while also providing information about its energy spectrum. We describe the monitor and present measurements which show a gradual 40% increase in count rate coinciding with the declining phase of the solar cycle. Modulations of up to 10% were also observed with periods of 13 and 26 days that are associated with co-rotating interaction regions and heliospheric current sheet crossings. These variations in the flux above the monitor detection threshold (approximately 70 MeV) are shown to be coherent with measurements made by the IREM monitor on-board the Earth orbiting INTEGRAL spacecraft. Finally we use the measured deposited energy spectra, in combination with a GEANT4 model, to estimate the galactic cosmic ray differential energy spectrum over the course of the mission.

We introduce the cross spectrum based FRB (Fast Radio Burst) search method for VLBI observation. This method optimizes the fringe fitting scheme in geodetic VLBI data post processing, which fully utilizes the cross spectrum fringe phase information and therefore maximizes the power of single pulse signals. Working with cross spectrum greatly reduces the effect of radio frequency interference (RFI) compared with using auto spectrum. Single pulse detection confidence increases by cross identifying detections from multiple baselines. By combining the power of multiple baselines, we may improve the detection sensitivity. Our method is similar to that of coherent beam forming, but without the computational expense to form a great number of beams to cover the whole field of view of our telescopes. The data processing pipeline designed for this method is easy to implement and parallelize, which can be deployed in various kinds of VLBI observations. In particular, we point out that VGOS observations are very suitable for FRB search.

Broad-lined type Ic supernovae (SNe Ic-BL) are a subclass of rare core collapse SNe whose energy source is debated in the literature. Recently a series of investigations on SNe Ic-BL with the magnetar (plus 56Ni) model were carried out. Evidence for magnetar formation was found for the well-observed SNe Ic-BL 1998bw and 2002ap. In this paper we systematically study a large sample of SNe Ic-BL not associated with gamma-ray bursts. We use photospheric velocity data determined in a homogeneous way. We find that the magnetar+56Ni model provides a good description of the light curves and velocity evolution of our sample of SNe Ic-BL, although some SNe (not all) can also be described by the pure-magnetar model or by the two-component pure-56Ni model (3 out of 12 are unlikely explained by two-component model). In the magnetar+56Ni model, the amount of 56Ni required to explain their luminosity is significantly reduced, and the derived initial explosion energy is, in general, in accordance with neutrino heating. Some correlations between different physical parameters are evaluated and their implications regarding magnetic field amplification and the total energy reservoir are discussed.

W UMa-type contact binaries belong to close binary systems whose components exactly overflow their Roche lobes and share a common convective envelope (CCE). In the last twenty years, the long-term variations of their orbital periods have been thought to depend on several mechanisms. Now, we suggest a new mechanism: CCE-dominated mechanism. The CCE-dominated mechanism is found based on our numerical result, especially at high mass ratios, that the orbital periods (P) of contact binaries change very much with their fill-out factors (f). Because f is taken as a measurement of the thickness of CCE, the physical cause for the variation of P is a mass transfer between CCE and components. Further, an f-dominated simplification model for this mechanism is introduced. According to it, P may change in a long-term oscillation way with a similar time scale of the thermal modulation, meanwhile q is decreasing slowly till the two components merge. It could be also applied to explain the presence of extremely short period, high mass ratio and deep contact binaries. Moreover, the CCE-dominated mechanism should always work due to mass transfer and mass loss both occurring via CCE. Therefor, the effect of CCE on the variations of orbital periods may have been underestimated before.

We carried out high-precision photometric observations of three eclipsing ultrashortperiod contact binaries (USPCBs). Theoretical models were fitted to the light-curves by means of the Wilson-Devinney code. The solutions suggest that the three targets have evolved to a contact phase. The photometric results are as follows: a) 1SWASP J030749.87-365201.7, q=0.439\pm0.003, f=0.0\pm3.6%; b) 1SWASP J213252.93-441822.6, q=0.560\pm0.003, f=14.2\pm1.9%; c) 1SWASP J200059.78+054408.9, q=0.436\pm0.008, f=58.4\pm1.8%. The light curves show O'Connell effects, which can be modeled by assumed cool spots. The cool spots models are strongly supported by the night-to-night variations in the I-band light curves of 1SWASP J030749.87-365201.7. For a comparative study, we collected the whole set of 28 well-studied USPCBs with P < 0.24 day. Thus, we found that most of them (17 of 28) are in shallow contact (i.e. fill-out factors f<20%). Only 4 USPCBs have deep fill-out factors (i.e. f>50%). Generally, contact binaries with deep fill-out factors are going to merge, but it is believed that USPCBs have just evolved to a contact phase. Hence, the deep USPCB 1SWASP J200059.78+054408.9 seems to be a contradiction, making it very interesting. Particularly, 1SWASP J030749.87-365201.7 is a zero contact binary within thermal equilibrium, implying that it should be a turn-off sample as predicted by the thermal relaxation oscillation (TRO) theory.

Very recently, the gravitational-wave (GW) event GW170817 was discovered to be associated with the short gamma-ray burst (GRB) 170817A. Multi-wavelength follow-up observations were carried out, and X-ray, optical, and radio counterparts to GW170817 were detected. The observations undoubtedly indicate that GRB 170817A originates from a binary neutron star merger. However, the GRB falls into the low-luminosity class that could have a higher statistical occurrence rate and detection probability than the normal (high-luminosity) class. This implies the possibility that GRB 170817A is intrinsically powerful, but we are off-axis and only observe its side emission. In this Letter, we provide a timely modeling of the multi-wavelength afterglow emission from this GRB and the associated kilonova signal from the merger ejecta, under the assumption of a structured jet, a two-component jet, and an intrinsically less-energetic quasi-isotropic fireball, respectively. Comparing the afterglow properties with the multi-wavelength follow-up observations, we can distinguish between these three models. Furthermore, a few model parameters (e.g., the ejecta mass and velocity) can be constrained.

We study the most luminous known supernova (SN) associated with a gamma-ray burst (GRB), SN 2011kl. The photospheric velocity of SN 2011kl around peak brightness is $21,000\pm7,000$ km s$^{-1}$. Owing to different assumptions related to the light-curve (LC) evolution (broken or unbroken power-law function) of the optical afterglow of GRB 111209A, different techniques for the LC decomposition, and different methods (with or without a near-infrared contribution), three groups derived three different bolometric LCs for SN 2011kl. Previous studies have shown that the LCs without an early-time excess preferred a magnetar model, a magnetar+$^{56}$Ni model, or a white dwarf tidal disruption event model rather than the radioactive heating model. On the other hand, the LC shows an early-time excess and dip that cannot be reproduced by the aforementioned models, and hence the blue-supergiant model was proposed to explain it. Here we reinvestigate the energy sources powering SN 2011kl. We find that the two LCs without the early-time excess of SN 2011kl can be explained by the magnetar+$^{56}$Ni model, and the LC showing the early excess can be explained by the magnetar+$^{56}$Ni model taking into account the cooling emission from the shock-heated envelope of the SN progenitor, demonstrating that this SN might primarily be powered by a nascent magnetar.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project we present an estimate of the mass of the supermassive black hole (SMBH) in the nearby fast-rotating early-type galaxy NGC4429, that is barred and has a boxy/peanut-shaped bulge. This estimate is based on Atacama Large Millimeter/submillimeter Array (ALMA) cycle-2 observations of the 12CO(3-2) emission line with a linear resolution of ~13 pc (0"18x0"14). NGC4429 has a relaxed, flocculent nuclear disc of molecular gas that is truncated at small radii, likely due to the combined effects of gas stability and tidal shear. The warm/dense 12CO(3-2) emitting gas is confined to the inner parts of this disc, likely again because the gas becomes more stable at larger radii, preventing star formation. The gas disc has a low velocity dispersion of 2.2$^{+0.68}_{-0.65}$ km/s. Despite the inner truncation of the gas disc, we are able to model the kinematics of the gas and estimate a mass of (1.5$\pm0.1^{+0.15}_{-0.35}$)$\times$10$^8$ Msun for the SMBH in NGC4429 (where the quoted uncertainties reflect the random and systematic uncertainties, respectively), consistent with a previous upper limit set using ionised gas kinematics. We confirm that the V-band mass-to-light ratio changes by ~30% within the inner 400 pc of NGC4429, as suggested by other authors. This SMBH mass measurement based on molecular gas kinematics, the sixth presented in the literature, once again demonstrates the power of ALMA to constrain SMBH masses.

We demonstrate a path to hitherto unachievable differential photometric precisions from the ground, both in the optical and near-infrared (NIR), using custom-fabricated beam-shaping diffusers produced using specialized nanofabrication techniques. Such diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. This PSF reshaping significantly increases the achievable dynamic range of our observations, increasing our observing efficiency and thus better averages over scintillation. Diffusers work in both collimated and converging beams. We present diffuser-assisted optical observations demonstrating $62^{+26}_{-16}$ppm precision in 30 minute bins on a nearby bright star 16-Cygni A (V=5.95) using the ARC 3.5m telescope---within a factor of $\sim$2 of Kepler's photometric precision on the same star. We also show a transit of WASP-85-Ab (V=11.2) and TRES-3b (V=12.4), where the residuals bin down to $180^{+66}_{-41}$ppm in 30 minute bins for WASP-85-Ab---a factor of $\sim$4 of the precision achieved by the K2 mission on this target---and to 101ppm for TRES-3b. In the NIR, where diffusers may provide even more significant improvements over the current state of the art, our preliminary tests have demonstrated $137^{+64}_{-36}$ppm precision for a $K_S =10.8$ star on the 200" Hale Telescope. These photometric precisions match or surpass the expected photometric precisions of TESS for the same magnitude range. This technology is inexpensive, scalable, easily adaptable, and can have an important and immediate impact on the observations of transits and secondary eclipses of exoplanets.

A binary neutron star (BNS) merger has been widely argued to be one of the progenitors of a short gamma-ray burst (SGRB). This central engine can be verified if its gravitational-wave (GW) event is detected simultaneously. Once confirmed, this kind of association will be a landmark in multi-messenger astronomy and will greatly enhance our understanding of the BNS merger processes. Due to the limited detection horizon of BNS mergers for the advanced LIGO/Virgo GW observatories, we are inclined to local SGRBs within few hundreds of mega-parsecs. Since normal SGRBs rarely fall into such a close range, to make it more observationally valuable, we have to focus on low-luminosity SGRBs which have a higher statistical occurrence rate and detection probability. However, there is a possibility that an observed low-luminosity SGRB is intrinsically powerful but we are off-axis and only observe its side emission. In this paper, we provide some theoretical predictions of both the off-axis afterglow emission from a nearby SGRB under the assumption of a structured jet and the macronova signal from the ejecta of this GW-detectable BNS merger. From the properties of the afterglow emission, we could distinguish an off-axis normal SGRB from an intrinsically low-energy quasi-isotropic class. Furthermore, with follow-up multi-wavelength observations, a few parameters for BNS mergers (e.g. the medium density and the ejecta mass and velocity) would be constrained.

We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H$_2$O data using a state-of-the-art 3D radiative transfer code which includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust- and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: A warm ($T_{\rm dust}\,\sim\,40-70\,{\rm K}$), dense ($n({\rm H})\,\sim\,10^5-10^6\,{\rm cm^{-3}}$) phase which dominates the emission of medium-excitation H$_2$O lines. This gas phase also dominates the FIR emission and the CO intensities for $J_{\rm up} > 8$. In addition a cold ($T_{\rm dust}\,\sim\,20-30\,{\rm K}$), dense ($n({\rm H})\sim\,10^4- 10^5\,{\rm cm^{-3}}$) more extended phase is present. It outputs the emission in the low-excitation H$_2$O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R$_s\,\le\,100$ pc) region is present which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with $E_{\rm up}\le300\,{\rm K}$ and $E_{\rm up}\le800\,{\rm K}$ in the warm and hot component, respectively. Higher energy levels are mainly excited by IR pumping.

The rest-frame UV-optical (i.e., $NUV-B$) color is sensitive to both low-level recent star formation (specific star formation rate - sSFR) and dust. In this Letter, we extend our previous work on the origins of $NUV-B$ color gradients in star-forming galaxies (SFGs) at $z\sim1$ to those at $z\sim2$. We use a sample of 1335 large (semi-major axis radius $R_{\rm SMA}>0.''18$) SFGs with extended UV emission out to $2R_{\rm SMA}$ in the mass range $M_{\ast} = 10^{9}-10^{11}M_{\odot}$ at $1.5<z<2.8$ in the CANDELS/GOODS-S and UDS fields. We show that these SFGs generally have negative $NUV-B$ color gradients (redder centres), and their color gradients strongly increase with galaxy mass. We also show that the global rest-frame $FUV-NUV$ color is approximately linear with $A_{\rm V}$, which is derived by modeling the observed integrated FUV to NIR spectral energy distributions of the galaxies. Applying this integrated calibration to our spatially-resolved data, we find a negative dust gradient (more dust extinguished in the centers), which steadily becomes steeper with galaxy mass. We further find that the $NUV-B$ color gradients become nearly zero after correcting for dust gradients regardless of galaxy mass. This indicates that the sSFR gradients are negligible and dust reddening is likely the principal cause of negative UV-optical color gradients in these SFGs. Our findings support that the buildup of the stellar mass in SFGs at the Cosmic Noon is self-similar inside $2R_{\rm SMA}$.

In this paper, we identified the magnetic source locations of 142 quasi- homologous (QH) coronal mass ejections (CMEs), of which 121 are from solar cycle (SC) 23 and 21 from SC 24. Among those CMEs, 63% originated from the same source location as their predecessor (defined as S-type), while 37% origi- nated from a different location within the same active region as their predecessor (defined as D-type). Their distinctly different waiting time distribution, peaking around 7.5 and 1.5 hours for S- and D-type CMEs, suggests that they might involve different physical mechanisms with different characteristic time scales. Through detailed analysis based on non-linear force free (NLFF) coronal mag- netic field modeling of two exemplary cases, we propose that the S-type QH CMES might involve a recurring energy release process from the same source location (by magnetic free energy replenishment), whereas the D-type QH CMEs can happen when a flux tube system disturbed by a nearby CME.

Recently, a hydrogen-poor superluminous supernova (SLSN) iPTF15esb at redshift $z=0.224$ was reported, whose light curve (LC) and spectrum show several unusual characteristics. Its late-time spectrum shows a strong, broad H$\alpha $ emission line and the bolometric LC exhibits two peaks and a post-peak plateau. Here we propose an ejecta-circumstellar medium (CSM) interaction model involving multiple shells and winds to explain this double-peak SN. We find that the theoretical LC reproduced by this model can well match the observations of iPTF15esb. Based on this result, we infer that the progenitor has undergone at least three violent mass-loss processes before the SN explosion. Furthermore, we find that the masses of the CSM wind and shells from the outermost shell (the first eruption) to the innermost wind (the final eruption) decrease sequentially but their densities increase. The variation trend of the inferred densities of the shells and wind is consistent with the stellar structure before an SN explosion. Further investigations for similar SLSNe would provide a probe for the mass-loss history of their progenitors.

In this paper we collect 19 hydrogen-deficient superluminous supernovae (SLSNe) and fit their light curves, temperature evolution, and velocity evolution based on the magnetar-powered model. To obtain the best-fitting parameters, we incorporate the Markov Chain Monte Carlo approach. We get rather good fits for 7 events ($\chi^2$/d.o.f = $0.24-0.96$) and good fits for other 7 events ($\chi^2$/d.o.f = $1.37-3.13$). We find that the initial periods ($P_{0}$) and magnetic strength ($B_{p}$) of the magnetars supposed to power these SLSNe are in the range of $\sim 1.2-8.3$ ms and $\sim (0.2-8.8)\times 10^{14}\,$G, respectively; the inferred masses of the ejecta of these SLSNe are between 1 and 27.6 $M_{\odot }$, and the values of the gamma-ray opacity $\kappa_{\gamma }$ are between 0.01 and 0.82 cm$^2$ ~g$^{-1}$. We also calculate the fraction of the initial rotational energy of the magnetars harbored in the centers of the remnants of these SLSNe that is converted to the kinetic energy of the ejecta and find that the fraction is $\sim 19-97$ for different values of $P_{0}$ and $B_{p}$, indicating that the acceleration effect cannot be neglected. Moreover, we find that the initial kinetic energies of most of these SLSNe are so small ($\lesssim 2\times 10^{51}$ erg) that they can be easily explained by the neutrino-driven mechanism. These results can help clarify some important issues related to the energy-source mechanisms and explosion mechanisms and reveal the nature of SLSNe.

We report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

Broad-lined type Ic supernovae (SNe Ic-BL) are peculiar stellar explosions that distinguish themselves from ordinary SNe. Some SNe Ic-BL are associated with long-duration (\gtrsim 2 s) gamma-ray bursts (GRBs). Black holes and magnetars are two types of compact objects that are hypothesized to be central engines of GRBs. In spite of decades of investigations, no direct evidence for the formation of black holes or magnetars has been found for GRBs so far. Here we report the finding that the early peak (t \lesssim 50 days) and late-time (t \gtrsim 300 days) slow decay displayed in the light curves of both SNe 1998bw (associated with GRB 980425) and 2002ap (not GRB-associated) can be attributed to magnetar spin-down with initial rotation period P0 \sim 20 ms, while the intermediate-time (50 \lesssim t \lesssim 300 days) linear decline is caused by radioactive decay of 56Ni. The connection between the early peak and late-time slow decline in the light curves is unexpected in alternative models. We thus suggest that GRB 980425 and SN 2002ap were powered by magnetars.

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

Solar active regions (ARs) are the major sources of two kinds of the most violent solar eruptions, namely flares and coronal mass ejections (CMEs). The largest AR in the past 24 years, NOAA AR 12192, crossed the visible disk from 2014 October 17 to 30, unusually produced more than one hundred flares, including 32 M-class and 6 X-class ones, but only one small CME. Flares and CMEs are believed to be two phenomena in the same eruptive process. Why is such a flare-rich AR so CME-poor? We compared this AR with other four ARs; two were productive in both and two were inert. The investigation of the photospheric parameters based on the SDO/HMI vector magnetogram reveals that the flare-rich AR 12192, as the other two productive ARs, has larger magnetic flux, current and free magnetic energy than the two inert ARs, but contrast to the two productive ARs, it has no strong, concentrated current helicity along both sides of the flaring neutral line, indicating the absence of a mature magnetic structure consisting of highly sheared or twisted field lines. Furthermore, the decay index above the AR 12192 is relatively low, showing strong constraint. These results suggest that productive ARs are always large and have enough current and free energy to power flares, but whether or not a flare is accompanied by a CME is seemingly related to (1) if there is mature sheared or twisted core field serving as the seed of the CME, (2) if the constraint of the overlying arcades is weak enough.

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are $\lesssim -21$ mag can be explained by the $^{56}$Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta-circumstellar medium (CSM) interaction model. Recently, iPTF13ehe challenges these energy-source models, because the spectral analysis shows that $\sim 2.5M_\odot$ of $^{56}$Ni have been synthesized but are inadequate to power the peak bolometric emission of iPTF13ehe, while the rebrightening of the late-time light-curve (LC) and the H$\alpha$ emission lines indicate that the ejecta-CSM interaction must play a key role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar together with some amount ($\lesssim 2.5M_\odot$) of $^{56}$Ni may power the early LC of iPTF13ehe while the late-time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we suggest that iPTF13ehe is a genuine core-collapse supernova rather than a pulsational pair-instability supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on their explosion and energy-source mechanisms.

(Abridged) We present Herschel/HIFI spectra of the H2O 1113 GHz and H2O+ 1115 GHz lines toward five nearby prototypical starburst/AGN systems, and OH+ 971 GHz spectra toward three of these. The beam size of 20" corresponds to resolutions between 0.35 and 7 kpc. The observed line profiles range from pure absorption (NGC 4945, M82) to P-Cygni indicating outflow (NGC 253, Arp 220) and inverse P-Cygni indicating infall (Cen A). The similarity of the H2O, OH+, and H2O+ profiles to each other and to HI indicates that diffuse and dense gas phases are well mixed. We estimate column densities assuming negligible excitation (for absorption features) and using a non-LTE model (for emission features), adopting calculated collision data for H2O and OH+, and rough estimates for H2O+. Column densities range from ~10^13 to ~10^15 cm^-2 for each species, and are similar between absorption and emission components, indicating that the nuclear region does not contribute much to the emission in these ground-state lines. The N(H2O)/N(H2O+) ratios of 1.4-5.6 indicate an origin of the lines in diffuse gas, and the N(OH+)/N(H2O+) ratios of 1.6-3.1 indicate a low H2 fraction (~11%) in the gas. Adopting recent Galactic values for the average gas density and the ionization efficiency, we find ionization rates for our sample galaxies of ~3x10^-16 s^-1 which are similar to the value for the Galactic disk, but ~10x below that of the Galactic Center and ~100x below estimates for AGN from excited-state H3O+ lines. We conclude that the ground-state lines of water and its associated ions probe primarily non-nuclear gas in the disks of these centrally active galaxies. Our data thus provide evidence for a decrease in ionization rate by a factor of ~10 from the nuclei to the disks of galaxies, as found before for the Milky Way.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

The study of short-duration gamma-ray bursts provides growing evidence that a good fraction of double neutron star mergers lead to the formation of stable millisecond magnetars. The launch of Poynting flux by the millisecond magnetars could leave distinct electromagnetic signatures that reveal the energy dissipation processes in the magnetar wind. In previous studies (Wang & Dai 2013b; Wang et al. 2015), we assume that the magnetar wind becomes completely lepton-dominated so that electrons/positrons in the magnetar wind are accelerated by a diffusive shock. However, theoretical modeling of pulsar wind nebulae shows that in many cases the magnetic field energy in the pulsar wind may be strong enough to suppress diffusive shock acceleration. In this paper, we investigate the reverse shock emission as well as the forward shock emission with an arbitrary magnetization parameter $\sigma$ of a magnetar wind. We find that the reverse shock emission strongly depends on $\sigma$, and in particular, $\sigma \sim 0.3$ leads to the strongest reverse shock emission. Future observations would be helpful to diagnose the composition of the magnetar wind.

There is growing evidence that a stable magnetar could be formed from the coalescence of double neutron stars. In previous papers, we investigated the signature of formation of stable millisecond magnetars in radio and optical/ultraviolet bands by assuming that the central rapidly rotating magnetar deposits its rotational energy in the form of a relativistic leptonized wind. We found that the optical transient PTF11agg could be the first evidence for the formation of post-merger millisecond magnetars. To enhance the probability of finding more evidence for the post-merger magnetar formation, it is better to extend the observational channel to other photon energy bands. In this paper we propose to search the signature of post-merger magnetar formation in X-ray and especially gamma-ray bands. We calculate the SSC emission of the reverse shock powered by post-merger millisecond magnetars. We find that the SSC component peaks at $1\,{\rm GeV}$ in the spectral energy distribution and extends to $\gtrsim 10\,{\rm TeV}$ for typical parameters. These energy bands are quite suitable for Fermi/LAT and CTA, which, with their current observational sensitivities, can detect the SSC emission powered by post-merger magnetars up to $1\,{\rm Gpc}$. NuSTAR, sensible in X-ray bands, can detect the formation of post-merger millisecond magnetars at redshift $z\sim 1$. Future improvement in sensitivity of CTA can also probe the birth of post-merger millisecond magnetars at redshift $z\sim 1$. However, because of the $\gamma$-$\gamma$ collisions, strong high-energy emission is clearly predicted only for ejecta masses lower than $10^{-3}M_\odot$.

Recent X-ray observations of star-forming galaxies such as M82 have shown the Ly beta/Ly alpha line ratio of Ne X to be in excess of predictions for thermal electron impact excitation. Here we demonstrate that the observed line ratio may be due to charge exchange and can be used to constrain the ion kinetic energy to be <500 eV/u. This is accomplished by computing spectra and line ratios via a range of theoretical methods and comparing these to experiments with He over astrophysically relevant collision energies. The charge exchange emission spectra calculations were performed for Ne[10+] +H and Ne[10+] +He using widely applied approaches including the atomic orbital close coupling, classical trajectory Monte Carlo, and multichannel Landau- Zener (MCLZ) methods. A comparison of the results from these methods indicates that for the considered energy range and neutrals (H, He) the so-called "low-energy l-distribution" MCLZ method provides the most likely reliable predictions.

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

Atmospheric optical turbulence seriously limits the performance of high angular resolution instruments. An 8-night campaign of measurements was carried out at the LAMOST site in 2011, to characterize the optical turbulence. Two instruments were set up during the campaign: a Differential Image Motion Monitor (DIMM) used to measure the total atmospheric seeing, and a Single Star Scidar (SSS) to measure the vertical profiles of the turbulence C_n^2(h) and the horizontal wind velocity V(h). The optical turbulence parameters are also calculated with the Weather Research and Forecasting (WRF) model coupled with the Trinquet-Vernin model, which describes optical effects of atmospheric turbulence by using the local meteorological parameters. This paper presents assessment of the optical parameters involved in high angular resolution astronomy. Its includes seeing, isoplanatic angle, coherence time, coherence etendue, vertical profiles of optical turbulence intensity _n^2(h)$ and horizontal wind speed V(h). The median seeing is respectively 1.01 arcsec, 1.17 arcsec and 1.07arcsec as measured with the DIMM, the SSS and predicted with WRF model. The history of seeing measurements at the LAMOST site are reviewed, and the turbulence measurements in this campaign are compared with other astronomical observatories in the world.

To better characterize the global star formation (SF) activity in a galaxy, one needs to know not only the star formation rate (SFR) but also the rest-frame, far-infrared (FIR) color (e.g., the 60-to-100 $\mu$m color, $C(60/100)$] of the dust emission. The latter probes the average intensity of the dust heating radiation field and scales statistically with the effective SFR surface density in star-forming galaxies including (ultra-)luminous infrared galaxies [(U)LIRGs]. To this end, we exploit here a new spectroscopic approach involving only two emission lines: CO\,(7$-$6) at 372 $\mu$m and [NII] at 205 $\mu$m. For local (U)LIRGs, the ratios of the CO (7$-$6) luminosity ($L_{\rm CO\,(7-6)}$) to the total infrared luminosity ($L_{\rm IR}$; 8$-$1000 $\mu$m) are fairly tightly distributed (to within $\sim$0.12 dex) and show little dependence on $C(60/100)$. This makes $L_{\rm CO\,(7-6)}$ a good SFR tracer, which is less contaminated by active galactic nuclei (AGN) than $L_{\rm IR}$ and may also be much less sensitive to metallicity than $L_{\rm CO\,(1-0)}$. Furthermore, the logarithmic [NII] 205 $\mu$m to CO (7$-$6) luminosity ratio is fairly steeply (at a slope of $\sim$$-1.4$) correlated with $C(60/100)$, with a modest scatter ($\sim$0.23 dex). This makes it a useful estimator on $C(60/100)$ with an implied uncertainty of $\sim$0.15 [or $\lesssim$4 K in the dust temperature ($T_{\rm dust}$) in the case of a graybody emission with $T_{\rm dust} \gtrsim 30$ K and a dust emissivity index $\beta \ge 1$]. Our locally calibrated SFR and $C(60/100)$ estimators are shown to be consistent with the published data of (U)LIRGs of $z$ up to $\sim$6.5.

We study the global SF law - the relation between gas and SFRs in a sample of 181 local galaxies with L_IR spanning almost five orders of magnitude, which includes 115 normal galaxies and 66 (U)LIRGs. We derive their atomic, molecular gas and dense molecular gas masses using newly available HI, CO and HCN data from the literature, and SFRs are determined both from total IR and 1.4 GHz radio continuum (RC) luminosities. In order to derive the disk-averaged surface densities of gas and SFRs, we have used high-resolution RC observations to measure the radio sizes for all galaxies. We find that dense molecular gas (as traced by HCN) has the tightest correlation with that of SFRs, and is linear in (N=1.01 +/- 0.02) across the full galaxy sample. The correlation between densities of molecular gas (traced by CO) and SFRs is sensitive to the adopted value of the alpha_CO used to infer molecular gas masses from CO luminosities. For a fixed value of alpha_CO, a slope of 1.14+/-0.02 is found. If instead we adopt values of 4.6 and 0.8 for disk galaxies and (U)LIRGs, respectively, we find the two distinct relations. If applying a continuously varying alpha_CO to our sample, we recover a single relation with slope of 1.60+/-0.03. The SFRs is a steeper function of total gas than that of molecular gas, and is tighter among low-luminosity galaxies. We find no correlation between SFRs and atomic gas.

The group velocity of light in material around the AGN jet is acquiescently one (c as a unit), but this is only a hypothesis. Here, we re-derive apparent superluminal and Doppler formulas for the general case (it is assumed that the group velocity of light in the uniform and isotropic medium around a jet (a beaming model) is not necessarily equal to one, e.g., Araudo et al. (2010) thought that there may be dense clouds around AGN jet base), and show that the group velocity of light close to one could seriously affect apparent superluminal phenomena and Doppler effect in the AGN jet (when the viewing angle and Lorentz factor take some appropriate values).

New multi-color light curves of the very short period K-type eclipsing binary V1799 Ori were obtained and analyzed with the W-D code. The photometric solutions reveal that the system is a W-type shallow-contact binary with a mass ratio of $q=1.335(\pm0.005)$ and a degree of contact about $f = 3.5(\pm1.1)\%$. In general, the results are in good agreement with which is reported by Samec. The remarkable O'Connell effects in the light curves are well explained by employing star spots on the binary surface, which confirms that the system is active at present. Several new times of light minimum were obtained. All the available times of light minimum were collected, along with the recalculated and new obtained. Applying a least-squares method to the constructed O-C diagram, a new ephemeris was derived for V1799 Ori. The orbital period is found to show a continuous weak increase at a rate of $1.8(\pm0.6)\times10^{-8}$ days$\cdot$yr$^{-1}$. The extreme-shallow contact, together with the period increase, suggests that the binary may be at a critical stage predicted by the TRO theory. \keywords{binaries : close -- binaries : eclipsing -- stars: individual (V1799 Ori)

Usually, in order to investigate the evolution of a theory, one may find the critical points of the system and then perform perturbations around these critical points to see whether they are stable or not. This local method is very useful when the initial values of the dynamical variables are not far away from the critical points. Essentially, the nonlinear effects are totally neglected in such kind of approach. Therefore, one can not tell whether the dynamical system will evolute to the stable critical points or not when the initial values of the variables do not close enough to these critical points. Furthermore, when there are two or more stable critical points in the system, local analysis can not provide the informations that which one the system will finally evolute to. In this paper, we have further developed the nullcline method to study the bifurcation phenomenon and global dynamical behaviour of the $f(T)$ theory. We overcome the shortcoming of local analysis. And it is very clear to see the evolution of the system under any initial conditions.

NGC 7789 is an intermediate-age open cluster with an age that is similar to the mean age of contact binary stars. V12 is a bright UMa-type binary star with an orbital period of 0.397 day. The first complete light curves of V12 in V, R, and I bands are presented and analyzed with the W-D (Wilson and Devinney) method. The results shows that V12 is an intermediate-contact binary (f = 43%) with a mass ratio of 3.848, and it is a W-type contact binary where the less massive component is slightly hotter than the more massive one. The asymmetries of the light curves are explained by the presence of a dark spot on the more massive component. The derived orbital inclination (i = 83.6o) indicates that it is a total eclipsing binary, which suggests that the determined parameters are reliable. The orbital period may show a long-term increase at a rate of P = +2.48 x 10-6 days/year that reveals a rapid mass transfer from the less massive component to the more massive one. However, more observations are needed to confirm this conclusion. The presence of an intermediate-contact binary in an intermediate-age open cluster may suggest that some contact binaries have a very short pre-contact timescale. The presence of a third body or/and stellar collision may help to shorten the pre-contact evolution.

As indicated by the observed X-ray flares, a great amount of energy could be intermediately released from the postburst central engine of gamma-ray bursts (GRBs). As a natural consequence, the GRB external shock could be energized over and over. With such a multiple energy injection model, we explore the unique X-ray afterglow light curve of GRB 050712, which exists four apparent shallow decay plateaus. Together with three early X-ray flares, the central engine of GRB 050712 is supposed to release energy at least seven times after the burst. Furthermore we find that the energy released during four plateaus are all on the same order of magnitude, but the luminosity decreases with time significantly. These results may provide some interesting implications for the GRB central engine.

A Single Star Scidar system(SSS) has been developed for remotely sensing atmospheric turbulence profiles. The SSS consists of computing the spatial auto/cross-correlation functions of short exposure images of the scintillation patterns produced by a single star, and provides the vertical profiles of optical turbulence intensity C2n(h) and wind speed V(h). The SSS needs only a 40 cm aperture telescope, so that can be portable and equipped easily to field candidate sites. Some experiments for the SSS have been made in Beijing last year, successfully retrieving atmospheric turbulence and wind profiles from the ground to 30 km. The SSS observations has recently been made at the Xinglong station of NAOC, characterizing atmospheric parameters at this station. We plan to automatize SSS instrument and run remote observation via internet; a more friendly auto-SSS system will be set up and make use at the candidate sites in Tibet and Dome A.

Why and how may some active regions (ARs) frequently produce coronal mass ejections (CMEs)? It is one of the key questions to deepen our understanding of the mechanisms and processes of energy accumulation and sudden release in ARs and to improve our capability of space weather prediction. Although some case studies have been made, the question is still far from fully answered. This issue is now being tried to address statistically through an investigation of waiting times of quasi-homologous CMEs from super ARs in solar cycle 23. It is found that the waiting times of quasi-homologous CMEs have a two-component distribution with a separation at about 18 hours. The first component is a Gaussian-like distribution with a peak at about 7 hours, which indicates a tight physical connection between these quasi-homologous CMEs. The likelihood of occurrences of two or more CMEs faster than 1200 km s-1 from the same AR within 18 hours is about 20%. Furthermore, the correlation analysis among CME waiting times, CME speeds and CME occurrence rates reveals that these quantities are independent to each other, suggesting that the perturbation by preceding CMEs rather than free energy input be the direct cause of quasi-homologous CMEs. The peak waiting time of 7 hours probably characterize the time scale of the growth of instabilities triggered by preceding CMEs. This study uncovers more clues from a statistical perspective for us to understand quasi-homologous CMEs as well as CME-rich ARs.

We plan to install an infrared telescope at the new site of Tibet, China. The primary mirror diameter is 50cm, and the focal ratio F8. The Xenics 640\times512 near infrared camera is employed, equipped with a dedicated high speed InGaAs detector array, working up to 1.7{\mu}m. The new site is located on 5100m mountain, near Gar town, Ali, where is an excellent site for both infrared and submillimeter observations. The telescope will be remotely controlled through internet. The goal of IRT is to make site testing, detect variable stars, and search for extrasolar planets.

We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.

By using six new determined mid-eclipse times together with those collected from the literature, we found that the Observed-Calculated (O-C) curve of RR Cae shows a cyclic change with a period of 11.9 years and an amplitude of 14.3s, while it undergoes an upward parabolic variation (revealing a long-term period increase at a rate of dP/dt =+4.18(+-0.20)x10^(-12). The cyclic change was analyzed for the light-travel time effect that arises from the gravitational influence of a third companion. The mass of the third body was determined to be M_3*sin i' = 4.2(+-0.4) M_{Jup} suggesting that it is a circumbinary giant planet when its orbital inclination is larger than 17.6 degree. The orbital separation of the circumbinary planet from the central eclipsing binary is about 5.3(+-0.6)AU. The period increase is opposite to the changes caused by angular momentum loss via magnetic braking or/and gravitational radiation, nor can it be explained by the mass transfer between both components because of its detached configuration. These indicate that the observed upward parabolic change is only a part of a long-period (longer than 26.3 years) cyclic variation, which may reveal the presence of another giant circumbinary planet in a wide orbit.

Light transport in graded index media follows a curved trajectory determined by the Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.

We study the global star formation law - the relation between the gas and star formation rate (SFR) in a sample of 130 local galaxies with infrared (IR) luminosities spanning over three orders of magnitude (10^9-10^12 Lsun), which includes 91 normal spiral galaxies and 39 (ultra)luminous IR galaxies [(U)LIRGs]. We derive their total (atomic and molecular) gas and dense molecular gas masses using newly available HI, CO and HCN data from the literature. The SFR of galaxies is determined from total IR (8-1000 um) and 1.4 GHz radio continuum (RC) luminosities. The galaxy disk sizes are defined by the de-convolved elliptical Gaussian FWHM of the RC maps. We derive the galaxy disk-averaged SFRs and various gas surface densities, and investigate their relationship. We find that the galaxy disk-averaged surface densities of dense molecular gas mass has the tightest correlation with that of SFR (scatter ~ 0.26 dex), and is linear in log-log space (power-law slope of N=1.03 +/- 0.02) across the full galaxy sample. The correlation between the total gas and SFR surface densities for the full sample has a somewhat larger scatter (~ 0.48 dex), and is best fit by a power-law with slope 1.45 +/- 0.02. However, the slope changes from ~ 1 when only normal spirals are considered, to ~ 1.5 when more and more (U)LIRGs are included in the fitting. When different CO-to-H2 conversion factors are used to infer molecular gas masses for normal galaxies and (U)LIRGs, the bi-modal relations claimed recently in CO observations of high-redshift galaxies appear to also exist in local populations of star-forming galaxies.

We introduce a probabilistic approach to the problem of counting dwarf satellites around host galaxies in databases with limited redshift information. This technique is used to investigate the occurrence of satellites with luminosities similar to the Magellanic Clouds around hosts with properties similar to the Milky Way in the object catalog of the Sloan Digital Sky Survey. Our analysis uses data from SDSS Data Release 7, selecting candidate Milky-Way-like hosts from the spectroscopic catalog and candidate analogs of the Magellanic Clouds from the photometric catalog. Our principal result is the probability for a Milky-Way-like galaxy to host N_{sat} close satellites with luminosities similar to the Magellanic Clouds. We find that 81 percent of galaxies like the Milky Way are have no such satellites within a radius of 150 kpc, 11 percent have one, and only 3.5 percent of hosts have two. The probabilities are robust to changes in host and satellite selection criteria, background-estimation technique, and survey depth. These results demonstrate that the Milky Way has significantly more satellites than a typical galaxy of its luminosity; this fact is useful for understanding the larger cosmological context of our home galaxy.

Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.

The new discovered short-period close binary star, XY LMi, was monitored photometrically since 2006. It is shown that the light curves are typical EW-type and show complete eclipses with an eclipse duration of about 80 minutes. By analyzing the complete B, V, R, and I light curves with the 2003 version of the W-D code, photometric solutions were determined. It is discovered that XY LMi is a high fill-out, extreme mass ratio overcontact binary system with a mass ratio of q=0.148 and a fill-out factor of f=74.1%, suggesting that it is on the late evolutionary stage of late-type tidal-locked binary stars. As observed in other overcontact binary stars, evidence for the presence of two dark spots on both components are given. Based on our 19 epoches of eclipse times, it is found that the orbital period of the overcontact binary is decreasing continuously at a rate of dP/dt=-1.67\times10^{-7}\,days/year, which may be caused by the mass transfer from the primary to the secondary or/and angular momentum loss via magnetic stellar wind. The decrease of the orbital period may result in the increase of the fill-out, and finally, it will evolve into a single rapid-rotation star when the fluid surface reaching the outer critical Roche Lobe.

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.

The Sun affects the Earth's physical phenomena in multiple ways, in particular the material in interplanetary space comes from coronal expansion in the form of inhomogeneous plasma flow (solar wind), which is the primary source of the interplanetary medium. Ground-based Interplanetary Scintillation (IPS) observations are an important and effective method for measuring solar wind speed and the structures of small diameter radio sources. We discuss one mode of ground-based single-station observations: Single-Station Single-Frequency (SSSF) mode. To realize the SSSF mode, a new system has been established at Urumqi Astronomical Observatory (UAO), China, and a series of experimental observations were carried out successfully from May to December, 2008.

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.

The instability of current-driven ion-acoustic waves in the collisionless magnetic-field-free space plasma is investigated by using a nonextensive approach. The ions and the electrons are thought of in the power-law distributions that can be described by the generalized q-Maxwellian velocity distribution and are considered with the different nonextensive q-parameters. The generalized q-wave frequency and the generalized instability q-growth rate for the ion-acoustic waves are derived. The numerical results show that the nonextensive effects on the ion-acoustic waves are not apparent when the electron temperature is much more than the ion temperature, but they are salient when the electron temperature is not much more than the ion temperature. As compared with the electrons, the ions play a dominant role in the nonextensive effects.

HS 0705+6700 is a short-period (P=2.3 hours), close binary containing a hot sdB-type primary and a fully convective secondary. We have monitored this eclipsing binary for more than 2 years and as a result, 32 times of light minimum were obtained. Based on our new eclipse times together with these compiled from the literature, it is discovered that the O-C curve of HS 0705+6700 shows a cyclic variation with a period of 7.15 years and a semiamplitude of 92.4 s. The periodic change was analyzed for the light-travel time effect that may be due to the presence of a tertiary companion. The mass of the third body is determined to be M3 sin i = 0.0377 (+/-0.0043) Msun when a total mass of 0.617 Msun for HS 0705+6700 is adopted. For orbital inclinations i >= 32.8, the mass of the tertiary component would be below the stable hydrogen-burning limit of M3~0.072 Msun, and thus it would be a brown dwarf. The third body is orbiting the sdB-type binary at a distance shorter than 3.6 astronomical units (AU). HS 0705+6700 was formed through the evolution of a common envelope after the primary becomes a red giant. The detection of a sub-stellar companion in HS 0705+6700 system at this distance from the binary could give some constraints on stellar evolution in such systems and the interactions between red giants and their companions.

In measuring the power spectrum of the distribution of large numbers of dark matter particles in simulations, or galaxies in observations, one has to use Fast Fourier Transforms (FFT) for calculational efficiency. However, because of the required mass assignment onto grid points in this method, the measured power spectrum $\la |\delta^f(k)|^2\ra$ obtained with an FFT is not the true power spectrum $P(k)$ but instead one that is convolved with a window function $|W(\vec k)|^2$ in Fourier space. In a recent paper, Jing (2005) proposed an elegant algorithm to deconvolve the sampling effects of the window function and to extract the true power spectrum, and tests using N-body simulations show that this algorithm works very well for the three most commonly used mass assignment functions, i.e., the Nearest Grid Point (NGP), the Cloud In Cell (CIC) and the Triangular Shaped Cloud (TSC) methods. In this paper, rather than trying to deconvolve the sampling effects of the window function, we propose to select a particular function in performing the mass assignment that can minimize these effects. An ideal window function should fulfill the following criteria: (i) compact top-hat like support in Fourier space to minimize the sampling effects; (ii) compact support in real space to allow a fast and computationally feasible mass assignment onto grids. We find that the scale functions of Daubechies wavelet transformations are good candidates for such a purpose. Our tests using data from the Millennium Simulation show that the true power spectrum of dark matter can be accurately measured at a level better than 2% up to $k=0.7k_N$, without applying any deconvolution processes. The new scheme is especially valuable for measurements of higher order statistics, e.g. the bi-spectrum,........

The electrostatic charging of the test mass in ASTROD I (Astrodynamical Space Test of Relativity using Optical Devices I) mission can affect the quality of the science data as a result of spurious Coulomb and Lorentz forces. To estimate the size of the resultant disturbances, credible predictions of charging rates and the charging noise are required. Using the GEANT4 software toolkit, we present a detailed Monte Carlo simulation of the ASTROD I test mass charging due to exposure of the spacecraft to galactic cosmic-ray (GCR) protons and alpha particles (3He, 4He) in the space environment. A positive charging rate of 33.3 e+/s at solar minimum is obtained. This figure reduces by 50% at solar maximum. Based on this charging rate and factoring in the contribution of minor cosmic-ray components, we calculate the acceleration noise and stiffness associated with charging. We conclude that the acceleration noise arising from Coulomb and Lorentz effects are well below the ASTROD I acceleration noise limit at 0.1 mHz both at solar minimum and maximum. The coherent Fourier components due to charging are investigated, it needs to be studied carefully in order to ensure that these do not compromise the quality of science data in the ASTROD I mission.

We use a multi-dimensional hydrodynamics code to study the gravitational interaction between an embedded planet and a protoplanetary disk with emphasis on the generation of vortensity (potential vorticity) through a Baroclinic Instability and subsequent development of Rossby-Wave-Instability (RWI). It is found that the generation of potential vorticity is very common and effective in non-barotropic disks through the Baroclinic Instability, especially within the coorbital region. Our results also complement previous studies by \citet{KLL03} that non-axisymmetric RWIs are likely to develop at local minima of potential vorticity distribution that are generated by the interaction between a planet and a inviscid barotropic disk. This second instability appears to be very common and robust, regardless of the equation of state, initial density distribution, and rotational law of the disk. The development of RWIs results in non-axisymmetric density blobs, which exert stronger torques onto the planet when they travel in the vicinity of the planet. As a result of that, large amplitude oscillations are introduced to the time behavior of the total torque acted on the planet by the disk. In our current simulations, RWIs do not change the overall picture of inward orbital migration but bring in a non-monotonic behavior to the migration speed. As a side effect, RWIs also introduce interesting structures into the disk. These structures may help the formation of Earth-like planets in the Habitable Zone or Hot Earths interior to a close-in giant planet.

As ASTROD I travels through space, its test mass will accrue charge due to galactic cosmic-rays and solar energetic particles incident on the spacecraft. This test mass charge will result in Coulomb forces between the test mass and the surrounding electrodes. In earlier work using the GEANT4 toolkit, we predicted a net charging rate of nearly 9.0 +e/s from cosmic-ray protons between 0.1 and 1000 GeV at solar maximum, and rising to 26.5 +e/s at solar minimum. Here we use GEANT4 to simulate the charging process due to solar energetic particle events and to estimate the magnitude of acceleration noise due to this charging. The predicted charging rates range from 2840 to 64300 +e/s, at peak intensity, for the 4 largest SEP events in September and October 1989. For the 2 larger events, the acceleration disturbances due to charging exceeds the ASTROD I acceleration noise budget requirement. Continuous discharge should be considered for suppressing this charging noise. The acceleration noise during the 2 small events is well below the design target, although during these events, the net charging rate will be dominated by these solar fluxes.

ASTROD I is planned as a single spacecraft mission. It will use interferometric and pulse ranging techniques between the spacecraft and ground stations, to make high precision measurements of the parameters that describe the solar system, and to test relativistic gravity with improved accuracy. At the heart of the spacecraft is a test mass, which the spacecraft will follow using a drag-free control system. The mission critically depends on maintaining the geodesic motion of the test mass. Charging of the test mass due to cosmic rays and solar particles will disturb its geodesic motion. We have modelled the charging process using the GEANT4 toolkit and a simplified, geometrical model and estimate that the ASTROD I test mass will charge positively, at a rate of 24 +/-7 e+/s, due to cosmic ray protons and alpha particles (3He and 4He) at solar minimum. We have used the results of this simulation to estimate the magnitude of disturbances associated with test mass charging, for the worst-case scenario, taking into account uncertainties in the model and potential charging contributions from minor cosmic-ray components.

Some nonlinear dynamical techniques, including state-space reconstruction and correlation integral, are used to analyze the light curve of 3C 273. The result is compared with a chaotic model. The similarity between them suggests that there is a low-dimensional chaotic attractor in the light curve of 3C 273.

(Abridged)We numerically investigated the dynamical architecture of 47 UMa with the planetary configuration of the best-fit orbital solutions by Fischer et al. We systematically studied the existence of Earth-like planets in the region 0.05 AU $\leq a \leq 2.0$ AU for 47 UMa with numerical simulations, and we also explored the packed planetary geometry and Trojan planets in the system. In the simulations, we found that "hot Earths" at 0.05 AU $\leq a < $ 0.4 AU can dynamically survive at least for 1 Myr. The Earth-like planets can eventually remain in the system for 10 Myr in areas involved in the mean motion resonances (MMR) (e.g., 3:2 MMR) with the inner companion. Moreover, we showed that the 2:1 and 3:1 resonances are on the fringe of stability, while the 5:2 MMR is unstable. Additionally, the 2:1 MMR marks out a remarkable boundary between chaotic and regular motions, inside, most of the orbits can survive, outside, they are mostly lost in the orbital evolution. In a dynamical sense, the most likely candidate for habitable environment is Earth-like planets with orbits in the ranges 0.8 AU $\leq a < 1.0$ AU and 1.0 AU $ < a < 1.30$ AU (except several unstable cases) with relatively low eccentricities. The Trojan planets with low eccentricities and inclinations can secularly last at the triangular equilibrium points of two massive planets. Hence, the 47 UMa planetary system may be a close analog to our solar system.

In the mean field approximation of the relativistic $\sigma$-$\omega$-$\rho$ model, the magnetic fields are incorporated, and its influence on the properties of n-p-e neutron star matter are studied. When the strength of the magnetic field is weaker than $\sim 10^{18}$G, the particle fractions and chemical potentials, matter energy density and pressure hardly change with the magnetic field; when the strength of the magnetic field is stronger than $\sim 10^{20}$G, the above quantities change with the magnetic field evidently. Furthermore, the pressure is studied in both thermodynamics and hydrodynamics. The difference between these two ways exits in the high density region, that is, the thermal self-consistency may not be satisfied in this region if the magnetic field is considered.

We investigated the apsidal motion for the multi-planet systems. In the simulations, we found that the two planets of HD 37124, HD 12661, 47 Uma and HD 82943 separately undergo apsidal alignment or antialignment. But the companions of GJ 876 and $\upsilon$ And are only in apsidal lock about $0^{\circ}$. Moreover, we obtained the criteria with Laplace-Lagrange secular theory to discern whether a pair of planets for a certain system are in libration or circulation.

To date, more than 100 giant Jupiter-like planets have been discovered in Doppler surveys of solar-type stars. In this paper, we perform simulations to investigate three systems: GJ 876, HD 82943 and 55 Cnc. The former two systems both have a pair of planets in the 2:1 Mean Motion Resonance (MMR), while the inner two companions of the later is close to 3:1 MMR. By integrating hundreds of the planetary orbits of three systems for million years, we find that for GJ 876 and HD 82943, the critical argument $\lambda_{1} - 2\lambda_{2} + \varpi_{1}$ and $\lambda_{1} - 2\lambda_{2} + \varpi_{2}$ librate about $0^{\circ}$ or $180^{\circ}$, indicating 2:1 MMR can play an important role in stabilizing the motion of the planets so that they are protected from frequent close encounters. As for 55 Cnc, we further show the three resonant arguments for 3:1 MMR execute librations for millions of years respectively, which reveals the evidence of the resonance for this system. Additionally, we should emphasize another vital mechanism is the apsidal phase-locking between a couple of planets for a certain system. For GJ 876 and HD 82943, we discover the relative apsidal longitudes $\varpi_{1} - \varpi_{2}$ move about $0^{\circ}$ or $180^{\circ}$, respectively; but for 55 Cnc, we find that there exists an asymmetric apsidal libration between two inner planets. Finally, we made a brief discussion about the Habitable Zones in the exoplanetary systems.

The influence of the rotation on the total masses and radii of the neutron stars are calculated by the Hartle's slow rotation formalism, while the equation of state is considered in a relativistic $\sigma-\omega$ model. Comparing with the observation, the calculating result shows that the double neutron star binaries are more like hyperon stars and the neutron stars of X-ray binaries are more like traditional neutron stars. As the changes of the mass and radius to a real neutron star caused by the rotation are very small comparing with the total mass and radius, one can see that Hartle's approximate method is rational to deal with the rotating neutron stars. If three property values: mass, radius and period are observed to the same neutron star, then the EOS of this neutron star could be decided entirely.