We propose a concept of multiplexing lobster-eye (MuLE) optics to achieve significant reductions in the number of focal plane imagers in lobster-eye (LE) wide-field X-ray monitors. In the MuLE configuration, an LE mirror is divided into several segments and the X-rays reflected on each of these segments are focused on a single image sensor in a multiplexed configuration. If each LE segment assumes a different rotation angle, the azimuthal rotation angle of a cross-like image reconstructed from a point source by the LE optics identifies the specific segment that focuses the X-rays on the imager. With a focal length of 30 cm and LE segments with areas of 10 x 10 cm^2, ~1 sr of the sky can be covered with 36 LE segments and only four imagers (with total areas of 10 x 10 cm^2). A ray tracing simulation was performed to evaluate the nine-segment MuLE configuration. The simulation showed that the flux (0.5 to 2 keV) associated with the 5-sigma detection limit was ~2 x 10^-10 erg cm^-2 s^-1 (10 mCrab) for a transient with a duration of 100 s. The simulation also showed that the direction of the transient for flux in the range of 14 to 17 mCrab at 0.6 keV was determined correctly with 99.7% confidence limit. We conclude that the MuLE configuration can become an effective on-board device for small satellites for future X-ray wide-field transient monitoring.

A TeV spectral break in the total flux of cosmic-ray electrons and positrons (CREs) at which the spectral power index softens from $\sim 3$ to $\sim 4$ has been observed by H.E.S.S. and recently confirmed by DAMPE with a high significance of $6.6~\sigma$. Such an observation is apparently inconsistent with the data from other experiments such as Fermi-LAT, AMS-02 and CALET. We perform a global analysis to the latest CRE data including Fermi-LAT, AMS-02, CALET, DAMPE and H.E.S.S. with energy scale uncertainties taken into account to improve the consistency between the data sets. The fit result strongly favors the existence of the break at $\sim 1$ TeV with an even higher statistical significance of $13.3~\sigma$. In view of the tentative CRE break, we revisit a number of models of nearby sources, such as a single generic Pulsar Wind Nebula (PWN), known multiple PWNe from the ATNF catalog, and their combinations with either an additional Dark Matter (DM) component or a Supernova Remnant (SNR). We show that the CRE break at $\sim 1$ TeV, together with the known CR positron excess points towards the possibility that the nearby sources should be highly charge asymmetric. Among the models under consideration, the one with a PWN plus SNR is most favored by the current data. The favoured distance and age of the PWN and SNR sources are both within $0.6$ kpc and around $10^{5}$ yr respectively. Possible candidate sources include PSR J0954-5430, Vela and Monogem ring, etc. We find that for the models under consideration, the additional DM component is either unnecessary, or predicts too much photons in tension with the H.E.S.S. data of $\gamma$-rays from the direction of the Galactic Center. We also show that the current measurement of the anisotropies in the arrival direction of the CRE can be useful in determining the property of the sources.

We use our recent electric dipole moment (EDM) measurement data to constrain the possibility that the HfF$^+$ EDM oscillates in time due to interactions with candidate dark matter axion-like particles (ALPs). We employ a Bayesian analysis method which accounts for both the look-elsewhere effect and the uncertainties associated with stochastic density fluctuations in the ALP field. We find no evidence of an oscillating EDM over a range spanning from 27 nHz to 400 mHz, and we use this result to constrain the ALP-gluon coupling over the mass range $10^{-22}-10^{-15}$ eV. This is the first laboratory constraint on the ALP-gluon coupling in the $10^{-17}-10^{-15}$ eV range, and the first laboratory constraint to properly account for the stochastic nature of the ALP field.

{We investigate the atmospheric composition of the long period ($P_{\rm orb}=$~10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets.} % methods heading (mandatory) { We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines to ensure a robust detection of the $3 \sigma$-water spectrum. High resolution measurements were taken with VLT/ESPRESSO in the optical.} % results heading (mandatory) {Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$~bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 \sigma$. We cannot constrain the temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.}

Proton decay is a smoking gun signature of Grand Unified Theories (GUTs). Searches by Super-Kamiokande have resulted in stringent limits on the GUT symmetry breaking scale. The large-scale multipurpose neutrino experiments DUNE, Hyper-Kamiokande and JUNO will either discover proton decay or further push the symmetry breaking scale above $10^{16}$ GeV. Another possible observational consequence of GUTs is the formation of a cosmic string network produced during the breaking of the GUT to the Standard Model gauge group. The evolution of such a string network in the expanding Universe produces a stochastic background of gravitational waves which will be tested by a number of gravitational wave detectors over a wide frequency range. We demonstrate the non-trivial complementary between the observation of proton decay and gravitational waves produced from cosmic strings in determining $SO(10)$ GUT breaking chains. We show that such observations could exclude $SO(10)$ breaking via flipped $SU(5)\times U(1)$ or standard $SU(5)$, while breaking via a Pati-Salam intermediate symmetry, or standard $SU(5)\times U(1)$, may be favoured if a large separation of energy scales associated with proton decay and cosmic strings is indicated.

The last decade has seen a rapid development in asteroseismology thanks to the CoRoT and Kepler missions. With more detailed asteroseismic observations available, it is becoming possible to infer exactly how oscillations are driven and dissipated in solar-type stars. We have carried out three-dimensional (3D) stellar atmosphere simulations together with one-dimensional (1D) stellar structural models of key benchmark turn-off and subgiant stars to study this problem from a theoretical perspective. Mode excitation and damping rates are extracted from 3D and 1D stellar models based on analytical expressions. Mode velocity amplitudes are determined by the balance between stochastic excitation and linear damping, which then allows the estimation of the frequency of maximum oscillation power, $\nu_{\max}$, for the first time based on ab initio and parameter-free modelling. We have made detailed comparisons between our numerical results and observational data and achieved very encouraging agreement for all of our target stars. This opens the exciting prospect of using such realistic 3D hydrodynamical stellar models to predict solar-like oscillations across the HR-diagram, thereby enabling accurate estimates of stellar properties such as mass, radius and age.

We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range $500 \leq \ell <3000$, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

The presence of mixed modes makes subgiants excellent targets for asteroseismology, providing a probe for the internal structure of stars. Here we study 36 \emph{Kepler} subgiants with solar-like oscillations and report their oscillation mode parameters. We performed a so-called peakbagging exercise, i.e. estimating oscillation mode frequencies, linewidths, and amplitudes with a power spectrum model, fitted in the Bayesian framework and sampled with a Markov Chain Monte Carlo algorithm. The uncertainties of the mode frequencies have a median value of 0.180 $\mu$Hz. We obtained seismic parameters from the peakbagging, analysed their correlation with stellar parameters, and examined against scaling relations. The behaviour of seismic parameters (e.g. $\Delta\nu$, $\nu_{\rm max}$, $\epsilon_p$) is in general consistent with theoretical predictions. We presented the observational p--g diagrams: $\gamma_1$--$\Delta\nu$ for early subgiants and $\Delta\Pi_1$--$\Delta\nu$ for late subgiants, and demonstrate their capability to estimate stellar mass. We also found a $\log g$ dependence on the linewidths and a mass dependence on the oscillation amplitudes and the widths of oscillation excess. This sample will be valuable constraints for modelling stars and studying mode physics such as excitation and damping.

We present temporal variations of the Si IV line profiles at the flare ribbons in three solar flares observed by the Interface Region Imaging Spectrograph (IRIS). In the M1.1 flare on 2014 September 6 and the X1.6 flare on 2014 September 10, the Si IV line profiles evolve from wholly redshifted to red-wing enhanced with the flare development. However, in the B1.8 flare on 2016 December 2, the Si IV line profiles are wholly redshifted throughout the flare evolution. We fit the wholly redshifted line profiles with a single Gaussian function but the red-asymmetric ones with a double Gaussian function to deduce the corresponding Doppler velocities. In addition, we find that hard X-ray emission above 25 keV shows up in the two large flares, implying a nonthermal electron beam heating. In the microflare, there only appears weak hard X-ray emission up to 12 keV, indicative of a thermal heating mostly. We interpret the redshifts or red asymmetries of the Si IV line at the ribbons in the three flares as spectral manifestations of chromospheric condensation. We propose that whether the line appears to be wholly redshifted or red-asymmetric depends on the heating mechanisms and also on the propagation of the condensation.

The radiative decay of neutral fermions has been studied for decades but $CP$ violation induced within such a paradigm has evaded attention. $CP$ violation in these processes can produce an asymmetry between circularly polarised directions of the radiated photons and produces an important source of net circular polarisation in particle and astroparticle physics observables. The results presented in this work outlines the general connection between $CP$ violation and circular polarisation for both Dirac and Majorana fermions and can be used for any class of models that produce such radiative decays. The total $CP$ violation is calculated based on a widely studied Yukawa interaction considered in both active and sterile neutrino radiative decay scenarios as well as searches for dark matter via direct detection and collider signatures. Finally, the phenomenological implications of the formalism on keV sterile neutrino decay, leptogenesis-induced right-handed neutrino radiative decay and IceCube-driven heavy dark matter decay are discussed.

The IRIS Mg II k line serves as a very good tool to diagnose the heating processes in solar flares. Recent studies have shown that apart from the usual red asymmetries which are interpreted as the result of condensation downflows, this line could also show a blue-wing enhancement. To investigate how such a blue asymmetry is formed, we perform a grid of radiative hydrodynamic simulations and calculate the corresponding line profiles. We find that such a spectral feature is likely to originate from the upward plasma motion in the upper chromosphere. However, the formation region that is responsible for the blue-wing enhancement could be located in an evaporation region, in an upward moving blob, and even an upward moving condensation region. We discuss how the electron beam parameters affect these different dynamics of the atmosphere.

Multi-wavelength time-resolved observations of rotationally modulated variability from brown dwarfs and giant exoplanets are the most effective method for constraining their heterogeneous atmospheric structures. In a companion paper (Bowler et al. 2020), we reported the discovery of strong near-infrared variability in HST/WFC3/G141 light curves of the very red L-dwarf companion VHS J1256-1257b. In this paper, we present a follow-up 36-hr Spitzer/IRAC Channel 2 light curve together with an in-depth analysis of the HST and the Spitzer data. The combined dataset provides time-resolved light curves of VHS1256b sampling 1.1 to 4.5 $\mu$m. The Spitzer light curve is best-fit with a single sine wave with a period of $22.04\pm0.05$ hr and a peak-to-peak amplitude of $5.76\pm0.04$%. Combining the period with a previously measured projected rotational velocity ($v\sin i$), we find that VHS1256b is most consistent with equatorial viewing geometry. The HST/G141+Spitzer spectral energy distribution favors a $T_{\rm eff}$ of 1000~K, low surface gravity model with disequilibrium chemistry. The spectral variability of VHS1256b is consistent with predictions from partly cloudy models, suggesting heterogeneous clouds are the dominant source of the observed modulations. We find evidence at the $3.3\sigma$ level for amplitude variations within the 1.67$\mu$m CH$_{4}$ band, which is the first such detection for a variable L-dwarf. We compare the HST/G141 time-resolved spectra of three red L-dwarfs with high-amplitude near-infrared rotational modulations and find that although their time-averaged spectra are similar, their spectroscopic variabilities exhibit notable differences. This diversity reinforces the advantage of time-resolved spectroscopic observations for understanding the atmospheres of brown dwarfs and directly-imaged exoplanets.

Rotationally-modulated variability of brown dwarfs and giant planets provides unique information about their surface brightness inhomogeneities, atmospheric circulation, cloud evolution, vertical atmospheric structure, and rotational angular momentum. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared time-series spectroscopic observations of three companions with masses in or near the planetary regime: VHS J125601.92-125723.9 b, GSC 6214-210 B, and ROXs 42 B b. VHS J1256-1257 b exhibits strong total intensity and spectral variability with a brightness difference of 19.3% between 1.1-1.7 $\mu$m over 8.5 hours and even higher variability at the 24.7% level at 1.27 $\mu$m. The light curve of VHS J1256-1257 b continues to rise at the end of the observing sequence so these values represent lower limits on the full variability amplitude at this epoch. This observed variability rivals (and may surpass) the most variable brown dwarf currently known, 2MASS J21392676+0220226. The implied rotation period of VHS J1256-1257 b is $\approx$21-24 hr assuming sinusoidal modulations, which is unusually long for substellar objects. No significant variability is evident in the light curves of GSC 6214-210 B ($<$1.2%) and ROXs 42 B b ($<$15.6%). With a spectral type of L7, an especially red spectrum, and a young age, VHS J1256-1257 b reinforces emerging patterns between high variability amplitude, low surface gravity, and evolutionary phase near the L/T transition.

The effective vacuum energy density contributed by the non-trivial contortion distribution and the bare vacuum energy density can be viewed as the energy density of the auxiliary quintessence field potential. We find that the negative bare vacuum energy density from string landscape leads to a monotonically decreasing quintessence potential while the positive one from swampland leads to the meta stable or stable de Sitter like potential. Moreover, the non-trivial Brans-Dicke like coupling between quintessence field and gravitation field is necessary in the latter case.

HD106906b is an ~11$M_{\mathrm{Jup}}$, ~15Myr old directly-imaged exoplanet orbiting at an extremely large distance from its host star. The wide separation (7.11 arcsec) between HD106906b and its host star greatly reduces the difficulty in direct-imaging observations, making it one of the most favorable directly-imaged exoplanets for detailed characterization. In this paper, we present HST/WFC3/IR time-resolved observations of HD106906b in the F127M, F139M, and F153M bands. We have achieved ~1% precision in the lightcurves in all three bands. The F127M lightcurve demonstrates marginally-detectable ($2.7\sigma$ significance) variability with a best-fitting period of 4 hr, while the lightcurves in the other two bands are consistent with flat lines. We construct primary-subtracted deep images and use these images to exclude additional companions to HD106906 that are more massive than 4$M_{\mathrm{Jup}}$ and locate at projected distances of more than ~500 au. We measure the astrometry of HD106906b in two HST/WFC3 epochs and achieve precisions better than 2.5 mas. The position angle and separation measurements do not deviate from those in the 2004 HST/ACS/HRC images for more than $1\sigma$ uncertainty. We provide the HST/WFC3 astrometric results for 25 background stars that can be used as reference sources in future precision astrometry studies. Our observations also provide the first 1.4-micron water band photometric measurement for HD106906b. HD106906b's spectral energy distribution and the best-fitting BT-Settl model have an inconsistency in the 1.4-micron water absorption band, which highlights the challenges in modeling atmospheres of young planetary-mass objects.

We introduce a model in which the genesis of dark matter and neutrino masses is associated with a first order phase transition of a scalar singlet field. During the phase transition a source right-handed neutrino acquires a spacetime-dependent mass dynamically, a small fraction of which is converted via resonant oscillations into a very weakly mixed dark right-handed neutrino with the observed dark matter relic abundance. Neutrino masses are generated via a traditional two right-handed neutrino type-I seesaw between a third right-handed neutrino and the source neutrino. The gravitational waves produced during the phase transition have a peak frequency that increases with the dark matter mass, and are detectable at future gravitational wave interferometers for dark matter masses in the 1 MeV - 20 GeV range. For source right-handed neutrinos heavier than a GeV, successful leptogenesis is also possible.

The polarisation of sunlight after scattering off the atmosphere was first described by Chandrasekhar using a geometrical description of Rayleigh interactions. Kosowsky later extended Chandrasekhar's formalism by using Quantum Field Theory (QFT) to describe the polarisation of the Cosmological Microwave Background radiation. Here we focus on a case that is rarely discussed in the literature, namely the polarisation of high energy radiation after scattering off particles. After demonstrating why the geometrical and low energy QFT approaches fail in this case, we establish the transport formalism that allows to describe the change of polarisation of high energy photons when they propagate through space or the atmosphere. We primarily focus on Compton interactions but our approach is general enough to describe e.g. the scattering of high energy photons off new particles or through new interactions. Finally we determine the conditions for a circularly polarised $\gamma$--ray signal to keep the same level of circular polarisation as it propagates through its environment.

Standard stellar evolution model predicts a severe depletion of lithium (Li) abundance during the first dredge up process (FDU). Yet a small fraction of giant stars are still found to preserve a considerable amount of Li in their atmospheres after the FDU. Those giants are usually identified as Li-rich by a widely used criterion, A(Li) > 1.5 dex. A large number of works dedicated to searching for investigating this minority of the giant family, and the amount of Li-rich giants, has been largely expanded on, especially in the era of big data. In this paper, we present a catalog of Li-rich giants found from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey with Li abundances derived from a template-matching method developed for LAMOST low-resolution spectra.

Among the greatest challenges in understanding ultra-cool brown dwarf and exoplanet atmospheres is the evolution of cloud structure as a function of temperature and gravity. In this study, we present the rotational modulations of GU Psc b -- a rare mid-T spectral type planetary-mass companion at the end of the L/T spectral type transition. Based on the HST/WFC3 1.1-1.67$\rm\, \mu m$ time-series spectra, we observe a quasi-sinusoidal light curve with a peak-to-trough flux variation of 2.7 % and a minimum period of eight hours. The rotation-modulated spectral variations are weakly wavelength-dependent, or largely gray between 1.1-1.67$\rm\,\mu$m. The gray modulations indicate that heterogeneous clouds are present in the photosphere of this low-gravity mid-T dwarf. We place the color and brightness variations of GU Psc b in the context of rotational modulations reported for mid-L to late-T dwarfs. Based on these observations, we report a tentative trend: mid-to-late T dwarfs become slightly redder in $J-H$ color with increasing $J$-band brightness, while L dwarfs become slightly bluer with increasing brightness. If this trend is verified with more T-dwarf samples, it suggests that in addition to the mostly gray modulations, there is a second-order spectral-type dependence on the nature of rotational modulations.

This paper was withdrawed from the ApJ after the comments from the referee, please be carefully.

We test the possible dipole anisotropy of a Finslerian cosmological model and other three dipole-modulated cosmological models, i.e., the dipole-modulated $\rm{\Lambda}$CDM, $w$CDM and Chevallier--Polarski--Linder (CPL) model by using the recently released Pantheon sample of SNe Ia. The Markov chain Monte Carlo (MCMC) method is used to explore the whole parameter space. We find that the dipole anisotropy is very weak in all cosmological models used. Although the dipole amplitudes of four cosmological models are consistent with zero within $1\sigma$ uncertainty, the dipole directions are close to the axial direction to the plane of SDSS subsample among Pantheon. It may imply that the weak dipole anisotropy in the Pantheon sample originates from the inhomogeneous distribution of the SDSS subsample. More homogeneous distribution of SNe Ia is necessary to constrain the cosmic anisotropy.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

The Medium Energy X-ray telescope (ME) is one of the three main telescopes on board the Insight Hard X-ray Modulation Telescope (Insight-HXMT) astronomy satellite. ME contains 1728 pixels of Si-PIN detectors sensitive in 5-30 keV with a total geometrical area of 952 cm2. Application Specific Integrated Circuit (ASIC) chips, VA32TA6, is used to achieve low power consumption and low readout noise. The collimators define three kinds of field of views (FOVs) for the telescope, 1{\deg}{\times}4{\deg}, 4{\deg}{\times}4{\deg}, and blocked ones. Combination of such FOVs can be used to estimate the in-orbit X-ray and particle background components. The energy resolution of ME is ~3 keV at 17.8 keV (FWHM) and the time resolution is 255 {\mu}s. In this paper, we introduce the design and performance of ME.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

Using the very recently reported mass $2.14^{+0.10}_{-0.09}M_\odot$ of PSR J0740+6620 together with the data of finite nuclei and the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities from flow data in heavy-ion collisions, we show that the symmetry energy $E_{\rm sym}(n)$ cannot be supersoft so that it becomes negative at suprasaturation densities in neutron stars (NSs) and thus may make the NS have a pure neutron matter core. This is in contrast to the fact that using the mass $2.01 \pm 0.04 M_\odot$ of PSR J0348+0432 as the NS maximum mass cannot rule out the supersoft high-density $E_{\rm sym}(n)$. Furthermore, we find the stiffer high-density $E_{\rm sym}(n)$ based on the existence of $2.14M_\odot$ NSs leads to a strong constraint of $\Lambda_{1.4} \ge 348^{+88}_{-51}$ for the dimensionless tidal deformability of the canonical $1.4M_\odot$ NS.

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.

We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs.

The amplitude of solar p-mode oscillations is governed by stochastic excitation and mode damping, both of which take place in the surface convection zone. However, the time-dependent, turbulent nature of convection makes it difficult to self-consistently study excitation and damping processes through the use of traditional one-dimensional hydrostatic models. To this end, we carried out \textit{ab initio} three-dimensional, hydrodynamical numerical simulations of the solar atmosphere to investigate how p-modes are driven and dissipated in the Sun. The description of surface convection in the simulations is free from the tuneable parameters typically adopted in traditional one-dimensional models. Mode excitation and damping rates are computed based on analytical expressions whose ingredients are evaluated directly from the three-dimensional model. With excitation and damping rates both available, we estimate the theoretical oscillation amplitude and frequency of maximum power, $\nu_{\max}$, for the Sun. We compare our numerical results with helioseismic observations, finding encouraging agreement between the two. The numerical method presented here provides a novel way to investigate the physical processes responsible for mode driving and damping, and should be valid for all solar-type oscillating stars.

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.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

We employ the hemisphere comparison (HC) method and the dipole fitting (DF) method to investigate the cosmic anisotropy in the recently released Pantheon sample of type Ia supernovae (SNe Ia) and five combinations among Pantheon. For the HC method, we find the maximum anisotropy level in the full Pantheon sample is $\mathrm{AL}_{max}=0.361\pm0.070$ and corresponding direction $(l,b)=({123.05^{\circ}}^{+11.25^{\circ}}_{-4.22^{\circ}}, {4.78^{\circ}}^{+1.80^{\circ}}_{-8.36^{\circ}})$. A robust check shows the statistical significance of maximum anisotropy level is about $2.1\sigma$. We also find that the Low-$z$ and SNLS subsamples have decisive impact on the overall anisotropy while other three subsamples have little impact. Moreover, the anisotropy level map significantly rely on the inhomogeneous distribution of SNe Ia in the sky. For the DF method, we find the dipole anisotropy in the Pantheon sample is very weak. The dipole magnitude is constrained to be less than $1.16\times10^{-3}$ at $95\%$ confidence level. However, the dipole direction is well inferred by MCMC method and it points towards $(l,b)=({306.00^{\circ}}^{+82.95^{\circ}}_{-125.01^{\circ}}, {-34.20^{\circ}}^{+16.82^{\circ}}_{-54.93^{\circ}})$. This direction is very close to the axial direction to the plane of SDSS subsample. It may imply that SDSS subsample is the decisive part to the dipole anisotropy in the full Pantheon sample. All these facts imply that the cosmic anisotropy found in Pantheon sample significantly rely on the inhomogeneous distribution of SNe Ia in the sky. More homogeneous distribution of SNe Ia is necessary to search for a more convincing cosmic anisotropy.

While the equation of state (EOS) of symmetric nuclear matter (SNM) at suprasaturation densities has been relatively well constrained from heavy-ion collisions, the EOS of high-density neutron-rich matter is still largely uncertain due to the poorly known high-density behavior of the symmetry energy. Using the constraints on the EOS of SNM at suprasaturation densities from heavy-ion collisions together with the data of finite nuclei and the existence of $2M_\odot$ neutron stars from electromagnetic (EM) observations, we show that the high-density symmetry energy cannot be too soft, which leads to lower bounds on dimensionless tidal deformability of $\Lambda_{1.4} \ge 193$ and radius of $R_{1.4} \ge 11.1$ km for $1.4M_\odot$ neutron star. Furthermore, we find that the recent constraint of $\Lambda_{1.4} \le 580$ from the gravitational wave signal GW170817 detected from the binary neutron star merger by the LIGO and Virgo Collaborations rules out too stiff high-density symmetry energy, leading to an upper limit of $R_{1.4} \le 13.3$ km. All these terrestrial nuclear experiments and astrophysical observations based on strong, EM and gravitational measurements together put stringent constraints on the high-density symmetry energy and the EOS of SNM, pure neutron matter and neutron star matter.

Correlation between the variations of quasar absorption lines and the ionizing continuum have been recently confirmed in systematic studies. However, no convincing individual case is reported. We present a statistical analysis of the variable C iv broad absorption line (BAL) in the quasar SDSS J141007.74+541203.3, which have been observed with 44 epochs by the Sloan Digital Sky Survey Data Release 14. Grier et al. (2015) has recently concluded that the most likely cause of the variability of the BAL in SDSS J141007.74+541203.3 is a rapid response to changes in the incident ionizing continuum. In this paper, we confirm the anticorrelation between the equivalent width of BALs and the flux of the continuum based on the spectra of this quasar that show significant variations, which serve as another independent evidence for the view of Grier et al. (2015).

The discovery of Li-rich giant has cast a new challenge for the standard stellar evolution models, and to resolve this issue, the number of this type object has been rapidly increased because of the development of worldwide surveys these days. Taking advantage of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, 44 newly Li-rich giants are reported, which are confirmed with high-resolution observations. Based on the high-resolution and high signal-to-noise spectra, we derived the atmospheric parameters and elemental abundances with the spectral synthesis method. We performed a detailed analysis of their evolutionary stages, infrared excess, projected rotational velocity (v sin i), and the stellar population. We find that (1) The Li-rich giants concentrate at the evolutionary status of the red giant branch bump, red clump, and asymptotic giant branch; (2) Three of them are fast rotators and none exhibit infrared excess. Our results imply that the origins of Li enrichment are most likely to be associated with the extra mixing in the stellar interior, and the external sources maybe only make a minor contribution. Moreover, various Li-rich episodes take place at different evolutionary stages.

Measurements of photometric variability at different wavelengths provide insights into the vertical cloud structure of brown dwarfs and planetary-mass objects. In seven Hubble Space Telescope consecutive orbits, spanning $\sim$10 h of observing time}, we obtained time-resolved spectroscopy of the planetary-mass T8-dwarf Ross 458C using the near-infrared Wide Field Camera 3. We found spectrophotometric variability with a peak-to-peak signal of 2.62$\pm$0.02 % (in the 1.10-1.60~$\mu$m white light curve). Using three different methods, we estimated a rotational period of 6.75$\pm$1.58~h for the white light curve, and similar periods for narrow $J$- and $H$- band light curves. Sine wave fits to the narrow $J$- and $H$-band light curves suggest a tentative phase shift between the light curves with wavelength when we allow different periods between both light curves. If confirmed, this phase shift may be similar to the phase shift detected earlier for the T6.5 spectral type 2MASS J22282889-310262. We find that, in contrast with 2M2228, the variability of Ross~458C shows evidence for a {color trend} within the narrow $J$-band, but gray variations in the narrow $H$-band. The spectral time-resolved variability of Ross 458C might be potentially due to heterogeneous sulfide clouds in the atmosphere of the object. Our discovery extends the study of spectral modulations of condensate clouds to the coolest T dwarfs, planetary-mass companions.

Milgrom's modified Newtonian dynamics (MOND) can explain well the mass discrepancy problem in galaxy without invoking dark matter. The MOND theory predicts a universal constant acceleration scale in galaxy, below which the Newtonian dynamics is no longer valid. Recently, Rodrigues et al. analyzed 193 high-quality disk galaxies by using Bayesian inference and they concluded that the probability of existence of a fundamental acceleration scale is essentially 0. In this paper, we use the same method to revisit the same question: Is there a fundamental acceleration scale in galaxies? We fit each galaxy rotation curve of 175 SPARC galaxies with Gaussian priors on galaxy parameters and a variable acceleration scale. After marginalizing over the nuisance parameters, we find that the marginalized posterior distributions of the acceleration scale become broad. The incompatibility between the global best fit and the marginalized posterior distributions of the acceleration scale is greatly reduced. However, there still exists evidence that rejects MOND as a fundamental theory. More accurate observations are needed to exclude or confirm the existence of a fundamental acceleration scale in galaxies.

Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope.

The Earths core formation process has decisive effect in the chemical differentiation between the Earths core and its mantle. Here, we propose a new core formation model which is caused by a special giant impact. This model suggests that the impactors core can be kept intact by its own sticky mantle under appropriate impacting conditions and let it merge into the targets core without contact with the targets mantle. We call this special giant impact that caused the new core formation mode as glue ball impact model (GBI). By simulating hundreds of giant impacts with the sizes from planetesimals to planets, the conditions that can lead to GBI have been found out. If with small impact angle (i.e., less than 20 degree), small impact velocity and small impactors mass but larger than 0.07 Mearth, there is a good chance to produce a GBI at the final stage of the Earths accretion. We find that it will be much easier to have GBIs at the late stage of the Earths accretion rather than at the early stage of it. The GBI model will pose a great challenge to many problems between the equilibrium of Earths core and mantle. It provides an additional source for the excess of highly siderophile elements in the Earths mantle and also brings excessive lithophile elements to the Earths core. The GBI model may shed light on the study of Moon-formation and chemical differentiations of the pro-Earth.

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

Bayesian atmospheric retrieval tools can place constraints on the properties of brown dwarfs and hot Jupiters atmospheres. To fully exploit these methods, high signal-to-noise spectral libraries with well-understood uncertainties are essential. We present a high signal-to-noise spectral library (1.10-1.69 microns) of the thermal emission of 76 brown dwarfs and hot Jupiters. All our spectra have been acquired with the Hubble Space Telescope's Wide Field Camera 3 instrument and its G141 grism. The near-infrared spectral types of these objects range from L4 to Y1. Eight of our targets have estimated masses below the deuterium-burning limit. We analyze the database to identify peculiar objects and/or multiple systems, concluding that this sample includes two very-low-surface-gravity objects and five intermediate-surface-gravity objects. In addition, spectral indices designed to search for composite atmosphere brown dwarfs, indicate that eight objects in our sample are strong candidates to have such atmospheres. None of these objects are overluminous, thus their composite atmospheres are unlikely a companion-induced artifact. Five of the eight confirmed candidates have been reported as photometrically variable, suggesting that composite atmospheric indices are useful in identifying brown dwarfs with strongly heterogeneous cloud covers. We compare hot Jupiters and brown dwarfs in a near-infrared color-magnitude diagram. We confirm that the coldest hot Jupiters in our sample have spectra similar to mid-L dwarfs, and the hottest hot Jupiters have spectra similar to those of M-dwarfs. Our sample provides a uniform dataset of a broad range of ultracool atmospheres, allowing large-scale, comparative studies, and providing a HST legacy spectral library.

Supercritical accretion onto compact objects powers a massive wind that is optically thick and Eddington-limited. If most of the hard X-rays from the central disk are obscured by the wind, the source will display a blackbody-like spectrum with a luminosity scaled with the mass of the compact object. From the Chandra archive of nearby galaxies, we selected a sample of luminous and very soft sources and excluded contaminations from foreground objects and supernova remnants. They are found to be preferentially associated with late-type galaxies. The majority of sources in our sample are either too hot or too luminous to be explained by nuclear burning on the surface of white dwarfs, and are argued to be powered by accretion. The most likely explanation is that they are due to emission from the photosphere of a wind driven by supercritical accretion onto compact objects. Their blackbody luminosity ranges from ~10^37 to nearly 10^40 erg/s, indicative of the presence of both neutron stars and stellar-mass black holes. The blackbody luminosity also shows a possible bimodal distribution, albeit at low significance, peaked around the Eddington limit for neutron stars and stellar-mass black holes, respectively. If this can be confirmed, it will be smoking gun evidence that supercritical accretion powers thick winds. Based on a wind model, the inferred mass accretion rate of these objects is around a few hundred times the Eddington rate, suggesting that they may be intermediate between the canonical ultraluminous X-ray sources and SS 433 in terms of the accretion rate.

Based on an extended Skyrme interaction that includes the terms in relative momenta up to sixth order, corresponding to the so-called Skyrme pseudopotential up to next-to-next-to-next-to leading order (N3LO), we derive the expressions of Hamiltonian density and single nucleon potential under general non-equilibrium conditions which can be applied in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei. While the conventional Skyrme interactions, which include the terms in relative momenta up to second order, predict an incorrect behavior as a function of energy for nucleon optical potential in nuclear matter, the present extended N3LO Skyrme interaction can give a nice description for the empirical nucleon optical potential. We also construct three interaction sets with different high-density behaviors of the symmetry energy, by fitting both the empirical nucleon optical potential up to energy of $1$ GeV and the empirical properties of isospin asymmetric nuclear matter. These extended N3LO Skyrme interactions will be useful in transport model simulations of heavy-ion collisions induced by neutron-rich nuclei at intermediate and high energies, and they can also be useful in nuclear structure studies within the mean-field model.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

Understanding morphological types of galaxies is a key parameter for studying their formation and evolution. Neural networks that have been used previously for galaxy morphology classification have some disadvantages, such as not being invariant under rotation. In this work, we studied the performance of Capsule Network, a recently introduced neural network architecture that is rotationally invariant and spatially aware, on the task of galaxy morphology classification. We designed two evaluation scenarios based on the answers from the question tree in the Galaxy Zoo project. In the first scenario, we used Capsule Network for regression and predicted probabilities for all of the questions. In the second scenario, we chose the answer to the first morphology question that had the highest user agreement as the class of the object and trained a Capsule Network classifier, where we also reconstructed galaxy images. We achieved promising results in both of these scenarios. Automated approaches such as the one introduced here will greatly decrease the workload of astronomers and will play a critical role in the upcoming large sky surveys.

About one percent of giants\upcite{Brown1989} are detected to have anomalously high lithium (Li) abundances in their atmospheres, conflicting directly with the prediction of the standard stellar evolution models\upcite{Iben1967}, and making the production and evolution of Li more intriguing, not only in the sense of the Big Bang nucleosynthesis\upcite{Cyburt2016,Spite1982} or the Galactic medium\upcite{Tajitsu2015}, but also the evolution of stars. Decades of efforts have been put into explaining why such outliers exist\upcite{Sackmann1999, Denissenkov2004, Charbonnel2010}, yet the origins of Li-rich giants are still being debated. Here we report the discovery of the most Li-rich giant known to date, with a super-high Li abundance of 4.51. This rare phenomenon was snapshotted together with another short-term event that the star is experiencing its luminosity bump on the red giant branch. Such high Li abundance indicates that the star might be at the very beginning of its Li-rich phase, which provides a great opportunity to investigate the origin and evolution of Li in the Galaxy. A detailed nuclear simulation is presented with up-to-date reaction rates to recreate the Li enriching process in this star. Our results provide tight constraints on both observational and theoretical points of view, suggesting that low-mass giants can produce Li inside themselves to a super high level via $^{7}$Be transportation during the red giant phase.

The seven approximately Earth-sized transiting planets in the \object{TRAPPIST-1} system provide a unique opportunity to explore habitable zone and non-habitable zone small planets within the same system. Its habitable zone exoplanets -- due to their favorable transit depths -- are also worlds for which atmospheric transmission spectroscopy is within reach with the Hubble Space Telescope (HST) and with the James Webb Space Telescope (JWST). We present here an independent reduction and analysis of two \textit{HST} Wide Field Camera 3 (WFC3) near-infrared transit spectroscopy datasets for six planets (b through g). Utilizing our physically-motivated detector charge trap correction and a custom cosmic ray correction routine, we confirm the general shape of the transmission spectra presented by \textbf{\citet{deWit2016, deWit2018}}. Our data reduction approach leads to a 25\% increase in the usable data and reduces the risk of confusing astrophysical brightness variations (e.g., flares) with instrumental systematics. No prominent absorption features are detected in any individual planet's transmission spectra; by contrast, the combined spectrum of the planets shows a suggestive decrease around 1.4\,$\micron$ similar to an inverted water absorption feature. Including transit depths from \textit{K2}, the SPECULOOS-South Observatory, and \textit{Spitzer}, we find that the complete transmission spectrum is fully consistent with stellar contamination owing to the transit light source effect. These spectra demonstrate how stellar contamination can overwhelm planetary absorption features in low-resolution exoplanet transit spectra obtained by \textit{HST} and \textit{JWST} and also highlight the challenges in combining multi epoch observations for planets around rapidly rotating spotted stars.

The eXTP (enhanced X-ray Timing and Polarimetry) mission is a major project of the Chinese Academy of Sciences (CAS) and China National Space Administration (CNSA) currently performing an extended phase A study and proposed for a launch by 2025 in a low-earth orbit. The eXTP scientific payload envisages a suite of instruments (Spectroscopy Focusing Array, Polarimetry Focusing Array, Large Area Detector and Wide Field Monitor) offering unprecedented simultaneous wide-band X-ray spectral, timing and polarimetry sensitivity. A large European consortium is contributing to the eXTP study and it is expected to provide key hardware elements, including a Large Area Detector (LAD). The LAD instrument for eXTP is based on the design originally proposed for the LOFT mission within the ESA context. The eXTP/LAD envisages a deployed 3.4 m2 effective area in the 2-30 keV energy range, achieved through the technology of the large-area Silicon Drift Detectors - offering a spectral resolution of up to 200 eV FWHM at 6 keV - and of capillary plate collimators - limiting the field of view to about 1 degree. In this paper we provide an overview of the LAD instrument design, including new elements with respect to the earlier LOFT configuration.

The Maunakea Spectroscopic Explorer (MSE) project will transform the CFHT 3.6m optical telescope to a 10m class dedicated multi-object spectroscopic facility, with an ability to measure thousands of objects with three spectral resolution modes respectively low resolution of R~3,000, moderate resolution of R~6,000 and high resolution of R~40,000. Two identical multi-object high resolution spectrographs are expected to simultaneously produce 1084 spectra with high resolution of 40,000 at Blue (401-416nm) and Green (472-489nm) channels, and 20,000 at Red (626-674nm) channel. At the Conceptual Design Phase (CoDP), different optical schemes were proposed to meet the challenging requirements, especially a unique design with a novel transmission image slicer array, and another conventional design with oversize Volume Phase Holographic (VPH) gratings. It became clear during the CoDP that both designs presented problems of complexity or feasibility of manufacture, especially high line density disperser (general name for all kinds of grating, grism, prism). At the present, a new design scheme is proposed for investigating the optimal way to reduce technical risk and get more reliable estimation of cost and timescale. It contains new dispersers, F/2 fast collimator and so on. Therein, the disperser takes advantage of a special grism and a prism to reduce line density on grating surface, keep wide opening angle of optical path, and get the similar spectrum layout in all three spectral channels. For the fast collimator, it carefully compares on-axis and off-axis designs in throughput, interface to fiber assembly and technical risks. The current progress is more competitive and credible than the previous design, but it also indicates more challenging work will be done to improve its accessibility in engineering.

Bi-site time-series photometric and high-resolution spectroscopic observations were made for the double-mode high-amplitude $\delta$ Scuti star VX\ Hya. The fundamental frequency $f_{0}=4.4763\ \rm{c\ days^{-1}}$, the first overtone $f_{1}=5.7897\ \rm{c\ days^{-1}}$ and 23 harmonics and linear combinations of $f_{0}$ and $f_{1}$ are detected by pulsation analysis. From the spectroscopic data, we get $\rm{[Fe/H] = -0.2\pm0.1\ dex}$. The period change rate of the fundamental mode is obtained by using the Fourier-phase diagram method, which gives the value of $(1/P_{0})(dP_{0}/dt)=(1.81\pm0.09) \times 10^{-7}\ \rm{yr}^{-1}$. With these results from the observations, we perform theoretical explorations with the stellar evolution code MESA, and constrain the models by fitting $f_{0}$, $f_{1}$, and $(1/P_{0})(dP_{0}/dt)$ within $3\sigma$ deviations. The results show that the period change of VX\ Hya could be ascribed to the evolutionary effect. The stellar parameters of VX\ Hya could be derived as: the mass $2.385\pm0.025\ M_{\odot}$, the luminosity $\log(L/L_{\odot})=1.93\pm0.02$ and the age $(4.43\pm0.13)\times 10^8$ years. VX\ Hya is found to locate at the post-main-sequence stage with a helium core and a hydrogen-burning shell on the H${-}$R diagram.

We propose a framework of Finsler space-time to explain the observed parity asymmetry and the power deficit in the low-$\ell$ ($2\leqslant \ell \leqslant 29$) multipole range of cosmic microwave background (CMB) temperature anisotropies. In the $3+1$ dimensional space-time, the three-dimensional space is described by a Randers-Finsler space, which is spatially irreversible, inducing the parity asymmetry in the CMB angular power spectrum. We estimate the constraints on the two parameters introduced by Finsler space-time via analyzing the low-$\ell$ angular power spectrum in PLANCK 2015 CMB temperature data. We see that the low-$\ell$ power suppression in the CMB temperature anisotropies can also be resolved in this scenario. Our study shows that the two low-$\ell$ anomalies, i.e., parity asymmetry and power deficit, may have a common origin.

Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory.

We report a possible dipole anisotropy on acceleration scale $g_{\dag}$ with 147 rotationally supported galaxies in local Universe. It is found that a monopole and dipole correction for the radial acceleration relation can better describe the SPARC data set. The monopole term is negligible but the dipole magnitude is significant. It is also found that the dipole correction is mostly induced by the anisotropy on the acceleration scale. The magnitude of $\hat{g}_{\dag}$-dipole reaches up to $0.25\pm0.04$, and its direction is aligned to $(l,b) = (171.30^{\circ}\pm7.18^{\circ}, -15.41^{\circ}\pm4.87^{\circ})$, which is very close to the maximum anisotropy direction from the hemisphere comparison method. Furthermore, robust check shows that the dipole anisotropy couldn't be reproduced by isotropic mock data set. However, it is still premature to claim that the Universe is anisotropic due to the small data samples and uncertainty in the current observations.

Helium is the second-most abundant element in the Universe after hydrogen and is one of the main constituents of gas-giant planets in our Solar System. Early theoretical models predicted helium to be among the most readily detectable species in the atmospheres of exoplanets, especially in extended and escaping atmospheres. Searches for helium, however, have hitherto been unsuccessful. Here we report observations of helium on an exoplanet, at a confidence level of 4.5 standard deviations. We measured the near- infrared transmission spectrum of the warm gas giant WASP-107b and identified the narrow absorption feature of excited metastable helium at 10,833 angstroms. The amplitude of the feature, in transit depth, is 0.049 +/- 0.011 per cent in a bandpass of 98 angstroms, which is more than five times greater than what could be caused by nominal stellar chromospheric activity. This large absorption signal suggests that WASP-107b has an extended atmosphere that is eroding at a total rate of 10^10 to 3 x 10^11 grams per second (0.1-4 per cent of its total mass per billion years), and may have a comet-like tail of gas shaped by radiation pressure.

The lithium abundances in a few percent of giants exceed the value predicted by the standard stellar evolution models, and the mechanisms of Li enhancement are still under debate. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey has obtained over six million spectra in the past five years, and thus provides a great opportunity to search these rare objects and to more clearly understand the mechanisms of Li enhancement. Based on the high-resolution spectrum we obtained the stellar parameters ($T_\mathrm{eff}$, $\log g$, [Fe/H]), and determined the elemental abundances of Li, C, N, $\alpha$, Fe-peak, r-process, s-process elements, and the projected rotational velocity. For a better understanding of the effect of mixing processes, we also derived the $^{12}\rm{C}$ to $^{13}\rm{C}$ ratio, and constrained the evolutionary status of TYC\,3251-581-1 based on the BaSTI stellar isochrones. The super Li-rich giant TYC\,3251-581-1 has $\rm{A(Li)} = 3.51$, the average abundance of two lithium lines at $\lambda = 6708$ \AA\ and 6104 \AA\ based on the non-local thermodynamic equilibrium (NLTE) analysis. The atmospheric parameters show that our target locates on the luminosity function bump. The low carbon isotopic ratio ($^{12}\rm{C}/^{13}\rm{C} = 9.0 $), a slow rotational velocity $v\sin i = 2.2\ \rm{km\,s}^{-1}$, and no sign of IR excess suggest that additional mixing after first dredge up (FDU) should occur to bring internal synthesized Li to the surface. The low carbon ($[\rm{C}/\rm{Fe}] \sim -0.34$ ) and enhanced nitrogen ($[\rm{N}/\rm{Fe}] \sim 0.33$) are also consistent with the sign of mixing. Given the evolutionary stage of TYC\,3251-581-1 with the relatively low $^{12}\rm{C}/^{13}\rm{C}$, the internal production which replenishes Li in the outer layer is the most likely origin of Li enhancement for this star.

It is well known that a gravitational wave (GW) experiences the damping effect when it propagates in a fluid with nonzero shear viscosity. In this paper, we propose a new method to constrain the GW damping rate and thus the fluid shear viscosity. By defining the effective distance which incorporates damping effects, we can transform the GW strain expression in a viscous Universe into the same form as that in a perfect fluid. Therefore, the constraints of the luminosity distances from the observed GW events by LIGO and Virgo can be directly applied to the effective distances in our formalism. We exploit the lognormal likelihoods for the available GW effective distances and a Gaussian likelihood for the luminosity distance inferred from the electromagnetic radiation observation of the binary neutron star merger event GW170817. Our fittings show no obvious damping effects in the current GW data, and the upper limit on the damping rate with the combined data is $6.75 \times 10^{-4}\,{\rm Mpc}^{-1}$ at 95\% confidence level. By assuming that the dark matter self-scatterings are efficient enough for the hydrodynamic description to be valid, we find that a GW event from its source at a luminosity distance $D\gtrsim 10^4\;\rm Mpc$ can be used to put a constraint on the dark matter self-interactions.

Nearby sources may contribute to cosmic-ray electron (CRE) structures at high energies. Recently, the first DAMPE results on the CRE flux hinted at a narrow excess at energy ~1.4 TeV. We show that in general a spectral structure with a narrow width appears in two scenarios: I) "Spectrum broadening" for the continuous sources with a delta-function-like injection spectrum. In this scenario, a finite width can develop after propagation through the Galaxy, which can reveal the distance of the source. Well-motivated sources include mini-spikes and subhalos formed by dark matter (DM) particles $\chi_{s}$ which annihilate directly into e+e- pairs. II) "Phase-space shrinking" for burst-like sources with a power-law-like injection spectrum. The spectrum after propagation can shrink at a cooling-related cutoff energy and form a sharp spectral peak. The peak can be more prominent due to the energy-dependent diffusion. In this scenario, the width of the excess constrains both the power index and the distance of the source. Possible such sources are pulsar wind nebulae (PWNe) and supernova remnants (SNRs). We analysis the DAMPE excess and find that the continuous DM sources should be fairly close within ~0.3 kpc, and the annihilation cross sections are close to the thermal value. For the burst-like source, the narrow width of the excess suggests that the injection spectrum must be hard with power index significantly less than two, the distance is within ~(3-4) kpc, and the age of the source is ~0.16 Myr. In both scenarios, large anisotropies in the CRE flux are predicted. We identify possible candidates of mini-spike (PWN) sources in the current Fermi-LAT 3FGL (ATNF) catalog. The diffuse gamma-rays from these sources can be well below the Galactic diffuse gamma-ray backgrounds and less constrained by the Ferm-LAT data, if they are located at the low Galactic latitude regions.

The substructures of the Galactic dark matter halo such as dark matter subhalos and dwarf galaxies have very low velocity dispersions, which makes them useful in constraining the scenario of Sommerfeld-enhanced dark matter annihilation. We calculate the velocity distribution of dark matter particles in dark matter halo substructures using the Eddington's formula with NFW density profile. We parametrize the effect of Sommerfeld enhancement of s-wave dark matter annihilation on the gamma-ray flux as the Sommerfeld-enhanced J-factors, and explicitly calculate their values for 15 known dwarf spheroidal galaxies. Using the Fermi-LAT 3FGL data on the unassociated point-sources and the N-body simulation results on the dark matter subhalo distribution, we derive upper limits on the dark matter annihilation cross sections with Sommerfeld enhacement. Similar upper limits are derived from the Ferm-LAT gamma-ray data on the dwarf spheroidal galaxies. We find that in a wide region of parameter space, the constraints can be a few orders of magnitude more stringent than that in the case without the Sommerfeld enhancement. For dark matter annihilation dominantly into $e^+e^-$, $\mu^+\mu^-$, $\tau^{+}\tau^{-}$ and $b\bar b$, those constraints can exclude thermal relic dark matter for the dark matter mass below about 1 TeV.

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms.

We propose a new mechanism to generate a lepton asymmetry based on the vacuum CP-violating phase transition (CPPT). This approach differs from classical thermal leptogenesis as a specific seesaw model, and its UV completion, need not be specified. The lepton asymmetry is generated via the dynamically realised coupling of the Weinberg operator during the phase transition. This mechanism provides a connection with low-energy neutrino observables.

Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared G141 taken in six consecutive orbits observations of HN Peg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the $1.1-1.7\mu$m broadband light curve has the amplitude of $1.206\pm0.025\%$ and period of $15.4\pm0.5$ hr. The modulation amplitude has no detectable wavelength dependence except in the 1.4 $\mu$m water absorption band, indicating that the characteristic condensate particle sizes are large ($>1\mu$m). We detect significantly ($4.4\sigma$) lower modulation amplitude in the 1.4$\mu$m water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3.

Study of primary stars lying in Sirius-like systems with various masses of WD companions and orbital separations is one of the key aspects to understand the origin and nature of Barium (Ba) stars. In this paper, based on high resolution and high S/N spectra, we present systematic analysis of photospheric abundances for 18 FGK primary stars of Sirius-like systems including six giants and 12 dwarfs. Atmospheric parameters, stellar masses, and abundances of 24 elements (C, Na, Mg, Al, Si, S, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Sr, Y, Zr, Ba, La, Ce and Nd) are determined homogeneously. The abundance patterns in these sample stars show that most of the elements in our sample follow the behavior of field stars with similar metallicity. As expected, s-process elements in four known Ba giants show overabundance. A weak correlation was found between anomalies of s-process elemental abundance and orbital separation, suggesting the orbital separation of the binaries could not be the main constraint to differentiate strong Ba stars from mild Ba stars. Our study shows that the large mass (>0.51 M ) of a WD companion in a binary system is not a sufficient condition to form a Ba star, even if the separation between the two components is small. Although not sufficient it seems to be a necessary condition since Ba stars with lower mass WDs in the observed sample were not found. Our results support that [s/Fe] and [hs/ls] ratios of Ba stars are anti-correlated with the metallicity. However, the different levels of s-process overabundance among Ba stars may not to be dominated mainly by the metallicity.

Recently, an instability due to the nonlinear coupling of p-modes to g-modes in tidally deformed neutron stars in coalescing binaries has been studied in some detail. The result is significant because it could influence the inspiral and leave an imprint on the gravitational wave signal that depends on the neutron star equation of state (EOS). Because of its potential importance, the details of the instability should be further elucidated and its sensitivity to the EOS should be investigated. To this end, we carry out a numerical analysis with six representative EOSs for both static and non-static tides. We confirm that the absence of the p-g instability under static tides, as well as its return under non-static tides, is generic across EOSs, and further reveal a new contribution to it that becomes important for moderately high-order p-g pairs (previous studies concentrated on very high order modes), whose associated coupling strength can vary by factors of ~10-100 depending on the EOS. We find that, for stars with stiffer EOSs and smaller buoyancy frequencies, the instability onsets earlier in the inspiral and the unstable modes grow faster. These results suggest that the instability's impact on the gravitational wave signal might be sensitive to the neutron star EOS. To fully assess this prospect, future studies will need to investigate its saturation as a function of the EOS and the binary parameters.

The shape of LOSVDs (line-of-sight velocity distributions) carries important information about the internal dynamics of galaxies. The skewness of LOSVDs represents their asymmetric deviation from a Gaussian profile. Correlations between the skewness parameter ($h_3$) and the mean velocity (\vm) of a Gauss-Hermite series reflect the underlying stellar orbital configurations of different morphological components. Using two self-consistent $N$-body simulations of disk galaxies with different bar strengths, we investigate \hv\ correlations at different inclination angles. Similar to previous studies, we find anticorrelations in the disk area, and positive correlations in the bar area when viewed edge-on. However, at intermediate inclinations, the outer parts of bars exhibit anticorrelations, while the core areas dominated by the boxy/peanut-shaped (B/PS) bulges still maintain weak positive correlations. When viewed edge-on, particles in the foreground/background disk (the wing region) in the bar area constitute the main velocity peak, whereas the particles in the bar contribute to the high-velocity tail, generating the \hv\ correlation. If we remove the wing particles, the LOSVDs of the particles in the outer part of the bar only exhibit a low-velocity tail, resulting in a negative \hv\ correlation, whereas the core areas in the central region still show weakly positive correlations. We discuss implications for IFU observations on bars, and show that the variation of the \hv\ correlation in the disk galaxy may be used as a kinematic indicator of the bar and the B/PS bulge.

We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method, the formation of current singularity has previously been confirmed in the Hahm--Kulsrud--Taylor problem in 2D. In this paper, we extend this problem to 3D line-tied geometry. The linear solution, which is singular in 2D, is found to be smooth for arbitrary system length. However, with finite amplitude, the linear solution can become pathological when the system is sufficiently long. The nonlinear solutions turn out to be smooth for short systems. Nonetheless, the scaling of peak current density versus system length suggests that the nonlinear solution may become singular at finite length. With the results in hand, we can neither confirm nor rule out this possibility conclusively, since we cannot obtain solutions with system length near the extrapolated critical value.

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.

Observations of rotational modulations of brown dwarfs and giant exoplanets allow the characterization of condensate cloud properties. As of now rotational spectral modulations have only been seen in three L-type brown dwarfs. We report here the discovery of rotational spectral modulations in LP261-75B, an L6-type intermediate surface gravity companion to an M4.5 star. As a part of the Cloud Atlas Treasury program we acquired time-resolved Wide Field Camera 3 grism spectroscopy (1.1--1.69~$\mu$m) of LP261-75B. We find gray spectral variations with the relative amplitude displaying only a weak wavelength dependence and no evidence for lower-amplitude modulations in the 1.4~$\mu$m water band than in the adjacent continuum. The likely rotational modulation period is 4.78$\pm$0.95 h, although the rotational phase is not well sampled. The minimum relative amplitude in the white light curve measured over the whole wavelength range is 2.41$\pm$0.14\%. We report an unusual light curve with seemingly three peaks approximately evenly distributed in rotational phase. The spectral modulations suggests that the upper atmosphere cloud properties in {LP261-75B} are similar to two other mid-L dwarfs of typical infrared colors, but differ from that of the extremely red L-dwarf WISE0047.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

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.

The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 micron silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8 +/- 6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources.

Coronal heating has been a long-standing conundrum in solar physics. Parker's conjecture that spontaneous current singularities lead to nanoflares that heat the corona has been controversial. In ideal magnetohydrodynamics (MHD), can genuine current singularities emerge from a smooth 3D line-tied magnetic field? To numerically resolve this issue, the schemes employed must preserve magnetic topology exactly to avoid artificial reconnection in the presence of singular current densities. Structure-preserving numerical methods are favorable for mitigating numerical dissipation, and variational integration is a powerful machinery for deriving them. In this thesis, we develop variational integrators for ideal MHD in Lagrangian labeling by discretizing Newcomb's Lagrangian on a moving mesh using discretized exterior calculus. With the built-in frozen-in equation, the schemes are free of artificial reconnection, hence optimal for studying current singularity formation. Using this method, we first study a fundamental prototype problem in 2D, the Hahm-Kulsrud-Taylor (HKT) problem. It considers the effect of boundary perturbations on a 2D plasma magnetized by a sheared field, and its linear solution is singular. We find that with increasing resolution, the nonlinear solution converges to one with a current singularity. We then extend the HKT problem to 3D line-tied geometry, which models the solar corona by anchoring the field lines in the boundaries. The linear solution, which is singular in 2D, is found to be smooth. The nonlinear solution turns out smooth for short systems. Nonetheless, the scaling of peak current density vs. system length suggests that the nonlinear solution may become singular at a finite length. With the results in hand, we cannot confirm or rule out this possibility conclusively, since we cannot obtain solutions with system lengths near the extrapolated critical value.

It is well known that the Milgrom's MOND (modified Newtonian dynamics) explains well the mass discrepancy problem in galaxy rotation curves. The MOND predicts a universal acceleration scale below which the Newtonian dynamics is invalid yet. The universal acceleration scale we got from the SPARC dataset is $g_{\dag}=1.02\times10^{-10} \rm m~s^{-2}$. Milgrom suggested that the acceleration scale may be a fingerprint of cosmology on local dynamics and related with the Hubble constant $g_{\dag}\sim cH_0$. In this paper, we use the hemisphere comparison method with the SPARC dataset to investigate the spatial anisotropy on the acceleration scale. We find that the hemisphere of the maximum acceleration scale is in the direction $(l,b) = ({175.5^\circ}^{+6^\circ}_{-10^\circ}, {-6.5^\circ}^{+8^\circ}_{-3^\circ})$ with $g_{\dag,max}=1.10\times10^{-10} \rm m~s^{-2}$, while the hemisphere of the minimum acceleration scale is in the opposite direction $(l,b) = ({355.5^\circ}^{+6^\circ}_{-10^\circ}, {6.5^\circ}^{+3^\circ}_{-8^\circ})$ with $g_{\dag,min}=0.76\times10^{-10} \rm m~s^{-2}$. The maximum anisotropy level reaches up to $0.37\pm0.04$. Robust tests present that such a level of anisotropy can't be reproduced by a statistically isotropic data. In addition, we show that the spatial anisotropy on the acceleration scale has little correlation with the non-uniform distribution of the SPARC data points in sky. We also find that the maximum anisotropy direction is close with other cosmological preferred directions, especially the direction of the "Australia dipole" for the fine structure constant.

The observation and research of the solar radio emission have unique scientific values in solar and space physics and related space weather forecasting applications, since the observed spectral structures may carry important information about energetic electrons and underlying physical mechanisms. In this study, we present the design of a novel dynamic spectrograph that is installed at the Chashan solar radio station operated by Laboratory for Radio Technologies, Institute of Space Sciences at Shandong University. The spectrograph is characterized by the real-time storage of digitized radio intensity data in the time domain and its capability to perform off-line spectral analysis of the radio spectra. The analog signals received via antennas and amplified with a low-noise amplifier are converted into digital data at a speed reaching up to 32 k data points per millisecond. The digital data are then saved into a high-speed electronic disk for further off-line spectral analysis. Using different word length (1 k - 32 k) and time cadence (5 ms - 10 s) for the off-line fast Fourier transform analysis, we can obtain the dynamic spectrum of a radio burst with different (user-defined) temporal (5 ms - 10 s) and spectral (3 kHz ~ 320 kHz) resolution. This brings a great flexibility and convenience to data analysis of solar radio bursts, especially when some specific fine spectral structures are under study.

The large-amplitude longitudinal oscillations of solar filaments have been observed and explored for more than ten years. Previous studies are mainly based on the one-dimensional rigid flux tube model with a single magnetic dip. However, it is noticed that there might be two magnetic dips, and hence two threads, along one magnetic field line. Following the previous work, we intend to investigate the kinematics of the filament longitudinal oscillations when two threads are magnetically connected, which is done by solving one-dimensional radiative hydrodynamic equations with the numerical code MPI-AMRVAC. Two different types of perturbations are considered, and the difference from previous works resulting from the filament thread-thread interaction is investigated. We find that even with the inclusion of the thread-thread interaction, the oscillation period is modified weakly, by at most 20% compared to the traditional pendulum model with one thread. However, the damping timescale is significantly affected by the thread-thread interaction. Hence, we should take it into account when applying the consistent seismology to the filaments where two threads are magnetically connected.

The Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) near-IR channel is extensively used in time-resolved observations, especially for transiting exoplanet spectroscopy and brown dwarf and directly imaged exoplanet rotational phase mapping. The ramp effect is the dominant source of systematics in the WFC3 for time-resolved observations, which limits its photometric precision. Current mitigation strategies are based on empirical fits and require additional orbits "to help the telescope reach a thermal equilibrium". We show that the ramp effect profiles can be explained and corrected with high fidelity using charge trapping theories. We also present a model for this process that can be used to predict and to correct charge trap systematics. Our model is based on a very small number of parameters that are intrinsic to the detector. We find that these parameters are very stable between the different datasets, and we provide best-fit values. Our model is tested with more than 120 orbits ($\sim40$ visits) of WFC3 observations and is proved to be able to provide near photon noise limited corrections for observations made with both staring and scanning modes of transiting exoplanets as well as for starting-mode observations of brown dwarfs. After our model correction, the light curve of the first orbit in each visit has the same photometric precision as subsequent orbits, so data from the first orbit need no longer be discarded. Near IR arrays with the same physical characteristics (e.g., JWST/NIRCam) may also benefit from the extension of this model, if similar systematic profiles are observed.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

It has been revealed that the magnetic topology in the solar atmosphere displays hemispheric preference, i.e., negative/positive helicity in the northern/southern hemisphere, respectively. However, the strength of the hemispheric rule and its cyclic variation are controversial. In this paper, we apply a new method based on filament drainage to 571 erupting filaments from 2010 May to 2015 December in order to determine the filament chirality and its hemispheric preference. It is found that $91.6\%$ of our sample of erupting filaments follow the hemispheric rule of helicity sign. It is found that the strength of the hemispheric preference of the quiescent filaments decreases slightly from $\sim$$97\%$ in the rising phase to $\sim$$85\%$ in the declining phase of solar cycle 24, whereas the strength of the intermediate filaments keeps a high value around $96\pm4\%$ all the time. Only the active-region filaments show significant variations. Their strength of the hemispheric rule rises from $\sim$$63\%$ to $\sim$$95\%$ in the rising phase, and keeps a high value $82\pm5\%$ during the declining phase. Furthermore, during a half-year period around the solar maximum, their hemispheric preference totally vanishes. Besides, we also diagnose the magnetic configurations of the filaments based on our indirect method, and found that in our sample of erupting events, $89\%$ are inverse-polarity filaments with a flux rope magnetic configuration, whereas $11\%$ are normal-polarity filaments with a sheared arcade configuration.

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-Ks=2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy we find a best-fit rotational period (13.20$\pm$0.14 hours) with a larger amplitude at 1.1 micron than at 1.7 micron. This is the third largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as haze particle, we successfully explain the wavelength dependent amplitude with submicron-sized haze particles sizes of around 0.4 {\mu}m. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed, the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition.

The current accelerated expansion of the universe has been one of the most important fields in physics and astronomy since 1998. Many cosmological models have been proposed in the literature to explain this mysterious phenomenon. Since the nature and cause of the cosmic acceleration are still unknown, model-independent approaches to study the evolution of the universe are welcome. One of the powerful model-independent approaches is the so-called cosmography. It only relies on the cosmological principle, without postulating any underlying theoretical model. However, there are several shortcomings in the usual cosmography. For instance, it is plagued with the problem of divergence (or an unacceptably large error), and it fails to predict the future evolution of the universe. In the present work, we try to overcome or at least alleviate these problems, and we propose two new generalizations of cosmography inspired by the Pad\'e approximant. One is to directly parameterize the luminosity distance based on the Pad\'e approximant, while the other is to generalize cosmography with respect to a so-called $y_\beta$-shift $y_\beta=z/(1+\beta z)$, which is also inspired by the Pad\'e approximant. Then, we confront them with the observational data with the help of the Markov chain Monte Carlo (MCMC) code emcee, and find that they work fairly well.

The first data from the LHC Run-2 have shown a possible excess in diphoton events with invariant mass $\sim 750$ GeV, suggesting the existence of a new resonance which may decay dominantly into dark matter (DM) particles. We show in a simple model that the reported diphoton excess at the LHC is consistent with another photon excess, the $2$ GeV excess in cosmic gamma-ray fluxes towards the Galactic Center observed by the Fermi-LAT. Both the excesses can be simultaneously explained by a $\sim 60$ GeV scalar DM particle annihilating dominantly into two gluons with a typical thermal annihilation cross section, which leads to the prediction of a large width to mass ratio $\Gamma/M\approx \mathcal{O}(10^{-2})$ of the resonance. The upper limit on the dijet search at LHC Run-1 leads to a $lower$ limit on the predicted cross section for DM annihilating into $\gamma\gamma$ final states $\langle\sigma v\rangle_{\gamma\gamma} \gtrsim\mathcal{O}(10^{-30})~\mbox{cm}^{3}\mbox{s}^{-1}$. Both the predictions can be tested by the LHC, Fermi-LAT and future experiments.

Recently a variational integrator for ideal magnetohydrodynamics in Lagrangian labeling has been developed. Its built-in frozen-in equation makes it optimal for studying current sheet formation. We use this scheme to study the Hahm-Kulsrud-Taylor problem, which considers the response of a 2D plasma magnetized by a sheared field under sinusoidal boundary forcing. We obtain an equilibrium solution that preserves the magnetic topology of the initial field exactly, with a fluid mapping that is non-differentiable. Unlike previous studies that examine the current density output, we identify a singular current sheet from the fluid mapping. These results are benchmarked with a constrained Grad-Shafranov solver. The same signature of current singularity can be found in other cases with more complex magnetic topologies.

Rotational modulations of brown dwarfs have recently provided powerful constraints on the properties of ultra-cool atmospheres, including longitudinal and vertical cloud structures and cloud evolution. Furthermore, periodic light curves directly probe the rotational periods of ultra-cool objects. We present here, for the first time, time-resolved high-precision photometric measurements of a planetary-mass companion, 2M1207b. We observed the binary system with HST/WFC3 in two bands and with two spacecraft roll angles. Using point spread function-based photometry, we reach a nearly photon-noise limited accuracy for both the primary and the secondary. While the primary is consistent with a flat light curve, the secondary shows modulations that are clearly detected in the combined light curve as well as in different subsets of the data. The amplitudes are 1.36% in the F125W and 0.78% in the F160W filters, respectively. By fitting sine waves to the light curves, we find a consistent period of $10.7^{+1.2}_{-0.6}$ hours and similar phases in both bands. The J- and H-band amplitude ratio of 2M1207b is very similar to a field brown dwarf that has identical spectral type but different J-H color. Importantly, our study also measures, for the first time, the rotation period for a directly imaged extra-solar planetary-mass companion.

Today, $f(T)$ theory has been one of the popular modified gravity theories to explain the accelerated expansion of the universe without invoking dark energy. In this work, we consider the so-called Hojman symmetry in $f(T)$ theory. Unlike Noether conservation theorem, the symmetry vectors and the corresponding conserved quantities in Hojman conservation theorem can be obtained by using directly the equations of motion, rather than Lagrangian or Hamiltonian. We find that Hojman symmetry can exist in $f(T)$ theory, and the corresponding exact cosmological solutions are obtained. We find that the functional form of $f(T)$ is restricted to be the power-law or hypergeometric type, while the universe experiences a power-law or hyperbolic expansion. These results are different from the ones obtained by using Noether symmetry in $f(T)$ theory. Therefore, it is reasonable to find exact cosmological solutions via Hojman symmetry.

The first measurement on the antiproton to proton ratio made by the AMS-02 collaboration agrees with the expection from conventional cosmic-ray secondaries in the kinetic energy range $\sim 10-100$ GeV, which can be turned into stringent upper limits on the dark matter (DM) annihilation cross sections above $\sim 300$ GeV. Using the GALPROP code, we derive the upper limits in various propagation models and DM profiles. We show that in the "conventional" propagation model, for the $q\bar q$, $b\bar b$, and $WW$ final states, the constraints can be more stringent than that derived from the recent Ferm-LAT gamma-ray data on the dwarf spheroidal satellite galaxies. Making use of the typical minimal, median and maximal models obtained from a previous global fit, we show that the variation of the upper limits is around a factor of five. The possibility of DM contribution to the high energy $\bar p/p$ data is discussed.

The AMS-02 experiment is measuring the high energy charged cosmic rays with unprecedented accuracy. We explore the possibility of determining the cosmic-ray propagation models using the AMS-02 data $alone$. A global Bayesian analysis of the constraints on the cosmic-ray propagation models from the latest AMS-02 data on the Boron to Carbon nuclei flux ratio and proton flux is performed, with the assumption that the primary nucleon source is a broken power law in rigidity. The ratio of the diffusion coefficient $D_{0}$ to the diffusive halo height $Z_{h}$ is found to be determined with high accuracy $D_{0}/Z_{h}\simeq 2.00\pm0.07\text{cm}^{2}\text{s}^{-1}\text{kpc}^{-1}$. The best-fit value of the halo width is $Z_{h}\simeq 3.3$ kpc with uncertainty less than $50\$. As a consequence, the typical uncertainties in the positron fraction is within a factor of two, and that in the antiproton flux is within an order of magnitude. Both of them are significantly smaller than that from the analyses prior to AMS-02. Taking into account all the uncertainties and correlations in the propagation parameters we derive conservative upper limits on the cross sections for DM annihilating into various standard model final states from the current PAMELA antiproton data. We also investigate the reconstruction capability of the future AMS-02 antiproton data on the DM properties. The result shows that for DM particles lighter than 100 GeV and with typical thermal annihilation cross section, the cross section can be well reconstructed with uncertainties about a factor of two for the AMS-02 three-year data taking.

We report and analyze the observational evidence of global kink oscillations in a solar filament as observed in H alpha by National Solar Observatory (NSO)/Global Oscillation Network Group (GONG) instrument. An M1.1-class flare in active region 11692 on 2013 March 15 induced a global kink mode in the filament lying in the south-west of AR11692.We find periods of about 61 - 67 minutes and damping times of 92 - 117 minutes at three vertical slice positions chosen in and around the filament apex. We find that the waves are damped. From the observed global kink mode period and damping time scale using the theory of resonant absorption we perform prominence seismology. We estimate a lower cut-off value for the inhomogeneity length-scale to be around 0.34 - 0.44 times the radius of the filament cross-section.

In an extended effective operator framework, we investigate in detail the effects of light mediators on the event spectra of dark matter (DM)-nucleus scatterings. The presence of light mediators changes the interpretation of the current experimental data, especially the determination of DM particle mass. We show by analytic and numerical illustrations that for all the operators relevant to spin-independent scatterings, the DM particle mass allowed by the data increases significantly when the mediator particle becomes lighter. For the data of CDMS-II-Si, the allowed DM particle mass can reach $\sim 50$ GeV, depending on the type of operators, and the increase of DM particle mass saturates when the mediator mass is below $\mathcal{O}(10)$ MeV. The upper limits from other experiments such as SuperCDMS, CDMSlite, CDEX, XENON10/100, LUX, PandaX etc. all tend to be weaker toward high DM mass regions. In a combined analysis, we show that the presence of light mediators can partially relax the tension in the current results of CDMS-II-Si, SuperCDMS and LUX.

The so-called Einstein-Aether theory is General Relativity coupled (at second derivative order) to a dynamical unit time-like vector field (the aether). It is a Lorentz-violating theory, and gained much attention in the recent years. In the present work, we study the cosmological evolution of Einstein-Aether models with power-law-like potential, by using the method of dynamical system. In the case without matter, there are two attractors which correspond to an inflationary universe in the early epoch, or a de Sitter universe in the late time. In the case with matter but there is no interaction between dark energy and matter, there are only two de Sitter attractors, and no scaling attractor exists. So, it is difficult to alleviate the cosmological coincidence problem. Therefore, we then allow the interaction between dark energy and matter. In this case, several scaling attractors can exist under some complicated conditions, and hence the cosmological coincidence problem could be alleviated.

In this paper, we introduce model-independent data analysis procedures for identifying inelastic WIMP-nucleus scattering as well as for reconstructing the mass and the mass splitting of inelastic WIMPs simultaneously and separately. Our simulations show that, with O(50) observed WIMP signals from one experiment, one could already distinguish the inelastic WIMP scattering scenarios from the elastic one. By combining two or more data sets with positive signals, the WIMP mass and the mass splitting could even be reconstructed with statistical uncertainties of less than a factor of two.

As is well known, in mathematics, any function could be approximated by the Pad\'e approximant. The Pad\'e approximant is the best approximation of a function by a rational function of given order. In fact, the Pad\'e approximant often gives better approximation of the function than truncating its Taylor series, and it may still work where the Taylor series does not converge. In the present work, we consider the Pad\'e approximant in two issues. First, we obtain the analytical approximation of the luminosity distance for the flat XCDM model, and find that the relative error is fairly small. Second, we propose several parameterizations for the equation-of-state parameter (EoS) of dark energy based on the Pad\'e approximant. They are well motivated from the mathematical and physical points of view. We confront these EoS parameterizations with the latest observational data, and find that they can work well. In these practices, we show that the Pad\'e approximant could be an useful tool in cosmology, and it deserves further investigation.

Measurements of accretion rates onto planetary mass objects may distinguish between different planet formation mechanisms, which predict different accretion histories. In this Letter, we use \HST/WFC3 UVIS optical photometry to measure accretion rates onto three accreting objects, GSC06214-00210 b, GQ Lup b, and DH Tau b, that are at the planet/brown dwarf boundary and are companions to solar mass stars. The excess optical emission in the excess accretion continuum yields mass accretion rates of $10^{-9}$ to $10^{-11}$ \Msol/yr for these three objects. Their accretion rates are an order of magnitude higher than expected from the correlation between mass and accretion rates measured from the UV excess, which is applicable if these wide planetary mass companions formed by protostellar core fragmentation. The high accretion rates and large separation from the central star demonstrate the presence of massive disks around these objects. Models for the formation and evolution of wide planetary mass companions should account for their large accretion rates. High ratios of H$\alpha$ luminosity over accretion luminosity for objects with low accretion rates suggest that searches for H$\alpha$ emission may be an efficient way to find accreting planets.

Hard X-ray emission in radio-quiet active galactic nuclei (AGNs) is believed to be produced via inverse Compton scattering by hot and compact coronae near the super massive black hole. However the origin and physical properties of the coronae, including geometry, kinematics and dynamics, yet remain poorly known. In this work, taking [OIV] 25.89um emission line as an isotropic indicator of AGN's intrinsic luminosity, we compare the intrinsic corona X-ray emission between Seyfert 1 and Compton-thin Seyfert 2 galaxies, which are viewed at different inclinations according to the unification scheme. We compile a sample of 130 Compton-thin Seyfert galaxies with both [OIV] 25.89um line luminosities measured with Spitzer-IRS and X-ray spectra observed by XMM-Newton, Chandra, Suzaku or Swift. Known radio-loud sources are excluded. We fit the X-ray spectra to obtain the absorption-corrected 2-10 keV continuum luminosities. We find that Seyfert 1 galaxies are intrinsically brighter in intrinsic 2-10 keV emission by a factor of 2.8+0.5-0.4 (2.2+0.9-0.3 in Swift-BAT 14-195 keV emission), comparing with Compton-thin Seyfert 2 galaxies. The Seyfert 1 and Compton-thin Seyfert 2 galaxies follow a statistically identical correlation between the absorption-corrected 2-10 keV luminosity and the 14-195 keV luminosity, indicating that our absorption correction to the 2-10 keV flux is sufficient. The difference in X-ray emission between the two populations is thus unlikely due to X-ray absorption, and instead implies an intrinsic anisotropy in the corona X-ray emission. This striking anisotropy of X-ray emission can be explained by a bipolar outflowing corona with a bulk velocity of ~0.3-0.5c. This would provide a natural link between the so-called coronae and weak jets in these systems. Other consequences of outflowing coronae are also discussed.

High-resolution H$\alpha$ observations indicate that filaments consist of an assembly of thin threads. In quiescent filaments, the threads are generally short, whereas in active region filaments, the threads are generally long. In order to explain these observational features, we performed one-dimensional radiative hydrodynamic simulations of filament formation along a dipped magnetic flux tube in the framework of the chromospheric evaporation-coronal condensation model. The geometry of a dipped magnetic flux tube is characterized by three parameters, i.e., the depth ($D$), the half-width ($w$), and the altitude ($h$) of the magnetic dip. The parameter survey in the numerical simulations shows that allowing the filament thread to grow in 5 days, the maximum length ($L_{th}$) of the filament thread increases linearly with $w$, and decreases linearly with $D$ and $h$. The dependence is fitted into a linear function $L_{th}=0.84w-0.88D-2.78h+17.31$ Mm. Such a relation can qualitatively explain why quiescent filaments have shorter threads and active region filaments have longer threads.

In light of the first measurement of the positron fraction by the AMS-02 experiment, we perform a detailed global analysis on the interpretation of the latest data of PAMELA, Fermi-LAT, and AMS-02 in terms of dark matter (DM) annihilation and decay in various propagation models. The allowed regions for the DM particle mass and annihilation cross section or decay life-time are obtained for channels with leptonic final states: $2e$, $2\mu$, $2\tau$, $4e$, $4\mu$ and $4\tau$. We show that for the conventional astrophysical background the AMS-02 positron fraction data alone favour a DM particle mass $\sim 500 \ (800)$ GeV if DM particles annihilate dominantly into $2\mu \ (4\mu)$ final states, which is significantly lower than that favoured by the Fermi-LAT data of the total flux of electrons and positrons. The allowed regions by the two experiments do not overlap at a high confidence level ($99.99999\%$C.L.). We consider a number of propagation models with different halo height $Z_{h}$, diffusion parameters $D_{0}$ and $\delta_{1/2}$, and power indices of primary nucleon sources $\gamma_{p1/p2}$. The normalization and the slope of the electron background are also allowed to vary. We find that the tension between the two experiments can be only slightly reduced in the propagation model with large $Z_{h}$ and $D_{0}$. The consistency of fit is improved for annihilation channels with $2\tau$ and $4\tau$ final states which favour TeV scale DM particle with large cross sections above $\sim 10^{-23} \text{cm}^3\text{s}^{-1}$. In all the considered leptonic channels, the current data favour the scenario of DM annihilation over DM decay. In the decay scenario, the charge asymmetric DM decay is slightly favoured.

We investigate the electron energy distributions (EEDs) and the corresponding acceleration processes in the jet of PKS 0447$-$439 and estimate its redshift through modeling its observed spectral energy distribution (SED) in the frame of a one-zone synchrotron-self Compton (SSC) model. Three EEDs formed in different acceleration scenarios are assumed: the power-law with exponential cut-off (PLC) EED (shock-acceleration scenario or the case of the EED approaching equilibrium in the stochastic-acceleration scenario), the log-parabolic (LP) EED (stochastic-acceleration scenario and the acceleration dominating) and the broken power law (BPL) EED (no acceleration scenario), and then the corresponding fluxes of both synchrotron and SSC are calculated. The model is applied to PKS 0447-439 and modeling SEDs are compared to the observed SED of this object by using the Markov Chain Monte Carlo (MCMC) method. Calculating results show that PLC model fails to fit the observed SED well, while the LP and BPL models give comparably good fits for the observed SED. The results indicate that it is possible that stochastic acceleration process acts in the emitting region of PKS 0447-439 and the EED is far from equilibrium (acceleration dominating) or no acceleration process works (in the emitting region). The redshift of PKS 0447-439 is also estimated in our fitting, and $z=0.16\pm0.05$ for LP case and $z=0.17\pm0.04$ for BPL case.

We consider extensions of the standard model with fourth generation fermions (SM4) in which extra symmetries are introduced such that the transitions between the fourth generation fermions and the ones in the first three generations are forbidden. In these models, the stringent lower bounds on the masses of fourth generation quarks from direct searches are relaxed, and the lightest fourth neutrino is allowed to be stable and light enough to trigger the Higgs boson invisible decay. In addition, the fourth Majorana neutrino can be a subdominant but highly detectable dark matter component. We perform a global analysis of the current LHC data on the Higgs production and decay in this type of SM4. The results show that the mass of the lightest fourth Majorana neutrino is confined in the range $\sim 41-59$ GeV. Within the allowed parameter space, the predicted effective cross-section for spin-independent DM-nucleus scattering is $\sim 3\times 10^{-48}-6\times 10^{-46} \text{cm}^{2}$, which is close to the current Xenon100 upper limit and is within the reach of the Xenon1T experiment in the near future. The predicted spin-dependent cross sections can also reach $\sim 8\times 10^{-40}\text{cm}^{2}$.

If the massive neutrinos are identified to be the Majorana particles via a convincing measurement of the neutrinoless double beta (0\nu\beta\beta) decay, how to determine the Majorana CP-violating phases in the 3 x 3 lepton flavor mixing matrix U will become a desirable experimental question. The answer to this question is to explore all the possible lepton-number-violating (LNV) processes in which the Majorana phases really matter. In this paper we carry out a systematic study of CP violation in neutrino-antineutrino oscillations, whose CP-conserving parts involve six independent 0\nu\beta\beta-like mass terms <m>_{\alpha\beta} and CP-violating parts are associated with nine independent Jarlskog-like parameters V^{ij}_{\alpha\beta} (for \alpha, \beta = e, \mu, \tau and i, j = 1, 2,3). With the help of current neutrino oscillation data, we analyze the sensitivities of |<m>_{\alpha\beta}| and V^{ij}_{\alpha\beta} to the three CP-violating phases of U, and illustrate the salient features of six independent CP-violating asymmetries between \nu_\alpha \to \bar{\nu}_\beta and \bar{\nu}_\alpha \to \nu_\beta oscillations. As a by-product, the effects of the CP-violating phases on the LNV decays of doubly- and singly-charged Higgs bosons are reexamined by taking account of the unsuppressed value of \theta_{13}. Such CP-conserving LNV processes can be complementary to the possible measurements of neutrino-antineutrino oscillations in the distant future.

The first AMS-02 measurement confirms the existence of an excess in the cosmic-ray positron fraction previously reported by the PAMELA and Fermi-LAT experiments. If interpreted in terms of thermal dark matter (DM) annihilation, the AMS-02 result still suggests that the DM annihilation cross section in the present day should be significantly larger than that at freeze out. The Sommerfeld enhancement of DM annihilation cross section as a possible explanation for the observed excesses is subject to constraint from the DM thermal relic density, mainly due to the annihilation of DM particles into force-carrier particles introduced by the mechanism. We show that the effect of the Sommerfeld enhancement and the relic density constraint depend significantly on the nature of the force-carrier. Three cases where the force-carrier is a vector boson, scalar and pseudoscalar particle are investigated and compared with each other. The results show that for the vector force-carrier, the Sommerfeld enhancement can marginally account for the AMS-02 data for $2\mu$ final states, while for scalar force-carrier the allowed Sommerfeld enhancement can be larger by a factor of two. For the case with a pesudoscalar force-carrier, the Sommerfeld enhancement can be large enough to accommodate the AMS-02 and Fermi-LAT result for a variety of DM annihilation final states.

The quasi-simultaneous multi-wavelength emission of TeV blazar 3C 66A is studied by using a one-zone multi-component leptonic jet model. It is found that the quasi-simultaneous spectral energy distribution (SED) of 3C 66A can be well reproduced, especially its {\it Fermi}-LAT first 3 months average spectrum can be well reproduced by the synchrotron-self Compton (SSC) component plus external Compton (EC) component of the broad line region (BLR). Clues on its redshift and gamma-ray emission location are obtained. The results indicate the following. (i) On the redshift; The theoretical intrinsic TeV spectra can be predicted by extrapolating the reproduced GeV spectra. Through comparing this extrapolated TeV spectra with the extragalactic background light (EBL) corrected observed TeV spectra, it is suggested that the redshift of 3C 66A could be between 0.1 and 0.3, the most likely value is $\sim$ 0.2. (ii) On the gamma-ray emission location; To well reproduce the GeV emission of 3C 66A under different assumptions on BLR, the gamma-ray emission region is always required to be beyond the inner zone of BLR. The BLR absorption effect on gamma-ray emission confirms this point.

In a simple extension of the standard electroweak theory where the phenomenon of lepton flavor mixing is described by a 3x3 unitary matrix V, the electric and magnetic dipole moments of three active neutrinos are suppressed not only by their tiny masses but also by the Glashow-Iliopoulos-Maiani (GIM) mechanism. We show that it is possible to lift the GIM suppression if the canonical seesaw mechanism of neutrino mass generation, which allows V to be slightly non-unitary, is taken into account. In view of current experimental constraints on the non-unitarity of V, we find that the effective electromagnetic transition dipole moments of three light Majorana neutrinos and the rates of their radiative decays can be maximally enhanced by a factor of O(10^2) and a factor of O(10^4), respectively. This important observation reveals an intrinsic and presumably significant correlation between the electromagnetic properties of massive neutrinos and the origin of their small masses.

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

We present XMM-Newton spectra of the Seyfert 2 Galaxy IRAS 00521-7054. A strong feature at ~6 keV (observer's frame) can be formally fitted with a strong (EW=1.3+-0.3 keV in the rest frame) and broad Fe Ka line, extending down to 3 keV. The underlying X-ray continuum could be fitted with an absorbed powerlaw (with Gamma = 1.8+-0.2 and N_H = 5.9 x 10^22 cm^-2) plus a soft component. If due to relativistically smeared reflection by an X-ray illuminated accretion disk, the spin of the supermassive black hole is constrained to be 0.97^+0.03_-0.13 (errors at 90% confidence level for one interesting parameter), and the accretion system is viewed at an inclination angle of 37+-4 degree. This would be the first type 2 AGN reported with strong red Fe Ka wing detected which demands a fast rotating SMBH. The unusually large EW would suggest that the light bending effect is strong in this source. Alternatively, the spectra could be fitted by a dual absorber model (though with a global chi^2 higher by ~ 6 for 283 d.o.f) with N_H1 = 7.0 x 10^22 cm^-2 covering 100% of the X-ray source, and N_H2 = 21.7 x 10^22 cm^-2 covering 71%, which does not require an extra broad Fe Ka line.

Using Chandra High Energy Grating (HEG) observations of 32 AGNs, we present a systematic study of the X-ray Baldwin effect (XBE, i.e. the anti-correlation between narrow Fe K{\alpha} line EW and X-ray continuum luminosity for AGN samples) with the highest spectral resolution currently available. We have previously reported an anti-correlation with EW\proptoL(2-10 keV)^-0.22 in a HEG sample, and the correlation is much weaker after averaging multiple observations of individual AGNs (EW\proptoL(2-10 keV)^-0.13). This indicates that rapid variation in X-ray continuum plays an important role in producing the XBE, and such an effect should also be visible in individual AGNs. In this paper, by normalizing the line EWs and continuum luminosities to the time-averaged values for each AGN in our sample with multiple HEG observations, we find a strong anti-correlation between EW and Lx (EW/<EW>\propto (L/<L>)^-0.82+/-0.10), consistent with the XBE expected in an individual AGN if the narrow line flux remains constant while the continuum varies. This is first observational evidence that the Fe K{\alpha} line flux in a large sample of AGNs lacks of a corresponding response to the continuum variation, supporting that the narrow Fe-K line emission originates from a region far from the nucleus. We then performed Monte-Carlo simulations to address whether the global XBE can be produced by X-ray continuum variation solely, and found that such interpretation of XBE cannot be ruled out statistically. One thus should be very cautious before reaching any scientific conclusion based on an observed XBE.

We have compiled a sample of 21 low redshift (z < 0.3), luminous active galactic nuclei (AGN) with large Balmer decrements (Halpha/Hbeta > 7) using the galaxy and QSO catalogs of the Sloan Digital Sky Survey Early Data Release. Using this sample we attempt to determine the fraction of quasars with large internal absorption. We find that these AGN have Strong [O III] and broad Halpha emission, and that starlight dominates the spectra in the blue band, suggesting that these objects are heavily reddened. Their narrow emission line ratios are similar to those of Seyfert 2 galaxies, yet the average [O III]5007 emission is ~10 times more luminous. Applying the empirical relation between the optical continuum and the Balmer line luminosity for blue quasars, we find that the intrinsic luminosities of these 21 objects are in the range for quasars. We propose that they are obscured, intermediate type quasars analogous to type 1.8 and 1.9 Seyfert galaxies. The ratio of these optically selected, intermediate type quasars to type 1 quasars are found to be around 1, similar to that for local Seyfert galaxies. Preliminary study indicates that most of these quasars are hosted in early type galaxies.

We present a detailed analysis of the iron K$\alpha$ line variabilities in NGC 4151 by using long ASCA observation data obtained in May 1995. Despite the relatively small amplitude variations in the continuum flux, the iron K$\alpha$ line flux and profile show dramatic variations. Particularly, the line profile changes from single peak to seeming double peaks and back in time scales of a few 10$^4$ sec. The seemingly double-peaked profiles can be well interpreted as line emission from a Keplerian ring around a massive black hole. An absorption line at around 5.9 keV is also marginnaly detected. We discussed current Fe K line models, but none of them can well explain the observed line and continuum variations.

Recent observations of Broad Absorption Line (BAL) quasars demonstrated that the soft X-ray emission of these objects is extremely weak and convincing evidence for very strong absorption by a high column density (~ 10^23.5 cm^-2) was obtained for PG 1411+442, even though it is one of the few BAL QSOs strongly detected in soft X-rays. This paper examines the ionization status and geometry of the X-ray absorber by combining the properties of the UV lines with the X-ray continuum absorption. We show that the gas has to have large column densities in ions of major UV absorption lines, such as CIV, NV, OVI and Ne VIII, in order to have sufficient opacity around 0.2 to 0.35 keV. The UV absorption lines have to be saturated if the X-ray absorber intersects the line of sight to the UV continuum emission region. A uniformly covering UV and X-ray absorption model can be constructed for PG 1411+442 but in some other soft X-ray detected BAL QSOs, such as PG 1001+054, the observed line optical depth is much lower than expected from the X-ray absorbing material. We propose a scheme in which a substantial fraction of the line of sight to the continuum source may be covered by either an optically thick flow or clouds in a narrow velocity range, but in which the total covering factor of either the whole flow or all clouds is close to unity.

Recent observations show strong evidence for the presence of Doppler-shifted emission lines in the spectrum of both black hole candidates and active galactic nuclei. These lines are likely to originate from relativistic outflows (or jets) in the vicinity of the central black hole. Consequently, the profile of such a line should be distorted by strong gravitational effects near the black hole, as well as special relativistic effects. In this paper, we present results from a detailed study on how each process affects the observed line profile. We found that the profile is sensitive to the intrinsic properties of the jets (Lorentz factor, velocity profile, and emissivity law), as well as to the spin of the black hole and the viewing angle (with respect to the axis of the jets). More specifically, in the case of approaching jets, an intrisically narrow line (blue-shifted) is seen as simply broadened at small viewing angles, but it shows a doubly peaked profile at large viewing angles for extreme Kerr black holes (due to the combination of gravitational focusing and Doppler effects); the profile is always singly peaked for Schwarzschild black holes. For receding jets, however, the line profile becomes quite complicated owing to complicated photon trajectories. To facilitate comparison with observations, we searched a large parameter space to derive representative line profiles. We show the results and discuss how to use emission lines as a potential tool for probing the inner region of a black hole jet system.

We present an analysis of the Fe K line profile of NGC 4151 by using long ASCA observation data obtained in May 1995. The unprecedented good data quality, which is much better in the energy band around 6.4 keV than that of the famous 4.2-day ASCA observation of MCG -6-30-15 in July 1994, offers a unique opportunity to study the details of Fe K line profile. Apart from those characteristics already noticed in earlier ASCA observations on this object (Yaqoob et al. 1995): a broad and skewed profile, with a strong peak at about 6.4 keV and a large red wing extending to $\sim$4 - 5 keV, which is remarkably similar to that of MCG -6-30-15, we also find a weak blue wing extending to about 8 keV, thanks to the good quality of the data. When fitted by a relativistic accretion disk line plus a narrow core at 6.4 keV, the data constrain the accretion disk to be nearly face-on, contrary to the edge-on geometry inferred from optical and UV observations. However, the extended blue wing can not be well fitted even after we include corresponding Fe K$\beta$ components. Ni K$\alpha$ line emission by an amount of 12% of Fe K$\alpha$ is statistically required. An alternative explanation is a model consisting of a narrow core and two disk lines with inclinations of 58$^{\rm o}$ and 0$^{\rm o}$, respectively. We suppose that the component with inclination of 58$^{\rm o}$ was observed directly, consistent with its edge-on geometry, and the component with inclination of 0$^{\rm o}$ was scattered into our line of sight by a Compton mirror, which might be the cool accretion disk corona proposed by Poutanen et al. (1996).

We present a detailed analysis on the variability of the Fe K emission line in NGC 4051 using ASCA data. Through simple Gaussian line fits, we find not only obvious Fe K line variability with no significant difference in the X-ray continuum flux between two ASCA observations which were separated by $\sim$ 440 days, but also rapid variability of Fe K line on time scales $\sim 10^4$ s within the second observation. During the second observation, the line is strong (EW = 733$^{+206}_{-219}$ eV) and broad ($\sigma = 0.96^{+0.49}_{-0.35}$ keV) when the source is brightest, and become weaker (EW = 165$^{+87}_{-86}$ eV) and narrower ($\sigma<0.09$ keV) whilst the source is weakest. The equivalent width of Fe K line correlates positively with the continuum flux, which shows an opposite trend with another Seyfert 1 galaxy MCG --6-30-15.

Anisotropic emission of the ionizing continuum is a general prediction of the accretion disk models. In this paper, we present the results of correlation analysis of the UV emission line and UV to X-ray continuum properties for a large sample of broad emission line AGNs observed with ROSAT, IUE and HST. We find strong correlations between the CIV/$Ly\alpha$ ratio, the equivalent width of CIV, and the UV to soft X-ray spectral slope. The results are in good agreement with the photoionization calculation, suggesting that the overall ionizing continuum can well match the observed UV to soft X-ray spectrum. These results are consistent with the assumption of isotropic ionizing continuum shape. Our analysis suggests a small range for the ``big blue bump'' cutoff energy for the objects in this sample, consistent with the similar results of Laor et al. (1997) and Walter & Fink (1993) based on the continuum properties . The mean UV-to-X-ray spectral slope is similar to the soft X-ray spectral slope. This similarity also holds for radio-loud and radio-quiet objects separately. This suggests that the two might be drawn from the same distribution. The two spectral slopes are only weakly correlated. The UV to X-ray spectral index is correlated with absolute optical magnitude. This result confirms the earlier suggestion that the ionizing continua are softer for higher luminosity objects.