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

We carry out a systematic study of the density structure of gas in the Central Molecular Zone (CMZ) in the Galactic center by extracting clumps from the APEX Telescope Large Area Survey of the Galaxy survey at 870 $\mu$m. We find that the clumps follow a scaling of $m = \rho_0 r^3$ which corresponds to a characteristic density of $n_{\rm H_2} = 1.6 \times 10^3\,\rm cm^{-3}$ ($\rho_0 =112\;M_{\odot}\;\rm pc^{-3}$) with a variation of $\approx 0.5\,\rm dex$, where we assumed a gas-to-dust mass ratio of 100. This characteristic density can be interpreted as the result of thermal pressure equilibrium between the molecular gas and the warm ambient interstellar medium. Such an equilibrium can plausibly be established since shear has approximately the same strength as self-gravity. Our findings may explain the fact that star formation in the CMZ is highly inefficient compared to the rest of the Milky Way disk. We also identify a population of clumps whose densities are two orders of magnitudes higher in the vicinity of the Sgr B2 region, which we propose are produced by collisions between the clumps of lower densities. For these collisions to occur, processes such as compressive tides probably have created the appropriate condition by assembling the clumps together.

GWAC will have been built an integrated FOV of 5,000 $degree^2$ and have already built 1,800 square $degree^2$. The limit magnitude of a 10-second exposure image in the moonless night is 16R. In each observation night, GWAC produces about 0.7TB of raw data, and the data processing pipeline generates millions of single frame alerts. We describe the GWAC Data Processing and Management System (GPMS), including hardware architecture, database, detection-filtering-validation of transient candidates, data archiving, and user interfaces for the check of transient and the monitor of the system. GPMS combines general technology and software in astronomy and computer field, and use some advanced technologies such as deep learning. Practical results show that GPMS can fully meet the scientific data processing requirement of GWAC. It can online accomplish the detection, filtering and validation of millions of transient candidates, and feedback the final results to the astronomer in real-time. During the observation from October of 2018 to December of 2019, we have already found 102 transients.

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young.

The initial stage of star formation is a complex area study because of its high density and low temperature. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the deuterated species could remain gaseous and are thus ideal tracers. We investigate the gas dynamics and NH$_2$D chemistry in eight massive pre/protocluster clumps. We present NH$_2$D 1$_{11}$-1$_{01}$ (at 85.926 GHz), NH$_3$ (1, 1) and (2, 2) observations in the eight clumps using the PdBI and the VLA, respectively. We find that the distribution between deuterium fractionation and kinetic temperature shows a number density peak at around $T_{\rm kin}=16.1$ K, and the NH$_2$D cores are mainly located at a temperature range of 13.0 to 22.0 K. We detect seven instances of extremely high deuterium fractionation of $1.0 \leqslant D_{\rm frac} \leqslant 1.41$. We find that the NH$_2$D emission does not appear to coincide exactly with either dust continuum or NH$_3$ peak positions, but often surrounds the star-formation active regions. This suggests that the NH$_{2}$D has been destroyed by the central young stellar object (YSO) due to its heating. The detected NH$_2$D lines are very narrow with a median width of $\rm 0.98\pm0.02 km/s$. The extracted NH$_2$D cores are gravitationally bound ($\alpha_{\rm vir} < 1$), are likely prestellar or starless, and can potentially form intermediate-mass or high-mass stars. Using NH$_3$ (1, 1) as a dynamical tracer, we find very complicated dynamical movement, which can be explained by a combined process with outflow, rotation, convergent flow, collision, large velocity gradient, and rotating toroids. High deuterium fractionation strongly depends on the temperature condition. NH$_2$D is a poor evolutionary indicator of high-mass star formation in evolved stages, but a useful tracer in the starless and prestellar cores.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

On September 10 2017, a fast coronal mass ejection (CME) erupted from the active region (AR) 12673, leading to a ground level enhancement (GLE) event at Earth. Using the 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we model the large solar energetic particle (SEP) event of 10 September 2017 observed at Earth, Mars and STEREO-A. Based on observational evidence,we assume that the CME-driven shock experienced a large lateral expansion shortly after the eruption, which is modelled by a double Gaussian velocity profile in this simulation. We use the in-situ shock arrival times and the observed CME speeds at multiple spacecraft near Earth and Mars as constraints to adjust the input model parameters. The modelled time intensity profiles and fluence for energetic protons are then compared with observations. Reasonable agreements with observations at Mars and STEREO-A are found. The simulated results at Earth differ from observations of GOES-15. Instead, the simulated results at a heliocentric longitude 20 degree west to Earth fit reasonably well with the GOES observation. This can be explained if the pre-event solar wind magnetic field at Earth is not described by a nominal Parker field. Our results suggest that a large lateral expansion of the CME-driven shock and a distorted interplanetary magnetic field due to previous events can be important in understanding this GLE event.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

Filament S242 is 25 pc long with massive clumps and YSO clusters concentrated in its end regions; it is considered a good example of edge collapse. We mapped this filament in the $^{12}$CO(1-0) and $^{13}$CO(1-0) lines. A large-scale velocity gradient along filament S242 has been detected; the relative velocity between the two end-clumps is $\sim$ 3 km s$^{-1}$, indicating an approaching motion between them. These signatures are consistent with the filament S242 being formed through the collapse of a single elongated entity, where an effect known as "gravitational focusing" drives the ends of the filament to collapse (edge collapse). Based on this picture, we estimate a collapse timescale of $\sim$ 4.2 Myr, which is the time needed for a finite and elongated entity evolving to the observed filament S242. For the whole filament, we find that increases in surface densities lead to increases in velocity dispersion, which can be consistently explained as the result of self-gravity. We also calculated the contribution of longitudinal collapse to the observed velocity dispersion and found it to be the dominant effect in driving the gas motion near the end-clumps. We propose that our filament S242 is formed through a two-stage collapse model, where the edge collapse of a truncated filament is followed by a stage of longitudinal accretion toward the dense end-clumps.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

Rotation and binary interaction play important roles in understanding the nature of massive stars (>8 Msun). Some interesting transients, such as the long-duration gamma-ray bursts, are thought to be originated from fast-rotating massive stars. Because the strong stellar wind can effectively spin down a metal-rich massive star with fast rotation, it is very hard to find single massive stars rotating critically in the Galaxy. In the present work reported is the discovery of the fastest rotator in the Galaxy, LAMOST J040643.69+542347.8, with a projected rotational velocity vsini ~540 km/s, which is ~ 100 km/s faster than that of the previous record holder HD 191423. The star has a spectral type of O6.5 Vnnn(f)p. Its He I 4471 absorption line is blueshifted and asymmetric, while its He II 4686 and Ha have central absorption reversals in their emissions. It is also a runaway star, which implies an origin in a close binary interaction. Compared to VFTS 285 and VFTS 102 (their vsini ~ 610 km/s) in the Large Magellanic Cloud, LAMOST J040643.69+542347.8 has its own peculiar spectral characteristics and earlier spectral type. Moreover, LAMOST J040643.69+542347.8 is bright (B ~ 13.9 mag) enough to allow future high-resolution spectroscopic follow-ups.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

The Moon's optical characteristics in visible and long-wavelength infrared (LWIR) have long been observed with our eyes or with instruments. What the mid-infrared (MIR) Moon looks like is still a mystery. For the first time we present detailed appearance of the MIR Moon observed by a high-resolution geostationary satellite and reveal the essence behind its appearance. The appearance of the MIR Moon is opposite to its normal visible appearance. In addition the MIR Moon shows limb darkening. Both the absolute and the relative brightness distribution of the MIR lunar disk changes with the solar incidence angle. The signatures of the MIR Moon are controlled by both the reflection and emission of the lunar surface. We also show first-ever brightness temperature maps of the lunar disk without needing a mosaic, which better show the temperature variation across the lunar disk. They reveal that the relationship between brightness temperature and solar incidence angle i is cos1/bi, and the power parameter is smaller than the Lambertian temperature model of cos1/4i observed for lunar orbit-based measurements. The slower decrease of the brightness temperature when moving away from the sub-solar point than the Lambertian model is due to topographic effects. The brightness temperature is dominated by albedo and the solar incidence angle and influenced by the topography. Our results indicate that the Moon in the MIR exhibits many interesting phenomena which were previously unknown, and contains abundant information about lunar reflection and thermal emission for future study.

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Ultra-short-period planets (USPs) provide important clues to planetary formation and migration. Recently, it is found that the mutual inclinations of the planetary systems are larger if the inner orbits are closer ($\lesssim 5R_*$) and if the planetary period ratios are larger ($P_2/P_1 \gtrsim 5$) (Dai et al. 2018). This suggests that the USPs experienced both inclination excitation and orbital shrinkage. Here we investigate the increase in the mutual inclination due to stellar oblateness. We find that the stellar oblateness (within $\sim 1$Gyr) is sufficient to enhance the mutual inclination to explain the observed signatures. This suggests that the USPs can migrate closer to the host star in a near coplanar configuration with their planetary companions (e.g., disk migration+tides or in-situ+tides), before mutual inclination gets excited due to stellar oblateness.

We present a large, homogeneous catalogue of molecular clouds within 4 kpc from the Sun at low Galactic latitudes ($|b|$ $<$ 10\degr) with unprecedented accurate distance determinations. Based on the three-dimensional dust reddening map and estimates of colour excesses and distances of over 32 million stars presented in Chen et al, we have identified 567 dust/molecular clouds with a hierarchical structure identification method and obtained their distance estimates by a dust model fitting algorithm. The typical distance uncertainty is less than 5 per cent. As far as we know, this is the first large catalogue of molecular clouds in the Galactic plane with distances derived in a direct manner. The clouds are seen to lie along the Sagittarius, Local and Perseus Arms. In addition to the known structures, we propose the existence of a possible {\it spur}, with a pitch angle of about 34\degr, connecting the Local and the Sagittarius Arms in the fourth quadrant. We have also derived the physical properties of those molecular clouds. The distribution of cloud properties in different parameter spaces agrees grossly with the previous results. Our cloud sample is an ideal starting point to study the concentration of dust and gas in the solar vicinity and their star formation activities.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

Planets that orbit only one of the stars in stellar binary systems (i.e., circumstellar) are dynamically constrained to a limited range of orbital parameters and thus understanding conditions on their stability is of great importance in exoplanet searches. We perform $\sim$700 million N-body simulations to identify how stability regions depend on properties of the binary, as well as, the starting planetary inclination and mean longitude relative to the binary orbit. Moreover, we provide grid interpolation maps and lookup tables for the community to use our results. Through Monte-Carlo methods we determine that planets with a semimajor axis $a_p$ $\lesssim$8\% of the binary semimajor axis $a_{bin}$ will likely be stable given the known distribution of binary star parameters. This estimate varies in the Lidov-Kozai regime or for retrograde orbits to 4\% or 10\% of $a_{bin}$, respectively. Our method to quickly determine the circumstellar stability limit is important for interpreting observations of binaries using direct imaging with JWST, photometry with TESS, or even astrometry with Gaia.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

Changes in planetary obliquity, or axial tilt, influence the climates on Earth-like planets. In the solar system, the Earth's obliquity is stabilized due to interactions with our moon and the resulting {small amplitude variations ($\sim$2.4\degree)} are beneficial for advanced life. Most Sun-like stars have at least one stellar companion and the habitability of circumstellar exoplanets is shaped by their stellar companion. We show that a stellar companion can dramatically change whether {Earth-like obliquity stability is} possible through planetary orbital precession relative to the binary orbit or resonant pumping of the obliquity through spin-orbit interactions. We present a new formalism for the planetary spin precession that accounts for orbital misalignments between the planet and binary. Using numerical modeling in $\alpha$ Centauri AB we show: a stark contrast between the planetary obliquity variations depending on the host star, planetary neighbors limit the possible spin states for {Earth-like obliquity stability}, and the presence of a moon can destabilize the obliquity, defying our Earth-based expectations. An Earth-like rotator orbiting the primary star would experience {small} obliquity variations for 87\%, 74\%, or 54\% of Solar type binaries, depending on the mass of the primary (0.8, 1.0, or 1.2 M$_\odot$, respectively). Thus, Earth-like planets likely experience much larger obliquity variations, with more extreme climates, unless they are in specific states, such as orbiting nearly planar with the binary and rotating retrograde (backwards) like Venus.

The ON stars are a rare subtype of O stars, of uncertain origin. We report two new, rapidly-rotating ON stars found in data acquired with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST. LS I +61 28 is an ON8.5 Vn dwarf with a projected equatorial rotational velocity of $v_{\textrm e}\sin{i} \simeq 298$ km$\;$s$^{-1}$, while HDE 236672 is an ON9 IVn subgiant with $v_{\textrm e}\sin{i} \simeq 253$ km$\;$s$^{-1}$ The former is the first rapidly-rotating ON dwarf to be found, and the latter is only the third ON subgiant. The luminosity classes of non-supergiant ON stars appear to be influenced by the axial inclination angle $i$: the rapidly-rotating giants are close to equator-on, while ON dwarfs with lower $v_{\textrm e}\sin{i}$ values are viewed more nearly pole-on. Combining parallaxes and proper motions from Gaia DR2 with radial-velocity measurements, we investigate the kinematics of non-supergiant ON stars, and infer that the dynamics, rapid rotation, and surface-nitrogen characteristics may all be consequences of binary interaction.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

We obtained 7 spectra of the Be star V423 Aur on Dec. 5th, 2017 using the LAMOST Medium-Resolution Spectrograph with exposures from 600 to 1200 seconds. These spectra show the irregular H{\alpha} emission line profile variations (LPVs). In the seven spectra, from the 4th to 7th, the left part of H{\alpha} profile even show excess. However, no variation can be seen from the follow-up observation of photometry by 1.26-m telescope and High-Resolution spectra by 2.16-m telescope. According to the High-Resolution spectra, we conclude that it is a B7V type star with E(B-V)= 0.709{\pm}0.036 and its vsini is {\sim} 221.8km/s. The short-term H{\alpha} LPVs could be explained as a result of the transient ejection of matter from rotating disk or shell around V423 Aur.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Mean motion resonances are commonly seen in planetary systems, e.g., in the formation of orbital structure of Jupiter's moons and the gaps in the rings of Saturn. In this work we study their effects in fully relativistic systems. We consider a model problem with two stellar mass black holes orbiting around a supermassive black hole. By adopting a two time-scale expansion technique and averaging over the fast varying orbital variables, we derive the effective Hamiltonian for the slowly varying dynamical variables. The formalism is illustrated with a n'_phi : n'_r : n_phi= 2:1:-2 resonance in Schwarzschild spacetime, which naturally becomes the 3:2 resonance widely studied in the Newtonian limit. We also derive the multi-body Hamiltonian in the post-Newtonian regime, where the radial and azimuthal frequencies are different because of the post-Newtonian precession. The capture and breaking conditions for these relativistic mean motion resonances are also discussed. In particular, pairs of stellar mass black holes surrounding the supermassive black hole could be locked into resonances as they enter the LISA band, and this would affect their gravitational wave waveforms.

N131 is a typical infrared dust bubble showing an expanding ringlike shell. We study what kinds of CO line ratios can be used to trace the interaction in the expanding bubble. We carry out new $\rm CO\,(3-2)$ observations towards bubble N131 using the JCMT 15-m telescope, and derive line ratios by combining with our previous $\rm CO\,(2-1)$ and $\rm CO\,(1-0)$ data from the IRAM 30-m telescope observations. To trace the interaction between the molecular gas and the ionized gas in the HII region, we use RADEX to model the dependence of CO line ratios on kinetic temperature and H$_2$ volume density, and examine the abnormal line ratios based on other simulations. We present $\rm CO\,(3-2)$, $\rm CO\,(2-1)$, and $\rm CO\,(1-0)$ integrated intensity maps convolved to the same angular resolution (22.5$"$). The three different CO transition maps show apparently similar morphology. The line ratios of $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}$ mostly range from 0.2 to 1.2 with a median of $0.54\pm0.12$, while the line ratios of $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}$ range from 0.5 to 1.6 with a median of $0.84\pm0.15$. The high CO line ratios $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}\gtrsim 0.8 $ and $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}\gtrsim 1.2$ are beyond the threshold predicted by numerical simulations based on the assumed density-temperature structure for the inner rims of ringlike shell, where are the compressed areas in bubble N131. These high CO integrated intensity ratios, such as $W_{\rm CO\,(3-2)}$/$W_{\rm CO\,(2-1)}\gtrsim0.8$ and $W_{\rm CO\,(2-1)}$/$W_{\rm CO\,(1-0)}\gtrsim1.2$, can be used as a tracer of gas compressed regions with a relatively high temperature and density. This further suggests that the non-Gaussian part of the line-ratio distribution can be used to trace the interaction between the molecular gas and the hot gas in the bubble.

We report our survey of $\gamma$\,Doradus stars from the 4-yr \textit{Kepler} mission. These stars pulsate mainly in g modes and r modes, showing period-spacing patterns in the amplitude spectra. The period-spacing patterns are sensitive to the chemical composition gradients and the near-core rotation, hence they are essential for understanding the stellar interior. We identified period-spacing patterns in 611 $\gamma$\,Dor stars. Almost every star pulsates in dipole g modes, while about 30\% of stars also show clear patterns for quadrupole g modes and 16\% of stars present r mode patterns. We measure periods, period spacings, and the gradient of the period spacings. These three observables guide the mode identifications and can be used to estimate the near-core rotation rate. We find many stars are hotter and show longer period-spacing patterns than theory. Using the Traditional Approximation of Rotation (TAR), we inferred the asymptotic spacings, the near-core rotation rates, and the radial orders of the g and r modes. Most stars have a near-core rotation rate around $1$\,$\mathrm{c/d}$ and an asymptotic spacing around 4000\,s. We also find that many stars rotate more slowly than predicted by theory for unclear reasons. 11 stars show rotational splittings with fast rotation rates. We compared the observed slope--rotation relation with the theory and find a large spread. We detected rotational modulations in 58 stars and used them to derive the core-to-surface rotation ratios. The interiors rotate faster than the cores in most stars, but by no more than 5\%.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

Quasi-biennial oscillation (QBO) of solar magnetic activities is intrinsic to dynamo mechanism, but still far from fully understood. In this work, the phase and amplitude asymmetry of solar QBO of Halpha flare activity in the northern and southern hemispheres is studied by the ensemble empirical mode decomposition, the cross-correlation analysis, and the wavelet transform technique. The following results are found: (1) solar QBO of Halpha flare index in the two hemispheres has a complicated phase relationship, but does not show any systematic regularity; (2) the solar cycle mode of solar Halpha flare index in the northern hemisphere generally leads that in the southern one by 9 months for the time interval from 1966 to 2014. The possible origin of these results is discussed.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

We present the first asteroseismic results for $\delta$ Scuti and $\gamma$ Doradus stars observed in Sectors 1 and 2 of the TESS mission. We utilise the 2-min cadence TESS data for a sample of 117 stars to classify their behaviour regarding variability and place them in the Hertzsprung-Russell diagram using Gaia DR2 data. Included within our sample are the eponymous members of two pulsator classes, $\gamma$ Doradus and SX Phoenicis. Our sample of pulsating intermediate-mass stars observed by TESS also allows us to confront theoretical models of pulsation driving in the classical instability strip for the first time and show that mixing processes in the outer envelope play an important role. We derive an empirical estimate of 74% for the relative amplitude suppression factor as a result of the redder TESS passband compared to the Kepler mission using a pulsating eclipsing binary system. Furthermore, our sample contains many high-frequency pulsators, allowing us to probe the frequency variability of hot young $\delta$ Scuti stars, which were lacking in the Kepler mission data set, and identify promising targets for future asteroseismic modelling. The TESS data also allow us to refine the stellar parameters of SX Phoenicis, which is believed to be a blue straggler.

We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.

In this work, BV RI light curves of 55 Type II supernovae (SNe II) from the Lick Observatory Supernova Search program obtained with the Katzman Automatic Imaging Telescope and the 1 m Nickel telescope from 2006 to 2018 are presented. Additionally, more than 150 spectra gathered with the 3 m Shane telescope are published. We conduct an analyse of the peak absolute magnitudes, decline rates, and time durations of different phases of the light and colour curves. Typically, our light curves are sampled with a median cadence of 5.5 days for a total of 5093 photometric points. In average V-band plateau declines with a rate of 1.29 mag (100 days)-1, which is consistent with previously published samples. For each band, the plateau slope correlates with the plateau length and the absolute peak magnitude: SNe II with steeper decline have shorter plateau duration and are brighter. A time-evolution analysis of spectral lines in term of velocities and pseudoequivalent widths is also presented in this paper. Our spectroscopic sample ranges between 1 and 200 days post-explosion and has a median ejecta expansion velocity at 50 days post-explosion of 6500 km/s (Halpha line) and a standard dispersion of 2000 km/s. Nebular spectra are in good agreement with theoretical models using a progenitor star having a mass <16 Msol. All the data are available to the community and will help to understand SN II diversity better, and therefore to improve their utility as cosmological distance indicators.

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

The hot massive luminous blue variables (LBVs) represent an important evolutionary phase of massive stars. Here, we report the discovery of a new LBV -- LAMOST J0037+4016 in the distant outskirt of the Andromeda galaxy. It is located in the south-western corner (a possible faint spiral arm) of M31 with an unexpectedly large projection distance of $\sim$ 22 kpc from the center. The optical light curve shows a 1.2 mag variation in $V$ band and its outburst and quiescence phases both last over several years. The observed spectra indicate an A-type supergiant at epoch close to the outburst phase and a hot B-type supergiant with weak [Fe II] emission lines at epoch of much dimmer brightness. The near-infrared color-color diagram further shows it follows the distribution of Galactic and M31 LBVs rather than B[e] supergiants. All the existing data strongly show that LAMOST J0037+4016 is an LBV. By spectral energy distribution fitting, we find it has a luminosity ($4.42 \pm 1.64$)$\times 10^5$ $L_{\odot}$ and an initial mass $\sim 30$ $M_{\odot}$, indicating its nature of less luminosity class of LBV.

The short-period (P $\approx1.7$ d), Algol-type eclipsing binary KIC 7385478 consists of an F-type primary star ($M_1 \approx 1.71M_{\odot}$) and an evolved K-type secondary ($M_2 \approx 0.37M_{\odot}$) (Ozdarcan \& Ali Dal 2017). We study the variability of the {\it Kepler} light curve and attribute many frequency peaks in the Fourier spectrum to the spot modulation. These frequencies are in the form of orbital harmonics and are highly variable in amplitude. They are most likely from the mass-accreting primary star. In addition, we identify a series prograde dipole g-modes from the primary star which show a quasi-linear period spacing pattern and are very stable in amplitude. The period spacing pattern reveals an asymptotic period spacing value in agreement with fundamental parameters of the primary star and also implies that the near-convective-core rotation rate is almost the same as the orbital period. Thus both the surface and the core of this Gamma Dor pulsator have synchronized with the binary orbit. We find that a lower stellar mass $\approx 1.50 M_{\odot}$ and higher effective temperature are needed in order to be compatible with the asteroseismic constraints from single star evolutionary models.

We present the characterization of KIC 4142768, an eclipsing binary with two evolved A-type stars in an eccentric orbit with a period of 14 days. We measure the fundamental parameters of the two components ($M_1=2.05M_{\odot}, R_1=2.96R_{\odot}$ and $M_2=2.05M_{\odot}, R_2=2.51R_{\odot}$) by combining {\it Kepler} photometry and spectra from {\it Keck} HIRES. The measured surface rotation rates are only one-fifth of the pseudo-synchronous rate of the eccentric orbit. Fourier spectrum of the light curve reveals hybrid pulsations of $\delta$ Scuti and $\gamma$ Doradus type, with pulsation frequencies at about $15-18$ day$^{-1}$ for p modes and about $0.2-1.2$ day$^{-1}$ for low-frequency g modes. Some of the g modes are exact orbital harmonics and are likely tidally excited. Their pulsation amplitudes and phases both agree with predictions from the linear tidal theory for $l=2, m=2$ prograde modes. We examine the period spacing patterns in the free oscillating g modes and identify them mostly as prograde sectoral dipole modes. The unstable frequency range and frequency spacing of p modes and the inferred asymptotic g-mode period spacings both agree with the stellar model for the primary star evolved to a late stage of the main sequence. The inferred rotation rate of the convective core boundary is very slow, similar to the small surface rotation rate inferred from the spectroscopy. The measured surface and near-core rotation rates provide constraints for testing the mechanism of angular momentum transfer and tidal synchronization in evolved eccentric binary star systems.

Pulsar timing arrays (PTAs) can be used to study the Solar-system ephemeris (SSE), the errors of which can lead to correlated timing residuals and significantly contribute to the PTA noise budget. Most Solar-system studies with PTAs assume the dominance of the term from the shift of the Solar-system barycentre (SSB). However, it is unclear to which extent this approximation can be valid, since the perturbations on the planetary orbits may become important as data precision keeps increasing. To better understand the effects of SSE uncertainties on pulsar timing, we develop the LINIMOSS dynamical model of the Solar system, based on the SSE of Guangyu Li. Using the same input parameters as DE435, the calculated planetary positions by LINIMOSS are compatible with DE435 at centimetre level over a 20-year timespan, which is sufficiently precise for pulsar-timing applications. We utilize LINIMOSS to investigate the effects of SSE errors on pulsar timing in a fully dynamical way, by perturbing one SSE parameter per trial and examining the induced timing residuals. For the outer planets, the timing residuals are dominated by the SSB shift, as assumed in previous work. For the inner planets, the variations in the orbit of the Earth are more prominent, making previously adopted assumptions insufficient. The power spectra of the timing residuals have complex structures, which may introduce false signals in the search of gravitational waves. We also study how to infer the SSE parameters using PTAs, and calculate the accuracy of parameter estimation.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time intensity profile and spectrum of SEPs accelerated at CME-driven shocks from stars of different ages traced by their rotation rates. We consider a solar-like (G-type) star with 6 different rotation rates varying from 0.5 to 3.0 times rotation speed of the sun. In all 6 cases, a fast CME is launched with the same speed of 1500 km/sec and the resulting time intensity profiles at 3 locations and and energy spectra at 5 locations at 1 AU are obtained. The maximum particle energy at the shock front as a function of r is also shown. Our results suggest that within 0.8 AU the maximum particle energy at the shock front increases with the rotation rate of the star. However, event integrated spectra for the five selected locations along the CME path show complicated patterns. This is because the Parker magnetic field for rapidly rotating stars is more tightly winded. Our results can be used in estimating the radiation environments of terrestrial-type exoplanets around solar-type stars.

The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds.

Massive clumps tend to fragment into clusters of cores and condensations, some of which form high-mass stars. In this work, we study the structure of massive clumps at different scales, analyze the fragmentation process, and investigate the possibility that star formation is triggered by nearby HII regions. We present a high angular resolution study of a sample of 8 massive proto-cluster clumps. Combining infrared data, we use few-arcsecond resolution radio- and millimeter interferometric data to study their fragmentation and evolution. Our sample is unique in the sense that all the clumps have neighboring HII regions. Taking advantage of that, we test triggered star formation using a novel method where we study the alignment of the centres of mass traced by dust emission at multiple scales. The eight massive clumps have masses ranging from 228 to 2279 $M_\odot$. The brightest compact structures within infrared bright clumps are typically associated with embedded compact radio continuum sources. The smaller scale structures of $R_{\rm eff}$ $\sim$ 0.02 pc observed within each clump are mostly gravitationally bound and massive enough to form at least a B3-B0 type star. Many condensations have masses larger than 8 $M_\odot$ at small scale of $R_{\rm eff}$ $\sim$ 0.02 pc. Although the clumps are mostly infrared quiet, the dynamical movements are active at clump scale ($\sim$ 1 pc). We studied the spatial distribution of the gas conditions detected at different scales. For some sources we find hints of external triggering, whereas for others we find no significant pattern that indicates triggering is dynamically unimportant. This probably indicates that the different clumps go through different evolutionary paths. In this respect, studies with larger samples are highly desired.

We use the Richardson-Lucy deconvolution algorithm to extract one dimensional (1D) spectra from LAMOST spectrum images. Compared with other deconvolution algorithms, this algorithm is much more fast. The practice on a real LAMOST image illustrates that the 1D resulting spectrum of this method has a higher SNR and resolution than those extracted by the LAMOST pipeline. Furthermore, our algorithm can effectively depress the ringings that are often shown in the 1D resulting spectra of other deconvolution methods.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

Kepler-1647b is the most recently discovered planet that transits two stars, i.e., a circumbinary planet (CBP). Due to its large orbital separation, Kepler-1647b stands out from the rest of the Kepler CBPs, which mostly reside on much tighter orbits near the stability limit. The large separation of Kepler-1647b challenges inward disk migration as a dominant formation pathway, suggested by the other Kepler CBPs. In this paper, we consider the possibility of an undetected planet near the stability limit by examining observational consequences of such a planet. We calculate the transit probability of the putative planet, transit timing variations (TTVs) of the known planet, and eclipsing timing variations (ETVs) of the host binary caused by the putative planet. We find the presence of a $\gtrsim$30M$_{\oplus}$ inner planet to be highly unlikely near the stability limit. In addition, we provide future TTV observation windows, which will further constrain possible undetected planets with lower masses.

We calculate the interplanetary magnetic field path lengths traveled by electrons in solar electron events detected by the WIND 3DP instrument from $1994$ to $2016$. The velocity dispersion analysis method is applied for electrons at energies of $\sim$ $27$ keV to $310$ keV. Previous velocity dispersion analyses employ the onset times, which are often affected by instrumental effects and the pre-existing background flux, leading to large uncertainties. We propose a new method here. Instead of using the peak or onset time, we apply the velocity dispersion analysis to the times that correspond to the rising phase of the fluxes that are a fraction, $\eta$, of the peak flux. We perform statistical analysis on selected events whose calculated path lengths have uncertainties smaller than $0.1$ AU. The mean and standard deviation, ($\mu$, $\sigma$), of the calculated path lengths corresponding to $\eta=$ $3/4$, $1/2$, and $1/3$ of the peak flux is ($1.17$ AU, $0.17$ AU), ($1.11$ AU, $0.14$ AU), and ($1.06$ AU, $0.15$ AU). The distribution of the calculated path lengths is also well fitted by a Gaussian distribution for the $\eta=3/4$ and $1/3$ cases. These results suggest that in these electron events the interplanetary magnetic field topology is close to the nominal Parker spiral with little field line meandering. Our results have important implications for particles' perpendicular diffusion.

Rossby modes are the oscillations in a rotating fluid, whose restoring force is the Coriolis force. They provide an additional diagnostic to understand the rotation of stars, which complicates asteroseismic modelling. We report 82 $\gamma$\,Doradus stars for which clear period spacing patterns of both gravity and Rossby modes have been detected. The period spacings of both show a quasi-linear relation with the pulsation period but the slope is negative for the gravity modes and positive for the Rossby modes. Most Rossby modes have $k=-2, m=-1$. For only one star a series of $k=-1,m=-1$ modes is seen. For each pattern, the mean pulsation period, the mean period spacing, and the slope are measured. We find that the slope correlates with the mean period for Rossby mode patterns. The traditional approximation of rotation is used to measure the near-core rotation rate, assuming the star rotates rigidly. We report the near-core rotation rates, the asymptotic period spacings, and the radial orders of excited modes of these 82 main-sequence stars. The near-core rotation rates lie between $0.6\,\mathrm{d^{-1}}$ and $2.3\,\mathrm{d^{-1}}$. Six stars show surface rotation modulations, among which only KIC\,3341457 shows differential rotation while the other five stars have uniform rotations. The radial orders of excited modes show different distributions for the dipole and quadrupole gravity modes versus the Rossby modes.

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline.

SWIFT J1658.2-4242 is an X-ray transient discovered recently in the Galactic plane, with severe X-ray absorption corresponding to an equivalent hydrogen column density of $N_{\rm H,abs}\sim2\times10^{23}$ cm$^{-2}$. Using new Chandra and XMM-Newton data, we discover a strong X-ray dust scattering halo around it. The halo profile can be well fitted by the scattering from at least three separated dust layers. During the persistent emission phase of SWIFT J1658.2-4242, the best-fit dust scattering $N_{\rm H,sca}$ based on the COMP-AC-S dust grain model is consistent with $N_{\rm H,abs}$. The best-fit halo models show that 85-90 percent of the intervening gas and dust along the line of sight of SWIFT J1658.2-4242 are located in the foreground ISM in the Galactic disk. The dust scattering halo also shows significant azimuthal asymmetry, which appears consistent with the inhomogeneous distribution of foreground molecular clouds. By matching the different dust layers to the distribution of molecular clouds along the line of sight, we estimate the source distance to be $\sim$10 kpc, which is also consistent with the results given by several other independent methods of distance estimation. The dust scattering opacity and the existence of a halo can introduce a significant spectral bias, the level of which depends on the shape of the instrumental point spread function and the source extraction region. We create the Xspec dscor model to correct for this spectral bias for different X-ray instruments. Our study reenforces the importance of considering the spectral effects of dust scattering in other absorbed X-ray sources.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Nuclear starburst discs (NSDs) are very compact star-forming regions in the centers of galaxies that have been studied as a possible origin for the absorbing gas around a central active galactic nucleus. NSDs may be most relevant at $z\sim 1$ when obscured accretion onto supermassive black holes (SMBHs) is common. This paper describes the characteristics of the stellar remnants of NSDs at $z=0.01$, taking into account the evolution from $z=1$. Using a stellar synthesis model, the colours, masses, and luminosities of the stellar remnants are computed for a suite of 192 two-dimensional NSD models. These properties are compared to observations of local nuclear star clusters (NSCs), and a good match is found between the predicted and observed properties. Dynamical effects will likely cause the final remnant to be a rotating, nearly spherical distribution. In addition, $\approx 20$% of the NSD remnants have half-light radii <~ 10 pc, consistent with NSCs hosted in both late-type and early-type galaxies, and all the remnants follow similar size-luminosity relationships as observed in nearby NSCs. NSDs require the presence of a central SMBH and the most massive and compact stellar remnants are associated with the most massive SMBHs, although stellar clusters with a variety of sizes can be produced by all considered SMBH masses. Overall, NSDs at $z\sim 1$ appear to be a promising origin for the $\gg 1$ Gyr NSC population in early- and late-type galaxies with large SMBHs.

We investigate the evolution of triple, non-hierarchical, black hole (BH) systems making use of $2.9\times10^4$ 3-body simulations. Varying the mutual orbital inclination, the three BH masses and the inner and outer eccentricities, we show that retrograde, nearly planar configurations lead to a significant shrinkage of the inner binary. We find an universal trend of triple systems, that they tend to evolve toward prograde configurations, Moreover, we demonstrate that the orbital flip, driven by the torque exerted on the inner BH binary (BHB) by the outer BH, leads in general to tighter inner orbits. In some cases, the resulting BHB undergoes coalescence within a Hubble time, releasing gravitational waves (GWs). Frequently, the inner BHB merger occurs after a component swap between one of its components and the outer BH. The mass spectrum of the BHBs that underwent the component exchange differs significantly from the case in which the BHB merge without any swap. A large fraction of merging BHBs with initial separation $1$ AU enter the $10^{-3}-10^{-1}$ Hz frequency band with large eccentricities, thus representing potential LISA sources. Mergers originating from initially tighter BHB ($a\sim 0.01$ AU), instead, have a large probability to have eccentricities above 0.7 in the $1$ Hz band. We find that the mergers' mass distribution in this astrophysical channel maps the original BH binary spectrum. This might have interesting consequences in light of the growing population of BH mergers detected by LIGO.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

Combining data from Gaia DR2, SDSS DR14 and LAMOST DR6, we update the fit to model of the properties of the stellar stream GD-1 and find that it has an age of $\sim 13$ Gyr, [Fe/H] of $-2.2\pm 0.12$, and a distance from the sun of $\sim 8$ kpc. We tabulate 6D phase-space fiducial points along the GD-1 stream orbit over a 90$^\circ$ arc. The fitted orbit shows that the stream has an eccentricity $e \sim 0.3$, perigalacticon of 14.2 kpc, apogalacticon of 27.0 kpc and inclination $i \sim 40^{\circ}$. There is evidence along the arc for 4 candidate stellar overdensities, one candidate gap, two candidate stellar underdensities and is cut off at $\phi_1 \sim 2^{\circ}$ (in the stream-aligned $(\phi_1,\phi_2)$ coordinate system of \citet{kop10}). The spur originating at $\phi_1 \sim -40^{\circ}$ implies stars were pulled away from the stream trace by an encounter (potentially a dark matter subhalo). The narrowest place ( FWHM $\sim 44.6$ pc) of the GD-1 trace is at $(\phi_1, \phi_2^c) \sim (-14^{\circ}, 0.15^{\circ})$, which is $\sim (178.18^{\circ}, 52.19^{\circ})$ in (R. A., Decl.), where the progenitor is possibly located. We also find six BHB and 10 BS spectroscopic stars in the GD-1 stream.

The existence of Planet Nine has been suggested to explain the pericenter clustering of extreme trans-Neptunian objects (TNOs). However, the underlying dynamics involving Planet Nine, test particles and Neptune is rich, and it remains unclear which dynamical processes lead to the alignment and how they depends on the properties of Planet Nine. Here, we investigate the secular interactions between an eccentric outer perturber and TNOs starting in a near-coplanar configuration. We find that a large number of TNOs could survive outside of mean motion resonances at 4Gyr, which differs from previous results obtained in the exact coplanar case with Neptune being treated as a quadrupole potential. In addition, secular dynamics leads to the orbital clustering seen in N-body simulations. We find that a near coplanar Planet Nine can flip TNO orbital planes, and when this happens, the geometrical longitudes of pericenter of the TNOs librate around $180^\circ$ during the flip. Orbital precession caused by the inner giant planets can suppress the flips while keeping the longitude of pericenter librating when $30\lesssim r_p \lesssim 80$ AU $\&$ $a\gtrsim 250$ AU. This results in the alignment of the pericenter of the low inclination TNOs ($i \lesssim 40^\circ$). We find the anti-aligned population and the flipped orbits could be produced by an eccentric ($e_9 \gtrsim 0.4$) outer planet of $\sim 10M_{\oplus}$ in a wide $a_9 \gtrsim 400 \sim 800$ AU orbit. Future surveys on the high inclination TNOs will help further constrain the properties of possible outer planets.

In $\gamma$ Doradus stars, the g-mode period spacing shows an approximately linear relation with period. The slope is a new asteroseismic diagnostic, related to the rotation rate and the azimuthal order $m$. We report two automated methods, the `moving-window Fourier transform' and the `cross-correlation', to detect and measure the period spacings based on four-year light curves from the \textit{Kepler} satellite. The results show that the cross-correlation method performs better at detecting the period spacings and their slopes. In this paper, we apply our method to 22 $\gamma$ Dor stars with g-mode multiplets split by rotation. The rotation periods are similar to the g-mode period spacings, causing the multiplets to overlap. To clarify the overlapping patterns, we use the \'echelle diagram and introduce a `copy-shift' diagram to discern and measure the splittings. The first observational relation between slopes and splittings is shown. The slope deviates from zero when the splitting increases, as the theory predicts. We found that what appears to be rotational splittings in two stars is in fact caused by two nearly-identical overlapping patterns from binaries.

The Serpens filament, as one of the nearest infrared dark clouds, is regarded as a pristine filament at a very early evolutionary stage of star formation. In order to study its molecular content and dynamical state, we mapped this filament in seven species. Among them, HCO$^{+}$, HNC, HCN, and CS show self-absorption, while C$^{18}$O is most sensitive to the filamentary structure. A kinematic analysis demonstrates that this filament forms a velocity-coherent (trans-)sonic structure, a large part of which is one of the most quiescent regions in the Serpens cloud. Widespread C$^{18}$O depletion is found throughout the Serpens filament. Based on the Herschel dust-derived H$_{2}$ column density map, the line mass of the filament is 36--41~M$_{\odot}$~pc$^{-1}$, and its full width at half maximum is 0.17$\pm$0.01~pc, while its length is ~1.6~pc. The inner radial column density profile of this filament can be well fitted with a Plummer profile with an exponent of 2.2$\pm$0.1, a scale radius of $0.018\pm 0.003$ pc, and a central density of $(4.0\pm 0.8)\times 10^{4}$~cm$^{-3}$. The Serpens filament appears to be slightly supercritical. The widespread blue-skewed HNC and CS line profiles and HCN hyperfine line anomalies across this filament indicate radial infall in parts of the Serpens filament. C$^{18}$O velocity gradients also indicate accretion flows along the filament. The velocity and density structures suggest that such accretion flows are likely due to a longitudinal collapse parallel to the filament's long axis. Both the radial infall rate and the longitudinal accretion rate along the Serpens filament are lower than all previously reported values in other filaments. This indicates that the Serpens filament lies at an early evolutionary stage when collapse has just begun, or that thermal and non-thermal support are effective in providing support against gravity.

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small.

Although the Galactic supernova rate is about 2 per century, only few supernova remnants are associated with historical records. There are a few ancient Chinese records of "guest stars" that are probably sightings of supernovae for which the associated supernova remnant is not established. Here we present an X-ray study of the supernova remnant G7.7-3.7, as observed by XMM-Newton, and discuss its probable association with the guest star of 386 CE. This guest star occurred in the ancient Chinese asterism Nan-Dou, which is part of Sagittarius. The X-ray morphology of G7.7-3.7 shows an arc-like feature in the SNR south, which is characterized by an under-ionized plasma with sub-solar abundances, a temperature of 0.4--0.8 keV, and a density of $\sim 0.5(d/4 kpc)^{-0.5}$ cm$^{-3}$. A small shock age of $1.2\pm 0.6 (d/4 kpc)^{0.5}$ kyr is inferred from the low ionization timescale of $2.4^{+1.1}_{-1.3}\times 10^{10}$ cm$^{-3}$ s of the X-ray arc. The low foreground absorption ($N_H=3.5\pm0.5\times 10^{21}$ cm$^{-2}$) of G7.7-3.7 made the supernova explosion visible to the naked eyes on the Earth. The position of G7.7-3.7 is consistent with the event of 386 CE, and the X-ray properties suggest that also its age is consistent. Interestingly, the association between G7.7-3.7 and guest star 386 would suggest the supernova to be a low-luminosity supernova, in order to explain the not very long visibility (2--4 months) of the guest star.

We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We investigated 64 pairs of interacting-CME events identified by the simultaneous observations of SOHO and STEREO spacecraft from 2010 January to 2014 August, to examine the relationship between the large SEP events in the energy of 25-60MeV and the properties of the interacting CMEs. We found that during CME interactions the large SEP events in this study were all generated by CMEs with the presence of enhanced type II radio bursts, which also have wider longitudinal distributions comparing to events with the absence of type II radio burst or its enhancement (almost associated with small SEP events). It seems that the signature of type II radio bursts enhancement is a good discriminator between large SEP and small or none SEP event producers during CME interactions. The type II radio burst enhancement is more likely to be generated by CME interactions, with the main CME having larger speed (v), angular width (WD), mass (m) and kinetic energy (E_k), that taking over the preceding CMEs which also have higher v, WD, m and E_k, than those preceding CMEs in CME pairs missing the type II radio bursts or enhancements. Our analysis also revealed that the intensities of associated SEP events correlate negatively with the intersection height of the two CMEs. Most of type-II-enhanced events and SEP events are coincidentally and almost always made by the fast and wide main CMEs that sweeping fully over a relatively slower and narrower preceding CMEs. We suggest that a fast CME with enough energy completely overtaking a relatively narrower preceding CME, especially in low height, can drive more energetic shock signified by the enhanced type II radio bursts. The shock may accelerate ambient particles and lead to large SEP event more easily.

Stars of spectral type Oe are very rare. To date, only 13 Oe stars have been identified within our Galaxy. In this paper, we present six new Oe stars and four new B0e stars found in LAMOST DR5. Repeated spectral observations of the same Oe stars show some emission line variability. The H$\beta$ emission of TYC 4801-17-1 shows rapid V/R variation. Phase lags in the V/R ratio of TYC 4801-17-1 spectra are also seen. We found the unusual O4.5 star RL 128 is an Oe star with variable H$\alpha$ intensity and its Ca II triplet emission appears when H$\alpha$ emission reaches maximum intensity. These newly identified early type Oe and B0e stars significantly increase the known sample.

We consider the Palatini formulation of $f(R,T)$ gravity theory, in which a nonminimal coupling between the Ricci scalar and the trace of the energy-momentum tensor is introduced, by considering the metric and the affine connection as independent field variables. The field equations and the equations of motion for massive test particles are derived, and we show that the independent connection can be expressed as the Levi-Civita connection of an auxiliary, energy-momentum trace dependent metric, related to the physical metric by a conformal transformation. Similarly to the metric case, the field equations impose the non-conservation of the energy-momentum tensor. We obtain the explicit form of the equations of motion for massive test particles in the case of a perfect fluid, and the expression of the extra-force, which is identical to the one obtained in the metric case. The thermodynamic interpretation of the theory is also briefly discussed. We investigate in detail the cosmological implications of the theory, and we obtain the generalized Friedmann equations of the $f(R,T)$ gravity in the Palatini formulation. Cosmological models with Lagrangians of the type $f=R-\alpha ^2/R+g(T)$ and $f=R+\alpha ^2R^2+g(T)$ are investigated. These models lead to evolution equations whose solutions describe accelerating Universes at late times.

Astrophysical systems, such as clumps that form star clusters share a density profile that is close to $\rho \sim r^{-2}$. We prove analytically this density profile is the result of the scale-free nature of the gravitational collapse. Therefore, it should emerge in many different situations as long as gravity is dominating the evolution for a period that is comparable or longer than the free-fall time, and this does not necessarily imply an isothermal model, as many have previously believed. To describe the collapse process, we construct a model called the turbulence-regulated gravitational collapse model, where turbulence is sustained by accretion and dissipates in roughly a crossing time. We demonstrate that a $\rho \sim r^{-2}$ profile emerges due to the scale-free nature the system. In this particular case, the rate of gravitational collapse is regulated by the rate at which turbulence dissipates the kinetic energy such that the infall speed can be $20 - 50 \%$ of the free-fall speed(which also depends on the interpretation of the crossing time based on simulations of driven turbulence). These predictions are consistent with existing observations, which suggests that these clumps are in the stage of turbulence-regulated gravitational collapse. Our analysis provides a unified description of gravitational collapse in different environments.

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

The convergence peaks, constructed from galaxy shape measurement in weak lensing, is a powerful probe of cosmology as the peaks can be connected with the underlined dark matter haloes. However the capability of convergence peak statistic is affected by the noise in galaxy shape measurement, signal to noise ratio as well as the contribution from the projected mass distribution from the large-scale structures along the line of sight (LOS). In this paper we use the ray-tracing simulation on a curved sky to investigate the correspondence between the convergence peak and the dark matter haloes at the LOS. We find that, in case of no noise and for source galaxies at $z_{\rm s}=1$, more than $65\%$ peaks with $\text{SNR} \geq 3$ (signal to noise ratio) are related to more than one massive haloes with mass larger than $10^{13} {\rm M}_{\odot}$. Those massive haloes contribute $87.2\%$ to high peaks ($\text{SNR} \geq 5$) with the remaining contributions are from the large-scale structures. On the other hand, the peaks distribution is skewed by the noise in galaxy shape measurement, especially for lower SNR peaks. In the noisy field where the shape noise is modelled as a Gaussian distribution, about $60\%$ high peaks ($\text{SNR} \geq 5$) are true peaks and the fraction decreases to $20\%$ for lower peaks ($ 3 \leq \text{SNR} < 5$). Furthermore, we find that high peaks ($\text{SNR} \geq 5$) are dominated by very massive haloes larger than $10^{14} {\rm M}_{\odot}$.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

Obliquity variability could play an important role in the climate and habitability of a planet. Orbital modulations caused by planetary companions and the planet's spin axis precession due to the torque from the host star may lead to resonant interactions and cause large-amplitude obliquity variability. Here we consider the spin axis dynamics of Kepler-62f and Kepler-186f, both of which reside in the habitable zone around their host stars. Using {\emph{N}}-body simulations and secular numerical integrations, we describe their obliquity evolution for particular realizations of the planetary systems. We then use a generalized analytic framework to characterize regions in parameter space where the obliquity is variable with large amplitude. We find that the locations of variability are fine-tuned over the planetary properties and system architecture in the lower-obliquity regimes ($\lesssim 40^\circ$). As an example, assuming a rotation period of 24 hr, the obliquities of both Kepler-62f and Kepler-186f are stable below $\sim 40^\circ$, whereas the high-obliquity regions ($60^\circ - 90^\circ$) allow moderate variabilities. However, for some other rotation periods of Kepler-62f or Kepler-186f, the lower-obliquity regions could become more variable owing to resonant interactions. Even small deviations from coplanarity (e.g. mutual inclinations $\sim 3^\circ$) could stir peak-to-peak obliquity variations up to $\sim 20^\circ$. Undetected planetary companions and/or the existence of a satellite could also destabilize the low-obliquity regions. In all cases, the high-obliquity region allows for moderate variations, and all obliquities corresponding to retrograde motion (i.e. $> 90^\circ$) are stable.

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

Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of the population of distant TNOs, similar to the observed sample. However, many aspects of the rich underlying dynamical processes induced by such a distant eccentric perturber have not been fully explored. We report the results of our investigation of the dynamics of coplanar test-particles that interact with a massive body on an circular orbit (Neptune) and a massive body on a more distant, highly eccentric orbit (the putative Planet Nine). We find that a detailed examination of our idealized simulations affords tremendous insight into the rich test-particle dynamics that are possible. In particular, we find that chaos and resonance overlap plays an important role in particles' dynamical evolution. We develop a simple mapping model that allows us to understand, in detail, the web of overlapped mean-motion resonances explored by chaotically evolving particles. We also demonstrate that gravitational interactions with Neptune can have profound effects on the orbital evolution of particles. Our results serve as a starting point for a better understanding of the dynamical behavior observed in more complicated simulations that can be used to constrain the mass and orbit of Planet 9.

The Fourier spectrum of the $\gamma$-Dor variable KIC 5608334 shows remarkable frequency groups at $\sim$3, $\sim$6, $\sim$9, and 11--12\,d$^{-1}$. We explain the four frequency groups as prograde sectoral g modes in a rapidly rotating star. Frequencies of intermediate-to-high radial order prograde sectoral g modes in a rapidly rotating star are proportional to $|m|$ (i.e., $\nu \propto |m|$) in the co-rotating frame as well as in the inertial frame. This property is consistent with the frequency groups of KIC 5608334 as well as the period vs. period-spacing relation present within each frequency group, if we assume a rotation frequency of $2.2$\,d$^{-1}$, and that each frequency group consists of prograde sectoral g modes of $|m| = 1, 2, 3,$ and 4, respectively. In addition, these modes naturally satisfy near-resonance conditions $\nu_i\approx\nu_j+\nu_k$ with $m_i=m_j+m_k$. We even find exact resonance frequency conditions (within the precise measurement uncertainties) in many cases, which correspond to combination frequencies.

A complete understanding of the onset and subsequent evolution of confined flares has not been achieved. Earlier studies mainly analyzed disk events so as to reveal their magnetic topology and cause of confinement. In this study, taking advantage of a tandem of instruments working at different wavelengths of X-rays, EUVs, and microwaves, we present dynamic details of a confined flare observed on the northwestern limb of the solar disk on July 24th, 2016. The entire dynamic evolutionary process starting from its onset is consistent with a loop-loop interaction scenario. The X-ray profiles manifest an intriguing double-peak feature. From spectral fitting, it is found that the first peak is non-thermally dominated while the second peak is mostly multi-thermal with a hot (~10 MK) and a super-hot (~30 MK) component. This double-peak feature is unique in that the two peaks are clearly separated by 4 minutes, and the second peak reaches up to 25-50 keV; in addition, at energy bands above 3 keV the X-ray fluxes decline significantly between the two peaks. This, together with other available imaging and spectral data, manifest a two-stage energy release process. A comprehensive analysis is carried out to investigate the nature of this two-stage process. We conclude that the second stage with the hot and super-hot sources mainly involves direct heating through loop-loop reconnection at a relatively high altitude in the corona. The uniqueness of the event characteristics and complete data set make the study a nice addition to present literature on solar flares.

Accurate reconstruction of the spatial distributions of the Point Spread Function (PSF) is crucial for high precision cosmic shear measurements. Nevertheless, current methods are not good at recovering the PSF fluctuations of high spatial frequencies. In general, the residual PSF fluctuations are spatially correlated, therefore can significantly contaminate the correlation functions of the weak lensing signals. We propose a method to correct for this contamination statistically, without any assumptions on the PSF and galaxy morphologies or their spatial distribution. We demonstrate our idea with the data from the W2 field of CFHTLenS.

In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.

Compact stellar binaries are expected to survive in the dense environment of the Galactic Center. The stable binaries may undergo Kozai-Lidov oscillations due to perturbations from the central supermassive black hole (Sgr A*), yet the General Relativistic precession can suppress the Kozai-Lidov oscillations and keep the stellar binaries from merging. However, it is challenging to resolve the binary sources and distinguish them from single stars. The close separation of the stable binaries allow higher eclipse probabilities. Here, we consider the massive star SO2 as an example and calculate the probability of detecting eclipses, assuming it is a binary. We find that the eclipse probability is ~ 30-50%, reaching higher values when the stellar binary is more eccentric or highly inclined relative to its orbit around Sgr A*.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

The interplay between gravity, turbulence and the magnetic field determines the evolution of the molecular ISM and the formation of the stars. In spite of growing interests, there remains a lack of understanding of the importance of magnetic field over multiple scales. We derive the magnetic energy spectrum -- a measure that constraints the multi-scale distribution of the magnetic energy, and compare it with the gravitational energy spectrum derived in Li \& Burkert (2016). In our formalism, the gravitational energy spectrum is purely determined by the surface density PDF, and the magnetic energy spectrum is determined by both the surface density PDF and the magnetic-field-density relation. If regions have density PDFs close to $P(\Sigma)\sim \Sigma^{-2}$ and a universal magnetic field-density relation $B\sim \rho^{1/2}$, we expect a multi-scale near equipartition between gravity and the magnetic fields. This equipartition is found to be true in NGC6334 where estimates of magnetic fields over multiple scales (from 0.1 pc to a few parsec) are available. However, the current observations are still limited in sample size. In the future, it is necessary to obtain multi-scale measurements of magnetic fields from different clouds with different surface density PDFs and apply our formalism to further study the gravity-magnetic field interplay.

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

Observations for the $\delta$ Scuti star AN Lyn have been made between 2008 and 2016 with the 85-cm telescope at Xinglong station of National Astronomical Observatories of China, the 84-cm telescope at SPM Observatory of Mexico and the Nanshan One meter Wide field Telescope of Xinjiang Observatory of China. Data in $V$ in 50 nights and in $R$ in 34 nights are obtained in total. The bi-site observations from both Xinglong Station and SPM Observatory in 2014 are analyzed with Fourier Decomposition to detect pulsation frequencies. Two independent frequencies are resolved, including one non-radial mode. A number of stellar model tracks are constructed with the MESA code and the fit of frequencies leads to the best-fit model with the stellar mass of $M = 1.70 \pm 0.05~\mathrm{M_{\odot}}$ , the metallicity abundance of $Z = 0.020 \pm 0.001$, the age of $1.33 \pm 0.01$ billion years and the period change rate $1/P\cdot \mathrm{d}P/\mathrm{d}t =1.06 \times 10^{-9} ~\mathrm{yr^{-1}} $, locating the star at the evolutionary stage close to the terminal age main sequence (TAMS). The O-C diagram provides the period change rate of $1/P \cdot \mathrm{d}P /\mathrm{d}t =4.5(8)\times10^{-7}~\mathrm{yr^{-1}}$. However, the period change rate calculated from the models is smaller in two orders than the one derived from the O-C diagram. Together with the sinusoidal function signature, the period variations are regarded to be dominated by the light-travel time effect of the orbital motion of a three-body system with two low-luminosity components, rather than the stellar evolutionary effect.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

Different from Milky-Way-like galaxies, discs of gas-rich galaxies are clumpy. It is believed that the clumps form because of gravitational instability. However, a necessary condition for gravitational instability to develop is that the disc must dissipate its kinetic energy effectively, this energy dissipation (also called cooling) is not well-understood. We propose that collisions (coagulation) between molecular clouds dissipate the kinetic energy of the discs, which leads to a dynamical cooling. The effectiveness of this dynamical cooling is quantified by the dissipation parameter $D$, which is the ratio between the free-fall time $t_{\rm ff}\approx 1/ \sqrt{G \rho_{\rm disc}}$ and the cooling time determined by the cloud collision process $t_{\rm cool}$. This ratio is related to the ratio between the mean surface density of the disc $\Sigma_{\rm disc}$ and the mean surface density of molecular clouds in the disc $\Sigma_{\rm cloud}$. When $D <1/3$ (which roughly corresponds to $\Sigma_{\rm disc} < 1/3 \Sigma_{\rm cloud}$), cloud collision cooling is inefficient, and fragmentation is suppressed. When $D > 1/3$ (which roughly corresponds to $\Sigma_{\rm disc} > 1/3 \Sigma_{\rm cloud}$), cloud-cloud collisions lead to a rapid cooling through which clumps form. On smaller scales, cloud-cloud collisions can drive molecular cloud turbulence. This dynamical cooling process can be taken into account in numerical simulations as a subgrid model to simulate the global evolution of disc galaxies.

LS I +61${^\circ}$303 is a $\gamma$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate these microflaring activity in LS I +61${^\circ}$303 we observed the source with the 100-m Effelsberg radio telescope at 4.85, 8.35, and 10.45 GHz and performed timing analysis on the obtained data. Radio oscillations of 15 hours time scales are detected at all three frequencies. We also compare the spectral index evolution of radio data to that of the photon index of GeV data observed by Fermi-LAT. We conclude that the observed QPO could result from multiple shocks in a jet.

We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.

We study the fragmentation of eight massive clumps using data from ATLASGAL 870 $\mu$m, SCUBA 850 and 450 $\mu$m, PdBI 1.3 and 3.5 mm, and probe the fragmentation from 1 pc to 0.01 pc scale. We find that the masses and the sizes of our objects follow $M \sim r^{1.68\pm0.05}$. The results are in agreements with the predictions of Li (2017) where $M \sim r^{5/3}$. Inside each object, the densest structures seem to be centrally condensed, with $\rho(r)\sim r^{-2}$. Our observational results support a scenario where molecular gas in the Milky Way is supported by a turbulence characterized by a constant energy dissipation rate, and gas fragments like clumps and cores are structures which are massive enough to be dynamically detached from the ambient medium.

Single-exposure spectra in large spectral surveys are valuable for time domain studies such as stellar variability, but there is no available method to eliminate cosmic rays for single-exposure, multi-fiber spectral images. In this paper, we describe a new method to detect and remove cosmic rays in multi-fiber spectroscopic single exposures. Through the use of two-dimensional profile fitting and a noise model that considers the position-dependent errors, we successfully detect as many as 80% of the cosmic rays and correct the cosmic ray polluted pixels to an average accuracy of 97.8%. Multiple tests and comparisons with both simulated data and real LAMOST data show that the method works properly in detection rate, false detection rate, and validity of cosmic ray correction.

Sky subtraction is the key technique in data reduction of multi-fiber spectra. Knowledge of the related instrument character is necessary to determine the method adopted in sky subtraction. In this study, we described the sky subtraction method designed for LAMOST(Large sky Area Multi-Object fiber Spectroscopic Telescope) survey. The method has been intergrated into LAMOST 2D Pipeline v2.6 and applied to data of LAMOST DR3 and later. For LAMOST, sky emission line calibration is used to alleviate the position-dependent (thus time-dependent) ~4% fiber throughput uncertainty and the small wavelength instability (0.1\AA ) during observation. PCA (Principal Component Analysis) sky subtraction further reduces 25% of the sky line residual of the OH lines in the red part of the LAMOST spectra after the mater sky spectrum, which is derived from a B-spline fit of 20 sky fibers in each spectrograph, is adjusted by sky emission line and subtracted from each fiber. Further analysis shows that our wavelength calibration accuracy is about 4.5km/s, and the average sky subtraction residuals are about 3% for sky emission lines and 3% for continuum region. The relative sky subtraction residuals vary with the moon light background brightness, could reach as low as 1.5% for the sky emission line regions in the dark night. Tests on the F stars of both similar sky emission line strength and similar object continuum intensity show that the sky emission line residual of LAMOST is smaller than those of SDSS survey.

The secular approximation of the hierarchical three body systems has been proven to be very useful in addressing many astrophysical systems, from planets, stars to black holes. In such a system two objects are on a tight orbit, and the tertiary is on a much wider orbit. Here we study the dynamics of a system by taking the tertiary mass to zero and solve the hierarchical three body system up to the octupole level of approximation. We find a rich dynamics that the outer orbit undergoes due to gravitational perturbations from the inner binary. The nominal result of the precession of the nodes is mostly limited for the lowest order of approximation, however, when the octupole-level of approximation is introduced the system becomes chaotic, as expected, and the tertiary oscillates below and above 90deg, similarly to the non-test particle flip behavior (e.g., Naoz 2016). We provide the Hamiltonian of the system and investigate the dynamics of the system from the quadrupole to the octupole level of approximations. We also analyze the chaotic and quasi-periodic orbital evolution by studying the surfaces of sections. Furthermore, including general relativity, we show case the long term evolution of individual debris disk particles under the influence of a far away interior eccentric planet. We show that this dynamics can naturally result in retrograde objects and a puffy disk after a long timescale evolution (few Gyr) for initially aligned configuration.

Reconstruction of the point spread function (PSF) is a critical process in weak lensing measurement. We develop a real-data based and galaxy-oriented pipeline to compare the performances of various PSF reconstruction schemes. Making use of a large amount of the CFHTLenS data, the performances of three classes of interpolating schemes - polynomial, Kriging, and Shepard - are evaluated. We find that polynomial interpolations with optimal orders and domains perform the best. We quantify the effect of the residual PSF reconstruction error on shear recovery in terms of the multiplicative and additive biases, and their spatial correlations using the shear measurement method of Zhang et al. (2015). We find that the impact of PSF reconstruction uncertainty on the shear-shear correlation can be significantly reduced by cross correlating the shear estimators from different exposures. It takes only 0.2 stars (SNR > 100) per square arcmin on each exposure to reach the best performance of PSF interpolation, a requirement that is generally satisfied in the CFHTlenS data.

The cold-dense plasma is occasionally detected in the solar wind with in situ data, but the source of the cold-dense plasma remains illusive. Interchange reconnections (IRs) between closed fields and nearby open fields are well known to contribute to the formation of solar winds. We present a confined filament eruption associated with a puff-like coronal mass ejection (CME) on 2014 December 24. The filament underwent successive activations and finally erupted, due to continuous magnetic flux cancellations and emergences. The confined erupting filament showed a clear untwist motion, and most of the filament material fell back. During the eruption, some tiny blobs escaped from the confined filament body, along newly-formed open field lines rooted around the south end of the filament, and some bright plasma flowed from the north end of the filament to remote sites at nearby open fields. The newly-formed open field lines shifted southward with multiple branches. The puff-like CME also showed multiple bright fronts and a clear southward shift. All the results indicate an intermittent IR existed between closed fields of the confined erupting filament and nearby open fields, which released a portion of filament material (blobs) to form the puff-like CME. We suggest that the IR provides a possible source of cold-dense plasma in the solar wind.

We present candidate members of the Pal 5, GD-1, Cetus Polar, and Orphan tidal stellar streams found in LAMOST DR3, SDSS DR9 and APOGEE catalogs. In LAMOST DR3, we find 20, 4, 24 high confidence candidates of tidal streams GD-1, Cetus Polar and Orphan respectively. We also list from the SDSS DR9 spectroscopic catalog 59, 118, 10 high confidence candidates of tidal streams Cetus Polar, Orphan and Pal 5, respectively. Furthermore, we find 7 high confidence candidates of the Pal 5 tidal stream in APOGEE data. Compared with SDSS, the new candidates from LAMOST DR3 are brighter, so that together, more of the color-magnitude diagram, including the giant branch can be explored. Analysis of SDSS data shows that there are 3 metallicity peaks of the Orphan stream and also shows some spatial separation. LAMOST data confirms multiple metallicities in this stream. The metallicity, given by the higher resolution APOGEE instrument, of the Pal 5 tidal stream is [Fe/H] $\sim -1.2$, higher than that given earlier by SDSS spectra. Many previously unidentified stream members are tabulated here for the first time, along with existing members, allowing future researchers to further constrain the orbits of these objects as they move within the Galaxy's dark matter potential.

Both General Relativistic (GR) precession and the Lidov-Kozai mechanism, separately, are known to play an important role in the orbital evolution of solar system bodies. Previous works have studied these two mechanisms independently in great detail. However, both these phenomena occurring at the same time in real solar system bodies have rarely been explored. In this work, we find a continuum connecting the GR precession dominant and Lidov-Kozai like mechanism dominant regimes, i.e. an intermediate regime where the competing effects of GR precession and Lidov-Kozai like oscillations co-exist simultaneously. We find some real examples in the solar system in this intermediate regime. Moreover we identify a rare example amongst them, comet 96P/Machholz 1, which shows significant changes in the rates of GR precession (an order of magnitude higher than Mercury's GR precession rate) due to sungrazing and sun colliding phases induced by Lidov-Kozai like oscillations. This comet's combination of orbital elements and initial conditions (at the present epoch) favour this measurable rapid change in GR precession (at some points peaking up to 60 times Mercury's GR precession rate) along with prograde-retrograde inclination flip (due to Lidov-Kozai like oscillations). Similar tests are performed for hundreds of bodies lying in the moderately low perihelion distance and moderately low semi-major axis phase space in the solar system, the present lowest perihelion distance asteroid 322P/SOHO 1, and further examples connected with 96P/Machholz 1 namely, the Marsden and Kracht families of sungrazing comets plus low perihelion meteoroid streams like Daytime Arietids (ARI) and Southern Delta Aquariids (SDA).

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

AST3-1 is the second-generation wide-field optical photometric telescope dedicated to time domain astronomy at Dome A, Antarctica. Here we present the results of $i$ band images survey from AST3-1 towards one Galactic disk field. Based on time-series photometry of 92,583 stars, 560 variable stars were detected with $i$ magnitude $\leq$ 16.5 mag during eight days of observations; 339 of these are previously unknown variables. We tentatively classify the 560 variables as 285 eclipsing binaries (EW, EB, EA), 27 pulsating variable stars ($\delta$~Scuti, $\gamma$~Doradus, $\delta$~Cephei variable and RR Lyrae stars) and 248 other types of variables (unclassified periodic, multi-periodic and aperiodic variable stars). Among the eclipsing binaries, 34 show O'Connell effects. One of the aperiodic variables shows a plateau light curve and another one shows a secondary maximum after peak brightness. We also detected a complex binary system with RS CVn-like light curve morphology; this object is being followed-up spectroscopically using the Gemini South telescope.

Context. Infrared dark clouds (IRDCs) are ubiquitous in the Milky Way, yet they play a crucial role in breeding newly-formed stars. Aims. With the aim of further understanding the dynamics, chemistry, and evolution of IRDCs, we carried out multi-wavelength observations on a small sample. Methods. We performed new observations with the IRAM 30 m and CSO 10.4 m telescopes, with tracers ${\rm HCO^+}$, HCN, ${\rm N_2H^+}$, ${\rm C^{18}O}$, DCO$^+$, SiO, and DCN toward six IRDCs G031.97+00.07, G033.69-00.01, G034.43+00.24, G035.39-00.33, G038.95-00.47, and G053.11+00.05. Results. We investigated 44 cores including 37 cores reported in previous work and seven newly-identified cores. Toward the dense cores, we detected 6 DCO$^+$, and 5 DCN lines. Using pixel-by-pixel spectral energy distribution (SED) fits of the $\textit{Herschel}$ 70 to 500 $\mu$m, we obtained dust temperature and column density distributions of the IRDCs. We found that ${\rm N_2H^+}$ emission has a strong correlation with the dust temperature and column density distributions, while ${\rm C^{18}O}$ showed the weakest correlation. It is suggested that ${\rm N_2H^+}$ is indeed a good tracer in very dense conditions, but ${\rm C^{18}O}$ is an unreliable one, as it has a relatively low critical density and is vulnerable to freezing-out onto the surface of cold dust grains. The dynamics within IRDCs are active, with infall, outflow, and collapse; the spectra are abundant especially in deuterium species. Conclusions. We observe many blueshifted and redshifted profiles, respectively, with ${\rm HCO^+}$ and ${\rm C^{18}O}$ toward the same core. This case can be well explained by model "envelope expansion with core collapse (EECC)".

We study the evolution of structures in turbulent, self-gravitating media, and present an analytical criterion $M_{\rm crit} \approx \epsilon_{\rm cascade}^{2/3} \eta^{-2/3} G^{-1} l^{5/3}$ (where $M_{\rm crit}$ is the critical mass, $l$ is the scale, $\epsilon_{\rm cascade}\approx \eta \sigma_{\rm v}^3 /l $ is the turbulence energy dissipation rate of the ambient medium, $G$ is the gravitational constant, $\sigma_{\rm v}$ is the velocity dispersion, $l$ is the scale and $\eta\approx 0.2$ is an efficiency parameter) for an object to undergo quasi-isolated gravitational collapse. The criterion also defines the critical scale ($l_{\rm crit} \approx \epsilon_{\rm cascade}^{1/2} \eta^{-1/2} G^{-3/4} \rho^{-3/4}$) for turbulent gravitational instability to develop. The analytical formalism explains the size dependence of the masses of the progenitors of star clusters ($M_{\rm cluster} \sim R_{\rm cluster}^{1.67}$) in our Galaxy.

The dynamical processes that control star formation in molecular clouds are not well understood, and in particular, it is unclear if rotation plays a major role in cloud evolution. We investigate the importance of rotation in cloud evolution by studying the kinematic structure of a spiral-shaped Galactic molecular cloud G052.24+00.74. The cloud belongs to a large filament, and is stretching over ~ 100 pc above the Galactic disk midplane. The spiral-shaped morphology of the cloud suggests that the cloud is rotating. We have analysed the kinematic structure of the cloud, and study the fragmentation and star formation. We find that the cloud exhibits a regular velocity pattern along west-east direction - a velocity shift of ~ 10 km/s at a scale of ~ 30 pc. The kinematic structure of the cloud can be reasonably explained by a model that assumes rotational support. Similarly to our Galaxy, the cloud rotates with a prograde motion. We use the formalism of Toomre (1964) to study the cloud's stability, and find that it is unstable and should fragment. The separation of clumps can be consistently reproduced assuming gravitational instability, suggesting that fragmentation is determined by the interplay between rotation and gravity. Star formation occurs in massive, gravitational bound clumps.

Gravity plays a determining role in the evolution of the molecular ISM. In \citet{2016arXiv160304342L}, we proposed a measure called gravitational energy spectrum to quantify the importance of gravity on multiple physical scales. In this work, using a wavelet-based decomposition technique, we derive the gravitational energy spectra of the Orion A and the Orion B molecular cloud from observational data. The gravitational energy spectra exhibit power-law-like behaviours. From a few pc down to $\sim 0.1 $ pc scale, the Orion A and Orion B molecular cloud have $E_{\rm p}(k)\sim k^{-1.88}$ and $E_{\rm p}(k)\sim k^{-2.09}$, respectively. These scaling exponents are close to the scaling exponents of the kinetic energy power spectrum of compressible turbulence (where $E\sim k^{-2}$), with a near-equipartition of turbulent versus gravitational energy on multiple scales. This provides a clear evidence that gravity is able to counteract effectively against turbulent motion for these length scales. The results confirm our earlier analytical estimates. For the Orion A molecular cloud, gravity inevitably dominates over turbulence inside the cloud. Our results provide a clear observational proof that gravity is playing a determining role in the evolution these molecular clouds from the cloud scale down to $\sim 0.1\;\rm pc$. { However, turbulence is likely to dominate in clouds like California. } The method is general and should be applicable to all the astrophysical problems where gravity plays a role.

In large SEP events, ions can be accelerated at CME-driven shocks to very high energies. Spectra of heavy ions in many large SEP events show features such as roll-overs or spectral breaks. In some events when the spectra are plotted in energy/nucleon they can be shifted relative to each other to make the spectral breaks align. The amount of shift is charge-to-mass ratio (Q/A) dependent and varies from event to event. This can be understood if the spectra of heavy ions are organized by the diffusion coefficients (Cohen et al., 2005). In the work of Li et al. (2009), the Q/A dependences of the scaling is related to shock geometry when the CME-driven shock is close to the Sun. For events where multiple in-situ spacecraft observations exist, one may expect that different spacecraft are connected to different portions of the CME-driven shock that have different shock geometries, therefore yielding different Q/A dependence. In this work, we examine one SEP event which occurred on 2013 November 4. We study the Q/A dependence of the energy scaling for heavy ion spectra using Helium, oxygen, and iron ions. Observations from STEREO-A, STEREO-B and ACE are examined. We find that the scalings are different for different spacecraft. We suggest that this is because ACE, STEREO-A and STEREO- B are connected to different parts of the shock that have different shock geometries. Our analysis indicates that studying the Q/A scaling of in-situ particle spectra can serve as a powerful tool to remotely examine the shock geometry for large SEP events.

The new discoveries of circumbinary planetary systems shed light on the understanding of planetary system formation. Learning the architectural properties of these systems is essential for constraining the different formation mechanisms. We first revisit the stability limit of circumbinary planets. Next, we focus on eclipsing stellar binaries and obtain an analytical expression for the transit probability in a realistic setting, where finite observation period and planetary orbital precession are included. Then, understanding of the architectural properties of the currently observed transiting systems is refined, based on Bayesian analysis and a series of hypothesis tests. We find 1) it is not a selection bias that the innermost planets reside near the stability limit for eight of the nine observed systems, and this is consistent with a log uniform distribution of the planetary semi-major axis; 2) it is not a selection bias that the planetary and stellar orbits are nearly coplanar ($\lesssim 3^\circ$), and this together with previous studies may imply an occurrence rate of circumbinary planets similar to that of single star systems; 3) the dominance of observed circumbinary systems with only one transiting planet may be caused by selection effects; 4) formation mechanisms involving Lidov-Kozai oscillations, which may produce misalignment and large separation between planet and stellar binaries, are consistent with the lack of transiting circumbinary planets around short-period stellar binaries, in agreement with previous studies. As a consequence of 4), eclipse timing variations may better suit the detection of planets in such configurations.

Type III and type-III-like radio bursts are produced by energetic electron beams guided along coronal magnetic fields. As a variant of type III bursts, Type N bursts appear as the letter "N" in the radio dynamic spectrum and reveal a magnetic mirror effect in coronal loops. Here, we report a well-observed N-shaped burst consisting of three successive branches at metric wavelength with both fundamental and harmonic components and a high brightness temperature ($>$10$^9$ K). We verify the burst as a true type N burst generated by the same electron beam from three aspects of the data. First, durations of the three branches at a given frequency increase gradually, may due to the dispersion of the beam along its path. Second, the flare site, as the only possible source of non-thermal electrons, is near the western feet of large-scale closed loops. Third, the first branch and the following two branches are localized at different legs of the loops with opposite sense of polarization. We also find that the sense of polarization of the radio burst is in contradiction to the O-mode and there exists a fairly large time delay ($\sim$3-5 s) between the fundamental and harmonic components. Possible explanations accounting for these observations are presented. Assuming the classical plasma emission mechanism, we can infer coronal parameters such as electron density and magnetic field near the radio source and make diagnostics on the magnetic mirror process.

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

In this work, we mainly study the magnification relations of quad lens models for cusp, fold and cross configurations. By dividing and ray-tracing in different image regions, we numerically derive the positions and magnifications of the four images for a point source lying inside of the astroid caustic. Then, based on the magnifications, we calculate the signed cusp and fold relations for the singular isothermal elliptical lenses. The signed fold relation map has positive and negative regions, and the positive region is usually larger than the negative region as has been confirmed before. It can also explain that for many observed fold image pairs, the fluxes of the Fermat minimum images are apt to be larger than those of the saddle images. We define a new quantity cross relation which describes the magnification discrepancy between two minimum images and two saddle images. Distance ratio is also defined as the ratio of the distance of two saddle images to that of two minimum images. We calculate the cross relations and distance ratios for nine observed Einstein crosses. In theory, for most of the quad lens models, the cross relations decrease as the distance ratios increase. In observation, the cross relations of the nine samples do not agree with the quad lens models very well, nevertheless, the cross relations of the nine samples do not give obvious evidence for anomalous flux ratio as the cusp and fold types do. Then, we discuss several reasons for the disagreement, and expect good consistencies for more precise observations and better lens models in the future.

(Abridged) We derive an analytical formula which provides estimates on multiscale gravitational energy distribution using the observed surface density PDF. Our analytical formalism also enables one to convert the observed column density PDF into an estimated volume density PDF, and to obtain average radial density profile $\rho(r)$. For a region with $N_{\rm col} \sim N^{-\gamma_{\rm N}}$, the gravitational energy spectra is $E_{\rm p}(k)\sim k^{-4(1 - 1/\gamma_{\rm N})}$. We apply the formula to observations of molecular clouds, and find that a scaling index of $-2$ of the surface density PDF implies that $\rho \sim r^{-2}$ and $E_{\rm p}(k) \sim k^{-2}$. The results are valid from the cloud scale (a few parsec) to around $\sim 0.1 \;\rm pc$. Because of the resemblance the scaling index of the gravitational energy spectrum and the that of the kinetic energy power spectrum of the Burgers turbulence (where $E\sim k^{-2}$), our result indicates that gravity can act effectively against turbulence over a multitude of physical scales. This is the critical scaling index which divides molecular clouds into two categories: clouds like Orion and Ophiuchus have shallower power laws, and the amount of gravitational energy is too large for turbulence to be effective inside the cloud. Because gravity dominates, we call this type of cloud g-type clouds. On the other hand, clouds like the California molecular cloud and the Pipe nebula have steeper power laws, and turbulence can overcome gravity if it can cascade effectively from the large scale. We call this type of cloud t-type clouds. The analytical formula can be used to determine if gravity is dominating cloud evolution when the column density probability distribution function (PDF) can be reliably determined.

[Abridged] Aims. We study the properties of filamentary structures from the ATLASGAL survey. Methods. We use the DisPerSE algorithm to identify spatially coherent structures located across the inner-Galaxy (300 < l < 60 and |b| < 1.5). Results. We have determined distances, masses and physical sizes for 241 of the filamentary structures. We find a median distance of 3.8 kpc, a mean mass of a few 10^3 m_sun, a mean length of ~6pc and a mass-to-length ratio of (M/L) ~200-2000M_sun/ pc. We also find that these filamentary structures are tightly correlated with the spiral arms in longitude and velocity, and that their semi-major axis is preferentially aligned parallel to the Galactic mid-plane and therefore with the direction of large-scale Galactic magnetic field. We find many examples where the dense filaments identified in ATLASGAL are associated with larger scale filamentary structures (~100 pc), and argue that this is likely to be common, and as such these may indicate a connection between large-scale Galactic dynamics and star formation. Conclusions. We have produced a large and Galaxy-wide catalogue of dense filamentary structures that are representative of a particular size and mass range not previously well studied in the literature. Analyses of the properties and distribution of these filaments reveals that they are correlated with the spiral arms and make a significant contribution to star formation in the Galaxy. Massive star formation is ongoing within ~20% of the filaments and is strongly correlated with the filaments with the largest mass-to- length ratios. The luminosity of the embedded sources has a similar distribution to the Galactic-wide samples of young massive stars and can therefore be considered to be representative.

We have derived the particle injections at the Sun for ten good electron/$^{3}$He-rich solar energetic particle (SEP) events, using a 1.2 AU particle path length (suggested by analysis of the velocity dispersion). The inferred solar injections of high-energy ($\sim$10 to 300 keV) electrons and of $\sim$MeV/nucleon ions (carbon and heavier) start with a delay of 17$\pm$3 minutes and 75$\pm$14 minutes, respectively, after the injection of low-energy ($\sim$0.4 to 9 keV) electrons. The injection duration (averaged over energy) ranges from $\sim$200 to 550 minutes for ions, from $\sim$90 to 160 minutes for low-energy electrons, and from $\sim$10 to 30 minutes for high-energy electrons. Most of the selected events have no reported H$\alpha$ flares or GOES SXR bursts, but all have type III radio bursts that typically start after the onset of a low-energy electron injection. All nine events with SOHO/LASCO coverage have a relatively fast ($>$570km/s), mostly narrow ($\lesssim$30$^{\circ}$), west-limb coronal mass ejection (CME) that launches near the start of the low-energy electron injection, and reaches an average altitude of $\sim$1.0 and 4.7 $R_{S}$, respectively, at the start of the high-energy electron injection and of the ion injection. The electron energy spectra show a continuous power law extending across the transition from low to high energies, suggesting that the low-energy electron injection may provide seed electrons for the delayed high-energy electron acceleration. The delayed ion injections and high ionization states may suggest an ion acceleration along the lower altitude flanks, rather than at the nose of the CMEs.

We fit the $\sim$0.1-500 MeV/nucleon H-Fe spectra in 46 large SEP events surveyed by Desai et al. (2016) with the double power-law Band function to obtain a normalization constant, low- and high-energy parameters $\gamma_a$ and $\gamma_b$; and break energy $E_B$. We also calculate the low-energy power-law spectral slope $\gamma_1$. We find that: 1) $\gamma_a$, $\gamma_1$, and $\gamma_b$ are species-independent within a given SEP event, and the spectra steepen with increasing energy; 2) $E_B$'s are well ordered by Q/M ratio, and decrease systematically with decreasing Q/M, scaling as (Q/M)$^\alpha$ with $\alpha$ varying between $\sim$0.2-3; 3) $\alpha$ is well correlated with Fe/O at $\sim$0.16-0.23 MeV/nucleon and CME speed; 4) In most events: $\alpha<$1.4, the spectra steepen significantly at higher energy with $\gamma_b$-$\gamma_a >$3; and 5) Seven out of 9 extreme SEP events (associated with faster CMEs and GLEs) are Fe-rich, have $\alpha >$1.4, have flatter spectra at low and high energies with $\gamma_b$-$\gamma_a <$3. The species-independence of $\gamma_a$, $\gamma_1$, and $\gamma_b$ and the systematic Q/M dependence of $E_B$ within an event, as well as the range of values for $\alpha$ suggest that the formation of double power-laws in SEP events occurs primarily due to diffusive acceleration at near-Sun CME shocks and not due to scattering in the interplanetary turbulence. In most events, the Q/M-dependence of $E_B$ is consistent with the equal diffusion coefficient condition while the event-to-event variations in $\alpha$ are probably driven by differences in the near-shock wave intensity spectra, which are flatter than the Kolmogorov turbulence spectrum but still weaker compared to that inferred for the extreme events.

Motivated by the report of a possible new planetary member of the Solar System, this work calculates cross sections for interactions between passing stars and this proposed Planet Nine. Evidence for the new planet is provided by the orbital alignment of Kuiper Belt objects, and other Solar System properties, which suggest a Neptune-mass object on an eccentric orbit with semimajor axis $a_9\approx400-1500$ AU. With such a wide orbit, Planet Nine has a large interaction cross section, and is susceptible to disruption by passing stars. Using a large ensemble of numerical simulations (several million), and Monte Carlo sampling, we calculate the cross sections for different classes of orbit-altering events: [A] scattering the planet into its proposed orbit from a smaller orbit, [B] ejecting it from the Solar System from its current orbit, [C] capturing the planet from another system, and [D] capturing a free-floating planet. Results are presented for a range of orbital elements with planetary mass $m_9=10M_{earth}$. Removing Planet Nine from the Solar System is the most likely outcome. Specifically, we obtain ejection cross sections $\sigma_{int}\sim5\times10^6$ AU$^2$ ($5\times10^4$ AU$^2$) for environments corresponding to the birth cluster (field). With these cross sections, Planet Nine is likely to be ejected if the Sun resides within its birth cluster longer than $\Delta{t} \gtrsim 100$ Myr. The probability of ejecting Planet Nine due to passing field stars is $\lesssim 3\%$ over the age of the Sun. Probabilities for producing the inferred Planet Nine orbit are low $(\lesssim 5\%)$.

The shear measurement is a crucial task in the current and the future weak lensing survey projects. And the reconstruction of the point spread function(PSF) is one of the essential steps. In this work, we present three different methods, including Gaussianlets, Moffatlets and EMPCA to quantify their efficiency on PSF reconstruction using four sets of simulated LSST star images. Gaussianlets and Moffatlets are two different sets of basis functions whose profiles are based on Gaussian and Moffat functions respectively. Expectation Maximization(EM) PCA is a statistical method performing iterative procedure to find principal components of an ensemble of star images. Our tests show that: 1) Moffatlets always perform better than Gaussianlets. 2) EMPCA is more compact and flexible, but the noise existing in the Principal Components (PCs) will contaminate the size and ellipticity of PSF while Moffatlets keeps them very well.

It has been recently suggested that large solar energetic particle (SEP) events are often caused by twin CMEs. In the twin-CME scenario, the preceding CME is to provide both an enhanced turbulence level and enhanced seed population at the main CME-driven shock. In this work, we study the effect of the preceding CMEs on the seed population. We examine event-integrated abundance of iron to oxygen ratio (Fe/O) at energies above 25 MeV/nuc for large SEP events in solar cycle 23. We find that the Fe/O ratio (normalized to the reference coronal value of $0.134$) $\leq2.0$ for almost all single-CME events and these events tend to have smaller peak intensities. In comparison, the Fe/O ratio of twin-CME events scatters in a larger range, reaching as high as $8$, suggesting the presence of flare material from perhaps preceding flares. For extremely large SEP events with peak intensity above $1000$ pfu, the Fe/O drop below $2$, indicating that in these extreme events the seed particles are dominated by coronal material than flare material. The Fe/O ratios of Ground level enhancement (GLE) events, all being twin-CME events, scatter in a broad range. For a given Fe/O ratio, GLE events tend to have larger peak intensities than non-GLE events. Using velocity dispersion analysis (VDA), we find that GLE events have lower solar particle release (SPR) heights than non-GLE events, \red{agreeing with earlier results by Reames 2009b.

Recently it has been suggested that the "twin-CME" scenario Li.etal2012 may be a very effective mechanism in causing extreme Solar Energetic Particle (SEP) events and in particular Ground Level Enhancement (GLE) events. Ding.etal2013 performed a statistical examination of the twin-CME scenario with a total of $126$ fast and wide western Coronal Mass Ejections (CMEs). They found that CMEs having a preceding CME with a speed $>$ 300 $km/s$ within $9$ hours from the same active region have larger probability of leading to large SEP events than CMEs that do not have preceding CMEs. The choice of $9$ hours being the time lag $\tau$ between the preceding CME and the main CME was based on some crude estimates of the decay time of the turbulence downstream of the shock driven by the preceding CME. In this work, we examine this choice. For the $126$ fast wide CMEs examined in Ding.etal2013, we vary the time lag $\tau$ from $1$ hour to $24$ hours with an increment of $1$ hour. By considering three quantities whose values depend on the choice of this time lag $\tau$, we show that the choice of $13$ hours for $\tau$ is more appropriate. Our study confirms our earlier result that twin CMEs are more likely to lead to large SEP events than single fast CMEs. The results shown here are of great relevance to space weather studies.

Magnetic clouds (MCs) are the interplanetary counterpart of coronal magnetic flux ropes. They can provide valuable information to reveal the flux rope characteristics at their eruption stage in the corona, which are unable to be explored in situ at present. In this paper, we make a comprehensive survey of the average iron charge state (<Q>Fe) distributions inside 96 MCs for solar cycle 23 using ACE (Advanced Composition Explorer) data. As the <Q>Fe in the solar wind are typically around 9+ to 11+, the Fe charge state is defined as high when the <Q>Fe is larger than 12+, which implies the existence of a considerable amount of Fe ions with high charge states (e.g., \geq 16+). The statistical results show that the <Q>Fe distributions of 92 (~ 96%) MCs can be classified into four groups with different characteristics. In group A (11 MCs), the <Q>Fe shows a bimodal distribution with both peaks higher than 12+. Group B (4 MCs) presents a unimodal distribution of <Q>Fe with its peak higher than 12+. In groups C (29 MCs) and D (48 MCs), the <Q>Fe remains higher and lower than 12+ throughout ACE passage through the MC, respectively. Possible explanations to these distributions are discussed.

Gravity plays important roles in the evolution of molecular clouds. We present an acceleration mapping method to estimate the acceleration induced by gravitational interactions in molecular clouds based on observational data. We find that the geometry of a region has a significant impact on the behavior of gravity. In the Pipe nebula which can be approximated as a gas filament, we find that gravitational acceleration can effectively compress the end of this filament, which may have triggered star formation. We identify this as the "gravitational focusing" effect proposed by Burkert & Hartman (2004). In the sheet-like IC348-B3 region, gravity can lead to collapse at its edge, while in the centrally condensed NGC1333 cluster-forming region gravity can drive accretion towards the center. In general, gravitational acceleration tends to be enhanced in the localized regions around the ends of the filaments and the edges of sheet-like structures. Neglecting magnetic fields, these "gravitational focusing" and "edge collapse" effects can promote the formation of dense gas in a timescale that is much shorter than the global dynamical time. Since the interstellar medium is in general structured, these edge effects should be prevalent.

Since the release of LAMOST (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope) catalog, we have the opportunity to use the LAMOST DR2 stellar catalog and \emph{WISE All-sky catalog} to search for 22 $\mu$m excess candidates. In this paper, we present 10 FGK candidates which show an excess in the infrared (IR) at 22 $\mu$m. The ten sources are all the newly identified 22 $\mu$m excess candidates. Of these 10 stars, 5 stars are F type and 5 stars are G type. The criterion for selecting candidates is $K_s-[22]_{\mu m}\geq0.387$. In addition, we present the spectral energy distributions (SEDs) covering wavelength from optical to mid-infrared band. Most of them show an obvious excess from 12 $\mu$m band and three candidates even show excess from 3.4 $\mu$m. To characterize the amount of dust, we also estimate the fractional luminosity of ten 22 $\mu$m excess candidates.

Corona structures and processes during the pre-impulsive stage of solar eruption are crucial to understanding the physics leading to the subsequent explosive energy release. Here we present the first microwave imaging study of a hot flux rope structure during the pre-impulsive stage of an eruptive M7.7 solar flare, with the Nobeyama Radioheliograph (NoRH) at 17 GHz. The flux rope is also observed by the SDO/AIA in its hot passbands of 94 and 131 \AA\. In the microwave data, it is revealed as an overall arcade-like structure consisting of several intensity enhancements bridged by generally weak emissions, with brightness temperatures ($T_B$) varying from ~10,000~K to ~20,000 K. Locations of microwave intensity enhancements along the structure remain relatively fixed at certain specific parts of the flux rope, indicating that the distribution of emitting electrons is affected by the large scale magnetic configuration of the twisted flux rope. Wavelet analysis shows a pronounced 2-min period of the microwave $T_B$ variation during the pre-impulsive stage of interest. The period agrees well with that reported for AIA sunward-contracting loops and upward ejective plasmoids (suggested to be reconnection outflows). This suggests that both periodicities are controlled by the same reconnection process that takes place intermittently at a 2-min time scale. We infer that at least a part of the emission is excited by non-thermal energetic electrons via the gyro-synchrotron mechanism. The study demonstrates the potential of microwave imaging in exploring the flux rope magnetic geometry and relevant reconnection process during the onset of solar eruption.

With a test-particle simulation, we investigate the effect of large-scale coronal magnetic fields on electron acceleration at an outward-propagating coronal shock with a circular front. The coronal field is approximated by an analytical solution with a streamer-like magnetic field featured by partially open magnetic field and a current sheet at the equator atop the closed region. We show that the large-scale shock-field configuration, especially the relative curvature of the shock and the magnetic field line across which the shock is sweeping, plays an important role in the efficiency of electron acceleration. At low shock altitudes, when the shock curvature is larger than that of magnetic field lines, the electrons are mainly accelerated at the shock flanks; at higher altitudes, when the shock curvature is smaller, the electrons are mainly accelerated at the shock nose around the top of closed field lines. The above process reveals the shift of efficient electron acceleration region along the shock front during its propagation. It is also found that in general the electron acceleration at the shock flank is not so efficient as that at the top of closed field since at the top a collapsing magnetic trap can be formed. In addition, we find that the energy spectra of electrons is power-law like, first hardening then softening with the spectral index varying in a range of -3 to -6. Physical interpretations of the results and implications on the study of solar radio bursts are discussed.

OB type stars have strong ionizing radiation, and drive energetic winds. The ultraviolet (UV) radiation from ionizing stars may heat dust and ionize gas to sweep up an expanding bubble shell. This shell may be the result of feedback leading to a new generation of stars. N131 is an infrared dust bubble residing in a molecular filament. We study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 $\mu$m, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and $^{13}$CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 $\mu$m and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of HII region, strong stellar winds, and fragments under self-gravity.

The obliquities of planet-hosting stars are clues about the formation of planetary systems. Previous observations led to the hypothesis that for close-in giant planets, spin-orbit alignment is enforced by tidal interactions. Here, we examine two problems with this hypothesis. First, Mazeh and coworkers recently used a new technique -- based on the amplitude of starspot-induced photometric variability -- to conclude that spin-orbit alignment is common even for relatively long-period planets, which would not be expected if tides were responsible. We re-examine the data and find a statistically significant correlation between photometric variability and planetary orbital period that is qualitatively consistent with tidal interactions. However it is still difficult to explain quantitatively, as it would require tides to be effective for periods as long as tens of days. Second, Rogers and Lin argued against a particular theory for tidal re-alignment by showing that initially retrograde systems would fail to be re-aligned, in contradiction with the observed prevalence of prograde systems. We investigate a simple model that overcomes this problem by taking into account the dissipation of inertial waves and the equilibrium tide, as well as magnetic braking. We identify a region of parameter space where re-alignment can be achieved, but it only works for close-in giant planets, and requires some fine tuning. Thus, while we find both problems to be more nuanced than they first appeared, the tidal model still has serious shortcomings.

Inspired by the close-proximity pair of planets in the Kepler-36 system, we consider two effects that may have important ramifications for the development of life in similar systems where a pair of planets may reside entirely in the habitable zone of the hosting star. Specifically, we run numerical simulations to determine whether strong, resonant (or non-resonant) planet-planet interactions can cause large variations in planet obliquity---thereby inducing large variations in climate. We also determine whether or not resonant interactions affect the rate of lithopanspermia between the planet pair---which could facilitate the growth and maintenance of life on both planets. We find that first-order resonances do not cause larger obliquity variations compared with non-resonant cases. We also find that resonant interactions are not a primary consideration in lithopanspermia. Lithopanspermia is enhanced significantly as the planet orbits come closer together---reaching nearly the same rate as ejected material falling back to the surface of the originating planet (assuming that the ejected material makes it out to the location of our initial conditions). Thus, in both cases our results indicate that close-proximity planet pairs in multihabitable systems are conducive to life in the system.

Dome A in the Antarctic plateau is likely one of the best astronomical observing sites on Earth. The first one of three Antarctic Survey Telescope (AST3-1), a 50/68 cm Schmidt-like equatorial-mount telescope, is the first trackable telescope of China operating in Antarctica and the biggest telescope located in Antarctic inland. AST3-1 obtained huge amounts of data in 2012 and we processed the time-series parts. Here we present light curves of 29 variable stars identified from ten-day observations in 2012 with AST3-1, including 22 newly discovered variable stars. 23 of them are eclipsing binaries and the others are pulsating stars. We present the properties of the 29 variable stars, including the classifications, periods and magnitude ranges in i band. For the 17 eclipsing binaries, the phased light curves are presented with the orbital period values well determined.

Recently the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) published its first data release (DR1), which is ranked as the largest stellar spectra dataset in the world so far. We combine the PASTEL catalog and SIMBAD radial velocities as a testing standard to validate the DR1 stellar parameters (effective temperature $T_{\mathrm{eff}}$, surface gravity $\log g$, metallicity $\mathrm{[Fe/H]}$ and radial velocity $V_{\mathrm{r}}$). Through cross-identification of the DR1 catalogs and the PASTEL catalog, we obtain a preliminary sample of 422 stars. After removal of stellar parameter measurements from problematic spectra and applying effective temperature constraints to the sample, we compare the stellar parameters from DR1 with those from PASTEL and SIMBAD to prove that the DR1 results are reliable in restricted $T_{\mathrm{eff}}$ ranges. We derive standard deviations of 110 K, 0.19 dex, 0.11 dex and 4.91 $\mathrm{km\,s^{-1}}$ , for $T_{\mathrm{eff}}$, $\log g$, $[\mathrm{Fe/H}]$ when $T_{\mathrm{eff}}<8000\,\mathrm{K}$, and for $V_{\mathrm{r}}$ when $T_{\mathrm{eff}}<10000\,\mathrm{K}$, respectively. Systematic errors are negligible except for that of $V_{\mathrm{r}}$. Besides, metallicities in DR1 are systematically higher than those in PASTEL, in the range of PASTEL $[\mathrm{Fe/H}]<-1.5$.

By means of fully kinetic simulations, we investigate electron acceleration during magnetic reconnection in a nonrelativistic proton--electron plasma with conditions similar to solar corona and flares. We demonstrate that reconnection leads to a nonthermally dominated electron acceleration with a power-law energy distribution in the nonrelativistic low-$\beta$ regime but not in the high-$\beta$ regime, where $\beta$ is the ratio of the plasma thermal pressure and the magnetic pressure. The accelerated electrons contain most of the dissipated magnetic energy in the low-$\beta$ regime. A guiding-center current description is used to reveal the role of electron drift motions during the bulk nonthermal energization. We find that the main acceleration mechanism is a \textit{Fermi}-type acceleration accomplished by the particle curvature drift motion along the electric field induced by the reconnection outflows. Although the acceleration mechanism is similar for different plasma $\beta$, low-$\beta$ reconnection drives fast acceleration on Alfv\'enic timescales and develops power laws out of thermal distribution. The nonthermally dominated acceleration resulting from magnetic reconnection in low-$\beta$ plasma may have strong implications for the highly efficient electron acceleration in solar flares and other astrophysical systems.

Band split of solar type II radio bursts, discovered several decades ago, is a fascinating phenomenon with the type-II lanes exhibiting two almost-parallel sub-bands with similar morphology. The underlying split mechanism remains elusive. One popular interpretation is that the splitting bands are emitted from the shock upstream and downstream, respectively, with their frequency ratio ({\gamma}) determined by the shock compression ratio. This interpretation has been taken as the physical basis for many published references. Here we report an observational analysis of type II events with nice split selected from the ground-based RSTN data from 2001 to 2014, in the metric-decametric wavelength. We investigate the temporal variation and distribution of {\gamma}, and conduct correlation analyses on the deduced spectral values. It is found that {\gamma} varies in a very narrow range with >80% of {\gamma} (one-minute averaged data) being between 1.15 to 1.25. For some well-observed and long-lasting events, {\gamma} does not show a systematic variation trend within observational uncertainties, from the onset to the termination of the splits. In addition, the parameters representing the propagation speed of the radio source (presumably the coronal shock) show a very weak or basically no correlation with {\gamma}. We suggest that these results do not favor the upstreamdownstream scenario of band splits.

We present CCD $UBVRI$ photometry of the field of the open cluster NGC 6866. Structural parameters of the cluster are determined utilizing the stellar density profile of the stars in the field. We calculate the probabilities of the stars being a physical member of the cluster using their astrometric data and perform further analyses using only the most probable members. The reddening and metallicity of the cluster were determined by independent methods. The LAMOST spectra and the ultraviolet excess of the F and G type main-sequence stars in the cluster indicate that the metallicity of the cluster is about the solar value. We estimated the reddening $E(B-V)=0.074 \pm 0.050$ mag using the $U-B$ vs $B-V$ two-colour diagram. The distance modula, the distance and the age of NGC 6866 were derived as $\mu = 10.60 \pm 0.10$ mag, $d=1189 \pm 75$ pc and $t = 813 \pm 50$ Myr, respectively, by fitting colour-magnitude diagrams of the cluster with the PARSEC isochrones. The Galactic orbit of NGC 6866 indicates that the cluster is orbiting in a slightly eccentric orbit with $e=0.12$. The mass function slope $x=1.35 \pm 0.08$ was derived by using the most probable members of the cluster.

We use about 200,000 FGK type main-sequence stars from the LAMOST DR1 data to map the local stellar kinematics. With the velocity de-projection technique, we are able to derive the averaged 3 dimensional velocity and velocity ellipsoids using only the line-of-sight velocity for the stars with various effective temperatures within $100 < |z| < 500$ pc. Using the mean velocities of the cool stars, we derive the solar motion of ($U_{\!\odot}$, $V_{\!\odot}$, $W_{\!\odot}$)=(9.58$\pm2.39$, 10.52$\pm1.96$, 7.01$\pm1.67$)$km\,s^{-1}$ with respect to the local standard of rest. Moreover, we find that the stars with ${T}_{\rm eff}>6000$K show a net asymmetric motion of $\sim3 km\,s^{-1}$ in $\langle W\rangle$ compared to the stars with ${T}_{\rm eff}<6000$K. And their azimuthal velocity increases when $|z|$ increases. This peculiar motion in the warmer stars is likely because they are young and not completely relaxed, although other reasons, such as the resonance induced by the central rotating bar or the spiral structures, and the perturbation of the merging dwarf galaxies, can not be ruled out. The derived velocity dispersions and cross terms for the data are approximately consistent with previous studies. We also find that the vertical gradients of $\sigma_{U}$ and $\sigma_V$ are larger than that of $\sigma_W$ . And the vertical gradient of $\sigma_U$ shows clear correlation with ${T}_{\rm eff}$, while the other two do not. Finally, our sample shows vertex deviation of about 11$^\circ$, at $300 < |z| < 500$pc, but roughly zero at $100 < |z| < 300$pc.

Hot channels (HCs), high temperature erupting structures in the lower corona of the Sun, have been proposed as a proxy of magnetic flux ropes (MFRs) since their initial discovery. However, it is difficult to make definitive proof given the fact that there is no direct measurement of magnetic field in the corona. An alternative way is to use the magnetic field measurement in the solar wind from in-situ instruments. On 2012 July 12, an HC was observed prior to and during a coronal mass ejection (CME) by the AIA high-temperature images. The HC is invisible in the EUVI low-temperature images, which only show the cooler leading front (LF). However, both the LF and an ejecta can be observed in the coronagraphic images. These are consistent with the high temperature and high density of the HC and support that the ejecta is the erupted HC. In the meanwhile, the associated CME shock was identified ahead of the ejecta and the sheath through the COR2 images, and the corresponding ICME was detected by \textit{ACE}, showing the shock, sheath and magnetic cloud (MC) sequentially, which agrees with the coronagraphic observations. Further, the MC contained a low-ionization-state center and a high-ionization-state shell, consistent with the pre-existing HC observation and its growth through magnetic reconnection. All of these observations support that the MC detected near the Earth is the counterpart of the erupted HC in the corona for this event. Therefore, our study provides strong observational evidence of the HC as an MFR.

Bolton and Schlegel presented a promising deconvolution method to extract 1D spectra from a 2D optical fiber spectral CCD image. The method could eliminate the PSF difference between fibers, extract spectra to the photo noise level, as well as improve the resolution. But the method is limited by its huge computation requirement and thus cannot be implemented in actual data reduction. In this article, we develop a practical computation method to solve the computation problem. The new computation method can deconvolve a 2D fiber spectral image of any size with actual PSFs, which may vary with positions. Our method does not require large amounts of memory and can extract a 4k multi 4k noise-free CCD image with 250 fibers in 2 hr. To make our method more practical, we further consider the influence of noise, which is thought to be an intrinsic illposed problem in deconvolution algorithms. We modify our method with a Tikhonov regularization item to depress the method induced noise. Compared with the results of traditional extraction methods, our method has the least residual and influence by cross talk and noise, even for extreme situation. Our method can convergent in 2 to 4 iterations, and the computation times are about 3.5 hr for the extreme fiber distance and about 2 hr for nonextreme cases. Finally, we apply our method to real LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope, a.k.a. Guo Shou Jing Telescope) data. We find that the 1D spectra extracted by our method have both higher signal-to-noise ratio and resolution than the traditional methods.

We perform a discrimination procedure with the spectral index diagram of TiO5 and CaH2+CaH3 to separate M giants from M dwarfs. Using the M giant spectra identified from the LAMOST DR1 with high signal-to-noise ratio (SNR), we have successfully assembled a set of M giant templates, which show more reliable spectral features. Combining with the M dwarf/subdwarf templates in Zhong et al. (2015), we present an extended M-type templates library which includes not only M dwarfs with well-defined temperature and metallicity grid but also M giants with subtype from M0 to M6. Then, the template-fit algorithm were used to automatically identify and classify M giant stars from the LAMOST DR1. The result of M giant stars catalog is cross-matched with 2MASS JHKs and WISE W1/W2 infrared photometry. In addition, we calculated the heliocentric radial velocity of all M giant stars by using the cross-correlation method with the template spectrum in a zero-velocity restframe. Using the relationship between the absolute infrared magnitude MJ and our classified spectroscopic subtype, we derived the spectroscopic distance of M giants with uncertainties of about 40%. A catalog of 8639 M giants is provided. As an additional search result, we also present 101690 M dwarfs/subdwarfs catalog which were classified by our classification pipeline.

An enhanced rate of stellar tidal disruption events (TDEs) may be an important characteristic of supermassive black hole (SMBH) binaries at close separations. Here we study the evolution of the distribution of stars around a SMBH binary due to the eccentric Kozai-Lidov (EKL) mechanism, including octupole effects and apsidal precession caused by the stellar mass distribution and general relativity. We identify a region around one of the SMBHs in the binary where the EKL mechanism drives stars to high eccentricities, which ultimately causes the stars to either scatter off the second SMBH or get disrupted. For SMBH masses 10^7 Msun and 10^8 Msun, the TDE rate can reach 10^{-2} yr and deplete a region of the stellar cusp around the secondary SMBH in ~0.5 Myr. As a result, the final geometry of the stellar distribution between 0.01 and 0.1 pc around the secondary SMBH is a torus. These effects may be even more prominent in nuclear stellar clusters hosting a supermassive and an intermediate mass black hole.

We present a method to estimate distances to stars with spectroscopically derived stellar parameters. The technique is a Bayesian approach with likelihood estimated via comparison of measured parameters to a grid of stellar isochrones, and returns a posterior probability density function for each star's absolute magnitude. This technique is tailored specifically to data from the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) survey. Because LAMOST obtains roughly 3000 stellar spectra simultaneously within each ~5-degree diameter "plate" that is observed, we can use the stellar parameters of the observed stars to account for the stellar luminosity function and target selection effects. This removes biasing assumptions about the underlying populations, both due to predictions of the luminosity function from stellar evolution modeling, and from Galactic models of stellar populations along each line of sight. Using calibration data of stars with known distances and stellar parameters, we show that our method recovers distances for most stars within ~20%, but with some systematic overestimation of distances to halo giants. We apply our code to the LAMOST database, and show that the current precision of LAMOST stellar parameters permits measurements of distances with ~40% error bars. This precision should improve as the LAMOST data pipelines continue to be refined.

Asteroseismology is one of the most accurate approaches to estimate the surface gravity of a star. However, most of the data from the current spectroscopic surveys do not have asteroseismic measurements, which is very expensive and time consuming. In order to improve the spectroscopic surface gravity estimates for a large amount of survey data with the help of the small subset of the data with seismic measurements, we set up a support vector regression model for the estimation of the surface gravity supervised by 1,374 LAMOST giant stars with Kepler seismic surface gravity. The new approach can reduce the uncertainty of the estimates down to about 0.1 dex, which is better than the LAMOST pipeline by at least a factor of 2, for the spectra with signal-to-noise ratio higher than 20. Compared with the logg estimated from the LAMOST pipeline, the revised logg values provide a significantly improved match to the expected distribution of red clump and RGB stars from stellar isochrones. Moreover, even the red bump stars, which extend to only about 0.1 dex in logg, can be discriminated from the new estimated surface gravity. The method is then applied to about 350,000 LAMOST metal-rich giant stars to provide improved surface gravity estimates. In general, the uncertainty of the distance estimate based on the SVR surface gravity can be reduced to about 12% for the LAMOST data.

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

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

Filament eruptions often lead to coronal mass ejections (CMEs), which can affect critical technological systems in space and on the ground when they interact with the geo-magnetosphere in high speeds. Therefore, it is an important issue to investigate the acceleration mechanisms of CMEs in solar/space physics. Based on observations and simulations, the resistive magnetic reconnection and the ideal instability of magnetic flux rope have been proposed to accelerate CMEs. However, it remains elusive whether both of them play a comparable role during a particular eruption. It has been extremely difficult to separate their contributions as they often work in a close time sequence during one fast acceleration phase. Here we report an intriguing filament eruption event, which shows two apparently separated fast acceleration phases and provides us an excellent opportunity to address the issue. Through analyzing the correlations between velocity (acceleration) and soft (hard) X-ray profiles, we suggest that the instability and magnetic reconnection make a major contribution during the first and second fast acceleration phases, respectively. Further, we find that both processes have a comparable contribution to accelerate the filament in this event.

A new line list including positions and absolute intensities (in the form of Einstein $A$ values and oscillator strengths) has been produced for the OH ground X\DP\ state rovibrational (Meinel system) and pure rotational transitions. All possible transitions are included with v$\primed$ and v$\Dprimed$ up to 13, and $J$ up to between 9.5 and 59.5, depending on the band. An updated fit to determine molecular constants has been performed, which includes some new rotational data and a simultaneous fitting of all molecular constants. The absolute line intensities are based on a new dipole moment function, which is a combination of two high level ab initio calculations. The calculations show good agreement with an experimental v=1 lifetime, experimental $\mu_\mathrm{v}$ values, and $\Delta$v=2 line intensity ratios from an observed spectrum. To achieve this good agreement, an alteration in the method of converting matrix elements from Hund's case (b) to (a) was made. Partitions sums have been calculated using the new energy levels, for the temperature range 5-6000 K, which extends the previously available (in HITRAN) 70-3000 K range. The resulting absolute intensities have been used to calculate O abundances in the Sun, Arcturus, and two red giants in the Galactic open and globular clusters M67 and M71. Literature data based mainly on [O I] lines are available for the Sun and Arcturus, and excellent agreement is found.

We present the G-virial method (available at http://gxli.github.io/G-virial/) which aims to quantify (1) the importance of gravity in molecular clouds in the position-position-velocity (PPV) space, and (2) properties of the gas condensations in molecular clouds. Different from previous approaches that calculate the virial parameter for different regions, our new method takes gravitational interactions between all the voxels in 3D PPV data cubes into account, and generates maps of the importance of gravity. This map can be combined with the original data cube to derive relations such as the mass-radius relation. Our method is important for several reasons. First, it offers the the ability to quantify the centrally condensed structures in the 3D PPV data cubes, and enables us to compare them in an uniform framework. Second, it allows us to understand the importance of gravity at different locations in the data cube, and provides a global picture of gravity in clouds. Third, it offers a robust approach to decomposing the data into different regions which are gravitationally coherent. To demonstrate the application of our method we identified regions from the Perseus and Ophiuchus molecular clouds, and analyzed their properties. We found an increase in the importance of gravity towards the centers of the individual molecular condensations. We also quantified the properties of the regions in terms of mass-radius and mass-velocity relations. Through evaluating the virial parameters based on the G-virial, we found that all our regions are almost gravitationally bound. Cluster-forming regions appear are more centrally condensed.

Modern large-scale surveys have allowed the identification of large numbers of white dwarfs. However, these surveys are subject to complicated target selection algorithms, which make it almost impossible to quantify to what extent the observational biases affect the observed populations. The LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope) Spectroscopic Survey of the Galactic anti-center (LSS-GAC) follows a well-defined set of criteria for selecting targets for observations. This advantage over previous surveys has been fully exploited here to identify a small yet well-characterised magnitude-limited sample of hydrogen-rich (DA) white dwarfs. We derive preliminary LSS-GAC DA white dwarf luminosity and mass functions. The space density and average formation rate of DA white dwarfs we derive are 0.83+/-0.16 x 10^{-3} pc^{-3} and 5.42 +/- 0.08 x 10^{-13} pc^{-3} yr^{-1}, respectively. Additionally, using an existing Monte Carlo population synthesis code we simulate the population of single DA white dwarfs in the Galactic anti-center, under various assumptions. The synthetic populations are passed through the LSS-GAC selection criteria, taking into account all possible observational biases. This allows us to perform a meaningful comparison of the observed and simulated distributions. We find that the LSS-GAC set of criteria is highly efficient in selecting white dwarfs for spectroscopic observations (80-85 per cent) and that, overall, our simulations reproduce well the observed luminosity function. However, they fail at reproducing an excess of massive white dwarfs present in the observed mass function. A plausible explanation for this is that a sizable fraction of massive white dwarfs in the Galaxy are the product of white dwarf-white dwarf mergers.

We review five often used quad lens models, each of which has analytical solutions and can produce four images at most. Each lens model has two parameters, including one that describes the intensity of non-dimensional mass density, and the other one that describes the deviation from the circular lens. In our recent work, we have found that the cusp and the fold summations are not equal to 0, when a point source infinitely approaches a cusp or a fold from inner side of the caustic. Based on the magnification invariant theory, which states that the sum of signed magnifications of the total images of a given source is a constant, we calculate the cusp summations for the five lens models. We find that the cusp summations are always larger than 0 for source on the major cusps, while can be larger or smaller than 0 for source on the minor cusps. We also find that if these lenses tend to the circular lens, the major and minor cusp summations will have infinite values, and with positive and negative signs respectively. The cusp summations do not change significantly if the sources are slightly deviated from the cusps. In addition, through the magnification invariants, we also derive the analytical signed cusp relations on the axes for three lens models. We find that both on the major and the minor axes the larger the lenses deviated from the circular lens, the larger the signed cusp relations. The major cusp relations are usually larger than the absolute minor cusp relations, but for some lens models with very large deviation from circular lens, the minor cusp relations can be larger than the major cusp relations.

Hot channel (HC) is a high temperature ($\sim$10 MK) structure in the inner corona revealed first by Atmospheric Imaging Assembly (AIA) on board \textit{Solar Dynamics Observatory}. Eruption of HC is often associated with flare and coronal mass ejection. Previous studies suggest that HC is a good proxy of magnetic flux rope (MFR) in the inner corona, in addition to another well-known MFR candidate, the prominence-cavity structure that is with a normal coronal temperature ($\sim$1-2 MK). In this paper, we report a high temperature structure (HTS, $\sim$1.5 MK) contained in an interplanetary coronal mass ejection induced by an HC eruption. According to the observations of bidirectional electrons, high temperature and density, strong magnetic field, and its association with the shock, sheath, and plasma pile-up region, we suggest that the HTS is the interplanetary counterpart of the HC. The scale of the measured HTS is around 14 R$_\odot$, and it maintained a much higher temperature than the background solar wind even at 1 AU. It is significantly different from the typical magnetic clouds (MCs), which usually have a much lower temperature. Our study suggests that the existence of a corotating interaction region ahead of the HC formed a magnetic container to inhibit the HC expansion and cooling down to a low temperature.

Most planetary systems are formed within stellar clusters, and these environments can shape their properties. This paper considers scattering encounters between solar systems and passing cluster members, and calculates the corresponding interaction cross sections. The target solar systems are generally assumed to have four giant planets, with a variety of starting states, including circular orbits with the semimajor axes of our planets, a more compact configuration, an ultra-compact state with multiple mean motion resonances, and systems with massive planets. We then consider the effects of varying the cluster velocity dispersion, the relative importance of binaries versus single stars, different stellar host masses, and finite starting eccentricities of the planetary orbits. For each state of the initial system, we perform an ensemble of numerical scattering experiments and determine the cross sections for eccentricity increase, inclination angle increase, planet ejection, and capture. This paper reports results from over 2 million individual scattering simulations. Using supporting analytic considerations, and fitting functions to the numerical results, we find a universal formula that gives the cross sections as a function of stellar host mass, cluster velocity dispersion, starting planetary orbital radius, and final eccentricity. The resulting cross sections can be used in a wide variety of applications. As one example, we revisit constraints on the birth aggregate of our Solar System due to dynamical scattering and find $N<10^4$ (consistent with previous estimates).

We use observed optical to near infrared spectral energy distributions (SEDs) of 266 galaxies in the COSMOS survey to derive the wavelength dependence of the dust attenuation at high redshift. All of the galaxies have spectroscopic redshifts in the range z = 2 to 6.5. The presence of the CIV absorption feature, indicating that the rest-frame UV-optical SED is dominated by OB stars, is used to select objects for which the intrinsic, unattenuated spectrum has a well-established shape. Comparison of this intrinsic spectrum with the observed broadband photometric SED then permits derivation of the wavelength dependence of the dust attenuation. The derived dust attenuation curve is similar in overall shape to the Calzetti curve for local starburst galaxies. We also see the 2175 \AA~bump feature which is present in the Milky Way and LMC extinction curves but not seen in the Calzetti curve. The bump feature is commonly attributed to graphite or PAHs. No significant dependence is seen with redshift between sub-samples at z = 2 - 4 and z = 4 - 6.5. The 'extinction' curve obtained here provides a firm basis for color and extinction corrections of high redshift galaxy photometry.

As a major component of the LAMOST Galactic surveys, the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC) aims to survey a significant volume of the Galactic thin/thick discs and halo for a contiguous sky area of over 3,400 deg$^2$ centred on the Galactic anti-centre ($|b| \leq 30^{\circ}$, $150 \leq l \leq 210^{\circ}$), and obtain $\lambda\lambda$3700 -- 9000 low resolution ($R \sim 1,800$) spectra for a statistically complete sample of $\sim 3$\,M stars of all colours down to a limiting magnitued of $r$ $\sim$ 17.8\,mag (to 18.5\,mag for limited fields). Together with Gaia, the LSS-GAC will yield a unique dataset to advance our understanding of the structure and assemblage history of the Galaxy, in particular its disk(s). In addition to the main survey, the LSS-GAC will also target hundreds of thousands objects in the vicinity fields of M\,31 and M\,33 and survey a significant fraction (over a million) of randomly selected very bright stars (VB; $r \le 14$ mag) in the northern hemisphere. During the Pilot and the first year Regular Surveys of LAMOST, a total of 1,042,586 [750,867] spectra of a signal to noise ratio S/N(7450\AA) $\ge$ 10 [S/N(4650\AA) $\ge$ 10] have been collected. In this paper, we present a detailed description of the target selection algorithm, survey design, observations and the first data release of value-added catalogues (including radial velocities, effective temperatures, surface gravities, metallicities, values of interstellar extinction, distances, proper motions and orbital parameters) of the LSS-GAC.

As a major component of the LAMOST Galactic surveys, the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC) aims to survey a significant volume of the Galactic thin/thick discs and halo for a contiguous sky area of over 3,400 deg$^2$ centred on the Galactic anti-centre ($|b| \leq 30^{\circ}$, $150 \leq l \leq 210^{\circ}$), and obtain $\lambda\lambda$3700 -- 9000 low resolution ($R \sim 1,800$) spectra for a statistically complete sample of $\sim 3$\,M stars of all colours down to a limiting magnitued of $r$ $\sim$ 17.8\,mag (to 18.5\,mag for limited fields). Together with Gaia, the LSS-GAC will yield a unique dataset to advance our understanding of the structure and assemblage history of the Galaxy, in particular its disk(s). In addition to the main survey, the LSS-GAC will also target hundreds of thousands objects in the vicinity fields of M\,31 and M\,33 and survey a significant fraction (over a million) of randomly selected very bright stars (VB; $r \le 14$ mag) in the northern hemisphere. During the Pilot and the first year Regular Surveys of LAMOST, a total of 1,042,586 [750,867] spectra of a signal to noise ratio S/N(7450\AA) $\ge$ 10 [S/N(4650\AA) $\ge$ 10] have been collected. In this paper, we present a detailed description of the target selection algorithm, survey design, observations and the first data release of value-added catalogues (including radial velocities, effective temperatures, surface gravities, metallicities, values of interstellar extinction, distances, proper motions and orbital parameters) of the LSS-GAC.

We introduce the LAMOST Stellar Parameter Pipeline at Peking University --- LSP3, developed and implemented for the determinations of radial velocity $V_{\rm r}$ and stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$, metallicity [Fe/H]) for the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC). We describe the algorithms of LSP3 and examine the accuracy of parameters yielded by it. The precision and accuracy of parameters yielded are investigated by comparing results of multi-epoch observations and of candidate members of open and globular clusters, with photometric calibration, as well as with independent determinations available from a number of external databases, including the PASTEL archive, the APOGEE, SDSS and RAVE surveys, as well as those released in the LAMOST DR1. The uncertainties of LSP3 parameters are characterized and quantified as a function of the spectral signal-to-noise ratio (SNR) and stellar atmospheric parameters. We conclude that the current implementation of LSP3 has achieved an accuracy of 5.0\,km\,s$^{-1}$, 150\,K, 0.25\,dex, 0.15\,dex for the radial velocity, effective temperature, surface gravity and metallicity, respectively, for LSS-GAC spectra of FGK stars of SNRs per pixel higher than 10. The LSP3 has been applied to over a million LSS-GAC spectra collected hitherto. Stellar parameters yielded by the LSP3 will be released to the general public following the data policy of LAMOST, together with estimates of the interstellar extinction $E(B-V)$ and stellar distances, deduced by combining spectroscopic and multi-band photometric measurements using a variety of techniques.

We have developed and implemented an iterative algorithm of flux calibration for the LAMOST Spectroscopic Survey of the Galactic anti-center (LSS-GAC). For a given LSS-GAC plate, the spectra are first processed with a set of nominal spectral response curves (SRCs) and used to derive initial stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$ and metallicity [Fe/H]) as well as dust reddening $E(B-V)$ of all targeted stars. For each of the sixteen spectrographs, several F-type stars of good signal-to-noise ratios (SNRs) are then selected as flux standard stars for further, iterative spectral flux calibration. Comparison of spectrophotometric colours, deduced from the flux-calibrated spectra, with the photometric measurements yields average differences of 0.02$\pm$0.07 and $-$0.04$\pm$0.09\,mag for the $(g-r)$ and $(g-i)$, respectively. The relatively large negative offset in $(g-i)$ is due to the fact that we have opted not to correct for the telluric bands, most notably the atmospheric A-band in the wavelength range of $i$-band. Comparison of LSS-GAC multi-epoch observations of duplicate targets indicates that the algorithm has achieved an accuracy of about 10 per cent in relative flux calibration for the wavelength range 4000 -- 9000\,\AA. The shapes of SRC deduced for the individual LAMOST spectrographs are found to vary by up to 30 per cent for a given night, and larger for different nights, indicating that the derivation of SRCs for the individual plates is essential in order to achieve accurate flux calibration for the LAMOST spectra.

We introduce the LAMOST Stellar Parameter Pipeline at Peking University --- LSP3, developed and implemented for the determinations of radial velocity $V_{\rm r}$ and stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$, metallicity [Fe/H]) for the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC). We describe the algorithms of LSP3 and examine the accuracy of parameters yielded by it. The precision and accuracy of parameters yielded are investigated by comparing results of multi-epoch observations and of candidate members of open and globular clusters, with photometric calibration, as well as with independent determinations available from a number of external databases, including the PASTEL archive, the APOGEE, SDSS and RAVE surveys, as well as those released in the LAMOST DR1. The uncertainties of LSP3 parameters are characterized and quantified as a function of the spectral signal-to-noise ratio (SNR) and stellar atmospheric parameters. We conclude that the current implementation of LSP3 has achieved an accuracy of 5.0\,km\,s$^{-1}$, 150\,K, 0.25\,dex, 0.15\,dex for the radial velocity, effective temperature, surface gravity and metallicity, respectively, for LSS-GAC spectra of FGK stars of SNRs per pixel higher than 10. The LSP3 has been applied to over a million LSS-GAC spectra collected hitherto. Stellar parameters yielded by the LSP3 will be released to the general public following the data policy of LAMOST, together with estimates of the interstellar extinction $E(B-V)$ and stellar distances, deduced by combining spectroscopic and multi-band photometric measurements using a variety of techniques.

We have developed and implemented an iterative algorithm of flux calibration for the LAMOST Spectroscopic Survey of the Galactic anti-center (LSS-GAC). For a given LSS-GAC plate, the spectra are first processed with a set of nominal spectral response curves (SRCs) and used to derive initial stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$ and metallicity [Fe/H]) as well as dust reddening $E(B-V)$ of all targeted stars. For each of the sixteen spectrographs, several F-type stars of good signal-to-noise ratios (SNRs) are then selected as flux standard stars for further, iterative spectral flux calibration. Comparison of spectrophotometric colours, deduced from the flux-calibrated spectra, with the photometric measurements yields average differences of 0.02$\pm$0.07 and $-$0.04$\pm$0.09\,mag for the $(g-r)$ and $(g-i)$, respectively. The relatively large negative offset in $(g-i)$ is due to the fact that we have opted not to correct for the telluric bands, most notably the atmospheric A-band in the wavelength range of $i$-band. Comparison of LSS-GAC multi-epoch observations of duplicate targets indicates that the algorithm has achieved an accuracy of about 10 per cent in relative flux calibration for the wavelength range 4000 -- 9000\,\AA. The shapes of SRC deduced for the individual LAMOST spectrographs are found to vary by up to 30 per cent for a given night, and larger for different nights, indicating that the derivation of SRCs for the individual plates is essential in order to achieve accurate flux calibration for the LAMOST spectra.

Magnetic reconnection is believed to be the dominant energy release mechanism in solar flares. The standard flare model predicts both downward and upward outflow plasmas with speeds close to the coronal Alfv\'{e}n speed. Yet, spectroscopic observations of such outflows, especially the downflows, are extremely rare. With observations of the newly launched Interface Region Imaging Spectrograph (IRIS), we report the detection of greatly redshifted ($\sim$125 km s$^{-1}$ along line of sight) Fe {\sc{xxi}} 1354.08\AA{} emission line with a $\sim$100 km s$^{-1}$ nonthermal width at the reconnection site of a flare. The redshifted Fe {\sc{xxi}} feature coincides spatially with the loop-top X-Ray source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We interpret this large redshift as the signature of downward-moving reconnection outflow/hot retracting loops. Imaging observations from both IRIS and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) also reveal the eruption and reconnection processes. Fast downward-propagating blobs along these loops are also found from cool emission lines (e.g., Si {\sc{iv}}, O {\sc{iv}}, C {\sc{ii}}, Mg {\sc{ii}}) and images of AIA and IRIS. Furthermore, the entire Fe {\sc{xxi}} line is blueshifted by $\sim$260 km s$^{-1}$ at the loop footpoints, where the cool lines mentioned above all exhibit obvious redshift, a result that is consistent with the scenario of chromospheric evaporation induced by downward-propagating nonthermal electrons from the reconnection site.

Two solar type II radio bursts, separated by ~24 hours in time, are examined together. Both events are associated with coronal mass ejections (CMEs) erupting from the same active region (NOAA 11176) beneath a well-observed helmet streamer. We find that the type II emissions in both events ended once the CME/shock fronts passed the white-light streamer tip, which is presumably the magnetic cusp of the streamer. This leads us to conjecture that the closed magnetic arcades of the streamer may play a role in electron acceleration and type II excitation at coronal shocks. To examine such a conjecture, we conduct a test-particle simulation for electron dynamics within a large-scale partially-closed streamer magnetic configuration swept by a coronal shock. We find that the closed field lines play the role of an electron trap, via which the electrons are sent back to the shock front for multiple times, and therefore accelerated to high energies by the shock. Electrons with an initial energy of 300eV can be accelerated to tens of keV concentrating at the loop apex close to the shock front with a counter-streaming distribution at most locations. These electrons are energetic enough to excite Langmuir waves and radio bursts. Considering the fact that most solar eruptions originate from closed field regions, we suggest that the scenario may be important to the generation of more metric type IIs. This study also provides an explanation to the general ending frequencies of metric type IIs at or above 20-30 MHz and the disconnection issue between metric and interplanetary type IIs.

The Earth's obliquity is stabilized by the Moon, which facilitates a rapid precession of the Earth's spin-axis, de-tuning the system away from resonance with orbital modulation. It is however, likely that the architecture of the Solar System underwent a dynamical instability-driven transformation, where the primordial configuration was more compact. Hence, the characteristic frequencies associated with orbital perturbations were likely faster in the past, potentially allowing for secular resonant encounters. In this work we examine if at any point in the Earth's evolutionary history, the obliquity varied significantly. Our calculations suggest that even though the orbital perturbations were different, the system nevertheless avoided resonant encounters throughout its evolution. This indicates that the Earth obtained its current obliquity during the formation of the Moon.

Kepler-56 is a multi-planet system containing two coplanar inner planets that are in orbits misaligned with respect to the spin axis of the host star, and an outer planet. Various mechanisms have been proposed to explain the broad distribution of spin-orbit angles among exoplanets, and these theories fall under two broad categories. The first is based on dynamical interactions in a multi-body system, while the other assumes that disk migration is the driving mechanism in planetary configuration and that the star (or disk) is titled with respect to the planetary plane. Here we show that the large observed obliquity of Kepler-56 system is consistent with a dynamical origin. In addition, we use observations by Huber et al. (2013) to derive the obliquity's probability distribution function, thus improving the constrained lower limit. The outer planet may be the cause of the inner planets' large obliquities, and we give the probability distribution function of its inclination, which depends on the initial orbital configuration of the planetary system. We show that even in the presence of precise measurement of the true obliquity, one cannot distinguish the initial configurations. Finally we consider the fate of the system as the star continues to evolve beyond the main sequence, and we find that the obliquity of the system will not undergo major variations as the star climbs the red giant branch. We follow the evolution of the system and find that the innermost planet will be engulfed in ~129 Myr. Furthermore we put an upper limit of ~155 Myr for the engulfment of the second planet. This corresponds to ~ 3% of the current age of the star.

In this Letter, we investigate the effects of dark energy on $P-V$ criticality of charged AdS black holes by considering the case of the RN-AdS black holes surrounded by quintessence. By treating the cosmological constant as thermodynamic pressure, we study its thermodynamics in the extended phase space. It is shown that quintessence dark energy does not affect the existence of small/large black hole phase transition. For the case $\omega_q=-2/3$ we derive analytic expressions of critical physical quantities, while for cases $\omega_q\neq-2/3$ we appeal to numerical method for help. It is shown that quintessence dark energy affects the critical physical quantities near the critical point. Critical exponents are also calculated. They are exactly the same as those obtained before for arbitrary other AdS black holes, which implies that quintessence dark energy does not change the critical exponents.

The Kozai-Lidov mechanism can be applied to a vast variety of astrophysical systems involving hierarchical three-body systems. Here, we study the Kozai-Lidov mechanism systematically in the test particle limit at the octupole level of approximation. We investigate the chaotic and quasiperiodic orbital evolution by studying surfaces of section and the Lyapunov exponents. We find that the resonances introduced by the octupole level of approximation cause orbits to flip from prograde to retrograde and back as well as cause significant eccentricity excitation, and the chaotic behaviors occur when the mutual inclination between the inner and the outer binary is high. We characterize the parameter space that allows large amplitude oscillations in eccentricity and inclination.

The variation of a planet's obliquity is influenced by the existence of satellites with a high mass ratio. For instance, the Earth's obliquity is stabilized by the Moon, and would undergo chaotic variations in the Moon's absence. In turn, such variations can lead to large-scale changes in the atmospheric circulation, rendering spin-axis dynamics a central issue for understanding climate. The relevant quantity for dynamically-forced climate change is the rate of chaotic diffusion. Accordingly, here we reexamine the spin-axis evolution of a Moonless Earth within the context of a simplified perturbative framework. We present analytical estimates of the characteristic Lyapunov coefficient as well as the chaotic diffusion rate and demonstrate that even in absence of the Moon, the stochastic change in the Earth's obliquity is sufficiently slow to not preclude long-term habitability. Our calculations are consistent with published numerical experiments and illustrate the putative system's underlying dynamical structure in a simple and intuitive manner.

Simultaneous radio and extreme ultraviolet (EUV)/white-light imaging data are examined for a solar type II radio burst occurring on 2010 March 18 to deduce its source location. Using a bow-shock model, we reconstruct the 3-dimensional EUV wave front (presumably the type-II emitting shock) based on the imaging data of the two STEREO spacecraft. It is then combined with the Nan\c{c}ay radio imaging data to infer the 3-dimensional position of the type II source. It is found that the type II source coincides with the interface between the CME EUV wave front and a nearby coronal ray structure, providing evidence that the type II emission is physically related to the CME-ray interaction. This result, consistent with those of previous studies, is based on simultaneous radio and EUV imaging data for the first time.

The dynamical evolution of a hierarchical three body system is well characterized by the eccentric Kozai-Lidov mechanism, where the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40-140 degree) can produce long timescale modulations that drives the eccentricity to extremely large value and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ~180 degree, rolling over its major axis. The ~180 degree flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exo-planetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems.

It is of great importance to track the solar wind back to its photospheric source region and identify the related current sheets; this will provide key information for investigating the origin and predictions of the solar wind. We report a statistical study relating the photospheric footpoint motion and in-situ observation of current sheets in the solar wind. We used the potential force-free source--surface (PFSS) model and the daily synoptic charts to trace the solar wind back from 1 AU, as observed by the Wind spacecraft, to the solar surface. As the footpoints move along the solar surface we obtain a time series of the jump times between different points. These jumps can be within a cell and between adjacent cells. We obtained the distribution of the jump times and the distribution for a subset of the jump times in which only jumps between adjacent cells were counted. For both cases, the distributions clearly show two populations. These distributions are compared with the distribution of in-situ current sheets reported in an earlier work of Miao Peng and Li (2011). Its implications on the origin of the current sheets are discussed.

We present the observation of a major solar eruption that is associated with fast sunspot rotation. The event includes a sigmoidal filament eruption, a coronal mass ejection, and a GOES X2.1 flare from NOAA active region 11283. The filament and some overlying arcades were partially rooted in a sunspot. The sunspot rotated at $\sim$10$^\circ$ per hour rate during a period of 6 hours prior to the eruption. In this period, the filament was found to rise gradually along with the sunspot rotation. Based on the HMI observation, for an area along the polarity inversion line underneath the filament, we found gradual pre-eruption decreases of both the mean strength of the photospheric horizontal field ($B_h$) and the mean inclination angle between the vector magnetic field and the local radial (or vertical) direction. These observations are consistent with the pre-eruption gradual rising of the filament-associated magnetic structure. In addition, according to the Non-Linear Force-Free-Field reconstruction of the coronal magnetic field, a pre-eruption magnetic flux rope structure is found to be in alignment with the filament, and a considerable amount of magnetic energy was transported to the corona during the period of sunspot rotation. Our study provides evidences that in this event sunspot rotation plays an important role in twisting, energizing, and destabilizing the coronal filament-flux rope system, and led to the eruption. We also propose that the pre-event evolution of $B_h$ may be used to discern the driving mechanism of eruptions.

We report observations of the acceleration and trapping of energetic ions and electrons between a pair of corotating interaction regions (CIRs). The event occurred in Carrington Rotation 2060. Observed at spacecraft STEREO-B, the two CIRs were separated by less than 5 days. In contrast to other CIR events, the fluxes of energetic ions and electrons in this event reached their maxima between the trailing-edge of the first CIR and the leading edge of the second CIR. The radial magnetic field (Br) reversed its sense and the anisotropy of the flux also changed from sunward to anti-sunward between the two CIRs. Furthermore, there was an extended period of counter-streaming suprathermal electrons between the two CIRs. Similar observations for this event were also obtained for ACE and STEREO-A. We conjecture that these observations were due to a "U-shape" large scale magnetic field topology connecting the reverse shock of the first CIR and the forward shock of the second CIR. Such a disconnected U-shaped magnetic field topology may have formed due to magnetic reconnection in the upper corona.

Formation of stars is now believed to be tightly linked to the dynamical evolution of interstellar filaments in which they form. In this paper we analyze the density structure and kinematics of a small network of infrared dark filaments, SDC13, observed in both dust continuum and molecular line emission with the IRAM 30m telescope. These observations reveal the presence of 18 compact sources amongst which the two most massive, MM1 and MM2, are located at the intersection point of the parsec-long filaments. The dense gas velocity and velocity dispersion observed along these filaments show smooth, strongly correlated, gradients. We discuss the origin of the SDC13 velocity field in the context of filament longitudinal collapse. We show that the collapse timescale of the SDC13 filaments (from 1 Myr to 4 Myr depending on the model parameters) is consistent with the presence of Class I sources in them, and argue that, on top of bringing more material to the centre of the system, collapse could generate additional kinematic support against local fragmentation, helping the formation of starless super-Jeans cores.

Star formation occurs in molecular gas. In previous studies, the structure of the molecular gas has been studied in terms of molecular clouds, but has been overlooked beyond the cloud scale. We present an observational study of the molecular gas at 49.5 degree <l<52.5 degree and -5.0 km/s <v_lsr <17.4 km/s. The molecular gas is found in the form of a huge (>= 500 pc) filamentary gas wisp. This has a large physical extent and a velocity dispersion of ~5 km/s. The eastern part of the filamentary gas wisp is located ~130 pc above the Galactic disk (which corresponds to 1.5-4 e-folding scale-heights), and the total mass of the gas wisp is >= 1 X 10^5 M_sun. It is composed of two molecular clouds and an expanding bubble. The velocity structure of the gas wisp can be explained as a smooth quiescent component disturbed by the expansion of a bubble. That the length of the gas wisp exceeds by much the thickness of the molecular disk of the Milky Way is consistent with the cloud-formation scenario in which the gas is cold prior to the formation of molecular clouds. Star formation in the filamentary gas wisp occurs at the edge of a bubble (G52L nebula), which is consistent with some models of triggered star formation.

Protostellar outflow is a prominent process that accompanies the formation of stars. It is generally agreed that wide-angled protostellar outflows come from the interaction between the wind from a forming star and the ambient gas. However, it is still unclear how the interaction takes place. In this work, we theoretically investigate the possibility that the outflow results from interaction between the wind and the ambient gas in the form of turbulent entrainment. In contrast to the previous models, turbulent motion of the ambient gas around the protostar is taken into account. In our model, the ram-pressure of the wind balances the turbulent ram-pressure of the ambient gas, and the outflow consists of the ambient gas entrained by the wind. The calculated outflow from our modelling exhibits a conical shape. The total mass of the outflow is determined by the turbulent velocity of the envelope as well as the outflow age, and the velocity of the outflow is several times higher than the velocity dispersion of the ambient gas. The outflow opening angle increases with the strength of the wind and decreases with the increasing ambient gas turbulence. The outflow exhibits a broad line width at every position. We propose that the turbulent entrainment process, which happens ubiquitously in nature, plays a universal role in shaping protostellar outflows.

A physical model of Geant4-based simulation of galactic cosmic ray (GCR) particles interaction with the lunar surface matter has been developed to investigate the production rate of cosmogenic nuclei. In this model the GCRs, mainly very high energy protons and $\alpha$ particles, bombard the surface of the Moon and produce many secondary particles such as protons and neutrons. The energies of proton and neutron at different depths are recorded and saved into ROOT files, and the analytical expressions for the differential proton and neutron fluxes are obtained through the best-fit procedure under the ROOT software. To test the validity of this model, we calculate the production rates of long-lived nuclei $^{10}$Be and $^{26}$Al in the Apollo 15 long drill core by combining the above differential fluxes and the newly evaluated spallation reaction cross sections. Numerical results show that the theoretical production rates agree quite well with the measured data. It means that this model works well. Therefore, it can be expected that this model can be used to investigate the cosmogenic nuclei in lunar samples returned by Chinese lunar exploration program and can be extended to study other objects, such as the meteorites and the Earth's atmosphere.

The Canadian Penning Trap mass spectrometer has made mass measurements of 33 neutron-rich nuclides provided by the new Californium Rare Isotope Breeder Upgrade (CARIBU) facility at Argonne National Laboratory. The studied region includes the 132Sn double shell closure and ranges in Z from In to Cs, with Sn isotopes measured out to A = 135, and the typical measurement precision is at the 100 ppb level or better. The region encompasses a possible major waiting point of the astrophysical r process, and the impact of the masses on the r process is shown through a series of simulations. These first-ever simulations with direct mass information on this waiting point show significant increases in waiting time at Sn and Sb in comparison with commonly used mass models, demonstrating the inadequacy of existing models for accurate r-process calculations.

Two major processes have been proposed to convert the coronal magnetic energy into the kinetic energy of a coronal mass ejection (CME): resistive magnetic reconnection and ideal macroscopic magnetohydrodynamic instability of magnetic flux rope. However, it remains elusive whether both processes play a comparable role or one of them prevails during a particular eruption. To shed light on this issue, we carefully studied energetic but flareless CMEs, \textit{i.e.}, fast CMEs not accompanied by any flares. Through searching the Coordinated Data Analysis Workshops (CDAW) database of CMEs observed in Solar Cycle 23, we found 13 such events with speeds larger than 1000 km s$^{-1}$. Other common observational features of these events are: (1) none of them originated in active regions; they were associated with eruptions of well-developed long filaments in quiet-Sun regions, (2) no apparent enhancement of flare emissions was present in soft X-ray, EUV and microwave data. Further studies of two events reveal that (1) the reconnection electric fields, as inferred from the product of the separation speed of post-eruption ribbons and the photospheric magnetic field measurement, were in general weak; (2) the period with a measurable reconnection electric field is considerably shorter than the total filament-CME acceleration time. These observations indicate that, for these fast CMEs, the magnetic energy was released mainly via the ideal flux rope instability through the work done by the large scale Lorentz force acting on the rope currents rather than via magnetic reconnections. We also suggest that reconnections play a less important role in accelerating CMEs in quiet Sun regions of weak magnetic field than those in active regions of strong magnetic field.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

Current sheets are ubiquitous in the solar wind.They are a major source of the solar wind MHD turbulence intermittency. They may result from non-linear interactions of the solar wind MHD turbulence or are the boundaries of flux tubes that originate from the solar surface. Some current sheets appear in pairs and are the boundaries of transient structures such as magnetic holes and reconnection exhausts, or the edges of pulsed Alfv\'{e}n waves. For an individual current sheet, discerning whether it is a flux tube boundary or due to non-linear interactions, or the boundary of a transient structure is difficult. In this work, using data from the {\sl Wind} spacecraft, we identify two three-current-sheet events. Detailed examination of these two events suggest that they are best explained by the flux tube crossing scenario. Our study provides a convincing evidence supporting the scenario that the solar wind consists of flux tubes where distinct plasmas reside.

It has been noted for a long time that the spectra of observed continuum emissions in many solar flares are consistent with double power laws with a hardening at energies $\sim > $ 300 keV. It is now largely believed that at least in electron-dominated events the hardening in photon spectrum reflects an intrinsic hardening in the source electron spectrum. In this paper, we point out that a power law spectrum of electron with a hardening at high energies can be explained by diffusive shock acceleration of electrons at a termination shock with a finite width. Our suggestion is based on an early analytical work by Drury et al., where the steady state transport equation at a shock with a tanh profile was solved for a $p$-independent diffusion coefficient. Numerical simulations with a $p$-dependent diffusion coefficient show hardenings in the accelerated electron spectrum which are comparable with observations. One necessary condition for our proposed scenario to work is that high energy electrons resonate with the inertial range of the MHD turbulence and low energy electrons resonate with the dissipation range of the MHD turbulence at the acceleration site, and the spectrum of the dissipation range $\sim k^{-2.7}$. A $\sim k^{-2.7}$ dissipation range spectrum is consistent with recent solar wind observations.

Some reports of supernova (SN) discoveries turn out not to be true core-collapse explosions. One such case was SN 2009ip, which was recognized to be a luminous blue variable (LBV) eruption. This source had a massive hot progenitor star identified in pre-explosion data, it had documented evidence of pre-outburst variability, and it was subsequently discovered to have a 2nd outburst in 2010. This same source rebrightened again in 2012, and early spectra showed the same narrow-line profiles as before, suggesting another LBV-like eruption. We present new photometry and spectroscopy of SN 2009ip, indicating that it has transitioned into a true SN. The most striking discovery in these data is that unlike previous reports, the spectrum exhibited Balmer lines with very broad P-Cygni profiles characteristic of normal Type II supernovae (SNe II), in addition to narrow emission lines seen in SNe IIn and LBVs. Emission components have FWHM 8000 km/s, while the P-Cygni absorption component has blue wings extending to -13,000 km/s. These velocities are typical of SNe II, but have never been associated with emission lines from a nonterminal LBV-like eruption. Initially, the peak absolute magnitude seemed fainter than that of normal SNe. However, after a brief period of fading, the source quickly brightened again to M_R=-17.5 mag over a couple days. The broad lines mostly disappeared, and the spectrum began to resemble the early optically thick phases of SNe IIn. Two weeks later the source leveled off near -18 mag, after which broad emission lines again developed in the spectrum as the source faded. We conclude that the 2012 outburst of SN 2009ip was the result of a true core-collapse SN IIn that occured when the progenitor star was in an LBV-like outburst phase, and where the SN was initially faint and then rapidly brightened due to interaction with circumstellar material (abridged).

In recent studies (Feng et al., 2012; Kong et al., 2012), we proposed that source properties of type II radio bursts can be inferred through a causal relationship between the special shape of the type II dynamic spectrum (e.g., bump or break) and simultaneous extreme ultraviolet (EUV)/white light imaging observations (e.g., CME-shock crossing streamer structures). As a further extension of these studies, in this paper we examine the CME event dated on December 31 2007 associated with a multiple type II radio burst. We identify the presence of two spectral bump features on the observed dynamic spectrum. By combining observational analyses of the radio spectral observations and the EUV-white light imaging data, we conclude that the two spectral bumps are resulted from CME-shock propagating across dense streamers on the southern and northern sides of the CME, respectively. It is inferred that the corresponding two type II emissions originate separately from the two CME-shock flanks where the shock geometries are likely quasi-perpendicular or oblique. Since the emission lanes are bumped as a whole within a relatively short time, it suggests that the type IIs with bumps of the study are emitted from spatially confined sources (with a projected lateral dimension smaller than 0.05-0.1 R$_\odot$ at a fundamental frequency level of 20-30 MHz).

An arbitrary surface mass density of gravitational lens can be decomposed into multipole components. We simulate the ray-tracing for the multipolar mass distribution of generalized SIS (Singular Isothermal Sphere) model, based on the deflection angles which are analytically calculated. The magnification patterns in the source plane are then derived from inverse shooting technique. As have been found, the caustics of odd mode lenses are composed of two overlapping layers for some lens models. When a point source traverses such kind of overlapping caustics, the image numbers change by \pm 4, rather than \pm 2. There are two kinds of images for the caustics. One is the critical curve and the other is the transition locus. It is found that the image number of the fold is exactly the average value of image numbers on two sides of the fold, while the image number of the cusp is equal to the smaller one. We also focus on the magnification patterns of the quadrupole (m = 2) lenses under the perturbations of m = 3, 4 and 5 mode components, and found that one, two, and three butterfly or swallowtail singularities can be produced respectively. With the increasing intensity of the high-order perturbations, the singularities grow up to bring sixfold image regions. If these perturbations are large enough to let two or three of the butterflies or swallowtails contact, eightfold or tenfold image regions can be produced as well. The possible astronomical applications are discussed.

We study the effects on the number counts of sub-millimetre galaxies due to gravitational lensing. We explore the effects on the magnification cross section due to halo density profiles, ellipticity and cosmological parameter (the power-spectrum normalisation $\sigma_8$). We show that the ellipticity does not strongly affect the magnification cross section in gravitational lensing while the halo radial profiles do. Since the baryonic cooling effect is stronger in galaxies than clusters, galactic haloes are more concentrated. In light of this, a new scenario of two halo population model is explored where galaxies are modeled as a singular isothermal sphere profile and clusters as a Navarro, Frenk and White (NFW) profile. We find the transition mass between the two has modest effects on the lensing probability. The cosmological parameter $\sigma_8$ alters the abundance of haloes and therefore affects our results. Compared with other methods, our model is simpler and more realistic. The conclusions of previous works is confirm that gravitational lensing is a natural explanation for the number count excess at the bright end.

Monitoring the orbits of stars around Sgr A* offers the possibility of detecting the precession of their orbital planes due to frame dragging, of measuring the spin and quadrupole moment of the black hole, and of testing the no-hair theorem. Here we investigate whether the deviations of stellar orbits from test-particle trajectories due to wind mass loss and tidal dissipation of the orbital energy compromise such measurements. We find that the effects of stellar winds are, in general, negligible. On the other hand, for the most eccentric orbits (e>0.96) for which an optical interferometer, such as GRAVITY, will detect orbital plane precession due to frame dragging, the tidal dissipation of orbital energy occurs at timescales comparable to the timescale of precession due to the quadrupole moment of the black hole. As a result, this non-conservative effect is a potential source of systematic uncertainty in testing the no-hair theorem with stellar orbits.

We consider the long-term tidal heating of a star by the supermassive black hole at the Galactic center, SgrA*. We show that gravitational interaction with background stars leads to a linear growth of the tidal excitation energy with the number of pericenter passages near SgrA*. The accumulated heat deposited by excitation of modes within the star over many pericenter passages can lead to a runaway disruption of the star at a pericenter distance that is 4-5 times farther than the standard tidal disruption radius. The accumulated heating may explain the lack of massive ($\gtrsim 10M_{\odot}$) S-stars closer than several tens of AU from SgrA*.

The morphologic properties of the magnetic networks during Carrington Rotations (CR) 1955 to 2091 (from 1999 to 2010) have been analyzed by applying the watershed algorithm to magnetograms observed by the Michelson Doppler Interferometer (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft. We find that the average area of magnetic cells on the solar surface at lower latitudes (within +-50 degree) are smaller than those at higher latitudes (beyond +-50 degree). Statistical analysis of these data indicates that the magnetic networks are of fractal in nature, and the average fractal dimension is D_f = 1.253+-0.011. We also find that both the fractal dimension and the size of the magnetic networks are anti-correlated with the sunspot area. This is perhaps because a strong magnetic field can suppress spatially modulated oscillation, compress the boundaries of network cells, leading to smoother cell boundaries. The fractal dimension of the cell deviates that predicted from an isobar of Kolmogorov homogeneous turbulence.

We present the results from the first public blind PSF reconstruction challenge, the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Star Challenge. Reconstruction of a spatially varying PSF, sparsely sampled by stars, at non-star positions is a critical part in the image analysis for weak lensing where inaccuracies in the modelled ellipticity and size-squared can impact the ability to measure the shapes of galaxies. This is of importance because weak lensing is a particularly sensitive probe of dark energy, and can be used to map the mass distribution of large scale structure. Participants in the challenge were presented with 27,500 stars over 1300 images subdivided into 26 sets, where in each set a category change was made in the type or spatial variation of the PSF. Thirty submissions were made by 9 teams. The best methods reconstructed the PSF with an accuracy of ~0.00025 in ellipticity and ~0.00074 in size squared. For a fixed pixel scale narrower PSFs were found to be more difficult to model than larger PSFs, and the PSF reconstruction was severely degraded with the inclusion of an atmospheric turbulence model (although this result is likely to be a strong function of the amplitude of the turbulence power spectrum).

Model fitting is frequently used to determine the shape of galaxies and the point spread function, for examples, in weak lensing analyses or morphology studies aiming at probing the evolution of galaxies. However, the number of parameters in the model, as well as the number of objects, are often so large as to limit the use of model fitting for future large surveys. In this article, we propose a set of algorithms to speed up the fitting process. Our approach is divided into three distinctive steps: centroiding, ellipticity measurement, and profile fitting. We demonstrate that we can derive the position and ellipticity of an object analytically in the first two steps and thus leave only a small number of parameters to be derived through model fitting. The position, ellipticity, and shape parameters can then used in constructing orthonomal basis functions such as s\'ersiclets for better galaxy image reconstruction. We assess the efficiency and accuracy of the algorithms with simulated images. We have not taken into account the deconvolution of the point spread function, which most weak lensing analyses do.

Seesaw mechanism provides a natural explanation of light neutrino masses through suppression of heavy seesaw scale. In inverse seesaw models the seesaw scale can be much lower than that in the usual seesaw models. If terms inducing seesaw masses are further induced by loop corrections, the seesaw scale can be lowered to be in the range probed by experiments at the LHC without fine tuning. In this paper we construct models in which inverse seesaw neutrino masses are generated at two loop level. These models also naturally have dark matter candidates. Although the recent data from Xenon100 put stringent constraint on the models, they can be consistent with data on neutrino masses, mixing, dark matter relic density and direct detection. These models also have some interesting experimental signatures for collider and flavor physics.

We investigate the electromagnetic (EM) counterpart of gravitational waves (GWs) emitted by a supermassive black hole binary (SMBHB) through the viscous dissipation of the GW energy in an accretion disk and stars surrounding the SMBHB. We account for the suppression of the heating rate if the forcing period is shorter than the turnover time of the largest turbulent eddies. We find that the viscous heating luminosity in 0.1 solar mass stars can be significantly higher than their intrinsic luminosity, but still too low to be detected for extragalactic sources. The relative brightening is small for accretion disks.

It has been suggested that type II radio bursts are due to energetic electrons accelerated at coronal shocks. Radio observations, however, have poor or no spatial resolutions to pinpoint the exact acceleration locations of these electrons. In this paper, we discuss a promising approach to infer the electron acceleration location by combining radio and white light observations. The key assumption is to relate specific morphological features (e.g. spectral bumps) of the dynamic spectra of type II radio bursts, to imaging features (e.g. CME going into a streamer) along the CME (and its driven shock) propagation. In this study, we examine the CME-streamer interaction for the solar eruption dated on 2003 November 1. The presence of spectral bump in the relevant type II radio burst is identified, which is interpreted as a natural result of the shock-radio emitting region entering the dense streamer structure. The study is useful for further determinations of the location of type II radio burst and the associated electron acceleration by CME-driven shock.

We discuss an intriguing type II radio burst that occurred on 2011 March 27. The dynamic spectrum was featured by a sudden break at about 43 MHz on the well-observed harmonic branch. Before the break, the spectrum drifted gradually with a mean rate of about -0.05 MHz/s. Following the break, the spectrum jumped to lower frequencies. The post-break emission lasted for about three minutes. It consisted of an overall slow drift which appeared to have a few fast drift sub-bands. Simultaneous observations from the Solar TErrestrial RElations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) were also available and are examined for this event. We suggest that the slow-drift period before the break was generated inside a streamer by a coronal eruption driven shock, and the spectral break as well as the relatively wide spectrum after the break is a consequence of the shock crossing the streamer boundary where density drops abruptly. It is suggested that this type of radio bursts can be taken as a unique diagnostic tool for inferring the coronal density structure, as well as the radio emitting source region.

Active reflector is one of the key technologies for constructing large telescopes, especially for the millimeter/sub-millimeter radio telescopes. This article introduces a new efficient laser angle metrology system for the active reflector antenna of the large radio telescopes, with a plenty of active reflector experiments mainly about the detecting precisions and the maintaining of the surface shape in real time, on the 65-meter radio telescope prototype constructed by Nanjing Institute of Astronomical Optics and Technology (NIAOT). The test results indicate that the accuracy of the surface shape segmenting and maintaining is up to micron dimension, and the time-response can be of the order of minutes. Therefore, it is proved to be workable for the sub-millimeter radio telescopes.

We investigate the evolution of the MBH-{\sigma} relation by examining the relationship between the intrinsic scatter in the MBH-{\sigma} relation and galaxy bolometric nuclear luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our sample is composed of galaxies with classical bulges when possible, of which 38 have dynamically measured BHs masses, and 17 have BHs masses measured by reverberation mapping. In order to obtain the bolometric nuclear luminosity for galaxies with low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra to bolometric luminosity. We find that the scatter in the MBH-{\sigma} relation is uncorrelated with nuclear luminosity over seven orders of magnitude in luminosity, with the high luminosity end approaching the Eddington luminosity. This suggests that at the present epoch galaxies evolve along the MBH-{\sigma} relation. This conclusion is consistent with the standard paradigm that BHs grow contemporaneously with their host stellar spheroids.

In this paper, with a survey through the Large Angle and Spectrometric Coronagraph (LASCO) data from 1996 to 2009, we present 11 events with plasma blobs flowing outwards sequentially along a bright coronal ray in the wake of a coronal mass ejection. The ray is believed to be associated with the current sheet structure that formed as a result of solar eruption, and the blobs are products of magnetic reconnection occurring along the current sheet. The ray morphology and blob dynamics are investigated statistically. It is found that the apparent angular widths of the rays at a fixed time vary in a range of 2.1-6.6 (2.0-4.4) degrees with an average of 3.5 (2.9) degrees at 3 (4) Rs, respectively, and the observed durations of the events vary from 12 h to a few days with an average of 27 h. It is also found, based on the analysis of blob motions, that 58% (26) of the blobs were accelerated, 20% (9) were decelerated, and 22% (10) moved with a nearly-constant speed. Comparing the dynamics of our blobs and those that are observed above the tip of a helmet streamer, we find that the speeds and accelerations of the blobs in these two cases differ significantly. It is suggested that these differences of the blob dynamics stem from the associated magnetic reconnection involving different magnetic field configurations and triggering processes.

We present line lists for ammonia (NH$_{3}$) at high temperatures obtained by recording Fourier transform infrared emission spectra. Calibrated line lists are presented for twelve temperatures (300 -- 1300$^{\circ}$C in 100$^{\circ}$C intervals and 1370$^{\circ}$C) and each line list covers the 740 -- 2100 cm$^{-1}$ range, which includes the majority of the $\nu_{2}$ umbrella bending mode region (11 $\mu$m) and the $\nu_{4}$ asymmetric bending mode region (6.2 $\mu$m). We also demonstrate the useful technique of obtaining empirical lower state energies from a set of spectra recorded at different sample temperatures. Using our NH$_{3}$ spectra, we have estimated lower state energies ($E_{Low}$ in cm$^{-1}$) and our values have been incorporated into the line lists along with calibrated wavenumbers ($\tilde{\nu}$ in cm$^{-1}$) and calibrated line intensities ($S'$ in cm/molecule) for each line. We expect our hot NH$_{3}$ line lists to find direct application in the modelling of planetary atmospheres and brown dwarfs.

Galaxy clusters are considered as excellent probes for cosmology. For that purpose, their mass needs to be measured and their structural properties needs to be understood. We propose a method for galaxy cluster mass reconstruction which combines information from strong lensing, weak lensing shear and flexion. We extend the weak lensing analysis to the inner parts of the cluster and, in particular, improve the resolution of substructure. We use simulations to show that the method recovers the mass density profiles of the cluster. We find that the weak lensing flexion is sensitive to substructure. After combining the flexion data into the joint weak and strong lensing analysis, we can resolve the cluster properties with substructures.

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are investigated in this work. The resulting emergent spectra are calculated using the relativistic ray-tracing method, in which we neglect the self-irradiation of the accretion disk. The angular dependence of the disk luminosity, the effects of the heat advection and the disk thickness on the estimation of the black hole spin are discussed. Compare to the previous works, our improvements are that we use the self-consistent disk equations and we consider the disk self-shadowing effect. We find that at the moderate accretion rate, the radiation trapped in the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to the emergent spectra. At the high accretion rate, for the large inclination and large black hole spin, both the disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these effects will influence the measurement of the black hole spin based on the spectra fitting and influence the angular dependence of the luminosity. For the disks around Kerr black holes with $a=0.98$, if the disk inclination is greater than $60^\circ$, and their luminosity is beyond 0.2 Eddington luminosity, the spectral model which is based on the relativistic standard accretion disk is no longer applicable for the spectra fitting. We also confirm that the effect of the self-shadowing is significantly enhanced by the light-bending, which implies that the non-relativistic treatment of the self-shadowing is inaccurate. According to our results, the observed luminosity dependence of the measured spin suggests that the disk self-shadowing significantly shapes the spectra of GRS 1915+105, which might lead to the underestimation of the black hole spin for the high luminosity states.

We use a high-resolution $N$-body simulation to investigate the influence of background galaxy properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out for 10 massive clusters at two redshifts, i.e. $z_{\rm l} \sim 0.2$ and 0.3. The virial mass ($M_{\rm vir}$) of the simulated lens clusters at $z\sim0.2$ ranges from $6.8\times10^{14} h^{-1} {M_{\odot}}$ to $1.1\times 10^{15} h^{-1} M_{\odot}$. The information of background galaxies brighter than 25 magnitude in the $I$-band is taken from Cosmological Evolution Survey (COSMOS) imaging data. Around $1.7\times 10^5$ strong lensing realizations with these images as background galaxies have been performed for each set. We find that the efficiency for forming giant arcs for $z_{\rm l}=0.2$ clusters is broadly consistent with observations. The efficiency of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our principal finding is that a small proportion ($\sim 1/3$) of galaxies with elongated shapes (e.g. ellipticity $\epsilon=1-b/a>0.5$) can boost the number of giant arcs substantially. Compared with recent studies where a uniform ellipticity distribution from 0 to 0.5 is used for the sources, the adoption of directly observed shape distribution increases the number of giant arcs by a factor of $\sim2$. Our results indicate that it is necessary to account for source information and survey parameters (such as point-spread-function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological parameters.(abridged)

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

Lyman-$\alpha$ spheres, i.e. regions around the first stars which are illuminated by Lyman-$\alpha$ photons and show 21cm absorption feature against the CMB, are smoking guns at the dawn of the reionization epoch. Though overwhelming radio foreground makes their detections extremely difficult, we pointed out that, strong gravitational lensing can significantly improve their observational feasibility. Since Lyman-$\alpha$ spheres have ~10" sizes, comparable to the caustic size of galaxy clusters, individual images of each strongly lensed Lyman-$\alpha$ sphere often merge together and form single structures in the 21cm sky with irregular shapes. Using high-resolution N-body LCDM simulations, we found that the lensing probability to have magnification bigger than 10 is ~10^{-5}. This results in $\ga 10^6$ strongly lensed Lyman-$\alpha$ spheres across the sky, which should be the primary targets for first detections of Lyman-$\alpha$ spheres. Although the required total radio array collecting area for their detection is large (~100 km^2), the design of long fixed cylindrical reflectors can significantly reduce the total cost of such array to the level of the square kilometer array (SKA) and makes the detection of these very first objects feasible.

We use high-resolution N-body numerical simulations to study the number of predicted large-separation multiply-imaged systems produced by clusters of galaxies in the SDSS photometric and spectroscopic quasar samples. We incorporate the condensation of baryons at the centre of clusters by (artificially) adding a brightest central galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cosmological models: a LCDM model with a matter density $\Omega_0=0.3$, and $\sigma_8=0.9$ (LCDM0), and a model favoured by the WMAP three-year data with $\Omega_0=0.238$, and $\sigma_8=0.74$ (WMAP3). We found that the predicted multiply-imaged quasars with separation >10" is about 6.2 and 2.6 for the SDSS photometric (with an effective area 8000 deg$^2$) and spectroscopic (with an effective area 5000 deg$^2$) quasar samples respectively in the LCDM0 model; the predicted numbers of large-separation lensed quasars agree well with the observations. These numbers are reduced by a factor of 7 or more in the WMAP3 model, and are consistent with data at < 8% level. The predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.3 and 1. We find that the BCG creates a central circular region, comparable to the Einstein ring of the BCG, where the central image disappears in the usual three-image and five-image configurations. If we include four image systems as an extreme case of five-image systems (with an infinitely demagnified central image), we find that 68% of the central images are fainter by a factor of 100 than the brightest image, and about 80% are within 1.5"of the BCG.

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong lensing events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, at constraining cosmological parameters from statistical lensing, or at determining the inner structure of galaxy clusters through gravitational arcs. We do so by creating two-dimensional gradually smoothed, differently elliptical and asymmetric versions of some numerical models. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~40%, ~10% and ~30% of the total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and asymmetric clusters, even if located at large distances from the cluster centers (~1Mpc/h). Conversely, their effect is less important in highly asymmetric lenses. The morphology, position and flux of individual arcs are strongly affected by the presence of substructures in the clusters. Removing substructures on spatial scales <~50kpc/h, roughly corresponding to mass scales <~5 10^{10}M_\odot/h, alters the image multiplicity of ~35% of the sources used in the simulations and causes position shifts larger than 5'' for ~40% of the arcs longer than 5''. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong gravitational lensing is potentially a powerfull tool to measure the level of substructures and asymmetries in clusters.

We use high-resolution $N$-body simulations to investigate the optical depth of giant arcs with length-to-width ratio larger than 7.5 and 10 in the `standard' $\LCDM$ model with $\sigma_8=0.9$ and $\Omega_{\rm m,0}=0.3$ and a model based on three-year Wilkinson Microwave Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the lensing probability in the three-year WMAP model (with $\sigma_8=0.74$ and $\Omega_{\rm m,0}=0.238)$ decreases by a factor of $\sim 6$ compared with that in the `standard' $\LCDM$ model. The effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue that baryons will only increase the the lensing cross-section by a moderate factor, $\sim 2$. We conclude that the low central value of $\sigma_8$ and $\Omega_{\rm m,0}$ preferred by the WMAP three-year data may be too low to be compatible with observations if conventional assumptions of the background source population are correct.

We use high-resolution N-body simulations to study the galaxy-cluster cross-sections and the abundance of giant arcs in the $\Lambda$CDM model. Clusters are selected from the simulations using the friends-of-friends method, and their cross-sections for forming giant arcs are analyzed. The background sources are assumed to follow a uniform ellipticity distribution from 0 to 0.5 and to have an area identical to a circular source with diameter $1\arcsec$. We find that the optical depth scales as the source redshift approximately as $\tau_{1''} = 2.25 \times 10^{-6}/[1+(\zs/3.14)^{-3.42}]$ ($0.6<\zs<7$). The amplitude is about 50% higher for an effective source diameter of $0.5\arcsec$. The optimal lens redshift for giant arcs with the length-to-width ratio ($L/W$) larger than 10 increases from 0.3 for $\zs=1$, to 0.5 for $\zs=2$, and to 0.7-0.8 for $\zs>3$. The optical depth is sensitive to the source redshift, in qualitative agreement with Wambsganss et al. (2004). However, our overall optical depth appears to be only $\sim$ 10% to 70% of those from previous studies. The differences can be mostly explained by different power spectrum normalizations ($\sigma_8$) used and different ways of determining the $L/W$ ratio. Finite source size and ellipticity have modest effects on the optical depth. We also found that the number of highly magnified (with magnification $|\mu|>10$) and ``undistorted'' images (with $L/W<3$) is comparable to the number of giant arcs with $|\mu|>10$ and $L/W>10$. We conclude that our predicted rate of giant arcs may be lower than the observed rate, although the precise `discrepancy' is still unclear due to uncertainties both in theory and observations.

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

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.

Core collapse supernovae are dominated by energy transport from neutrinos. Therefore, some supernova properties could depend on symetries and features of the standard model weak interactions. The cross section for neutrino capture is larger than that for antineutrino capture by one term of order the neutrino energy over the nucleon mass. This reduces the ratio of neutrons to protons in the $\nu$-driven wind above a protoneutron star by approximately 20 % and may significantly hinder r-process nucleosynthesis.

The well known charge conjugation violating interactions in the Standard Model increase neutrino- and decrease anti-neutrino- nucleon cross sections. This impacts neutrino transport in core collapse supernovae through "recoil" corrections of order the neutrino energy $k$ over the nucleon mass $M$. All $k/M$ corrections to neutrino transport deep inside a protoneutron star are calculated from angular integrals of the Boltzmann equation. We find these corrections significantly modify neutrino currents at high temperatures. This produces a large mu and tau number for the protoneutron star and can change the ratio of neutrons to protons. In addition, the relative size of neutrino mean free paths changes. At high temperatures, the electron anti-neutrino mean free path becomes {\it longer} than that for mu or tau neutrinos.

Core collapse supernovae are dominated by energy transport from neutrinos. Therefore, some supernova properties could depend on symetries and features of the standard model weak interactions. The cross section for neutrino capture is larger than that for antineutrino capture by one term of order the neutrino energy over the nucleon mass. This reduces the ratio of neutrons to protons in the $\nu$-driven wind above a protoneutron star by approximately 20 % and may significantly hinder r-process nucleosynthesis.

The neutrino-nucleon cross section is slightly larger than that for $\bar\nu-N$. Therefore, $\bar\nu$ will escape more quickly from core collapse supernovae leaving the star $\nu$ rich. A diffusion formalism is used to calculate the time evolution of the mu and tau lepton number densities. These quickly reach steady state equilibria. We estimate that a protoneutron star with a maximum temperature near 50 MeV will contain over 50 % more $\nu_\mu$ and $\nu_\tau$ than $\bar\nu_\mu$ and $\bar\nu_\tau$. Supernovae may be the only known systems with large mu and or tau lepton numbers.

The recent past witnesses the growing interdependence between the physics of hadrons, the physics of relativistic heavy-ion collisions, and the physics of compact objects in astrophysics. A notable example is the kaon which plays special roles in all the three fields. In this paper, we first review the various theoretical investigations of kaon properties in nuclear medium, focusing on possible uncertainties in each model. We then present a detailed transport model study of kaon production in heavy-ion collisions at SIS energies. We shall discuss especially the elementary kaon and antikaon production cross sections in hadron-hadron interactions, that represent one of the most serious uncertainties in the transport model study of particle production in heavy-ion collisions. The main purpose of such a study is to constrain kaon in-medium properties from the heavy-ion data. This can provide useful guidances for the development of theoretical models of the kaon in medium. In the last part of the paper, we apply the kaon in-medium properties extracted from heavy-ion data to the study of neutron star properties. Based on a conventional equation of state of nuclear matter that can be considered as one of the best constrained by available experimental data on finite nuclei, we find that the maximum mass of neutron stars is about 2$M_\odot$, which is reduced to about 1.5$M_\odot$ once kaon condensation as constrained by heavy-ion data is introduced.

Supernovae provide a unique opportunity for large scale parity violation because they are dominated by neutrinos. We calculate the parity violating asymmetry A of neutrino emission in a strong magnetic field. We assume the neutrinos elastically scatter many times from slightly polarized neutrons. Because of multple interactions, A grows with the optical thickness of the proto-neutron star and may be much larger than previous estimates. As a result, the neutron star could recoil at a significant velocity.