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.

One of the very common in situ signatures of interplanetary coronal mass ejections (ICMEs), as well as other interplanetary transients, are Forbush decreases (FDs), i.e. short-term reductions in the galactic cosmic ray (GCR) flux. A two-step FD is often regarded as a textbook example, which presumably owes its specific morphology to the fact that the measuring instrument passed through the ICME head-on, encountering first the shock front (if developed), then the sheath and finally the CME magnetic structure. The interaction of GCRs and the shock/sheath region, as well as the CME magnetic structure, occurs all the way from Sun to Earth, therefore, FDs are expected to reflect the evolutionary properties of CMEs and their sheaths. We apply modelling to different ICME regions in order to obtain a generic two-step FD profile, which qualitatively agrees with our current observation-based understanding of FDs. We next adapt the models for energy dependence to enable comparison with different GCR measurement instruments (as they measure in different particle energy ranges). We test these modelling efforts against a set of multi-spacecraft observations of the same event, using the Forbush decrease model for the expanding flux rope (ForbMod). We find a reasonable agreement of the ForbMod model for the GCR depression in the CME magnetic structure with multi-spacecraft measurements, indicating that modelled FDs reflect well the CME evolution.

RV variable stars are important in astrophysics. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic survey has provided ~ 6.5 million stellar spectra in its Data Release 4 (DR4). During the survey, ~ 4.7 million unique sources were targeted and ~ 1 million stars observed repeatedly. The probabilities of stars being RV variables are estimated by comparing the observed radial velocity variations with the simulated ones. We build a catalog of 80,702 RV variable candidates with probability greater than 0.60 by analyzing the duplicate-observed multi-epoch sources covered by the LAMOST DR4. Simulations and cross-identifications show that the purity of the catalog is higher than 80%. The catalog consists of 77% binary systems and 7% pulsating stars as well as 16% pollution by single stars. 3,138 RV variables are classified through cross-identifications with published results in literatures. By using the 3,138 sources common to both LAMOST and a collection of published RV variable catalogs we are able to analyze LAMOST's RV variable detection rate. The efficiency of the method adopted in this work relies not only on the sampling frequency of observations but also periods and amplitudes of RV variables. With the progress of LAMOST, Gaia and other surveys, more and more RV variables would will be confirmed and classified. This catalog is valuable for other large-scale surveys, especially for RV variable searches. The catalog will be released according to the LAMOST Data Policy via http://dr4.lamost.org.

The concept of the Solar Ring mission was gradually formed from L5/L4 mission concept, and the proposal of its pre-phase study was funded by the National Natural Science Foundation of China in November 2018 and then by the Strategic Priority Program of Chinese Academy of Sciences in space sciences in May 2019. Solar Ring mission will be the first attempt to routinely monitor and study the Sun and inner heliosphere from a full 360-degree perspective in the ecliptic plane. The current preliminary design of the Solar Ring mission is to deploy six spacecraft, grouped in three pairs, on a sub-AU orbit around the Sun. The two spacecraft in each group are separated by about 30 degrees and every two groups by about 120 degrees. This configuration with necessary science payloads will allow us to establish three unprecedented capabilities: (1) determine the photospheric vector magnetic field with unambiguity, (2) provide 360-degree maps of the Sun and the inner heliosphere routinely, and (3) resolve the solar wind structures at multiple scales and multiple longitudes. With these capabilities, the Solar Ring mission aims to address the origin of solar cycle, the origin of solar eruptions, the origin of solar wind structures and the origin of severe space weather events. The successful accomplishment of the mission will advance our understanding of the star and the space environment that hold our life and enhance our capability of expanding the next new territory of human.

The portable adaptive optics (PAO) is a high-contrast instrument that is optimized for the direct imaging of exoplanets and sub-stellar companions. The PAO has a compact physical size, high flexibility and high performance. As a visiting instrument, PAO has been used for the imaging observation of exoplanets with current 3-4 meter class telescopes. In these observations, it has delivered a diffraction limited image in H band. Combined with our unique image rotation and subtraction (IRS) technique and the optimized IRS algorithm, our PAO has successfully recovered the known exoplanets of kappa And b, in recent observation at 3.5-meter ARC telescope at Apache Point Observatory. The observation contrast has reached 10E-5 at an angular distance of 1 arcsec. We have performed the associated astrometry and photometry analysis of the recovered kappa And b planet, which gives a project separation of 1.044 +/- 0.016 arcsec, a position angle of 53.5 +/- 0.5 degrees, and a mass of 10.15 (-1.255) (+2.19) MJup. It has fully demonstrated that our portable high-contrast imaging system can be used for the exoplanets direct imaging observation with middle-class telescopes.

Space radiation is a major risk for humans, especially on long-duration missions to outer space, e.g., a manned mission to Mars. Galactic cosmic rays (GCR) contribute a predictable radiation background, the main risk is due to the highly variable and currently unpredictable flux of solar energetic particles (SEPs). Such sporadic SEP events may induce acute health effects and are thus considered a critical mission risk for future human exploration of Mars. Therefore, it is of utmost importance to study, model, and predict the surface radiation environment during such events. It is well known that the deep-space SEP differential energy spectrum at high energies is often given by a power law. We use a measurement-validated particle transport code to show that, for large SEP events with proton energy extending above ~ 500 MeV with a power-law distribution, it is sufficient to measure the SEP flux at a pivot energy of ~ 300 MeV above the Martian atmosphere to predict the dose rate on the Martian surface. In conjunction with a validation by in-situ measurements from the Martian surface, this remarkable simplification and elegant quantification could enable instant predictions of the radiation environment on the surface of Mars upon the onset of large SEP events.

The Radiation Assessment Detector (RAD), onboard the Mars Science Laboratory (MSL) rover Curiosity, has been measuring the energetic charged and neutral particles and the radiation dose rate on the surface of Mars since the landing of the rover in August 2012. In contrast to charged particles, neutral particles (neutrons and gamma-rays) are measured indirectly: the energy deposition spectra produced by neutral particles are complex convolutions of the incident particle spectra with the detector response functions. An inversion technique has been developed and applied to jointly unfold the deposited energy spectra measured in two scintillators of different types (CsI for high gamma detection efficiency, and plastic for neutrons) to obtain the neutron and gamma-ray spectra. This result is important for determining the biological impact of the Martian surface radiation contributed by neutrons, which interact with materials differently from the charged particles. These first in-situ measurements on Mars provide (1) an important reference for assessing the radiation-associated health risks for future manned missions to the red planet and (2) an experimental input for validating the particle transport codes used to model the radiation environments within spacecraft or on the surface of planets. Here we present neutral particle spectra as well as the corresponding dose and dose equivalent rates derived from RAD measurement during a period (November 15, 2015 to January 15, 2016) for which the surface particle spectra have been simulated via different transport models.

The observation of cosmic sub-PeV gamma-rays from the Crab Nebula opens up the possibility of testing cosmic ray photon transparency at the multi-hundred TeV scale. Assuming no deviation from a source gamma-ray emission due to accelerated electron inverse-Compton scattering, higher event energies can extend constraints on the effects of new physics; We consider oscillation between gamma-rays and axions/dark photons, plus attenuation effects from gamma-ray absorption in the case of dark photon dark matter. Combining the recent AS$\gamma$ and HAWC sub-PeV data with earlier MAGIC and HEGRA data, axion-like particles are most constrained in the $10^{-7}-10^{-6}$ eV mass range, where the coupling $g_{a\gamma}$ is constrained to be below ${1.8}\times 10^{-10}$ GeV$^{-1}$. In comparison, gamma ray flux attenuation due to oscillation with a dark photon leads to a very weak constraint on the mixing parameter; $\epsilon\lesssim$ 0.2 for dark photon mass between $10^{-7}$ and $10^{-6}$ eV. Direct scattering from dark photon dark matter limits $\epsilon\lesssim$ 0.01 for masses between $6$ and $400$ eV.

We investigate the axion like particle (ALP)-photon oscillation effect in the high energy $\gamma$-ray spectra of PG 1553+113 and PKS 2155-304 measured by Fermi-LAT and H.E.S.S.. The choice of extragalactic background light (EBL) model, which induces the attenuate effect in observed $\gamma$-ray spectra, would affect the ALP implication. For the ordinary EBL model that prefers a null hypothesis, we set constraint on the ALP-photon coupling constant at 95% C.L. as $g_{a\gamma}\lesssim 5\times 10^{-11} ~\rm{GeV}^{-1}$ for the ALP mass $\sim 10$ neV. We also consider the CIBER observation of the cosmic infrared radiation, which shows an excess at the wave wavelength of $\sim 1~\mu$m after the substraction of foregrounds. The high energy gamma-rays from extragalactic sources at high redshifts would suffer from a more significant attenuate effect caused by this excess. In this case, we find that the ALP-photon oscillation would improve the fit to the observed spectra of PKS 2155-304 and PG 1553+113 and find a favored parameter region at 95% C.L..

For filter-based magnetographs, the linear calibration method under the weak-field assumption is usually adopted; this leads to magnetic saturation effect in the regions with strong magnetic field. This article explores a new method to overcome the above disadvantage using a multilayer perceptron network, which we call MagMLP, based on a back-propagation algorithm with one input layer, five hidden layers, and one output layer. We use the data from the \textit{Spectropolarimeter} (SP) on board \textit{Hinode} to simulate single-wavelength observations for the model training, and take into account the influence of the Doppler velocity field and the filling factor. The training results show that the linear fitting coefficient (LFC) of the transverse field reaches above 0.91, and that of the longitudinal field is above 0.98. The generalization of the models is good because the corresponding LFCs are above 0.9 for the test subsets. Compared with the linear calibration method, the MagMLP is much more effective on dealing with the magnetic saturation effect. Analyzing an active region, the results of the linear calibration present an evident magnetic saturation effect in the umbra regions; the corresponding systematic error reaches values greater than 1000 G in most areas, or even exceeds 2000 G at some pixels. However, the results of MagMLP at these locations are very close to the inversion results, and the systematic errors are basically within 300 G. In addition, we find that there are many "bright spots" and "dark spots" on the inclination angle images from the inversion results of \textit{Hinode}/SP with values of 180 and 0 degrees, respectively, where the inversion is not reliable and does not produce a good result; the MagMLP handles these points well.

Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von K\'arm\'an crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content.

In this paper, we report the multiwavelength observations of a very long filament in active region (AR) 11112 on 2010 October 18. The filament was composed of two parts, the eastern part (EP) and western part (WP). We focus on longitudinal oscillations of the EP, which were enhanced by two homologous C-class flares in the same AR. The C1.3 flare was confined without a CME. Both EP and WP of the filament were slightly disturbed and survived the flare. After 5 hrs, eruption of the WP generated a C2.6 flare and a narrow, jet-like CME. Three oscillating threads (thd$_a$, thd$_b$, thd$_c$) are obviously identified in the EP and their oscillations are naturally divided into three phases by the two flares. The initial amplitude ranges from 1.6 to 30 Mm with a mean value of $\sim$14 Mm. The period ranges from 34 to 73 minutes with a mean value of $\sim$53 minutes. The curvature radii of the magnetic dips are estimated to be 29 to 133 Mm with a mean value of $\sim$74 Mm. The damping times ranges from $\sim$62 to $\sim$96 minutes with a mean value of $\sim$82 minutes. The value of $\tau/P$ is between 1.2 and 1.8. For thd$_a$ in the EP, the amplitudes were enhanced by the two flares from 6.1 Mm to 6.8 Mm after the C1.3 flare and further to 21.4 Mm after the C2.6 flare. The period variation as a result of perturbation from the flares was within 20\%. The attenuation became faster after the C2.6 flare. To the best of our knowledge, this is the first report of large-amplitude, longitudinal filament oscillations enhanced by flares. Numerical simulations reproduce the oscillations of thd$_a$ very well. The simulated amplitudes and periods are close to the observed values, while the damping time in the last phase is longer, implying additional mechanisms should be taken into account apart from radiative loss.

Coronal mass ejections (CMEs) are the primary sources of intense disturbances at Earth, where their geo-effectiveness is largely determined by their dynamic pressure and internal magnetic field, which can be significantly altered during interactions with other CMEs in interplanetary space. We analyse three successive CMEs that erupted from the Sun during September 4-6, 2017, investigating the role of CME-CME interactions as source of the associated intense geomagnetic storm (Dst_min=-142 nT on September 7). To quantify the impact of interactions on the (geo-)effectiveness of individual CMEs, we perform global heliospheric simulations with the EUHFORIA model, using observation-based initial parameters with the additional purpose of validating the predictive capabilities of the model for complex CME events. The simulations show that around 0.45 AU, the shock driven by the September 6 CME started compressing a preceding magnetic ejecta formed by the merging of two CMEs launched on September 4, significantly amplifying its Bz until a maximum factor of 2.8 around 0.9 AU. The following gradual conversion of magnetic energy into kinetic and thermal components reduced the Bz amplification until its almost complete disappearance around 1.8 AU. We conclude that a key factor at the origin of the intense storm triggered by the September 4-6, 2017 CMEs was their arrival at Earth during the phase of maximum Bz amplification. Our analysis highlights how the amplification of the magnetic field of individual CMEs in space-time due to interaction processes can be characterised by a growth, a maximum, and a decay phase, suggesting that the time interval between the CME eruptions and their relative speeds are critical factors in determining the resulting impact of complex CMEs at various heliocentric distances (helio-effectiveness).

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

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

Since the mass loss rates are the function of the mean density of the planet and the stellar irradiation, we calculated about 450 models covering planets with different densities and stellar irradiation. Our results show that the mass loss rates are dependent on the stellar irradiation and the mean density. However, the mass loss rates predicted by the energy-limited equation are higher than that of hydrodynamic model when the XUV integrated flux is higher than $\sim$2$\times$10$^{4}$ erg/cm$^{2}$/s. The overestimation can be revised if the kinetic and thermal energy of the escaping atmosphere is included in the energy-limited equation. We found that the heating efficiencies are proportional to the product of the gravitational potential of the planet and the stellar irradiation. The mean absorption radii of stellar irradiation are 1.1-1.2 R$_{p}$ for the Jupiter-like planets while they vary in the range of 1.1-1.7 R$_{p}$ for the planets with smaller sizes. We evaluated the absorption of stellar Ly$\alpha$ by planetary atmosphere and found that the deeper Ly$\alpha$ absorptions tend to locate in the high stellar irradiation and low planetary mean density regions, and vice versa. Moreover, planets with mass loss rates higher than 10$^{11} g/s$ are likely to exhibit obvious absorptions. Finally, we suggested that the absorption levels are related to the inherent properties of the exoplanets. The planets with larger sizes (or lower mean density) show strong Ly$\alpha$ absorptions. Neptune-like and Earth-like planets tend to have weak Ly$\alpha$ absorptions because of their small sizes (or high densities).

On July 25 2017 a multi-step Forbush decrease (FD) with the remarkable total amplitude of more than 15\% was observed by MSL/RAD at Mars. We find that these particle signatures are related to very pronounced plasma and magnetic field signatures detected in situ by STEREO-A on July 24 2017, with a higher than average total magnetic field strength reaching more than 60 nT. In the observed time period STEREO-A was at a relatively small longitudinal separation (46 degrees) to Mars and both were located at the back side of the Sun as viewed from Earth. We analyse a number of multi-spacecraft and multi-instrument (both in situ and remote-sensing) observations, and employ modelling to understand these signatures. We find that the solar sources are two CMEs which erupted on July 23 2017 from the same source region on the back side of the Sun as viewed from Earth. Moreover, we find that the two CMEs interact non-uniformly, inhibiting the expansion of one of the CMEs in STEREO-A direction, whereas allowing it to expand more freely in the Mars direction. The interaction of the two CMEs with the ambient solar wind adds up to the complexity of the event, resulting in a long, sub-structured interplanetary disturbance at Mars, where different sub-structures correspond to different steps of the FD, adding-up to a globally large-amplitude FD.

We propose the thick-disc model of Gu et al. 2016 to interpret the transition between soft ultraluminous state (SUL) and supersoft ultraluminous (SSUL) state in NGC 247. As accretion rate increases, the inner disc will puff up and act as shield to block the innermost X-ray emission regions and absorb both soft and hard X-ray photons. The absorbed X-ray emission will be re-radiated as a much softer blackbody X-ray spectrum. Hence NGC 247 shows flux dips in the hard X-ray band and transits from the SUL state to the SSUL state. The $\sim 200$s transition timescale can be explained by the viscous timescale. According to our model, the inner disc in the super-soft state is thicker and has smaller viscous timescale than in the soft state. X-ray flux variability, which is assumed to be driven by accretion rate fluctuations, might be viscous time-scale invariant. Therefore, in the SSUL state, NGC 247 is more variable. The bolometric luminosity is saturated in the thick disc; the observed radius-temperature relation can therefore be naturally explained.

X-ray emission is an important indicator of stellar activity. In this paper, we study stellar X-ray activity using the XMM-Newton and LAMOST data for different types of stars. We provide a sample including 1259 X-ray emitting stars, of which 1090 have accurate stellar parameter estimations. Our sample size is much larger than those in previous works. We find a bimodal distribution of X-ray to optical flux ratio (log(fX/fV)) for G and K stars. We interpret that this bimodality is due to two subpopulations with different coronal heating rates. Furthermore, using the full widths at half maxima calculated from H{\alpha} and H\b{eta} lines, we show that these stars in the inactive peaks have smaller rotational velocities. This is consistent with the magnetic dynamo theory that stars with low rotational velocities have low levels of stellar activity. We also examine the correlation between log(fX/fV) and luminosity of the excess emission in the H{\alpha} line, and find a tight relation between the coronal and chromospheric activity indicators.

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

We probed stellar X-ray activity over a wide range of stellar parameters, using {\it Chandra} and LAMOST data. We measured the X-ray-to-bolometric luminosity ratio ($R_X = L_X/L_{\rm bol}$) for 484 main-sequence stars, and found a bimodal distribution for G and K types. We interpret this bimodality as evidence of two sub-populations with different coronal temperatures, which are caused by different coronal heating rates. Using the metallicity and velocity information, we find both of the two sub-populations are mostly located in the thin disk. We find no trend of $R_X$ with stellar age for stars older than $\sim$ 4 Gyr; there is a trough in the $R_X$ vs age distribution, with the lowest range of $R_X$ appearing at ages around 2 Gyr. We then examined the correlation between $R_X$ and $R_{\rm H\alpha}$ (proxy of chromospheric activity): we find that the two quantities are well correlated, as found in many earlier studies. Finally, we selected a sample of twelve stars with X-ray flares, and studied the light-curve morphology of the flares. The variety of flare profiles and timescales observed in our sample suggests the contribution of different processes of energy release.

The X8.2 event of 10 September 2017 provides unique observations to study the genesis, magnetic morphology and impulsive dynamics of a very fast CME. Combining GOES-16/SUVI and SDO/AIA EUV imagery, we identify a hot ($T\approx 10-15$ MK) bright rim around a quickly expanding cavity, embedded inside a much larger CME shell ($T\approx 1-2$ MK). The CME shell develops from a dense set of large AR loops ($\gtrsim$0.5 $R_s$), and seamlessly evolves into the CME front observed in LASCO C2. The strong lateral overexpansion of the CME shell acts as a piston initiating the fast EUV wave. The hot cavity rim is demonstrated to be a manifestation of the dominantly poloidal flux and frozen-in plasma added to the rising flux rope by magnetic reconnection in the current sheet beneath. The same structure is later observed as the core of the white light CME, challenging the traditional interpretation of the CME three-part morphology. The large amount of added magnetic flux suggested by these observations explains the extreme accelerations of the radial and lateral expansion of the CME shell and cavity, all reaching values of $5 - 10$ km s$^{-2}$. The acceleration peaks occur simultaneously with the first RHESSI $100-300$ keV hard X-ray burst of the associated flare, further underlining the importance of the reconnection process for the impulsive CME evolution. Finally, the much higher radial propagation speed of the flux rope in relation to the CME shell causes a distinct deformation of the white light CME front and shock.

This paper reports on the re-analysis of solar flares in which the hard X-rays (HXRs) come predominantly from the corona rather than from the more usual chromospheric footpoints. All of the 26 previously analyzed event time intervals, over 13 flares, are re-examined for consistency with a flare model in which electrons are accelerated near the top of a magnetic loop that has a sufficiently high density to stop most of the electrons by Coulomb collisions before they can reach the footpoints. Of particular importance in the previous analysis was the finding that the length of the coronal HXR source increased with energy in the 20 - 30 keV range. However, after allowing for the possibility that footpoint emission at the higher energies affects the inferred length of the coronal HXR source, and using analysis techniques that suppress the possible influence of such footpoint emission, we conclude that there is no longer evidence that the length of the HXR coronal sources increase with increasing energy. In fact, for the 6 flares and 12 time intervals that satisfied our selection criteria, the loop lengths decreased on average by 1.0 +/- 0.2 arcsec between 20 and 30 keV, with a standard deviation of 3.5 arcsec. We find strong evidence that the peak of the coronal HXR source increases in altitude with increasing energy. For the thermal component of the emission, this is consistent with the standard CHSKP flare model in which magnetic reconnection in a coronal current sheet results in new hot loops being formed at progressively higher altitudes. The explanation for the nonthermal emission is not so clear.

On 2017-09-10, solar energetic particles (SEPs) originating from the active region 12673 were registered as a ground level enhancement (GLE) at Earth and the biggest GLE on the surface of Mars as observed by the Radiation Assessment Detector (RAD) since the landing of the Curiosity rover in August 2012. Based on multi-point coronagraph images, we identify the initial 3D kinematics of an extremely fast CME and its shock front as well as another 2 CMEs launched hours earlier (with moderate speeds) using the Graduated Cylindrical Shell (GCS) model. These three CMEs interacted as they propagated outwards into the heliosphere and merged into a complex interplanetary CME (ICME). The arrival of the shock and ICME at Mars caused a very significant Forbush Decrease (FD) seen by RAD only a few hours later than that at Earth which is about 0.5 AU closer to the Sun. We investigate the propagation of the three CMEs and the consequent ICME together with the shock using the Drag Based Model (DBM) and the WSA-ENLIL plus cone model constrained by the in-situ SEP and FD/shock onset timing. The synergistic modeling of the ICME and SEP arrivals at Earth and Mars suggests that in order to better predict potentially hazardous space weather impacts at Earth and other heliospheric locations for human exploration missions, it is essential to analyze 1) the CME kinematics, especially during their interactions and 2) the spatially and temporally varying heliospheric conditions, such as the evolution and propagation of the stream interaction regions.

Mixing in the convective core is quite uncertain in core helium burning stars. In order to explore the overshooting mixing beyond the convective core, we incorporated the $k$-$\omega$ proposed by Li (2012, 2017) into MESA, and investigated the overshooting mixing in evolution of sdB models. We found that the development of the convective core can be divided into three stages. When the radiative temperature gradient $\nabla_{\rm rad}$ monotonically decreases from the stellar center, the overshooting mixing presents a behaviour of exponential decay similar with Herwig (2000), and the overshooting distance is to make $\nabla_{\rm rad} \simeq \nabla_{\rm ad}$ at the boundary of the convective core, in agreement with the prediction of the self-driving mechanism Castellani, Giannone \& Renzini (1971a). When the radiative temperature gradient $\nabla_{\rm rad}$ shows a minimum value in the convective core, the convective core may be divided into two convection zones if the minimum value of $\nabla_{\rm rad}$ is smaller than the adiabatic temperature gradient $\nabla_{\rm ad}$. For the single-zone case, the overshooting mixing still shows an exponential decay behaviour, but the overshooting distance is much smaller than in the initial stage. For the double-zone case, the overshooting mixing is similar to that of the single-case beyond the convective core, while it almost stops on both sides of the above convective shell. Our overshooting mixing scheme is similar to the maximal overshoot scheme of Constantino et al. (2015). In the final stage, continuous injection of fresh helium into the convective core by the overshooting mixing happens, which is similar to the "core breathing pulses" (Sweigart \& Demarque 1973; Castellani et al. 1985).

Magnetic flux rope (MFR) is the core structure of the greatest eruptions, i.e., the coronal mass ejections (CMEs), on the Sun, and magnetic clouds are post-eruption MFRs in interplanetary space. There is a strong debate about whether or not a MFR exists prior to a CME and how the MFR forms/grows through magnetic reconnection during the eruption. Here we report a rare event, in which a magnetic cloud was observed sequentially by four spacecraft near Mercury, Venus, Earth and Mars, respectively. With the aids of a uniform-twist flux rope model and a newly developed method that can recover a shock-compressed structure, we find that the axial magnetic flux and helicity of the magnetic cloud decreased when it propagated outward but the twist increased. Our analysis suggests that the `pancaking' effect and `erosion' effect may jointly cause such variations. The significance of the `pancaking' effect is difficult to be estimated, but the signature of the erosion can be found as the imbalance of the azimuthal flux of the cloud. The latter implies that the magnetic cloud was eroded significantly leaving its inner core exposed to the solar wind at far distance. The increase of the twist together with the presence of the erosion effect suggests that the post-eruption MFR may have a high-twist core enveloped by a less-twisted outer shell. These results pose a great challenge to the current understanding on the solar eruptions as well as the formation and instability of MFRs.

Stellar systems composed of single, double, triple or high-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multi-epoch observations, we analyze the binary fraction of dwarf stars observed with the LAMOST. Effects of different model assumptions such as orbital period distributions on the estimate of binary fractions, are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with $T_{\rm eff}$ and decreases with [Fe/H]. We first investigate the relation between $\alpha$-elements and binary fraction in such a large sample as the LAMOST. The old stars with high [$\alpha$/Fe] dominate higher binary fraction than young stars with low [$\alpha$/Fe]. At the same mass, former forming stars possess a higher binary fraction than newly forming ones, which may be related with the evolution of the Galaxy.

We present the ultraviolet magnitudes for over three million stars in the LAMOST survey, in which 2,202,116 stars are detected by $GALEX$. For 889,235 undetected stars, we develop a method to estimate their upper limit magnitudes. The distribution of (FUV $-$ NUV) shows that the color declines with increasing effective temperature for stars hotter than 7000 K in our sample, while the trend disappears for the cooler stars due to upper atmosphere emission from the regions higher than their photospheres. For stars with valid stellar parameters, we calculate the UV excesses with synthetic model spectra, and find that the (FUV $-$ NUV) vs. $R'_{\mathrm{FUV}}$ can be fitted with a linear relation and late-type dwarfs tend to have high UV excesses. There are 87,178 and 1,498,103 stars detected more than once in the visit exposures of $GALEX$ in the FUV and NUV, respectively. We make use of the quantified photometric errors to determine statistical properties of the UV variation, including intrinsic variability and the structure function on the timescale of days. The overall occurrence of possible false positives is below 1.3\% in our sample. UV absolute magnitudes are calculated for stars with valid parallaxes, which could serve as a possible reference frame in the NUV. We conclude that the colors related to UV provide good criteria to distinguish between M giants and M dwarfs, and the variability of RR Lyrae stars in our sample is stronger than that of other A and F stars.

The injection, propagation and arrival of solar energetic particles (SEPs) during eruptive solar events is an important and current research topic of heliospheric physics. During the largest solar events, particles may have energies up to a few GeVs and sometimes even trigger ground-level enhancements (GLEs) at Earth. We study the first GLE-event of solar cycle 24, from 17th May 2012, using data from multiple spacecraft (SOHO, GOES, MSL, STEREO-A, STEREO-B and MESSENGER). These spacecraft are located throughout the inner heliosphere, at heliocentric distances between 0.34 and 1.5 astronomical units (au), covering nearly the whole range of heliospheric longitudes. We present and investigate sub-GeV proton time profiles for the event at several energy channels, obtained via different instruments aboard the above spacecraft. We investigate issues due to magnetic connectivity, and present results of three-dimensional SEP propagation simulations. We gather virtual time profiles and perform qualitative and quantitative comparisons with observations, assessing longitudinal injection and transport effects as well as peak intensities. We distinguish different time profile shapes for well-connected and weakly connected observers, and find our onset time analysis to agree with this distinction. At select observers, we identify an additional low-energy component of Energetic Storm Particles (ESPs). Using well-connected observers for normalisation, our simulations are able to accurately recreate both time profile shapes and peak intensities at multiple observer locations. This synergetic approach combining numerical modeling with multi-spacecraft observations is crucial for understanding the propagation of SEPs within the interplanetary magnetic field. Our novel analysis provides valuable proof of the ability to simulate SEP propagation throughout the inner heliosphere, at a wide range of longitudes.

We consider the models of vacuum energy interacting with cold dark matter in this study, in which the coupling can change sigh during the cosmological evolution. We parameterize the running coupling $b$ by the form $b(a)=b_0a+b_e(1-a)$, where at the early-time the coupling is given by a constant $b_{e}$ and today the coupling is described by another constant $b_{0}$. We explore six specific models with (i) $Q(a)=b(a)H_0\rho_0$, (ii) $Q(a)=b(a)H_0\rho_{\rm de}$, (iii) $Q(a)=b(a)H_0\rho_{\rm c}$, (iv) $Q(a)=b(a)H\rho_0$, (v) $Q(a)=b(a)H\rho_{\rm de}$, and (vi) $Q(a)=b(a)H\rho_{\rm c}$. The current observational data sets we use to constrain the models include the JLA compilation of type Ia supernova data, the Planck 2015 distance priors data of cosmic microwave background observation, the baryon acoustic oscillations measurements, and the Hubble constant direct measurement. We find that, for all the models, we have $b_0<0$ and $b_e>0$ at around the 1$\sigma$ level, and $b_0$ and $b_e$ are in extremely strong anti-correlation. Our results show that the coupling changes sign during the evolution at about the 1$\sigma$ level, i.e., the energy transfer is from dark matter to dark energy when dark matter dominates the universe and the energy transfer is from dark energy to dark matter when dark energy dominates the universe.

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC) is a photometric sky survey that covers nearly 6 000 deg^2 towards Galactic anti-center in g r i bands. Half of its survey field locates on the Galactic Anti-center disk, which makes XSTPS-GAC highly suitable for searching new open clusters in the GAC region. In this paper, we report new open cluster candidates discovered in this survey, as well as properties of these open cluster candidates, such as age, distance and reddening, derived by isochrone fitting in the color-magnitude diagram (CMD). These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. (2008). Each candidate is inspected on its true color image composed from XSTPS-GAC three band images. Then its CMD is checked, in order to identify whether the central region stars have a clear isochrone-like trend differing from the background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC. Meanwhile, when these new candidates are able to be seen clearly on 2MASS, their parameters are also derived based on 2MASS (J-H, J) CMD. Finally, there are 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Further more, the parameters of these new candidates, as well as another 11 known recovered open clusters, are properly determined for the first time.

The propagation of 15 interplanetary coronal mass ejections (ICMEs) from Earth's orbit (1 AU) to Mars (~ 1.5 AU) has been studied with their propagation speed estimated from both measurements and simulations. The enhancement of magnetic fields related to ICMEs and their shock fronts cause the so-called Forbush decrease, which can be de- tected as a reduction of galactic cosmic rays measured on-ground. We have used galactic cosmic ray (GCR) data from in-situ measurements at Earth, from both STEREO A and B as well as GCR measurements by the Radiation Assessment Detector (RAD) instrument onboard Mars Science Laboratory (MSL) on the surface of Mars. A set of ICME events has been selected during the periods when Earth (or STEREO A or B) and Mars locations were nearly aligned on the same side of the Sun in the ecliptic plane (so-called opposition phase). Such lineups allow us to estimate the ICMEs' transit times between 1 and 1.5 AU by estimating the delay time of the corresponding Forbush decreases measured at each location. We investigate the evolution of their propagation speeds before and after passing Earth's orbit and find that the deceleration of ICMEs due to their interaction with the ambient solar wind may continue beyond 1 AU. We also find a substantial variance of the speed evolution among different events revealing the dynamic and diverse nature of eruptive solar events. Furthermore, the results are compared to simulation data obtained from two CME propagation models, namely the Drag-Based Model and ENLIL plus cone model.

The energetic particle environment on the Martian surface is influenced by solar and heliospheric modulation and changes in the local atmospheric pressure (or column depth). The Radiation Assessment Detector (RAD) on board the Mars Science Laboratory rover Curiosity on the surface of Mars has been measuring this effect for over four Earth years (about two Martian years). The anticorrelation between the recorded surface Galactic Cosmic Ray-induced dose rates and pressure changes has been investigated by Rafkin et al. (2014) and the long-term solar modulation has also been empirically analyzed and modeled by Guo et al. (2015). This paper employs the newly updated HZETRN2015 code to model the Martian atmospheric shielding effect on the accumulated dose rates and the change of this effect under different solar modulation and atmospheric conditions. The modeled results are compared with the most up-to-date (from 14 August 2012 to 29 June 2016) observations of the RAD instrument on the surface of Mars. Both model and measurements agree reasonably well and show the atmospheric shielding effect under weak solar modulation conditions and the decline of this effect as solar modulation becomes stronger. This result is important for better risk estimations of future human explorations to Mars under different heliospheric and Martian atmospheric conditions.

The Radiation Assessment Detector (RAD), on board Mars Science Laboratory's (MSL) Curiosity rover, has been measuring ground level particle fluxes along with the radiation dose rate at the surface of Mars since August 2012. Similar to neutron monitors at Earth, RAD sees many Forbush decreases (FDs) in the galactic cosmic ray (GCR) induced surface fluxes and dose rates. These FDs are associated with coronal mass ejections (CMEs) and/or stream/corotating interaction regions (SIRs/CIRs). Orbiting above the Martian atmosphere, the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft has also been monitoring space weather conditions at Mars since September 2014. The penetrating particle flux channels in the Solar Energetic Particle (SEP) instrument onboard MAVEN can also be employed to detect FDs. For the first time, we study the statistics and properties of a list of FDs observed in-situ at Mars, seen both on the surface by MSL/RAD and in orbit detected by the MAVEN/SEP instrument. Such a list of FDs can be used for studying interplanetary CME (ICME) propagation and SIR evolution through the inner heliosphere. The magnitudes of different FDs can be well-fitted by a power-law distribution. The systematic difference between the magnitudes of the FDs within and outside the Martian atmosphere may be mostly attributed to the energy-dependent modulation of the GCR particles by both the pass-by ICMEs/SIRs and the Martian atmosphere.

For future human missions to Mars, it is important to study the surface radiation environment during extreme and elevated conditions. In the long term, it is mainly Galactic Cosmic Rays (GCRs) modulated by solar activity that contributes to the radiation on the surface of Mars, but intense solar energetic particle (SEP) events may induce acute health effects. Such events may enhance the radiation level significantly and should be detected as immediately as possible to prevent severe damage to humans and equipment. However, the energetic particle environment on the Martian surface is significantly different from that in deep space due to the influence of the Martian atmosphere. Depending on the intensity and shape of the original solar particle spectra as well as particle types, the surface spectra may induce entirely different radiation effects. In order to give immediate and accurate alerts while avoiding unnecessary ones, it is important to model and well understand the atmospheric effect on the incoming SEPs including both protons and helium ions. In this paper, we have developed a generalized approach to quickly model the surface response of any given incoming proton/helium ion spectra and have applied it to a set of historical large solar events thus providing insights into the possible variety of surface radiation environments that may be induced during SEP events. Based on the statistical study of more than 30 significant solar events, we have obtained an empirical model for estimating the surface dose rate directly from the intensities of a power-law SEP spectra.

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

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

Distant luminous quasars provide important information on the growth of the first supermassive black holes, their host galaxies and the epoch of reionization. The identification of quasars is usually performed through detection of their Lyman-$\alpha$ line redshifted to $\sim$ 0.9 microns at z>6.5. Here, we report the discovery of a very Lyman-$\alpha$ luminous quasar, PSO J006.1240+39.2219 at redshift z=6.618, selected based on its red colour and multi-epoch detection of the Lyman-$\alpha$ emission in a single near-infrared band. The Lyman-$\alpha$-line luminosity of PSO J006.1240+39.2219 is unusually high and estimated to be 0.8$\times$10$^{12}$ Solar luminosities (about 3% of the total quasar luminosity). The Lyman-$\alpha$ emission of PSO J006.1240+39.2219 shows fast variability on timescales of days in the quasar rest frame, which has never been detected in any of the known high-redshift quasars. The high luminosity of the Lyman-$\alpha$ line, its narrow width and fast variability resemble properties of local Narrow-Line Seyfert 1 galaxies which suggests that the quasar is likely at the active phase of the black hole growth accreting close or even beyond the Eddington limit.

We analyze sunspots rotation and magnetic transients in NOAA AR 11429 during two X-class (X5.4 and X1.3) flares using the data from the Helioseismic and Magnetic Imager on board the \emph{Solar Dynamics Observatory}. A large leading sunspot with positive magnetic polarity rotated counterclockwise. As expected, the rotation was significantly affected by the two flares. The magnetic transients induced by the flares were clearly evident in the sunspots with negative polarity. They were moving across the sunspots with speed of order $3-7\ \rm km \ s^{-1}$. Furthermore, the trend of magnetic flux evolution of these sunspots exhibited changes associated with the flares. These results may shed light on the understanding of the evolution of sunspots.

We present a new tool for colour-magnitude diagram (CMD) studies, $Powerful~CMD$. This tool is built on the basis of the advanced stellar population synthesis (ASPS) model, in which single stars, binary stars, rotating stars, and star formation history have been taken into account. Via $Powerful~CMD$, the distance modulus, colour excess, metallicity, age, binary fraction, rotating star fraction, and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC6362, NGC6652, NGC6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those with the data of the Hubble Space Telescope (HST). Moreover, we find that the inclusion of binaries in theoretical stellar population models may lead to smaller colour excess compared to the case of single star population models.

A systematic study is performed on radio properties of H$_{2}$O megamaser host Seyfert 2 galaxies, through multi-band radio continuum observations (at 11cm, 6.0 cm, 3.6 cm, 2.0 cm and 1.3 cm) with the Effelsberg 100-m radio telescope within a total time duration of four days. For comparison, a control Seyfert 2 galaxy sample without detected maser emission was also observed. Spectral indices were determined for those sources for which measurements exist at two adjacent bands assuming a power-law dependence S$_\nu \propto \nu^{-\alpha}$, where S is the flux density and $\nu$ is the frequency. Comparisons of the radio continuum properties between megamaser and non-masing Seyfert 2s show no difference in spectral indices. However, a difference in radio luminosity is statistically significant, i.e. the maser galaxies tend to have higher radio luminosities by a factor of 2 to 3 than the non-masing ones, commonly reaching values above a critical threshold of 10$^{29}$ ergs$^{-1}$Hz$^{-1}$. This result confirms an earlier conclusion by Zhang et al. (2012), but is based on superior data with respect to the time interval within which the data were obtained, with respect to the observational facility (only one telescope used), the number of frequency bands.

Based on regularities in rotational splittings, we seek possible multiplets for the observed frequencies of CoRoT 102749568. Twenty-one sets of multiplets are identified, including four sets of multiplets with $l=1$, nine sets of multiplets with $l=2$, and eight sets of multiplets with $l=3$. In particular, there are three complete triplets ($f_{10}$, $f_{12}$, $f_{14}$), ($f_{31}$, $f_{34}$, $f_{35}$), and ($f_{41}$, $f_{43}$, $f_{44}$). The rotational period of CoRoT 102749568 is estimated to be $1.34^{+0.04}_{-0.05}$ days. When doing model fittings, three $l=1$ modes ($f_{12}$, $f_{34}$, and $f_{43}$) and the radial first overtone $f_{13}$ are used. Our results shows that the three nonradial modes ($f_{12}$, $f_{34}$, and $f_{43}$) are mixed modes, which mainly provide constraints on the helium core. The radial first overtone $f_{13}$ mainly provides constraint on the stellar envelope. Hence the size of the helium core of CoRoT 102749568 is determined to be $M_{\rm He}$ = 0.148 $\pm$ 0.003 $M_{\odot}$ and $R_{\rm He}$ = 0.0581 $\pm$ 0.0007 $R_{\odot}$. The fundamental parameters of CoRoT 102749568 are determined to be $M$ = 1.54 $\pm$ 0.03 $M_{\odot}$, $Z=$ 0.006, $f_{\rm ov}$ = 0.004 $\pm$ 0.002, $\log g$ = 3.696 $\pm$ 0.003, $T_{\rm eff}$ = 6886 $\pm$ 70 K, $R$ = 2.916 $\pm$ 0.039 $R_{\odot}$, and $L$ = 17.12 $\pm$ 1.13 $L_{\odot}$.

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

We present spectra of the extreme polar AR Ursae Majoris (AR UMa) which display a clear Al I absorption doublet, alongside spectra taken less than a year earlier in which that feature is not present. Re-examination of earlier SDSS spectra indicates that the Al I absorption doublet was also present $\approx$8 years before our first non-detection. We conclude that this absorbing material is unlikely to be on the surface of either the white dwarf (WD) or the donor star. We suggest that this Al I absorption feature arises in circumstellar material, perhaps produced by the evaporation of asteroids as they approach the hot WD. The presence of any remaining reservoir of rocky material in AR UMa might help to constrain the prior evolution of this unusual binary system. We also apply spectral decomposition to find the stellar parameters of the M dwarf companion, and attempt to dynamically measure the mass of the WD in AR UMa by considering both the radial velocity curves of the H$_\beta$ emission line and the Na I absorption line. Thereby we infer a mass range for the WD in AR UMa of 0.91 $M_{\odot}$ $<$ $M_{\mathrm{WD}}$ $<$ 1.24 $M_{\odot}$.

The $Chandra$ archival data is a valuable resource for various studies on different topics of X-ray astronomy. In this paper, we utilize this wealth and present a uniformly processed data set, which can be used to address a wide range of scientific questions. The data analysis procedures are applied to 10,029 ACIS observations, which produces 363,530 source detections, belonging to 217,828 distinct X-ray sources. This number is twice the size of the $Chandra$ Source Catalog (Version 1.1). The catalogs in this paper provide abundant estimates of the detected X-ray source properties, including source positions, counts, colors, fluxes, luminosities, variability statistics, etc. Cross-correlation of these objects with galaxies shows 17,828 sources are located within the $D_{25}$ isophotes of 1110 galaxies, and 7504 sources are located between the $D_{25}$ and 2$D_{25}$ isophotes of 910 galaxies. Contamination analysis with the log$N$--log$S$ relation indicates that 51.3\% of objects within 2$D_{25}$ isophotes are truly relevant to galaxies, and the "net" source fraction increases to 58.9\%, 67.3\%, and 69.1\% for sources with luminosities above $10^{37}$, $10^{38}$, and $10^{39}$ erg s$^{-1}$. Among the possible scientific uses of this catalog, we discuss the possibility to study intra-observation variability, inter-observation variability, and supersoft sources.

By varying the profiles of stellar extreme ultraviolet (EUV) spectral energy distribution (SED), we tested the influences of stellar EUV SEDs on the physical and chemical properties of the escaping atmosphere. We apply our model to study four exoplanets, HD\,189733b, HD\,209458b, GJ \,436b, and Kepler-11b. We found that the total mass loss rates of an exoplanet, which are determined mainly by the integrated fluxes, are moderately affected by the profiles of the EUV SED, but the composition and species distributions in the atmosphere can be dramatically modified by the different profiles of the EUV SED. For exoplanets with a high hydrodynamic escape parameter ($\lambda$), the amount of atomic hydrogen produced by photoionization at different altitudes can vary by one to two orders of magnitude with the variation of stellar EUV SEDs. The effect of photoionization of H is prominent when the EUV SED is dominated by the low-energy spectral region (400-900${\AA}$), which pushes the transition of H/H$^{+}$ to low altitudes. On the contrary, the transitions of H/H$^{+}$ move to higher altitudes when most photons concentrate in the high-energy spectral region (50-400${\AA}$). For exoplanets with a low $\lambda$, the lower temperatures of the atmosphere make many chemical reactions so important that photoionization alone can no longer determine the composition of the escaping atmosphere. For HD 189733b, it is possible to explain the time variability of \lya\ between 2010 and 2011 by a change in the EUV SED of the host K star, yet invoking only thermal H\,I in the atmosphere.

The formation of relativistic jets by an accreting compact object is one of the fundamental mysteries of astrophysics. While the theory is poorly understood, observations of relativistic jets from systems known as microquasars have led to a well-established phenomenology. Relativistic jets are not expected from sources with soft or supersoft X-ray spectra, although two such systems are known to produce relatively low-velocity bipolar outflows. Here we report optical spectra of an ultraluminous supersoft X-ray source (ULS) in the nearby galaxy M81 (M81 ULS-1) showing blueshifted broad H\alpha\ emission lines, characteristic of baryonic jets with relativistic speeds. The time variable jets have projected velocities ~17 per cent of the speed of light, and seem similar to those in the prototype microquasar SS 433. Such relativistic jets are not expected to be launched from white dwarfs, but an origin from a black hole or neutron star in M81 ULS-1 is hard to reconcile with its constant soft X-rays. The completely unexpected presence of relativistic jets in a ULS challenges the canonical theories for jet formation, but may possibly be explained by a long speculated super-critically accreting black hole with optically thick outflows

Optical counterparts can provide significant constraints on the physical nature of ultraluminous X-ray sources (ULXs). In this letter, we identify six point sources in the error circle of a ULX in M82, namely M82 X-1, by registering Chandra positions onto Hubble Space Telescope images. Two objects are considered as optical counterpart candidates of M82 X-1, which show F658N flux excess compared to the optical continuum that may suggest the existence of an accretion disk. The spectral energy distributions of the two candidates match well with the spectra for supergiants, with stellar types as F5-G0 and B5-G0, respectively. Deep spatially resolved spectroscopic follow-up and detailed studies are needed to identify the true companion and confirm the properties of this BH system.

The Radiation Assessment Detector (RAD), on board Mars Science Laboratory's (MSL) rover Curiosity, measures the {energy spectra} of both energetic charged and neutral particles along with the radiation dose rate at the surface of Mars. With these first-ever measurements on the Martian surface, RAD observed several effects influencing the galactic cosmic ray (GCR) induced surface radiation dose concurrently: [a] short-term diurnal variations of the Martian atmospheric pressure caused by daily thermal tides, [b] long-term seasonal pressure changes in the Martian atmosphere, and [c] the modulation of the primary GCR flux by the heliospheric magnetic field, which correlates with long-term solar activity and the rotation of the Sun. The RAD surface dose measurements, along with the surface pressure data and the solar modulation factor, are analysed and fitted to empirical models which quantitatively demonstrate} how the long-term influences ([b] and [c]) are related to the measured dose rates. {Correspondingly we can estimate dose rate and dose equivalents under different solar modulations and different atmospheric conditions, thus allowing empirical predictions of the Martian surface radiation environment.

M82 X-1 is the brightest ultraluminous X-ray source in starburst galaxy M82 and is one of the best intermediate mass black hole candidates. Previous studies based on the Rossi X-ray Timing Explorer/Proportional Counter Array (RXTE/PCA) reported a regular X-ray flux modulation of M82 with a period of 62 days, and attributed this periodic modulation to M82 X-1. However, this modulation is not necessarily from M82 X-1 because RXTE/PCA has a very poor spatial resolution of ~1 degree. In this work, we analyzed 1000 days of monitoring data of M82 X-1 from the Swift/X-ray telescope (XRT), which has a much better spatial resolution than RXTE/PCA. The periodicity distribution map of M82 reveals that the 62-day periodicity is most likely not from M82 X-1, but from the summed contributions of several periodic X-ray sources 4 arcsec southeast of M82 X-1. However, Swift/XRT is not able to resolve those periodic sources and locate the precise origin of the periodicity of M82. Thus, more long-term observations with higher spatial resolution are required.

As part of our ongoing effort to investigate transit timing variations (TTVs) of known exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs due to the known properties of their host systems based on the respective discovery papers. During the past three years 46 transit observations were carried out, mostly using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to refine the systems orbital parameters. In all cases we found no hints for significant TTVs, or changes in the system parameters inclination, fractional stellar radius and planet to star radius ratio. However, comparing our results with those available in the literature shows that we can confirm the already published values.

Aims: To predict the cruise radiation environment related to future human missions to Mars, the correlation between solar modulation potential and the dose rate measured by the Radiation Assessment Detector (RAD) has been analyzed and empirical models have been employed to quantify this correlation. Methods: The instrument RAD, onboard Mars Science Laboratory's (MSL) rover Curiosity, measures a broad spectrum of energetic particles along with the radiation dose rate during the 253-day cruise phase as well as on the surface of Mars. With these first ever measurements inside a spacecraft from Earth to Mars, RAD observed the impulsive enhancement of dose rate during solar particle events as well as a gradual evolution of the galactic cosmic ray (GCR) induced radiation dose rate due to the modulation of the primary GCR flux by the solar magnetic field, which correlates with long-term solar activities and heliospheric rotation. Results: We analyzed the dependence of the dose rate measured by RAD on solar modulation potentials and estimated the dose rate and dose equivalent under different solar modulation conditions. These estimations help us to have approximate predictions of the cruise radiation environment, such as the accumulated dose equivalent associated with future human missions to Mars. Conclusions: The predicted dose equivalent rate during solar maximum conditions could be as low as one-fourth of the current RAD cruise measurement. However, future measurements during solar maximum and minimum periods are essential to validate our estimations.

In this paper, we revisit the solar system tests of f(R) gravity. When the Sun sits in a vacuum, the field f' is light, which leads to a metric different from the observations. We reobtain this result in a simpler way by directly focusing on the equations of motion for f(R) gravity in the Jordan frame. The discrepancy between the metric in the f(R) gravity and the observations can be alleviated by the chameleon mechanism. The implications from the chameleon mechanism on the functional form f(R) are discussed. Considering the analogy of the solar system tests to the false vacuum decay problem, the effective potentials in different cases are also explored. The combination of analytic and numerical approaches enables us to ascertain whether an f(R) model can pass the solar system tests or not.

The cosmic age is an important physical quantity in cosmology. Based on the radiometric method, a reliable lower limit of the cosmic age is derived to be $15.68\pm 1.95$ Gyr by using the $r$-process abundances inferred for the solar system and observations in metal-poor stars. This value is larger than the latest cosmic age $13.813\pm 0.058$ Gyr from Planck 2013 results, while they still agree with each other within the uncertainties. The uncertainty of $1.95$ Gyr mainly originates from the error on thorium abundance observed in metal-poor star CS 22892-052, so future high-precision abundance observations on CS 22892-052 are needed to understand this age deviation.

In this paper, we study the cosmological viability conditions, the phase-space dynamics, and the cosmological evolution of f(R) gravity. In contrast to most previous works in the literature, which analyzed the background dynamics of f(R) gravity by means of a dynamical system, we proceed by focusing on the equivalent scalar field description of the theory, which we believe is a more intuitive way of treating the problem. In order to study how the physical solutions evolve in f(R) cosmology, we explore the cosmological dynamics of a range of f(R) models, including models that yield a large hierarchy of scales and are singularity free. We present generic features of the phase-space dynamics in f(R) cosmology. We study the global structure of the phase space in f(R) gravity by compactifying the infinite phase space into a finite space via the Poincare transformation. On the expansion branch of the phase space, the constraint surface has a repeller and a de Sitter attractor; while on the contraction branch, the constraint surface has an attractor and a de Sitter repeller. Generally, the phase currents originate from the repeller and terminate at the corresponding attractor in each space. The trajectories between the repeller and the attractor in the presence of matter density are different from those in the vacuum case. The phase analysis techniques developed in this paper are very general, and can be applied to other similar dynamical systems.

The radial basis function (RBF) approach is applied in predicting nuclear masses for 8 widely used nuclear mass models, ranging from macroscopic-microscopic to microscopic types. A significantly improved accuracy in computing nuclear masses is obtained, and the corresponding rms deviations with respect to the known masses is reduced by up to 78%. Moreover, strong correlations are found between a target nucleus and the reference nuclei within about three unit in distance, which play critical roles in improving nuclear mass predictions. Based on the latest Weizs\"{a}cker-Skyrme mass model, the RBF approach can achieve an accuracy comparable with the extrapolation method used in atomic mass evaluation. In addition, the necessity of new high-precision experimental data to improve the mass predictions with the RBF approach is emphasized as well.

In this paper, we study how to construct f(R) gravity models. Cosmological observations and local gravity tests imply that a viable f(R) model should be very close to the LCDM model. We create procedures to construct viable LCDM-like f(R) models, and present multiple models of three types. The connections between some of these models are discussed. We also study the cosmological evolution of the LCDM-like f(R) gravity. Exact numerical integration can generate accurate cosmological evolution, but the numerical simulation in the early-universe stage is slow due to the oscillations of the field phi(= f') near the minimum of the effective potential. To avoid this problem, we take the minimum of the effective potential as an approximate solution for phi, and obtain the cosmological evolution. This approximate method describes the cosmological evolution well except in the late universe. Therefore, we use the approximate method in the early-universe evolution, and use the exact method in the late-universe one.

Self-consistent proton-neutron quasiparticle random phase approximation based on the spherical nonlinear point-coupling relativistic Hartree-Bogoliubov theory is established and used to investigate the $\beta^+$/EC-decay half-lives of neutron-deficient Ar, Ca, Ti, Fe, Ni, Zn, Cd, and Sn isotopes. The isoscalar proton-neutron pairing is found to play an important role in reducing the decay half-lives, which is consistent with the same mechanism in the $\beta$ decays of neutron-rich nuclei. The experimental $\beta^+$/EC-decay half-lives can be well reproduced by a universal isoscalar proton-neutron pairing strength.

We examine the electron-ion streaming instabilities driven by drift velocities of the order of the electron thermal velocity in a nonmagnetized plasma by one-dimensional electrostatic particle-in-cell code which adopts an ion-to-electron mass ratio of 1600. An initial state is set up where the ion bulk speed is zero while the electrons drift relative to ions, and where electrons are much hotter. We examine in detail four runs where drift velocity is systematically varied from lower than to larger than the electron thermal velocity. In all runs the Langmuir waves with Doppler-shifted frequencies dominate early on when streaming instabilities are too weak to discern. And then intense ion-acoustic waves or Buneman instabilities appear, which tend to be accompanied by localized electron and ion beams. Ion-acoustic modes and Buneman modes co-exist in the system when the initial drift velocity is just over the electron thermal speed. Beam modes are excited when the localized beams with large enough velocities appear. In the developed stage of instabilities, the direction in which density depressions propagate is always opposite to that of the localized ion beams. When the initial drift velocity is close to the electron thermal speed, categorizing the relevant instabilities is not easy, and one needs to examine in detail the wave dispersion diagrams at various stages of the evolution of the system.

We analyze electron flux maps based on RHESSI hard X-ray imaging spectroscopy data for a number of extended coronal loop flare events. For each event, we determine the variation of the characteristic loop length $L$ with electron energy $E$, and we fit this observed behavior with models that incorporate an extended acceleration region and an exterior "propagation" region, and which may include collisional modification of the accelerated electron spectrum inside the acceleration region. The models are characterized by two parameters: the plasma density $n$ in, and the longitudinal extent $L_0$ of, the acceleration region. Determination of the best-fit values of these parameters permits inference of the volume that encompasses the acceleration region and of the total number of particles within it. It is then straightforward to compute values for the emission filling factor and for the {\it specific acceleration rate} (electrons s$^{-1}$ per ambient electron above a chosen reference energy). For the 24 events studied, the range of inferred filling factors is consistent with a value of unity. The inferred mean value of the specific acceleration rate above $E_0=20$ keV is $\sim10^{-2}$ s$^{-1}$, with a 1$\sigma$ spread of about a half-order-of-magnitude above and below this value. We compare these values with the predictions of several models, including acceleration by large-scale, weak (sub-Dreicer) fields, by strong (super-Dreicer) electric fields in a reconnecting current sheet, and by stochastic acceleration processes.

We constructed a multi-fluid two-dimensional hydrodynamic model with detailed radiative transfer to depict the escape of particles. We found that the tidal forces supply significant accelerations and result in anisotropic winds. An important effect of the tidal forces is that it severely depresses the outflow of particles near the polar regions. As a consequence, most particles escape the surface of the planet from the regions of low-latitude. Comparing the tidal and non-tidal locking cases, we found that their optical depths are very different so that the flows also emerge with a different pattern. In the case of the non-tidal locking, the radial velocities at the base of the wind are higher than the meridional velocities. However, in the case of tidal locking, the meridional velocities dominate the flow at the base of the wind, and they can transfer effectively mass and energy from the day sides to the night sides. Further, we also found that the differences of the winds show middle extent at large radii. It expresses that the structure of the wind at the base can be changed by the two-dimensional radiative transfer due to large optical depths, but the extent is reduced with an increase in radius. Because the escape is depressed in the polar regions, the mass loss rate predicted by the non-tidal locking model, in the order of magnitude of $10^{10}$ g s$^{-1}$, is a factor of 2 lower than that predicted by one-dimensional hydrodynamic model. The results of tidal locking show that the mass loss rate is decreased a order of magnitude, only 4.3$\times$ $10^{9}$ g s$^{-1}$, due to large optical depths on the night side. We also found that the distributions of hydrogen atoms show clear variations from the day side to night side, thus the origin of the excess absorption in Ly $\alpha$ should be reexamined by using multi-dimensional hydrodynamic models.

In this paper, we explore an idea of having Newton's constant change its value depending on the curvature scale involved. Such modification leads to a particular scalar-tensor gravity theory, with the Lagrangian derived from renormalization group (RG) flow arguments. Several of the well-known f(R) modified gravity models have remarkably simple description in terms of the infrared renormalization group, but not the "designer" types in general. We find that de Sitter-like accelerated expansion can be generated even in the absence of cosmological constant term, entirely due to running of the Newton's constant. In hopes of tackling the problem of cosmological constant's smallness, we explore the flows which are capable of generating exponential hierarchy between infrared and ultraviolet scales, and investigate cosmological evolution in the models thus derived.

The Magnetic cloud boundary layer (BL) is a dynamic region formed by the interaction of the magnetic cloud (MC) and the ambient solar wind. In the present study, we comparatively investigate the proton and electron mean flux variations in the BL, in the interplanetary reconnection exhaust (RE) and across the MC-driven shock by using the Wind 3DP and MFI data from 1995 to 2006. In general, the proton flux has higher increments at lower energy bands compared with the ambient solar wind. Inside the BL, the core electron flux increases quasi-isotropically and the increments decrease monotonously with energy from ~30% (at 18 eV) to ~10% (at 70 eV); the suprathermal electron flux usually increases in either parallel or antiparallel direction; the correlation coefficient of electron flux variations in parallel and antiparallel directions changes sharply from ~0.8 below 70 eV to ~0 above 70 eV. Similar results are also found for RE. However, different phenomena are found across the shock where the electron flux variations first increase and then decrease with a peak increment (>200%) near 100 eV. The correlation coefficient of electron flux variations in parallel and antiparallel directions is always around 0.8. The similar behavior of flux variations in BL and RE suggests that reconnection may commonly occur in BL. Our work also implies that the strong energy dependence and direction selectivity of electron flux variations, which are previously thought to have no enough relevance to magnetic reconnection, could be considered as an important signature of solar wind reconnection in the statistical point of view.

Using {\em RHESSI} hard X-ray imaging spectroscopy observations, we analyze electron flux maps for a number of extended coronal loop flares. For each event, we fit a collisional model with an extended acceleration region to the observed variation of loop length with electron energy $E$, resulting in estimates of the plasma density in, and longitudinal extent of, the acceleration region. These quantities in turn allow inference of the number of particles within the acceleration region and hence the filling factor $f$ -- the ratio of the emitting volume to the volume that encompasses the emitting region(s). We obtain values of $f$ that lie mostly between 0.1 and 1.0; the (geometric) mean value is $f = 0.20 \times \div 3.9$, somewhat less than, but nevertheless consistent with, unity. Further, coupling information on the number of particles in the acceleration region with information on the total rate of acceleration of particles above a certain reference energy (obtained from spatially-integrated hard X-ray data) also allows inference of the specific acceleration rate (electron s$^{-1}$ per ambient electron above the chosen reference energy). We obtain a (geometric) mean value of the specific acceleration rate $\eta(20$ keV) $ = (6.0 \times / \div 3.4) \times 10^{-3}$ electrons s$^{-1}$ per ambient electron; this value has implications both for the global electrodynamics associated with replenishment of the acceleration region and for the nature of the particle acceleration process.

A statistical study is carried out on the photospheric magnetic nonpotentiality in solar active regions and its relationship with associated flares. We select 2173 photospheric vector magnetograms from 1106 active regions observed by the Solar Magnetic Field Telescope at Huairou Solar Observing Station, National Astronomical Observatories of China, in the period of 1988-2008, which covers most of the 22nd and 23rd solar cycles. We have computed the mean planar magnetic shear angle (\bar{\Delta\phi}), mean shear angle of the vector magnetic field (\bar{\Delta\psi}), mean absolute vertical current density (\bar{|J_{z}|}), mean absolute current helicity density (\bar{|h_{c}|}), absolute twist parameter (|\alpha_{av}|), mean free magnetic energy density (\bar{\rho_{free}}), effective distance of the longitudinal magnetic field (d_{E}), and modified effective distance (d_{Em}) of each photospheric vector magnetogram. Parameters \bar{|h_{c}|}, \bar{\rho_{free}}, and d_{Em} show higher correlation with the evolution of the solar cycle. The Pearson linear correlation coefficients between these three parameters and the yearly mean sunspot number are all larger than 0.59. Parameters \bar{\Delta\phi}, \bar{\Delta\psi}, \bar{|J_{z}|}, |\alpha_{av}|, and d_{E} show only weak correlations with the solar cycle, though the nonpotentiality and the complexity of active regions are greater in the activity maximum periods than in the minimum periods. All of the eight parameters show positive correlations with the flare productivity of active regions, and the combination of different nonpotentiality parameters may be effective in predicting the flaring probability of active regions.

In order to study the acceleration and propagation of bremsstrahlung-producing electrons in solar flares, we analyze the evolution of the flare loop size with respect to energy at a variety of times. A GOES M3.7 loop-structured flare starting around 23:55 on 2002 April 14 is studied in detail using \textit{Ramaty High Energy Solar Spectroscopic Imager} (\textit{RHESSI}) observations. We construct photon and mean-electron-flux maps in 2-keV energy bins by processing observationally-deduced photon and electron visibilities, respectively, through several image-processing methods: a visibility-based forward-fit (FWD) algorithm, a maximum entropy (MEM) procedure and the uv-smooth (UVS) approach. We estimate the sizes of elongated flares (i.e., the length and width of flaring loops) by calculating the second normalized moments of the intensity in any given map. Employing a collisional model with an extended acceleration region, we fit the loop lengths as a function of energy in both the photon and electron domains. The resulting fitting parameters allow us to estimate the extent of the acceleration region which is between $\sim 13 \rm{arcsec}$ and $\sim 19 \rm{arcsec}$. Both forward-fit and uv-smooth algorithms provide substantially similar results with a systematically better fit in the electron domain.The consistency of the estimates from these methods provides strong support that the model can reliably determine geometric parameters of the acceleration region. The acceleration region is estimated to be a substantial fraction ($\sim 1/2$) of the loop extent, indicating that this dense flaring loop incorporates both acceleration and transport of electrons, with concurrent thick-target bremsstrahlung emission.

GM Cep in the young (~4 Myr) open cluster Trumpler 37 has been known to be an abrupt variable and to have a circumstellar disk with very active accretion. Our monitoring observations in 2009-2011 revealed the star to show sporadic flare events, each with brightening of < 0.5 mag lasting for days. These brightening events, associated with a color change toward the blue, should originate from an increased accretion activity. Moreover, the star also underwent a brightness drop of ~1 mag lasting for about a month, during which the star became bluer when fainter. Such brightness drops seem to have a recurrence time scale of a year, as evidenced in our data and the photometric behavior of GM Cep over a century. Between consecutive drops, the star brightened gradually by about 1 mag and became blue at peak luminosity. We propose that the drop is caused by obscuration of the central star by an orbiting dust concentration. The UX Orionis type of activity in GM Cep therefore exemplifies the disk inhomogeneity process in transition between grain coagulation and planetesimal formation in a young circumstellar disk.

This paper on ASTROD I is based on our 2010 proposal submitted for the ESA call for class-M mission proposals, and is a sequel and an update to our previous paper [Experimental Astronomy 23 (2009) 491-527; designated as Paper I] which was based on our last proposal submitted for the 2007 ESA call. In this paper, we present our orbit selection with one Venus swing-by together with orbit simulation. In Paper I, our orbit choice is with two Venus swing-bys. The present choice takes shorter time (about 250 days) to reach the opposite side of the Sun. We also present a preliminary design of the optical bench, and elaborate on the solar physics goals with the radiation monitor payload. We discuss telescope size, trade-offs of drag-free sensitivities, thermal issues and present an outlook. ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system; and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals.

A multi-fluid model for an atomic hydrogen-proton mixture in the upper atmosphere of extrosolar planet is presented when the continuity and momentum equations of each component have been already solved with an energy equation. The particle number density, the temperature distribution and the structure of velocity can be found by means of the model. We chose two special objects, HD 209458b and HD 189733b, as discussion samples and the conclusion is that their mass loss rates predicted by the model are in accordance with those of observation. The most important physical process in coupling each component is charge exchange which tightly couples atomic hydrogen with protons. Most of the hydrogen escaping from hot Jupiters is protons, especially in young star-planet system. We found that the single-fluid model can describe the escape of particles when the mass loss rate is higher than a few times $10^{9}$ g/s while below $10^{9}$ g/s the multi-fluid model is more suitable for it due to the decoupling of particles. We found that the predicted mass loss rates of HD 189733b with the assumption of energy-limit are a factor of 10 larger than that calculated by our models due to the high ionization degree. For the ionized wind which is almost compose of protons, the assumption of energy-limit is no longer effective. We fitted the mass loss rates of the ionized wind as a function of $F_{UV}$ by calculating the variation of the mass loss rates with UV fluxes.

X-ray observations of solar flares routinely reveal an impulsive high-energy and a gradual low-energy emission component, whose relationship is one of the key issues of solar flare study. The gradual and impulsive emission components are believed to be associated with, respectively, the thermal and nonthermal components identified in spectral fitting. In this paper, a prominent about 50 second hard X-ray (HXR) pulse of a simple GOES class C7.5 flare on 20 February 2002 is used to study the association between high energy, non-thermal and impulsive evolution, and low energy, thermal and gradual evolution. We use regularized methods to obtain time derivatives of photon fluxes to quantify the time evolution as a function of photon energy, obtaining a break energy between impulsive and gradual behavior. These break energies are consistent with a constant value of about 11 keV in agreement with those found spectroscopically between thermal and non-thermal components, but the relative errors of the former are greater than 15% and much greater than the a few percent errors found from the spectral fitting. These errors only weakly depend on assuming an underlying spectral model for the photons, pointing to the current data being inadequate to reduce the uncertainties rather than there being a problem associated with an assumed model. The time derivative method is used to test for the presence of a 'pivot energy' in this flare. Although these pivot energies are marginally consistent with a constant value of about 9 keV, its values in the HXR rise phase appear to be lower than those in the decay phase.

Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t~1Gyr, and have been even included in the volutionary population synthesis (EPS) models by Han et al. (2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (2007) and those of the Yunnan group (Zhang et al. 2004, 005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is vailable upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0-Sd and Irr types of galaxies by using the package of Bruzual & Charlot (2003, BC03). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (2007), we find that the photo-zs fluctuate with the SDSS spectroscopic redshifts, while the Sa-Sc galaxies in the catalogue of Fukugita et al. (2007) are classified earlier as Burst-E galaxies.

With more and more exoplanets being detected, it is paid closer attention to whether there are lives outside solar system. We try to obtain habitable zones and the probability distribution of terrestrial planets in habitable zones around host stars. Using Eggleton's code, we calculate the evolution of stars with masses less than 4.00 \mo. We also use the fitting formulae of stellar luminosity and radius, the boundary flux of habitable zones, the distribution of semimajor axis and mass of planets and the initial mass function of stars. We obtain the luminosity and radius of stars with masses from 0.08 to 4.00 \mo, and calculate the habitable zones of host stars, affected by stellar effective temperature. We achieve the probability distribution of terrestrial planets in habitable zones around host stars. We also calculate that the number of terrestrial planets in habitable zones of host stars is 45.5 billion, and the number of terrestrial planets in habitable zones around K type stars is the most, in the Milky Way.

Ultraviolet radiation is a double-edged sword to life. If it is too strong, the terrestrial biological systems will be damaged. And if it is too weak, the synthesis of many biochemical compounds can not go along. We try to obtain the continuous ultraviolet habitable zones, and compare the ultraviolet habitable zones with the habitable zones of host stars. Using the boundary ultraviolet radiation of ultraviolet habitable zone, we calculate the ultraviolet habitable zones of host stars with masses from 0.08 to 4.00 \mo. For the host stars with effective temperatures lower than 4,600 K, the ultraviolet habitable zones are closer than the habitable zones. For the host stars with effective temperatures higher than 7,137 K, the ultraviolet habitable zones are farther than the habitable zones. For hot subdwarf as a host star, the distance of the ultraviolet habitable zone is about ten times more than that of the habitable zone, which is not suitable for life existence.

A star will become brighter and brighter with stellar evolution, and the distance of its habitable zone will become farther and farther. Some planets outside the habitable zone of a host star during the main sequence phase may enter the habitable zone of the host star during other evolutionary phases. A terrestrial planet within the habitable zone of its host star is generally thought to be suited to life existence. Furthermore, a rocky moon around a giant planet may be also suited to life survive, provided that the planet-moon system is within the habitable zone of its host star. Using Eggleton's code and the boundary flux of habitable zone, we calculate the habitable zone of our Solar after the main sequence phase. It is found that Mars' orbit and Jupiter's orbit will enter the habitable zone of Solar during the subgiant branch phase and the red giant branch phase, respectively. And the orbit of Saturn will enter the habitable zone of Solar during the He-burning phase for about 137 million years. Life is unlikely at any time on Saturn, as it is a giant gaseous planet. However, Titan, the rocky moon of Saturn, may be suitable for biological evolution and become another Earth during that time. For low-mass stars, there are similar habitable zones during the He-burning phase as our Solar, because there are similar core masses and luminosities for these stars during that phase.