Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

By using Data Analysis Pipeline (DAP) products of Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), which are publicly available from the SDSS Data Release 15, we analyze the local properties at the SN explosion sites and global properties of different types of SN host galaxies to explore the explosion environments of different types of SNe. In our sample, there are 67 SN host galaxies in the field of view of MaNGA, including 32 Type Ia, 29 CCSNe, 1 super-luminous SN (SLSN), 1 Type I and 4 unclassified type of SNe, with which we can perform the K-S test for analysis and derive statistically robust results. Due to the limited sample size, we couldn't remove the mass dependence in this work, which is likely the true driver of the trends for the properties presented in this work. The global star formation rate (SFR) and EW(H$\alpha$) for SN Ia hosts is slightly lower than that for CCSN hosts on average. SN Ia host galaxies are $\sim$ 0.3 dex more massive than CCSN hosts, which implies that the number ratio of CCSNe to Type Ia SNe will decrease with the increasing of stellar mass of host galaxies. The stellar population age of SN Ia host galaxies is older than that of CCSN hosts on average. There is no significant difference between different types of SN hosts for some properties, including local SFR density ($\Sigma \rm SFR$), local and global gas-phase oxygen abundance. For most galaxies in our sample, the global gas-phase oxygen abundance estimated from the integrated spectra of SN hosts can represent the local gas-phase oxygen abundance at the SN explosion sites with small bias.

Although the physical origin of prompt emission in gamma-ray bursts (GRBs) remains inconclusive, previous studies have considered the synchrotron radiation of relativistic electrons as a promising mechanism. These works usually adopted a invariable injection rate of electrons ($Q$) which may be discordant with that in a Poynting-flux dominated jet. In a Poynting-flux dominated jet (e.g., ICMART model, Zhang & Yan 2011), the number of magnetic reconnections occurred simultaneously may grow rapidly with time and results in an increase of $Q$ with time. This paper is dedicated to study the synchrotron radiation spectrum in this scenario. It is found that the radiation spectrum would obviously get harder if an increasing $Q$ is adopted and a Band-like radiation spectrum can be obtained if the increase of $Q$ is fast enough. The latter is related to the fact that a bump-shape rather than a power-law spectrum appears in the low-energy regime of the obtained electron spectrum. This effect can strongly harden the low-energy radiation spectrum. It indicates that an increasing $Q$ can help to alleviate the "fast-cooling problem" of synchrotron radiation for GRBs. Our studies also reveal that a Poynting-flux dominated jet with a large emission radius, a small length of the magnetic reconnection region, or a low-minimum energy of injected electron would prefer to form a Band-like radiation spectrum. We suggest that the Band spectrum found in GRBs may be the synchrotron emission of the electrons with a bump-shape distribution in its low-energy regime.

The goal of Project GAUSS is to return samples from the dwarf planet Ceres. Ceres is the most accessible ocean world candidate and the largest reservoir of water in the inner solar system. It shows active cryovolcanism and hydrothermal activities in recent history that resulted in minerals not found in any other planets to date except for Earth's upper crust. The possible occurrence of recent subsurface ocean on Ceres and the complex geochemistry suggest possible past habitability and even the potential for ongoing habitability. Aiming to answer a broad spectrum of questions about the origin and evolution of Ceres and its potential habitability, GAUSS will return samples from this possible ocean world for the first time. The project will address the following top-level scientific questions: 1) What is the origin of Ceres and the origin and transfer of water and other volatiles in the inner solar system? 2) What are the physical properties and internal structure of Ceres? What do they tell us about the evolutionary and aqueous alteration history of icy dwarf planets? 3) What are the astrobiological implications of Ceres? Was it habitable in the past and is it still today? 4) What are the mineralogical connections between Ceres and our current collections of primitive meteorites? GAUSS will first perform a high-resolution global remote sensing investigation, characterizing the geophysical and geochemical properties of Ceres. Candidate sampling sites will then be identified, and observation campaigns will be run for an in-depth assessment of the candidate sites. Once the sampling site is selected, a lander will be deployed on the surface to collect samples and return them to Earth in cryogenic conditions that preserves the volatile and organic composition as well as the original physical status as much as possible.

Metals are thought to have profound effects on the internal structures of molecular clouds in which stars are born. The absence of metals is expected to prevent gas from efficient cooling and fragmentation in theory. However, this effect has not yet been observed at low metallicity environments, such as in the early Universe and local dwarf galaxies, because of the lack of high spatial resolution maps of gas. We carried out ALMA observations of the carbon monoxide (CO) J=2-1 emission line at 1.4-parsec resolutions of a molecular cloud in DDO 70 at 7% solar metallicity, the most metal-poor galaxy currently known with a CO detection. In total, five clumps have been identified and they are found to follow more or less the Larson's law. Since the CO emission exists in regions with visual extinction A_V around 1.0, we converted this A_V to the gas mass surface density using a gas-to-dust ratio of 4,594+-2,848 for DDO 70. We found that the CO clumps in DDO 70 exhibit significantly larger (on average four times) sizes than those at the same gas mass surface densities in massive star-formation regions of the Milky Way. The existence of such large clumps appears to be consistent with theoretical expectations that gas fragmentation in low metallicity clouds is suppressed. While our observation is only for one cloud in the galaxy, if it is representative, the above result implies suppressed gas fragmentation during the cloud collapse and star formation in the early Universe.

Many asteroids in the main belt are continuously pushed by Yarkovsky effect into regions of different mean motion resonances (MMRs) and then ejected out. They are considered as the principal source of near-Earth objects. We investigate in this paper the effects of the 7/3 MMR with Jupiter (J7/3 MMR) on the transportation of asteroids from Koronis and Eos families that reside respectively on the inner and outer side of the resonance. The fraction of asteroids that make successful crossing through the resonance and the escaping rate from the resonance are found to depend on the Yarkovsky drifting rate, the initial inclination and the migrating direction. The excitation of eccentricity and inclination due to the combined influence from both the resonance and Yarkovsky effect is discussed. Only the eccentricity can be pumped up considerably, and it is attributed mainly to the resonance. In the observational data, family members are also found in the resonance and on the opposite side of the resonance with respect to the corresponding family centre. The existence of these family members is explained using our results of numerical simulations. Finally, the replenishment of asteroids in the J7/3 MMR and the transportation of asteroids by it are discussed.

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $\Delta a$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $\Delta a$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $\Delta a$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $\omega_{\rm rot}$ of asteroid depends on its size $R$ by $\omega_{\rm rot}\propto R^{-1}$.

The high-inclination population in the 4:7 mean motion resonance (MMR) with Neptune has also substantial eccentricities ($e\gtrsim0.1$), with more inclined objects tending to occupy more eccentric orbits. For this high-order resonance, there are two different resonant modes. The principal one is the eccentricity-type mode, and we find that libration is permissible for orbits with $e\ge e_c^0$, where the critical eccentricity $e_c^0$ increases as a function of increasing inclination $i$. Correspondingly, we introduce a limiting curve $e_c^0(i)$, which puts constraints on the $(e, i)$ distribution of possible 4:7 resonators. We then perform numerical simulations on the sweep-up capture and long-term stability of the 4:7 MMR, and the results show that the simulated resonators are well-constrained by this theoretical limiting curve. The other 4:7 resonant mode is the mixed-$(e, i)$-type, and we show that stable resonators should exist at $i\gtrsim20^{\circ}$. We predict that the intrinsic number of these mixed-$(e, i)$-type resonators may provide a new clue into the Solar system's evolution, but, so far, only one real object has been observed resonating in this mode.

We aim to locate the stability region for Uranus Trojans (UT hereafter) and find out the dynamical mechanisms responsible for the structures in the phase space. Using the spectral number as the stability indicator, we construct the dynamical maps on the (a0, i0) plane. The proper frequencies of UTs are determined precisely so that we can depict the resonance web via a semi-analytical method. Two main stability regions are found, one each for the low-inclination (0-14deg) and high-inclination regime (32-59deg). There is also an instability strip in each of them, at 9deg and 51deg respectively. All stability regions are in the tadpole regime and no stable horseshoe orbits exist for UTs. The lack of moderate-inclined UTs is caused by the nu5 and nu7 secular resonances. The fine structures in the dynamical maps are shaped by high-degree secular resonances and secondary resonances. During the planetary migration, about 36.3% and 0.4% of the pre-formed orbits survive the fast and slow migrations (with migrating time scales of 1 and 10Myr) respectively, most of which are in high inclination. Since the low-inclined UTs are more likely to survive the age of the solar system, they make up 77% of all such long-life orbits by the end of the migration, making a total fraction up to 4.06E-3 and 9.07E-5 of the original population for the fast and slow migrations, respectively. About 3.81% UTs are able to survive the age of the solar system, among which 95.5% are on low-inclined orbits with i0<7.5deg. However, the depletion of the planetary migration seems to prevent a large fraction of such orbits, especially for the slow migration model.

We aim to locate the stability region for Uranus Trojans (UT hereafter) and find out the dynamical mechanisms responsible for the structures in the phase space. Using the spectral number as the stability indicator, we construct the dynamical maps on the (a0, i0) plane. The proper frequencies of UTs are determined precisely so that we can depict the resonance web via a semi-analytical method. Two main stability regions are found, one each for the low-inclination (0-14deg) and high-inclination regime (32-59deg). There is also an instability strip in each of them, at 9deg and 51deg respectively. All stability regions are in the tadpole regime and no stable horseshoe orbits exist for UTs. The lack of moderate-inclined UTs is caused by the nu5 and nu7 secular resonances. The fine structures in the dynamical maps are shaped by high-degree secular resonances and secondary resonances. During the planetary migration, about 36.3% and 0.4% of the pre-formed orbits survive the fast and slow migrations (with migrating time scales of 1 and 10Myr) respectively, most of which are in high inclination. Since the low-inclined UTs are more likely to survive the age of the solar system, they make up 77% of all such long-life orbits by the end of the migration, making a total fraction up to 4.06E-3 and 9.07E-5 of the original population for the fast and slow migrations, respectively. About 3.81% UTs are able to survive the age of the solar system, among which 95.5% are on low-inclined orbits with i0<7.5deg. However, the depletion of the planetary migration seems to prevent a large fraction of such orbits, especially for the slow migration model.

In this work, using the Gaussian Process, we explore the potentiality of future gravitational wave (GW) measurement to probe cosmic opacity through comparing its opacity-free luminosity distance (LD) with the opacity-dependent one from type Ia supernovae (SNIa). GW data points are simulated from the third generation Einstein Telescope, and SNIa data are taken from the Joint Light Analysis (JLA) or Pantheon compilation. The advantages of using Gaussian Process are that one may match SNIa data with GW data at the same redshift and use all available data to probe cosmic opacity. We obtain that the error bar of the constraint on cosmic opacity can be reduced to $\sigma_{\epsilon}\sim 0.011$ and $0.006$ at $1\sigma$ confidence level (CL) for JLA and Pantheon respectively in a cosmological-independent way. Thus, the future GW measurements can give competitive results on the cosmic opacity test. Furthermore, we propose a method to probe the spatial homogeneity of the cosmic transparency through comparing the reconstructed LD from the mock GW with the reconstructed one from SNIa data in a flat $\Lambda$CDM with the Gaussian Process. The result shows that a transparent universe is favored at $1\sigma$ CL, although the best-fit value of cosmic opacity is redshift-dependent.

We aim to understand the properties at the locations of supernova (SN) explosion in their host galaxies and compare with the global properties of the host galaxies. We use the integral field spectrograph (IFS) of Mapping Nearby Galaxies (MaNGA) at Apache Point Observatory (APO) to get the 2D maps of the parameter properties for eleven SN host galaxies. The sample galaxies are analyzed one by one in details on their properties of velocity field, star formation rate, oxygen abundance and stellar mass etc. This sample of SN host galaxies have redshifts around $z$ $\sim$ 0.03, which is higher than those of the previous related works. The higher redshift distribution allows us to obtain the properties of more distant SN host galaxies. Metallicity (gas-phase oxygen abundance) estimated from integrated spectra could represent the local metallicity at SN explosion sites with small bias. All the host galaxies in our sample are metal-rich galaxies (12+log(O/H) $>$ 8.5) except for NGC 6387, which means supernovae (SNe) may be more inclined to explode in rich-metallicity galaxies. There is a positive relation between global gas-phase oxygen abundance and the stellar mass of host galaxies. We also try to compare the differences of the host galaxies between SN Ia and SN II. In our sample, both SNe Ia and SNe II could explode in normal galaxies, while SNe II also could explode in an interactive or merger system, which has star formation in the galaxy.

We aim to determine the relative angle between the total angular momentum of the minor planets and that of the Sun-planets system, and to improve the orientation of the invariable plane of the solar system. By utilizing physical parameters available in public domain archives, we assigned reasonable masses to 718041 minor planets throughout the solar system, including near-Earth objects, main belt asteroids, Jupiter trojans, trans-Neptunian objects, scattered-disk objects, and centaurs. Then we combined the orbital data to calibrate the angular momenta of these small bodies, and evaluated the specific contribution of the massive dwarf planets. The effects of uncertainties on the mass determination and the observational incompleteness were also estimated. We determine the total angular momentum of the known minor planets to be $1.7817\times10^{46}$ g $\cdot$ cm$^2$ $\cdot$ s$^{-1}$. The relative angle $\alpha$ between this vector and the total angular momentum of the Sun-planets system is calculated to be 14.74 deg.. By excluding the dwarf planets Eris, Pluto, and Haumea, which have peculiar angular momentum directions, $\alpha$ drops sharply to 1.76 deg.; a similar result applies to each individual minor planet group. This suggests that, without these three most massive bodies, the plane perpendicular to the total angular momentum of the minor planets would be close to the invariable plane of the solar system. On the other hand, the inclusion of Eris, Haumea, and Makemake can produce a difference of 1254 mas in the inclination of the invariable plane, which is much larger than the difference of 9 mas induced by Ceres, Vesta, and Pallas as found previously. By taking into account the angular momentum contributions from all minor planets, including unseen ones, the orientation improvement of the invariable plane is larger than 1000 mas in inclination with a $1\sigma$ error of 50-140 mas.

The presented work investigates the possible formation of terrestrial planets in the habitable zone (HZ) of the exoplanetary system HD 141399. In this system the HZ is located approximately between the planets c (a = 0.7 au) and d (a = 2.1 au). Extensive numerical integrations of the equations of motion in the pure Newtonian framework of small bodies with different initial conditions in the HZ are performed. Our investigations included several steps starting with 500 massless bodies distributed between planets c and d in order to model the development of the disk of small bodies. It turns out that after some 10^6 years a belt-like structure analogue to the main belt inside Jupiter in our Solar System appears. We then proceed with giving the small bodies masses (Moon-mass) and take into account the gravitational interaction between these planetesimal-like objects. The growing of the objects - with certain percentage of water - due to collisions is computed in order to look for the formation of terrestrial planets. We observe that planets form in regions connected to mean motion resonances (MMR). So far there is no observational evidence of terrestrial planets in the system of HD 141399 but from our results we can conclude that the formation of terrestrial planets - even with an appropriate amount of water necessary for being habitable - in the HZ would have been possible.

The presented work investigates the possible formation of terrestrial planets in the habitable zone (HZ) of the exoplanetary system HD 141399. In this system the HZ is located approximately between the planets c (a = 0.7 au) and d (a = 2.1 au). Extensive numerical integrations of the equations of motion in the pure Newtonian framework of small bodies with different initial conditions in the HZ are performed. Our investigations included several steps starting with 500 massless bodies distributed between planets c and d in order to model the development of the disk of small bodies. It turns out that after some 10^6 years a belt-like structure analogue to the main belt inside Jupiter in our Solar System appears. We then proceed with giving the small bodies masses (Moon-mass) and take into account the gravitational interaction between these planetesimal-like objects. The growing of the objects - with certain percentage of water - due to collisions is computed in order to look for the formation of terrestrial planets. We observe that planets form in regions connected to mean motion resonances (MMR). So far there is no observational evidence of terrestrial planets in the system of HD 141399 but from our results we can conclude that the formation of terrestrial planets - even with an appropriate amount of water necessary for being habitable - in the HZ would have been possible.

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

To study the disc central surface brightness ($\mu_0$) distribution in optical and near-infrared bands, we select 708 disc-dominated galaxies within a fixed distance of 57 Mpc from SDSS DR7 and UKIDSS DR10. Then we fit $\mu_0$ distribution by using single and double Gaussian profiles with an optimal bin size for the final sample of 538 galaxies in optical $griz$ bands and near-infrared $YJHK$ bands. Among the 8 bands, we find that $\mu_{0}$ distribution in optical bands can not be much better fitted with double Gaussian profiles. However, for all the near-infrared bands, the evidence of being better fitted by using double Gaussian profiles is positive. Especially for $K$ band, the evidence of a double Gaussian profile being better than a single Gaussian profile for $\mu_{0}$ distribution is very strong, the reliability of which can be approved by 1000 times test for our sample. No dust extinction correction is applied. The difference of $\mu_{0}$ distribution between optical and near-infrared bands could be caused by the effect of dust extinction in optical bands. Due to the sample selection criteria, our sample is not absolutely complete. However, the sample incompleteness does not change the double Gaussian distribution of $\mu_{0}$ in $K$ band. Furthermore, we discuss some possible reasons for the fitting results of $\mu_{0}$ distribution in $K$ band. Conclusively, the double Gaussian distribution of $\mu_{0}$ in $K$ band for our sample may depend on bulge-to-disk ratio, color and disk scalelength, rather than the inclination of sample galaxies, bin size and statistical fluctuations.

The only discovery of Earth Trojan 2010 TK$_7$ and the subsequent launch of OSIRIS-REx motive us to investigate the stability around the triangular Lagrange points $L_4$ and $L_5$ of the Earth. In this paper we present detailed dynamical maps on the $(a_0,i_0)$ plane with the spectral number (SN) indicating the stability. Two main stability regions, separated by a chaotic region arising from the $\nu_3$ and $\nu_4$ secular resonances, are found at low ($i_0\leq 15^\circ$) and moderate ($24^\circ\leq {i_0}\leq 37^\circ$) inclinations respectively. The most stable orbits reside below $i_0=10^\circ$ and they can survive the age of the Solar System. The nodal secular resonance $\nu_{13}$ could vary the inclinations from $0^\circ$ to $\sim 10^\circ$ according to their initial values while $\nu_{14}$ could pump up the inclinations to $\sim 20^\circ$ and upwards. The fine structures in the dynamical maps are related to higher-degree secular resonances, of which different types dominate different areas. The dynamical behaviour of the tadpole and horseshoe orbits, reflected in their secular precession, show great differences in the frequency space. The secular resonances involving the tadpole orbits are more sensitive to the frequency drift of the inner planets, thus the instabilities could sweep across the phase space, leading to the clearance of tadpole orbits. We are more likely to find terrestrial companions on horseshoe orbits. The Yarkovsky effect could destabilize Earth Trojans in varying degrees. We numerically obtain the formula describing the stabilities affected by the Yarkovsky effect and find the asymmetry between the prograde and retrograde rotating Earth Trojans. The existence of small primordial Earth Trojans that avoid being detected but survive the Yarkovsky effect for 4.5\,Gyr is substantially ruled out.

The only discovery of Earth Trojan 2010 TK$_7$ and the subsequent launch of OSIRIS-REx motive us to investigate the stability around the triangular Lagrange points $L_4$ and $L_5$ of the Earth. In this paper we present detailed dynamical maps on the $(a_0,i_0)$ plane with the spectral number (SN) indicating the stability. Two main stability regions, separated by a chaotic region arising from the $\nu_3$ and $\nu_4$ secular resonances, are found at low ($i_0\leq 15^\circ$) and moderate ($24^\circ\leq {i_0}\leq 37^\circ$) inclinations respectively. The most stable orbits reside below $i_0=10^\circ$ and they can survive the age of the Solar System. The nodal secular resonance $\nu_{13}$ could vary the inclinations from $0^\circ$ to $\sim 10^\circ$ according to their initial values while $\nu_{14}$ could pump up the inclinations to $\sim 20^\circ$ and upwards. The fine structures in the dynamical maps are related to higher-degree secular resonances, of which different types dominate different areas. The dynamical behaviour of the tadpole and horseshoe orbits, reflected in their secular precession, show great differences in the frequency space. The secular resonances involving the tadpole orbits are more sensitive to the frequency drift of the inner planets, thus the instabilities could sweep across the phase space, leading to the clearance of tadpole orbits. We are more likely to find terrestrial companions on horseshoe orbits. The Yarkovsky effect could destabilize Earth Trojans in varying degrees. We numerically obtain the formula describing the stabilities affected by the Yarkovsky effect and find the asymmetry between the prograde and retrograde rotating Earth Trojans. The existence of small primordial Earth Trojans that avoid being detected but survive the Yarkovsky effect for 4.5\,Gyr is substantially ruled out.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

Close encounters (CEs) between celestial objects may exert significant influence on their orbits. The influence will be even enhanced when two groups of celestial objects are confined in stable orbital configurations, e.g. in adjacent mean motion resonances (MMRs). Plutinos and Neptune Trojans, trapped in the 2:3 and 1:1 MMRs with Neptune respectively, are such examples. %Meanwhile, many objects among these two groups have peculiar orbits, seemingly as the vestige of CEs. As the first part of our investigation, this paper provides a detailed description of CEs between Plutinos and Trojans and their potential influences on the Trojans' orbits. Statistical analyses of CE data from numerical simulations reveal the randomness lying in the CEs between the two planetesimals. The closest positions of CEs distribute symmetrically inside the given CE region and no particular bias is found between the positive and negative effects on the orbital elements of Trojans. Based on the Gaussian approximation on the distribution of the velocity orientation of Plutino, and the integral derivatives of Gaussian perturbation equations, a theoretical method is built to estimate the CE effects. To further verify the randomness of CEs, a Monte Carlo approach is applied, and it generates distribution features consistent with the numerical results. In summary, CEs brought by realistic Plutinos exert impartial effects and tiny total influence on the orbital elements of Trojans. However, driven by the random walk mechanism, tiny effects may accumulate to a prominent variation given sufficient CEs, which will be discussed in the accompanying paper.

We use high-resolution continuum images obtained at 870microns with the Atacama Large Millimeter Array (ALMA) to probe the surface density of star-formation in z~2 galaxies and study the different physical properties between galaxies within and above the star-formation main sequence of galaxies. This sample of eight star-forming galaxies at z~2 selected among the most massive Herschel galaxies in the GOODS-South field is supplemented with eleven galaxies from the public data of the 1.3 mm survey of the Hubble Ultra-Deep Field. ALMA reveals systematically dense concentrations of dusty star-formation close to the center of the stellar component of the galaxies. We identify two different starburst regimes: (i) the classical population of starbursts located above the SFR-M* main sequence, with enhanced gas fractions and short depletion times and (ii) a sub-population of galaxies located within the scatter of the main sequence that experience compact star formation with depletion timescales typical of starbursts of ~150 Myr. In both starburst populations, the far infrared and UV are distributed in distinct regions and dust-corrected star formation rates estimated using UV-optical-NIR data alone underestimate the total star formation rate. Starbursts hidden in the main sequence show instead the lowest gas fractions of our sample and could represent the last stage of star-formation before they become passive. Being Herschel-selected, these main sequence galaxies are located in the high-mass end of the main sequence, hence we do not know whether these "starbursts hidden in the main sequence" also exist below 10^11 Msun. Active galactic nuclei are found to be ubiquitous in these compact starbursts, suggesting that the triggering mechanism also feeds the central black hole or that the active nucleus triggers star formation.

We analyze in this paper the existence of the "out-of-plane" equilibrium points in the restricted three-body problem with oblateness. From the series expansion of the potential function of an oblate asteroid, we show analytically all equilibrium points locate on the orbital plane of primaries and how artificial equilibrium points may arise due to an inappropriate application of the potential function. Using the closed form of the potential of a triaxial ellipsoid, we analytically demonstrate that the gravitational acceleration in $z$-direction is always pointing toward the equatorial plane, thus it could not be balanced out at any value of $z\neq 0$ and the out-of-plane equilibrium points cannot exist. The out-of-plane equilibrium points appear only when additional acceleration other than the gravitation from primaries is taken into account. We suggest that special attention must be paid to the application of the spherical harmonics expansion of potential to find the equilibrium points, especially when these points may be very close to the celestial body.

We revisit the proposed extended Schmidt law (Shi et al. 2011) which points that the star formation efficiency in galaxies depends on the stellar mass surface density, by investigating spatially-resolved star formation rates (SFRs), gas masses and stellar masses of star formation regions in a vast range of galactic environments, from the outer disks of dwarf galaxies to spiral disks and to merging galaxies as well as individual molecular clouds in M33. We find that these regions are distributed in a tight power-law as Sigma_SFR ~(Sigma_star^0.5 Sigma_gas )^1.09, which is also valid for the integrated measurements of disk and merging galaxies at high-z. Interestingly, we show that star formation regions in the outer disks of dwarf galaxies with Sigma_SFR down to 10^(-5) Msun/yr/kpc^2, which are outliers of both Kennicutt-Schmidt and Silk-Elmegreen law, also follow the extended Schmidt law. Other outliers in the Kennicutt-Schmidt law, such as extremely-metal poor star-formation regions, also show significantly reduced deviations from the extended Schmidt law. These results suggest an important role for existing stars in helping to regulate star formation through the effect of their gravity on the mid-plane pressure in a wide range of galactic environments.

We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the dependence on the metallicity of the IR-radio relationship at 24 um, 70 um, 100 um and 160 um bands as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on average lower qIR than metal-rich ones with larger offsets at longer IR wavelengths, from -0.06 dex in q24um to -0.6 dex in q160um. The qIR of all galaxies as a whole at 160 um show positive trends with the metallicity and IR-to-FUV ratio, and negative trends with the IR color, while those at lower IR wavelengths show weaker correlations. We proposed a mechanism that invokes combined effects of low obscured-SFR/total-SFR fraction and warm dust temperature at low metallicity to interpret the above behavior of qIR, with the former reducing the IR radiation and the latter further reducing the IR emission at longer IR wavelength. Other mechanisms that are related to the radio emission including the enhanced magnetic field strength and increased thermal radio contribution are unable to reconcile the IR-wavelength-dependent differences of qIR between metal-poor and metal- rich galaxies. In contrast to qIR, the mean total-SFR/radio ratio of metal-poor galaxies is the same as the metal-rich one, indicating the 1.4 GHz radio emission is still an effective tracer of SFRs at low metallicity.

We investigate the energy sources of random turbulent motions of ionised gas from H$\alpha$ emission in eight local star-forming galaxies from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. These galaxies satisfy strict pure star-forming selection criteria to avoid contamination from active galactic nuclei (AGN) or strong shocks/outflows. Using the relatively high spatial and spectral resolution of SAMI, we find that -- on sub-kpc scales our galaxies display a flat distribution of ionised gas velocity dispersion as a function of star formation rate (SFR) surface density. A major fraction of our SAMI galaxies shows higher velocity dispersion than predictions by feedback-driven models, especially at the low SFR surface density end. Our results suggest that additional sources beyond star formation feedback contribute to driving random motions of the interstellar medium (ISM) in star-forming galaxies. We speculate that gravity, galactic shear, and/or magnetorotational instability (MRI) may be additional driving sources of turbulence in these galaxies.

In a previous paper, we proposed a heterotic dark energy model, called $\Lambda$HDE, in which dark energy is composed of two components: cosmological constant (CC) and holographic dark energy (HDE). The aim of this work is to give a more comprehensive and systematic investigation on the cosmological implications of the $\Lambda$HDE model. Firstly, we make use of the current observations to constrain the $\Lambda$HDE model, and compare its cosmology-fit results with the results of the $\Lambda$CDM and the HDE model. Then, by combining a qualitative theoretical analysis with a quantitative numerical study, we discuss the impact of considering curvature on the cosmic evolutions of fractional HDE density $\Omega_{hde}$ and fractional CC density $\Omega_{\Lambda}$, as well as on the ultimate cosmic fate. Finally, we explore the effects of adopting different types of observational data. We find that: (1) the current observational data cannot distinguish the $\Lambda$HDE model from the $\Lambda$CDM and the HDE model; this indicates that DE may contain multiple components. (2) the asymptotic solution of $\Omega_{hde}$ and the corresponding cosmic fate in a flat universe can be extended to the case of a non-flat universe; moreover, compared with the case of a flat universe, considering curvature will make HDE closer to a phantom dark energy. (3) compared with JLA dataset, SNLS3 data more favor a phantom type HDE; in contrast, using other types of observational data have no significant impact on the cosmic evolutions of the $\Lambda$HDE model.

In order to constrain the evolutionary state of the red supergiant Betelgeuse, we have produced a suite of models with ZAMS masses from 15 to 25 Msun in intervals of 1 Msun including the effects of rotation. The models were computed with the stellar evolutionary code MESA. For non-rotating models we find results that are similar to other work. It is somewhat difficult to find models that agree within 1 sigma of the observed values of R, Teff and L, but modestly easy within 3 sigma uncertainty. Incorporating the nominal observed rotational velocity, ~15 km/s, yields significantly different, and challenging, constraints. This velocity constraint is only matched when the models first approach the base of the red supergiant branch (RSB), having crossed the Hertzsprung gap, but not yet having ascended the RSB and most violate even generous error bars on R, Teff and L. Models at the tip of the RSB typically rotate at only ~0.1 km/s, independent of any reasonable choice of initial rotation. We discuss the possible uncertainties in our modeling and the observations, including the distance to Betelgeuse, the rotation velocity, and model parameters. We summarize various options to account for the rotational velocity and suggest that one possibility is that Betelgeuse merged with a companion star of about 1 Msun as it ascended the RSB, in the process producing the ring structure observed at about 7' away. A past coalescence would complicate attempts to understand the evolutionary history and future of Betelgeuse.

We present a study of nuclear activities in nearby S0 galaxies. After cross-matching the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) with the Third Reference Catalog of Bright Galaxies (RC3) and visually checking the SDSS images, we derive a sample of 583 S0 galaxies with the central spectrophotometric information. In order to separate nebular emission lines from the underlying stellar contribution, we fit the stellar population model to the SDSS spectra of these S0 galaxies. According to the BPT diagram, we find that $8\%$ of S0 galaxies show central star-forming activity, while the fractions of Seyfert, Composite and LINERs are 2\%, 8\%, and 21.4\%, respectively. We also find that star-forming S0s have the lowest stellar masses, over one magnitude lower than the others, and that the active S0s locate mainly in the sparse environment, while the normal S0s in the dense environment, which might suggest that the environment plays an important role in quenching star formation and/or AGN activity in S0 galaxies. By performing bulge-disk decomposition of 45 star-forming S0s in {\em g}- and {\em r}-bands with the 2D fitting software \textsc{Galfit}, as well as exploiting the catalog of 2D photometric decompositions of \citet{Meert2015}, we find that the bulges of $\sim$ 1/3 star-forming S0 galaxies (16/45) are bluer than their disks, while for other types of S0s, the bulge and disk components show similar color distributions. Besides, the S\'{e}rsic index of most star-forming S0s bulges is less than 2, while for normal S0s, it is between 2 and 6.

In this work, by applying the redshift tomography method to Joint Light-curve Analysis (JLA) supernova sample, we explore the possible redshift-dependence of stretch-luminosity parameter $\alpha$ and color-luminosity parameter $\beta$. The basic idea is to divide the JLA sample into different redshift bins, assuming that $\alpha$ and $\beta$ are piecewise constants. Then, by constraining the $\Lambda$CDM model, we check the consistency of cosmology-fit results given by the SN sample of each redshift bin. We also adopt the same technique to explore the possible evolution of $\beta$ in various subsamples of JLA. Using the full JLA data, we find that $\alpha$ is always consistent with a constant. In contrast, at high redshift $\beta$ has a significant trend of decreasing, at $\sim 3.5\sigma$ confidence level (CL). Moreover, we find that low-$z$ subsample favors a constant $\beta$; in contrast, SDSS and SNLS subsamples favor a decreasing $\beta$ at 2$\sigma$ and $3.3\sigma$ CL, respectively. Besides, by using a binned parameterization of $\beta$, we study the impacts of $\beta$'s evolution on parameter estimation. We find that compared with a constant $\beta$, a varying $\beta$ yields a larger best-fit value of fractional matter density $\Omega_{m0}$, which slightly deviates from the best-fit result given by other cosmological observations. However, for both the varying $\beta$ and the constant $\beta$ cases, the $1\sigma$ regions of $\Omega_{m0}$ are still consistent with the result given by other observations.

We present infrared (IR) spectral energy distributions (SEDs) of individual star-forming regions in four extremely metal poor (EMP) galaxies with metallicity Z around Zsun/10 as observed by the Herschel Space Observatory. With the good wavelength coverage of the SED, it is found that these EMP star-forming regions show distinct SED shapes as compared to those of grand design Spirals and higher metallicity dwarfs: they have on average much higher f70um/f160um ratios at a given f160um/f250um ratio; single modified black-body (MBB) fittings to the SED at \lambda >= 100 um still reveal higher dust temperatures and lower emissivity indices compared to that of Spirals, while two MBB fittings to the full SED with a fixed emissivity index (beta = 2) show that even at 100 um about half of the emission comes from warm (50 K) dust, in contrast to the cold (~20 K) dust component. Our spatially resolved images further reveal that the far-IR colors including f70um/f160um, f160um/f250um and f250um/f350um are all related to the surface densities of young stars as traced by far-UV, 24 um and SFRs, but not to the stellar mass surface densities. This suggests that the dust emitting at wavelengths from 70 um to 350 um is primarily heated by radiation from young stars.

Recently, a new dark energy model called $\Lambda$HDE was proposed. In this model, dark energy consists of two parts: cosmological constant $\Lambda$ and holographic dark energy (HDE). Two key parameters of this model are the fractional density of cosmological constant $\Omega_{\Lambda0}$, and the dimensionless HDE parameter $c$. Since these two parameters determine the dynamical properties of DE and the destiny of universe, it is important to study the impacts of different values of $\Omega_{\Lambda0}$ and $c$ on the $\Lambda$HDE model. In this paper, we apply various DE diagnostic tools to diagnose $\Lambda$HDE models with different values of $\Omega_{\Lambda0}$ and $c$; these tools include statefinder hierarchy \{$S_3^{(1)}, S_4^{(1)}$\}, fractional growth parameter $\epsilon$, and composite null diagnostic (CND), which is a combination of \{$S_3^{(1)}, S_4^{(1)}$\} and $\epsilon$. We find that: (1) adopting different values of $\Omega_{\Lambda0}$ only has quantitative impacts on the evolution of the $\Lambda$HDE model, while adopting different $c$ has qualitative impacts; (2) compared with $S_3^{(1)}$, $S_4^{(1)}$ can give larger differences among the cosmic evolutions of the $\Lambda$HDE model associated with different $\Omega_{\Lambda0}$ or different $c$; (3) compared with the case of using a single diagnostic, adopting a CND pair has much stronger ability to diagnose the $\Lambda$HDE model.

The Tsinghua University-National Astronomical Observatories of China (NAOC) Transient Survey (TNTS) is an automatic survey for a systematic exploration of optical transients (OTs), conducted with a 60/90 cm Schmidt telescope at Xinglong station of NAOC. This survey repeatedly covers ~ 1000 square degrees of the north sky with a cadence of 3-4 days. With an exposure of 60 s, the survey reaches a limited unfiltered magnitude of about 19.5 mag. This enables us to discover supernovae at their relatively young stages. In this paper, we describe the overall performance of our survey during the first year and present some preliminary results.

We compare the host galaxies of 902 supernovae, including SNe Ia, SNe II and SNe Ibc, which are selected by cross-matching the Asiago Supernova Catalog with the SDSS Data Release 7. We further selected 213 galaxies by requiring the light fraction of spectral observations $>$15%, which could represent well the global properties of the galaxies. Among them, 135 galaxies appear on the Baldwin-Phillips-Terlevich diagram, which allows us to compare the hosts in terms of star-forming, AGNs (including composites, LINERs and Seyfert 2s) and "Absorp" (their related emission-lines are weak or non-existence) galaxies. The diagrams related to parameters D$_n$(4000), H$\delta_A$, stellar masses, SFRs and specific SFRs for the SNe hosts show that almost all SNe II and most of SNe Ibc occur in SF galaxies, which have a wide range of stellar mass and low D$_n$(4000). The SNe Ia hosts as SF galaxies follow similar trends. A significant fraction of SNe Ia occurs in AGNs and Absorp galaxies, which are massive and have high D$_n$(4000). The stellar population analysis from spectral synthesis fitting shows that the hosts of SNe II have a younger stellar population than hosts of SNe Ia. These results are compared with those of the 689 comparison galaxies where the SDSS fiber captures less than 15% of the total light. These comparison galaxies appear biased towards higher 12+log(O/H) ($\sim$0.1dex) at a given stellar mass. Therefore, we believe the aperture effect should be kept in mind when the properties of the hosts for different types of SNe are discussed.

We select 947 star-forming galaxies from SDSS-DR7 with [O~{\sc iii}]$\lambda$4363 emission lines detected at a signal-to-noise {ratio }larger than 5$\sigma$. Their electron temperatures and direct oxygen abundances are {then }determined. {W}e compare the results from different methods. $t_2${, the} electron temperature in {the }low ionization region{,} estimated from $t_3${, that} in {the }high ionization region{,} {is} compared {using} three analysis relations between $t_2-t_3${. These} show obvious differences, which result in some different ionic oxygen abundances. The results of $t_3$, $t_2$, {$\rm O^{++}$/$\rm H^+$} and {$\rm O^{+}$/$\rm H^+$} derived by using methods from IRAF and literature are also compared. The ionic abundances $\rm O^{++}$/$\rm H^+$ {are} higher than $\rm O^{+}$/$\rm H^+$ for most cases. The{ different} oxygen abundances derived from $T_{\rm e}$ and the strong-line ratios show {a }clear discrepancy, which is more obvious following increasing stellar mass and strong-line ratio $R_{23}$. The sample{ of} galaxies from SDSS {with} detected [O~{\sc iii}]$\lambda$4363 have lower metallicites and higher {star formation rates}, {so} they may not be typical representatives of the whole{ population of} galaxies. Adopting data objects from {Andrews \& Martini}, {Liang et al.} and {Lee et al.} data, we derive new relations of stellar mass and metallicity for star-forming galaxies in a much wider stellar mass range: from $10^6\,M_\odot$ to $10^{11}\,M_\odot$.

As the second part of our study, in this paper we proceed to explore the dynamics of the high-inclination Twotinos in the 1:2 Neptune mean motion resonance (NMMR). Depending on the inclination $i$, we show the existence of two critical eccentricities $e_a(i)$ and $e_c(i)$, which are lower limits of the eccentricity $e$ for the resonant angle $\sigma$ to exhibit libration and asymmetric libration, respectively. Accordingly, we have determined the libration centres $\sigma_0$ for inclined orbits, which are strongly dependent on $i$. With initial $\sigma=\sigma_0$ on a fine grid of $(e, i)$, the stability of orbits in the 1:2 NMMR is probed by 4-Gyr integrations. It is shown that symmetric librators are totally unstable for $i\ge30^{\circ}$; while stable asymmetric librators exist for $i$ up to $90^{\circ}$. We further investigate the 1:2 NMMR capture and retention of planetesimals with initial inclinations $i_0\le90^{\circ}$ in the planet migration model using a time-scale of $2\times10^7$ yr. We find that: (1) the capture efficiency of the 1:2 NMMR decreases drastically with the increase of $i_0$, and it goes to 0 when $i_0\ge60^{\circ}$; (2) the probability of discovering Twotinos with $i>25^{\circ}$, beyond observed values, is roughly estimated to be $\le0.1$ per cent; (3) more particles are captured into the leading rather than the trailing asymmetric resonance for $i_0\le10^{\circ}$, but this number difference appears to be the opposite at $i_0=20^{\circ}$ and is continuously varying for even larger $i_0$; (4) captured Twotinos residing in the trailing resonance or having $i>15^{\circ}$ are practically outside the Kozai mechanism, like currently observed samples.

The dynamics of the high-inclination Plutinos is systematically studied. We first present the peculiar features of the 2:3 Neptune mean motion resonance (NMMR) for inclined orbits, especially for the correlation of resonant amplitude A_{\sigma} with inclination i. Using the numerical integrations for the age of the Solar system, the dynamical structure of the 2:3 NMMR is mapped out on the plane of semi-major axis versus i for different eccentricities. We have shown that i of stable resonant orbits could be as high as 90 deg; and the stable region is roughly surrounded by the contours of A_{\sigma} = 120 deg. These new findings allow us to further explore the 2:3 NMMR capture and retention of planetesimals with initial inclinations i0 =< 90 deg in the frame of the planet migration model. We find that the outward transportation of Plutinos is possible for any inclined or even perpendicular orbits. The role of i0 in the formation of Plutinos during Neptune's migration is highlighted and interesting results are obtained: (1) The capture efficiency of the 2:3 NMMR decreases drastically first with the increase of i0, but it then raises instead when i0 exceeds ~ 50 deg; (2) The magnitude of i-variation is limited to less than 5 deg for any i0, and moreover, for Plutinos with i > 48 deg, their i are forced to decrease throughout the outward migration; (3) Plutinos with i > 48 deg are certainly outside the Kozai mechanism, since an inclination increase is prohibited by the migrating 2:3 NMMR; (4) The 7:11 inclination-type NMMR could be responsible for nearly-circular Plutinos, and a minimum i0 ~ 15 deg is required to intrigue this mechanism.

Due to the angular momentum exchange with planetesimals, Neptune might have migrated outward to the current position, and captured many Kuiper belt objects (KBOs) into resonances. We set up a semi-analytic model to simulate the outward migration of Neptune, and the processes of resonant captures. Our model can naturally explain Neptune's currently observed semi-major axis and eccentricity. The results show that the current population ratio between 3:2 and 2:1 is mainly due to the original density distribution of KBOs, which might be related to drag-induced inward migrations of proto-KBOs.

The population of Neptune Trojans is believed to be bigger than that of Jupiter Trojans and that of asteroids in the main belt, although only eight members of this far distant asteroid swarm have been observed up to now. Six leading Neptune Trojans around the Lagrange point L4 discovered earlier have been studied in detail, but two trailing ones found recently around the L5 point, 2004 KV18 and 2008 LC18, have not been investigated yet. In this paper, we report our investigations on the dynamical behaviors of these two new Neptune Trojans. Our calculations show that the asteroid 2004 KV18 is a temporary Neptune Trojan. Most probably, it was captured into the trailing Trojan cloud no earlier than 203kyr ago, and it will not keep this identity no later than 165kyr in future. Based on the statistics on our orbital simulations, we argue that this object is more like a scattered Kuiper belt object. On the contrary, the orbit of asteroid 2008 LC18 is much more stable. Among the clone orbits spread within the orbital uncertainties, a considerable portion of clones may survive on the L5 tadpole orbits for 4Gyr. The strong dependence of the stability on the semimajor axis and resonant angle suggests that further observations are badly needed to confine the orbit in the stable region. We also discuss the implications of the existence and dynamics of these two trailing Trojans on the Solar system history.

Recently the first Earth Trojan has been observed (Mainzer et al., ApJ 731) and found to be on an interesting orbit close to the Lagrange point L4 (Connors et al., Nature 475). In the present study we therefore perform a detailed investigation on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our results are derived using different approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the phase space geometry, including the indirect influence of the additional planets of our Solar system. b) We use precise numerical methods to integrate the orbit forward and backward in time in different dynamical models. Based on a set of 400 clone orbits we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive numerical investigation of the stability region of the Earth's Lagrangian points. We present a detailed parameter study in the regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of stability.

The luminescent properties of CsI(Na) crystals are studied in this report. By using a TDS3054C oscilloscope with a sampling frequency of 5 GS/s, we find out that nuclear recoil signals are dominated by very fast light pulse with a decay time of ~20 ns, while {\gamma}-ray signals have a decay time of ~600 ns. The wavelength of nuclear recoil and {\gamma}-ray signals are also different. The study of n/{\gamma} separation shows that the rejection factor can reach an order of 10-7 with signal efficiency more than 80% at an equivalent electron recoil energy of 20 keV or more. Such a property makes CsI(Na) an ideal candidate for dark matter searches.

In a previous paper, we have presented a global view of the stability of Neptune Trojan (NT hereafter) on inclined orbit. We discuss in this paper the dependence of stability of NT orbits on the eccentricity. High-resolution dynamical maps are constructed using the results of extensive numerical integrations of orbits initialized on the fine grids of initial semimajor axis (a0) versus eccentricity (e0). The extensions of regions of stable orbits on the (a0, e0) plane at different inclinations are shown. The maximum eccentricities of stable orbits in three most stable regions at low (0, 12deg.), medium (22,36deg.) and high (51, 59deg.) inclination, are found to be 0.10, 0.12 and 0.04, respectively. The fine structures in the dynamical maps are described. Via the frequency analysis method, the mechanisms that portray the dynamical maps are revealed. The secondary resonances, concerning the frequency of the librating resonant angle and the frequency of the quasi 2:1 mean motion resonance between Neptune and Uranus, are found deeply involved in the motion of NTs. Secular resonances are detected and they also contribute significantly to the triggering of chaos in the motion. Particularly, the effects of the secular resonance v8, v18 are clarified. We also investigate the orbital stabilities of six observed NTs by checking the orbits of hundreds clones of them generated within the observing error bars. We conclude that four of them, except 2001 QR322 and 2005 TO74, are deeply inside the stable region. The 2001 QR322 is in the close vicinity of the most significant secondary resonance. The 2005 TO74 locates close to the boundary separating stable orbits from unstable ones, and it may be influenced by a secular resonance.

The area of stable motion for fictitious Trojan asteroids around Uranus' equilateral equilibrium points is investigated with respect to the inclination of the asteroid's orbit to determine the size of the regions and their shape. For this task we used the results of extensive numerical integrations of orbits for a grid of initial conditions around the points L4 and L5, and analyzed the stability of the individual orbits. Our basic dynamical model was the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected orbits up to the age of the Solar system of 5 billion years. One experiment has been undertaken for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor axes were varied. The extension of the stable region with respect to the initial semimajor axis lies between 19.05 < a < 19.3 AU but depends on the initial inclination. In another run the inclination of the asteroids' orbit was varied in the range 0 < i < 60 and the semimajor axes were fixed. It turned out that only four 'windows' of stable orbits survive: these are the orbits for the initial inclinations 0 < i < 7, 9 < i < 13, 31 < i < 36 and 38 < i < 50. We postulate the existence of at least some Trojans around the Uranus Lagrange points for the stability window at small and also high inclinations.

The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the FFT technique, we construct high resolution dynamical maps on the plane of initial semimajor axis $a_0$ versus inclination $i_0$. These maps show three most stable regions, with $i_0$ in the range of $(0^\circ,12^\circ), (22^\circ,36^\circ)$ and $(51^\circ,59^\circ)$ respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points $L_4$ and $L_5$ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the $\nu_8$ secular resonance around $i_0\sim44^\circ$ pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination ($>60^\circ$) and an unstable gap around $44^\circ$ on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of $(a_0,i_0)$. The fine structures in the dynamical maps can be explained by these secular resonances.

The dynamics of artificial asteroids on the Trojan-like orbits around Neptune is investigated in this paper. We describe the dependence of the orbital stability on the initial semimajor axis a and inclination i by constructing a dynamical map on the (a,i)-plane. Rich details are revealed in the dynamical map, especially a unstable gap at i=45 deg is determined and the mechanism triggering chaos in this region is figured out. Our investigation can be used to guide the observations.

We report in this paper the numerical simulations of the capture into the 3:1 mean-motion resonance between the planet b and c in the 55 Cancri system. The results show that this resonance can be obtained by a differential planetary migration. The moderate initial eccentricities, relatively slower migration and suitable eccentricity damping rate increase significantly the probability of being trapped in this resonance. Otherwise, the system crosses the 3:1 commensurability avoiding resonance capture, to be eventually captured into a 2:1 resonance or some other higher-order resonances. After the resonance capture, the system could jump from one orbital configuration to another one if the migration continues, making a large region of the configuration space accessible for a resonance system. These investigations help us understand the diversity of resonance configurations and put some constrains on the early dynamical evolution of orbits in the extra-solar planetary systems.

The inner two planets around the 55 Cancri were found to be trapped in the 3:1 mean motion resonance. In this paper, we study the dynamics of this extra-solar planetary system. Our numerical investigation confirms the existence of the 3:1 resonance and implies a complex orbital motion. Different stable motion types, with and without the apsidal corotation, are found. Due to the high eccentricities in this system, we apply a semi-analytical method based on a new expansion of the Hamiltonian of the planar three-body problem in the discussion. We analyse the occurrence of the apsidal corotation in this mean motion resonance and its influence on the stability of the system.