The accretion-induced collapse (AIC) scenario was proposed 40 years ago as an evolutionary end state of oxygen-neon white-dwarfs (ONe WDs), linking them to the formation of neutron star (NS) systems. However, there has been no direct detection of any AIC event so far, even though there exists a lot of indirect observational evidence. Meanwhile, the evolutionary pathways resulting in NS formation through AIC are still not well investigated. In this article, we review recent studies on the two classic progenitor models of AIC events, i.e., the single-degenerate model (including the ONe WD+MS/RG/He star channels and the CO WD+He star channel) and the double-degenerate model (including the double CO/ONe WD channels and the ONe WD+CO WD channel). Recent progress on these progenitor models is reviewed, including the evolutionary scenarios, the initial parameter space and the related objects. For the single-degenerate model, the pre-AIC systems could potentially be identified as supersoft X-ray sources, symbiotics and cataclysmic variables, whereas the post-AIC systems could be identified as low-/intermediate-mass X-ray binaries and the resulting low-/intermediate-mass binary pulsars, most notably millisecond pulsars. For the double-degenerate model, the pre-AIC systems are close double WDs, whereas the post-AIC systems are isolated NSs with peculiar properties. We also review the predicted rates of AIC events, the mass distribution of NSs, and the gravitational wave (GW) signals from double WDs that are potential GW sources in the Galaxy in the context of future space-based GW detectors. Recent theoretical and observational constraints on the detection of AIC events are summarized. In order to confirm the existence of the AIC process, and resolve this long-term issue presented by current stellar evolution theories, more numerical simulations and observational identifications are required.

The Baryon acoustic oscillations from Integrated Neutral Gas Observations (BINGO) telescope is a 40-m~class radio telescope under construction that has been designed to measure the large-angular-scale intensity of HI emission at 980--1260 MHz and hence to constrain dark energy parameters. A large focal plane array comprising of 1.7-metre diameter, 4.3-metre length corrugated feed horns is required in order to optimally illuminate the telescope. Additionally, very clean beams with low sidelobes across a broad frequency range are required, in order to facilitate the separation of the faint HI emission from bright Galactic foreground emission. Using novel construction methods, a full-sized prototype horn has been assembled. It has an average insertion loss of around 0.15 dB across the band, with a return loss around -25 dB. The main beam is Gaussian with the first sidelobe at around $-25 dB. A septum polariser to separate the signal into the two hands of circular polarization has also been designed, built and tested.

Streamer blowout (SBO) coronal mass ejections (CMEs) represent a particular class of CMEs that are characterized by a gradual swelling of the overlying streamer and a slow CME containing a flux-rope structure. SBO CMEs arising from the streamer flank fall into a special category of SBO CMEs involving three lower arches under the higher streamer arcade. However, the initiation mechanism for this special category of SBO CMEs remains elusive, due to the observational limitations. Here we report critical observations of a SBO CME associated with the eruption of a polar crown filament that originated from the streamer flank. The filament slowly rose toward the solar equator with the writhing motion, and underwent a sudden acceleration before its eruption. Interestingly, during the rising, the filament fields experienced gradual external reconnections, which is evidenced by the dip-shaped bottom of the enveloping flux-rope structure changing from a smooth concave, the slow inflows ($\sim$1.8 km s$^{-1}$) from both the filament fields and the coronal loops beneath, and the persistent brightenings around the interface between the filament fields and the coronal loops beneath. The newly formed lower loops at the filament source and the Y-shaped structure in the stretched tail fields indicate the internal reconnections for the filament eruption. The clear signatures of the external and internal reconnections shed light on the initiation mechanisms of SBO CMEs.

In this brief report, we investigate the existence of 4-dimensional static spherically symmetric black holes (BHs) in the Einstein-complex-scalar-Gauss-Bonnet (EcsGB) gravity with an arbitrary potential $V(\phi)$ and a coupling $f(\phi)$ between the scalar field $\phi$ and the Gauss-Bonnet (GB) term. We find that static regular BH solutions with complex scalar hairs do not exist. This conclusion does not depend on the coupling between the GB term and the scalar field, nor on the scalar potential $V(\phi)$ and the presence of a cosmological constant $\Lambda$ (which can be either positive or negative), as longer as the scalar field remains complex and is regular across the horizon.

The discovery of gravitational waves, which are ripples of space-time itself, opened a new window to test general relativity, because it predicts that there are only plus and cross polarizations for gravitational waves. For alternative theories of gravity, there may be up to six polarizations. The measurement of the polarization is one of the major scientific goals for future gravitational wave detectors. To evaluate the capability of the detector, we need to use the frequency dependent response functions averaged over the source direction and polarization angle. We derive the full analytical formulas of the averaged response functions for all six possible polarizations and present their asymptotic behaviors based on these analytical formulas. Compared with the numerical simulation, the full analytical formulas are more efficient and valid for any equal-arm interferometric gravitational wave detector without optical cavities in the arms and for a time-delay-interferometry Michelson combination.

Applying the distance sum rule in strong gravitational lensing (SGL) and type Ia supernova (SN Ia) observations, one can provide an interesting cosmological model-independent method to determine the cosmic curvature parameter $\Omega_k$. In this paper, with the newly compiled data sets including 161 galactic-scale SGL systems and 1048 SN Ia data, we place constraints on $\Omega_k$ within the framework of three types of lens models extensively used in SGL studies. Moreover, to investigate the effect of different mass lens samples on the results, we divide the SGL sample into three sub-samples based on the center velocity dispersion of intervening galaxies. In the singular isothermal sphere (SIS) and extended power-law lens models, a flat universe is supported with the uncertainty about 0.2, while a closed universe is preferred in the power-law lens model. We find that the choice of lens models and the classification of SGL data actually can influence the constraints on $\Omega_k$ significantly.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

Using the galaxy catalog built from ELUCID N-body simulation and the semi-analytical galaxy formation model, we have built a mock HI intensity mapping map. We have implemented the Finger-of-God (FoG) effect in the map by considering the galaxy HI gas velocity dispersion. By comparing the HI power spectrum in the redshift space with the measurement from IllustrisTNG simulation, we have found that such FoG effect can explain the discrepancy between current mock map built from N-body simulation and Illustris TNG simulation. Then we built a parameter-free FoG model and a shot-noise model to calculate the HI power spectrum. We found that our model can accurately fit both the monopole and quadrupole moments of the HI matter power spectrum. Our method of building the mock HI intensity map and the parameter-free FoG model will be widely useful for the up-coming 21cm intensity mapping experiments, such as CHIME, Tianlai, BINGO, FAST and SKA. It is also crucial for us to study the non-linear effects in 21cm intensity mapping.

We study the link between the kinematic-morphology of galaxies, as inferred from integral-field stellar kinematics, and their relation between mass and star formation rate (SFR). Our sample consists of $\sim 3200$ galaxies with integral-field spectroscopic data from the MaNGA survey with available determinations of their effective stellar angular momentum within the half-light radius $\lambda_{R_e}$. We find that for star-forming galaxies, namely along the star formation main sequence (SFMS), the $\lambda_{R_e}$ values remain large and almost unchanged over about two orders of magnitude in stellar mass, with the exception of the lowest masses $\mathcal{M}_{\star}\lesssim2\times10^{9} \mathcal{M}_{\odot}$, where $\lambda_{R_e}$ slightly decreases. The SFMS is dominated by spiral galaxies with small bulges. Below the SFMS, but above the characteristic stellar mass $\mathcal{M}_{\rm crit}\approx2\times10^{11} \mathcal{M}_{\odot}$, there is a sharp decrease in $\lambda_{R_e}$ with decreasing star formation rate: massive galaxies well below the SFMS are mainly slow-rotator early-type galaxies, namely genuinely spheroidal galaxies without disks. Below the SFMS and below $\mathcal{M}_{\rm crit}$ the decrease of $\lambda_{R_e}$ with decreasing SFR becomes modest or nearly absent: low-mass galaxies well below the SFMS, are fast-rotator early-type galaxies, and contain fast-rotating stellar disks like their star-forming counterparts. We also find a small but clear environmental dependence for the massive galaxies: in the mass range $10^{10.9}-10^{11.5} \mathcal{M}_{\odot}$, galaxies in rich groups or denser regions or classified as central galaxies have lower values of $\lambda_{R_e}$. While no environmental dependence is found for galaxies of lower mass. We discuss how our results can be understood as due to the different star formation and mass assembly histories of galaxies with varying mass.

The merging of double white dwarfs (WDs) may produce the events of accretion-induced collapse (AIC) and form single neutron stars (NSs). Meanwhile, it is also notable that the recently proposed WD+He subgiant scenario has a significant contribution to the production of massive double WDs, in which the primary WD grows in mass by accreting He-rich material from a He subgiant companion. In this work, we aim to study the binary population synthesis (BPS) properties of AIC events from the double WD mergers by considering the classical scenarios and also the contribution of theWD+He subgiant scenario to the formation of double WDs. First, we provided a dense and large model grid of WD+He star systems for producing AIC events through the double WD merger scenario. Secondly, we performed several sets of BPS calculations to obtain the rates and single NS number in our Galaxy. We found that the rates of AIC events from the double WD mergers in the Galaxy are in the range of 1.4-8.9*10^-3 yr^-1 for all ONe/CO WD+ONe/CO WD mergers, and in the range of 0.3-3.8*10^-3 yr^-1 when double COWD mergers are not considered.We also found that the number of single NSs from AIC events in our Galaxy may range from 0.328*10^7 to 1.072*10^8. The chirp mass of double WDs for producing AIC events distribute in the range of 0.55-1.25 Msun. We estimated that more than half of doubleWDs for producing AIC events are capable to be observed by the future space-based gravitational wave detectors.

The carbon-oxygen white dwarf (CO WD) + He star channel has been thought to be one of the promising scnarios to produce young type Ia supernovae (SNe Ia). Previous studies found that if the mass-accretion rate is greater than a critical value, the He-accreting CO WD will undergo inwardly propagating (off-centre) carbon ignition when it increases its mass close to the Chandrasekhar limit. The inwardly propagating carbon flame was supposed to reach the centre by previous works, leading to the production of an oxygen-neon (ONe) WD that may collapse into a neutron star but not an SN Ia. However, it is still uncertain how the carbon flame propagates under the effect of mixing mechanisms. In the present work, we aim to investigate the off-centre carbon burning of the He-accreting CO WDs by considering the effect of convective mixing. We found that the temperature of the flame is high enough to burn the carbon into silicon-group elements in the outer part of the CO core even if the convective overshooting is considered, but the flame would quench somewhere inside the WD, resulting in the formation of a C-O-Si WD. Owing to the inefficiency of thermohaline mixing, the C-O-Si WD may explode as an SN Ia if it continues to grow in mass. Our radiation transfer simulations show that the SN ejecta with the silicon-rich outer layer will form high-velocity absorption lines in Si II, leading to some similarities to a class of the high-velocity SNe Ia in the spectral evolution. We estimate that the birthrate of SNe Ia with Si-rich envelope is ~ 10^(-4)/yr in our galaxy.

The physics of the dark sector has remained one of the controversial areas of modern cosmology at present and hence it naturally attracts massive attention to the scientific community. With the developments of the astronomical data, the physics of the dark sector is becoming much more transparent than it was some twenty years back. The detection of gravitational waves (GWs) has now opened a cluster of possibilities in the cosmological regime. Being motivated by the detection of GWs and its possible impact on the physics of dark matter and dark energy, in this work we focus on the interacting dark energy models. Assuming the simplest possibility in which the vacuum energy with equation-of-state $w_x =-1$ is allowed to interact with the pressureless dark matter, we have extracted the constraints of the cosmological parameters. We find that the addition of the GWs data to the CMB measurements significantly improves up to a factor 4 of the parameters space, and up to a factor 2 for the full combination of current cosmological datasets, namely CMB+BAO+Pantheon+RSD+R16+CC+WL. The most affected parameters by the inclusion of the GWs are $\Omega_ch^2$, $\theta_{MC}$, $\xi$, and the derived parameters $\Omega_{m0}$, $\sigma_8$ and $H_0$.

We investigate what role the SKA neutral hydrogen (HI) intensity mapping (IM) sky survey observation will play in weighing neutrinos in cosmology. We use the simulated data of the baryon acoustic oscillation (BAO) measurements from the HI IM survey based on SKA1 and SKA2 to do the analysis. For the current observations, we use the Planck 2015 cosmic microwave background (CMB) anisotropies observation, the optical BAO measurements, the type Ia supernovae (SN) observation (Pantheon compilation), and the latest $H_0$ measurement. We consider three mass ordering cases for massive neutrinos, i.e., the normal hierarchy (NH), inverted hierarchy (IH), and degenerate hierarchy (DH) cases. It is found that the SKA observation can significantly improve the constraints on $\Omega_{\rm m}$ and $H_0$. Compared to the current observation, the SKA1 data can improve the constraints on $\Omega_{\rm m}$ by about 33%, and on $H_0$ by about 36%; the SKA2 data can improve the constraints on $\Omega_{\rm m}$ by about 58%, and on $H_0$ by about 66%. It is also found that the SKA observation can only slightly improve the constraints on $\sum m_\nu$. Compared to the current observation, the SKA1 data can improve the constraints on $\sum m_\nu$ by about 4%, 3%, and 10%, for the NH, IH, and DH cases, respectively; the SKA2 data can improve the constraints on $\sum m_\nu$ by about 7%, 7%, and 16%, for the NH, IH, and DH cases, respectively.

Galactic outflows play a major role in the evolution of galaxies, but the underlying physical processes are poorly understood. This is mainly because we have little information about the outflow structure, especially on large scales. In this paper, we probe the structure of galactic outflows in low-$z$ starburst by using a combination of ultra-violet spectroscopy and imaging of the fluorescence emission lines (associated with transitions to excited fine-structure levels) and spectroscopy of the corresponding strongly blue-shifted resonance absorption lines. We find that in the majority of cases the observed fluorescence emission lines are much weaker and narrower than the absorption lines, originating in the star-forming interstellar medium and/or the slowest-moving part of the inner outflow. In a minority of cases, the outflowing absorbing material does make a significant contribution to the fluorescence emission. These latter systems are characterized by both strong Ly$\alpha$ emission lines and weak low-ionization absorption lines (both known to be empirical signs of Lyman-continuum leakage). We argue that the observed weakness of emission from the outflow seen in the majority of cases is due to the missing emission arising on scales larger than those encompassed by the aperture of the {\it{Hubble Space Telescope}}. This implies shallow radial density profiles in these outflows, and suggests that most of the observed absorbing material must be created/injected at radii much larger than that of the starburst. This has important implications for our understanding of both the physics of galactic outflows and for our estimation of their principal properties.

Extreme ultraviolet (EUV) waves are impressive coronal propagating disturbances. They are closely associated with various eruptions, and can used for the global coronal seismology and the acceleration of solar energetic particles. Hence, the study of EUV waves plays an important role in solar eruptions and Space Weather. Here we present an EUV wave associated with a filament activation that did not evolve into any eruption. Due to the continuous magnetic flux emergence and cancellation around its one end, the filament rose with untwisting motion, and the filament mass flowed towards another end along the rising fields. Intriguingly, following the filament activation, an EUV wave formed with a fast constant speed ($\sim$500 km s$^{-1}$) ahead of the mass flow, and the overlying coronal loops expanded both in lateral and radial directions. Excluding the possibility of a remote flare and an absent coronal mass ejection, we suggest that the EUV wave was only closely associated with the filament activation. Furthermore, their intimate spacial and temporal relationship indicates that the EUV wave was likely directly triggered by the lateral expansion of overlying loops. We propose that the EUV wave can be interpreted as linear fast-mode wave, and the most vital key for the successful generation of the EUV wave is the impulsive early-phase lateral expansion of overlying loops that was driven by the activated filament mass flow without any eruption.

Strong gravitational lensing by galaxies provides us with a unique opportunity to understand the nature of gravity on galactic and extra-galactic scales. In this paper, we propose a new multimessenger approach using data from both gravitational wave (GW) and the corresponding electromagnetic (EM) counterpart to infer the constraint of the modified gravity (MG) theory denoted by the scale dependent phenomenological parameter. To demonstrate the robustness of this approach, we calculate the time-delay predictions by choosing various values of the phenomenological parameters and then compare them with that from the general relativity (GR). For the third generation ground-based GW observatory, with one typical strongly lensed GW+EM event, and assuming that the dominated error from the stellar velocity dispersions is 5\%, the GW time-delay data can distinguish an 18\% MG effect on a scale of tens of kiloparsecs with a $68\%$ confidence level. Assuming GR and a Singular Isothermal Sphere mass model, there exists a simplified consistency relationship between time-delay and imaging data. This relationship does not require for the velocity dispersion measurement, and hence can avoid major uncertainties. By using this relationship, the multimessenger approach is able to distinguish an $8\%$ MG effect. Our results show that the GW multimessenger approach can play an important role in revealing the nature of gravity on the galactic and extra-galactic scales.

While it is well known that mass transfer in binaries can pollute the surfaces of the accretors, it is still unclear whether this mechanism can reproduce the observed chemical inhomogeneities in globular clusters. We study the surface abundances of the accretors in low-mass binaries, as a first step towards understanding whether mass transfer in low-mass binaries is one of the potential origins of the aforementioned abundance anomalies in globular clusters. We use the Modules for Experiments in Stellar Astrophysics code to calculate binary evolutionary models with different initial donor masses between 0.9 and 1.9 $ \rm {M}_\odot$ for an initial metallicity of Z = 0.0034. The results show that in some low-mass binary systems, the accretors exhibit peculiar chemical patterns when they are still unevolved stars, e.g. C and O depletion; Na and N enhancement; constant Mg, Al and C+N+O. The abundance patterns of the accretors are significantly different from their initial abundances (or that of normal single stars), and can match the observed populations. These abundance patterns strongly depend not only on the initial parameters of binaries (donor mass, mass ratio, and orbital period), but also on the assumptions regarding mass-transfer efficiency and angular momentum loss. These results support the hypothesis that mass transfer in low-mass binaries is at least partly responsible for the unevolved anomalous stars in globular clusters. More work on binary evolutionary models and binary population synthesis is required to fully evaluate the contribution of this scenario.

If a fraction $f_{\rm dcdm}$ of the Dark Matter decays into invisible and massless particles (so-called "dark radiation") with the decay rate (or inverse lifetime) $\Gamma_{\rm dcdm}$, such decay will leave distinctive imprints on cosmological observables. With a full consideration of the Boltzmann hierarchy, we calculate the decay-induced impacts not only on the CMB but also on the redshift distortion and the kinetic Sunyaev-Zel'dovich effect, while providing detailed physical interpretations based on evaluating the evolution of gravitational potential. By using the current cosmological data with a combination of Planck 2015, Baryon Acoustic Oscillation and redshift distortion measurements which can improve the constraints, we update the $1\sigma$ bound on the fraction of decaying DM from $f_{\rm dcdm}\lesssim5.26\%$ to $f_{\rm dcdm}\lesssim2.73\%$ for the short-lived DM (assuming $\Gamma_{\rm dcdm}/H_0\gtrsim10^4$). However, no constraints are improved from RSD data ($f_{\rm dcdm}\lesssim0.94\%$) for the long-lived DM (i.e., $\Gamma_{\rm dcdm}/H_0\lesssim10^4$). We also find the fractional DM decay can only slightly reduce the $H_0$ and $\sigma_8$ tensions, which is consistent with other previous works. Furthermore, our calculations show that the kSZ effect in future would provide a further constraining power on the decaying DM.

In this paper, we present the application of a new method measuring Hubble parameter $H(z)$ by using the anisotropy of luminosity distance($d_{L}$) of the gravitational wave(GW) standard sirens of neutron star(NS) binary system. The method has never been put into practice so far due to the lack of the ability of detecting GW. However, LIGO's success in detecting GW of black hole(BH) binary system merger announced the potential possibility of this new method. We apply this method to several GW detecting projects, including Advanced LIGO(aLIGO), Einstein Telescope(ET) and DECIGO, and evaluate its constraint ability on cosmological parameters of $H(z)$. It turns out that the $H(z)$ by aLIGO and ET is of bad accuracy, while the $H(z)$ by DECIGO shows a good one. We simulate $H(z)$ data at every 0.1 redshift span using the error information of $H(z)$ by DECIGO, and put the mock data into the forecasting of cosmological parameters. Compared with the previous data and method, we get an obviously tighter constraint on cosmological parameters by mock data, and a concomitantly higher value of Figure of Merit(FoM, the reciprocal of the area enclosed by the $2\sigma$ confidence region). For a 3-year-observation by standard sirens of DECIGO, the FoM value is as high as 170.82. If a 10-year-observation is launched, the FoM could reach 569.42. For comparison, the FoM of 38 actual observed $H(z)$ data(OHD) is 9.3. We also investigate the undulant universe, which shows a comparable improvement on the constraint of cosmological parameters. These improvement indicates that the new method has great potential in further cosmological constraints.

So far most studies on the structure of coronal mass ejections (CMEs) are conducted through white-light coronagraphs, which demonstrate about one third of CMEs exhibit the typical three-part structure in the high corona (e.g., beyond 2 Rs), i.e., the bright front, the dark cavity and the bright core. In this paper, we address the CME structure in the low corona (e.g., below 1.3 Rs) through extreme-ultraviolet (EUV) passbands and find that the three-part CMEs in the white-light images can possess a similar three-part appearance in the EUV images, i.e., a leading edge, a low-density zone, and a filament or hot channel. The analyses identify that the leading edge and the filament or hot channel in the EUV passbands evolve into the front and the core later within several solar radii in the white-light passbands, respectively. What's more, we find that the CMEs without obvious cavity in the white-light images can also exhibit the clear three-part appearance in the EUV images, which means that the low-density zone in the EUV images (observed as the cavity in white-light images) can be compressed and/or transformed gradually by the expansion of the bright core and/or the reconnection of magnetic field surrounding the core during the CME propagation outward. Our study suggests that more CMEs can possess the clear three-part structure in their early eruption stage. The nature of the low-density zone between the leading edge and the filament or hot channel is discussed.

In this paper, we study the gravitational waveforms, polarizations and radiation powers of the first relativistic triple systems PSR J0337 + 1715, observed in 2014, by using the post-Newtonian approximations to their lowest order. Although they cannot be observed either by current or next generation of the detectors, they do provide useful information to test different theories of gravity. In particular, we carry out the studies in three different theories, general relativity (GR), Brans-Dicke (BD) gravity, and Einstein-aether ($\ae$) theory. The tensor modes $h_{+}$ and $h_{\times}$ exist in all three theories and have almost equal amplitudes. Their frequencies are all peaked at two locations, $ f^{(+, \times)}_1 \simeq 0.0686656 \mu$Hz and $ f^{(+, \times)}_2 \simeq 14.2138 \mu$Hz, which are about twice of the outer and inner orbital frequencies of the triple system, as predicted in GR. In $\ae$-theory, all the six polarization modes are different from zero, but the breathing $h_b$ and longitudinal $h_L$ modes are not independent and also peaked at two frequencies, but at the frequencies, $ f^{(b, L, \ae)}_{1} \simeq 0.0457771 \mu$Hz and $f^{(b, L, \ae)}_{ 2} \simeq 7.09545 \mu$Hz. A similar phenomenon is also observed in BD gravity, in which only the three modes $h_{+},\; h_{\times}$ and $h_{b}$ exit. We also study the radiation powers, and find that the quadrupole emission in each of the three theories has almost the same amplitude, but the dipole emission can be as big as the quadrupole emission in $\ae$-theory. This can provide a very promising window to obtain severe constraints on $\ae$-theory by the multi-band gravitational wave astronomy.

The source responsible for the reionization of the Universe is believed to be the population of star-forming galaxies at $z\sim6$ to 12. The biggest uncertainty concerns the fraction of Lyman-continuum photons that actually escape from the galaxies. In recent years, several relatively small samples of "leaky" galaxies have been uncovered, and clues have begun to emerge as to both the indirect signposts of leakiness and of the conditions/processes that enable the escape of ionizing radiation. In this paper we present the results of a pilot program aimed to test a new technique for finding leaky galaxies---using the weakness of the [SII] nebular emission-lines relative to typical star-forming galaxies as evidence that the interstellar medium is optically-thin to the Lyman continuum. We use the Cosmic Origins Spectrograph on the Hubble Space Telescope to detect significant emerging flux below the Lyman edge in two out of three [SII]-weak star-forming galaxies at $z\sim0.3$. We show that these galaxies differ markedly in their properties from the class of leaky "Green-Pea" galaxies at similar redshifts: our sample galaxies are more massive, more metal-rich, and less extreme in terms of their stellar population and the ionization state of the interstellar medium. Like the Green Peas, they have exceptionally high star-formation rates per unit area. They also share some properties with the known leaky galaxies at $z\sim3$, but are significantly dustier. Our results validate a new way to identify local laboratories for exploring the processes that made it possible for galaxies to reionize the Universe.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

We consider a quintessence model with Yukawa interaction between dark energy and dark matter and constrain this model by employing the recent cosmological data including the updated cosmic microwave background (CMB) measurements from Planck 2015, the weak gravitational lensing measurements from Kilo Degree Survey (KiDS) and redshift-space distortions. We find that an interaction in the dark sector is compatible with observations. The updated Planck data can significantly improve the constraints compared with the previous results from Planck 2013, while the KiDS data has less constraining power than Planck. The Yukawa interaction model is found to be moderately favored by Planck and able to alleviate the discordance between weak lensing measurements and CMB measurements as previously inferred from the standard Lambda cold dark matter model. N-body simulations for Yukawa interaction model is also performed. We find that using the halo density profile is plausible to improve the constraints significantly in the future.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We investigate predicted circumstellar properties of Type Ia supernova progenitor systems with non-degenerate helium star donors. It has been suggested that systems consisting of a carbon+oxygen white dwarf and a helium star can lead to Type Ia supernova explosions. Binary evolution calculations for the helium star donor channel predict that such a progenitor system is in either a stable helium-shell burning phase or a weak helium-shell flash phase at the time of the Type Ia supernova explosion. By taking the binary evolution models from our previous study, we show that a large fraction of the progenitor systems with a helium star donor have low enough density to explain the current non-detection of radio emission from Type Ia supernovae. Most of the progenitor systems in the weak helium-shell flash phase at the time of the Type Ia supernova explosions, which may dominate the prompt (short delay time) Type Ia supernova population, have both low circumstellar density and a faint helium star donor to account for the non-detection of radio emission and a pre-explosion companion star in SN 2011fe and SN 2014J. We also find some progenitor systems that are consistent with the properties of the companion star candidate identified at the explosion location of Type Iax SN 2012Z.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

It has been suggested that the Z-mode instability driven by energetic electrons with a loss-cone type velocity distribution is one candidate process behind the continuum and zebra pattern of solar type-IV radio bursts. Both the temperature of background plasma ($T_0$) and the energy of energetic electrons ($v_e$) are considered to be important to the variation of the maximum growth rate ($\gamma_{max}$). Here we present a detailed parameter study on the effect of $T_0$ and $v_e$, within a regime of the frequency ratio ($10 \leq \frac{\omega_{pe}}{\Omega_{ce}} \leq 30$). In addition to $\gamma_{max}$, we also analyze the effect on the corresponding wave frequency ($\omega^r_{max}$) and propagation angle ($\theta_{max}$). We find that (1) $\gamma_{max}$ in-general decreases with increasing $v_e$, while its variation with $T_0$ is more complex depending on the exact value of $v_e$; (2) with increasing $T_0$ and $v_e$, $\omega^r_{max}$ presents step-wise profiles with jumps separated by gradual or very-weak variations, and due to the warm-plasma effect on the wave dispersion relation $\omega^r_{max}$ can vary within the hybrid band (the harmonic band containing the upper hybrid frequency) and the band higher; (3) the propagation is either perpendicular or quasi-perpendicular, and $\theta_{max}$ presents variations in line with those of $\omega^r_{max}$, as constrained by the resonance condition. We also examine the profiles of $\gamma_{max}$ with $\frac{\omega_{pe}}{\Omega_{ce}}$ for different combinations of $T_0$ and $v_e$ to clarify some earlier calculations which show inconsistent results.

Super-Chandrasekhar remnants of double white dwarf mergers could sometimes collapse into a rapidly rotating neutron star (NS), accompanying with a mass ejection of a few times $0.01M_{\odot}$. Bright optical transient emission can be produced by the ejecta due to heating by radioactivities and particularly by energy injection from the NS. Since the merger remnants before collapse resemble a star evolving from the asymptotic giant branch phase to the planetary nebula phase, an intense dusty wind is considered to be driven about several thousand years ago before the collapse and surround the remnant at large radii. Therefore, the optical transient emission can be somewhat absorbed and scattered by the dusty wind, which can suppress the peak emission and cause a scattering plateau in optical light curves. Several years later, as the ejecta finally catches up with the wind material, the shock interaction between them can further give rise to a detectable radio transient emission on a timescale of tens of days. Discovery of and observations to such dust-affected optical transients and shock-driven radio transients can help to explore the nature of super-Chandrasekhar merger remnants and as well as the density and type ratios of double white dwarf systems, which is beneficial in assessing their gravitational wave contributions.

The Cosmology Large Angular Scale Surveyor (CLASS) is mapping the polarization of the Cosmic Microwave Background (CMB) at large angular scales ($2<\ell\lesssim200$) in search of a primordial gravitational wave B-mode signal down to a tensor-to-scalar ratio of $r \approx 0.01$. The same data set will provide a near sample-variance-limited measurement of the optical depth to reionization. Between June 2016 and March 2018, CLASS completed the largest ground-based Q-band CMB survey to date, covering over 31 000~square-degrees (75% of the sky), with an instantaneous array noise-equivalent temperature (NET) sensitivity of $32~\mu \mbox{K}_{cmb}\sqrt{\mbox{s}}$. We demonstrate that the detector optical loading ($1.6~\mbox{pW}$) and noise-equivalent power ($19~\mbox{aW}\sqrt{\mbox{s}}$) match the expected noise model dominated by photon bunching noise. We derive a $13.1\pm0.3~\mbox{K/pW}$ calibration to antenna temperature based on Moon observations, which translates to an optical efficiency of $0.48\pm0.04$ and a $27~\mbox{K}$ system noise temperature. Finally, we report a Tau A flux density of $308\pm11~\mbox{Jy}$ at $38.4\pm0.2~\mbox{GHz}$, consistent with the WMAP Tau A time-dependent spectral flux density model.

The current generation of radio and millimeter telescopes, particularly the Atacama Large Millimeter Array (ALMA), offers enormous advances in observing capabilities. While these advances represent an unprecedented opportunity to facilitate scientific understanding, the increased complexity in the spatial and spectral structure of these ALMA data cubes lead to challenges in their interpretation. In this paper, we perform a feasibility study for applying topological data analysis and visualization techniques never before tested by the ALMA community. Through techniques based on contour trees, we seek to improve upon existing analysis and visualization workflows of ALMA data cubes, in terms of accuracy and speed in feature extraction. We review our application development process in building effective analysis and visualization capabilities for the astrophysicists. We also summarize effective design practices by identifying domain-specific needs of simplicity, integrability, and reproducibility, in order to best target and service the large astrophysics community.

We combine constraints from linear and nonlinear scales, for the first time, to study the interaction between dark matter and dark energy. We devise a novel N-body simulation pipeline for cosmological models beyond $\Lambda$CDM. This pipeline is fully self-consistent and opens a new window to study the nonlinear structure formation in general phenomenological interacting dark energy models. By comparing our simulation results with the SDSS galaxy-galaxy weak lensing measurements, we are able to constrain the strength of interaction between dark energy and dark matter. Compared with the previous studies using linear examinations, we point to plausible improvements on the constraints of interaction strength by using small scale information from weak lensing. This improvement is mostly due to the sensitivity of weak lensing measurements on nonlinear structure formation at low redshift. With this new pipeline, it is possible to look for smoking gun signatures of dark matter-dark energy interaction.

The carbon-oxygen white dwarf (CO WD) + He star channel is one of the promising ways for producing type Ia supernovae (SNe Ia) with short delay times. Recent studies found that carbon under the He-shell can be ignited if the mass-accretion rate of CO WD is higher than a critical rate (about 2 * 10^-6 Msun/yr), triggering an inwardly propagating carbon flame. Previous studies usually supposed that the off-centre carbon flame would reach the centre, resulting in the formation of an oxygen-neon (ONe) WD that will collapse into a neutron star. However, the process of off-centre carbon burning is not well studied. This may result in some uncertainties on the final fates of CO WDs. By employing MESA, we simulated the long-term evolution of off-centre carbon burning in He-accreting CO WDs. We found that the inwardly propagating carbon flame transforms the CO WDs into OSi cores directly but not ONe cores owing to the high temperature of the burning front. We suggest that the final fates of the CO WDs may be OSi WDs under the conditions of off-centre carbon burning, or explode as iron-core-collapse SNe if the mass-accretion continues. We also found that the mass-fractions of silicon in the OSi cores are sensitive to the mass-accretion rates.

We study the time-dependent modulation effect and derive the local interstellar spectra (LIS) for the cosmic ray (CR) proton, helium, boron and carbon. A two-dimensional modulation model including the variation of the interplanetary environment with time is adopted to describe modulation process. The propagation equation of CRs in the heliosphere is numerically solved by the package Solarprop. We derive the LIS by fitting the latest results of several experiments, including Voyager 1, PAMELA, BESS-POLARII and ACE, during low solar activity periods. We further study the modulation in the polarity reversal periods with the PAMELA proton data. We find that the rigidity dependence of the diffusion coefficient is critical to explain the modulation effect during reversal periods. Our results also indicate a power law relation between the diffusion coefficient and the magnitude of the heliospheric magnetic field (HMF) at the Earth.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

Gravitationally bound hierarchies containing three or more components are very common in our Universe. In this paper we study {\em periodic} gravitational wave (GW) form, their polarizations, response function, its Fourier transform, and energy loss rate of a triple system through three different channels of radiation, the scalar, vector and tensor modes, in Einstein-aether theory of gravity. In the weak-field approximations and with the recently obtained constraints of the theory, we first analyze the energy loss rate of a binary system, and find that the dipole contributions from the scalar and vector modes could be of the order of ${\cal{O}}\left(c_{14}\right){\cal{O}}\left(G_Nm/d\right)^2$, where $c_{14} \; (\equiv c_{1} + c_{4})$ is constrained to $c_{14} \lesssim {\cal{O}}\left(10^{-5}\right)$ by current observations, where $c_i$'s are the four coupling constants of the theory. On the other hand, the "strong-field" effects for a binary system of neutron stars are about six orders lower than that of GR. So, in this paper we ignore these "strong-field" effects and first develop the general formulas to the lowest post-Newtonian order, by taking the coupling of the aether field with matter into account. Within this approximation, we find that the scalar breather mode and the scalar longitudinal mode are all suppressed by a factor of ${\cal{O}}\left(c_{14}\right)$ with respect to the transverse-traceless modes ($h_{+}$ and $h_{\times}$), while the vectorial modes $(h_{X}$ and $h_{Y}$) are suppressed by a factor of $c_{13} \lesssim {\cal{O}}\left(10^{-15}\right)$. Applying the general formulas to a triple system with periodic orbits, we find that the corresponding GW form, response function, and its Fourier transform depend sensitively on both the configuration of the triple system and their orientations with respect to the detectors.

We consider the reionization process in a cosmological model in which dark matter interacts with dark energy. Using a semi-analytical reionization model, we compute the evolution of the ionized fraction in terms of its spatial average and linear perturbations. We show that certain types of interactions between dark matter and dark energy can significantly affect the reionization history. We calculate the 21 cm signals in the interaction models, and compare the results with the predictions of the ${\rm \Lambda CDM}$ model.

The recently announced result by EDGES points an unexpected excess in the 21 cm global brightness temperature from cosmic dawn at $z\sim 17$, potentially indicating new phenomena beyond the $\Lambda$CDM model. A generic cosmological model which allows conformal and disformal couplings between dark matter and dark energy is employed to investigate the impact on the 21cm absorption signal and understand the EDGES anomaly. After exploring a wide range of parameter space for couplings, we find that the coupling effects can lead to a moderate change in the Hubble parameter while a negligible change in the spin temperature in the early Universe. Consequently, the decrease of the Hubble parameter from the mixed conformal and disformal couplings can reproduce the 21cm absorption approximately in consistent with the EDGES result at $z=17.5$. However, there is still tension in corresponding parameter space between EDGES and other cosmological observations for this model.

The morphologies of the wavefronts of extreme ultraviolet (EUV) waves can shed light on their physical nature and driving mechanism that are still strongly debated. In reality, the wavefronts always deform after interacting with ambient coronal structures during their propagation. Here, we focus on the initial wavefront morphologies of four selected EUV waves that are closely associated with jets or flux rope eruptions, using the high spatio-temporal resolution observations and different perspectives from the Solar Dynamics Observatory and the Solar-Terrestrial Relations Observatory. For the jet-driven waves, the jets originated from one end of the overlying closed loops, and the arc-shaped wavefront formed around the other far end of the expanding loops. The extrapolated field lines of the Potential Field Source Surface model show the close relationships between the jets, the wavefronts, and the overlying closed loops. For the flux-rope-driven waves, the flux ropes (sigmoids) lifted off beneath the overlying loops, and the circular wavefronts had an intimate spatio-temporal relation with the expanding loops. All the results suggest that the configuration of the overlying loops and their locations relative to the erupting cores are very important for the formation and morphology of the wavefronts, and both two jet-driven waves and two flux-rope-driven waves are likely triggered by the sudden expansion of the overlying closed loops. We also propose that the wavefront of EUV wave is possibly integrated by a chain of wave components triggered by a series of separated expanding loops.

The double-degenerate model, involving the merger of double carbon-oxygen white dwarfs (CO WDs), is one of the two classic models for the progenitors of type Ia supernovae (SNe Ia). Previous studies suggested that off-centre carbon burning would occur if the mass-accretion rate (Macc) is relatively high during the merging process, leading to the formation of oxygen-neon (ONe) cores that may collapse into neutron stars. However, the off-centre carbon burning is still incompletely understood, especially when the inwardly propagating burning wave reaches the centre. In this paper, we aim to investigate the propagating characteristics of burning waves and the subsequently evolutionary outcomes of these CO cores. We simulated the long-term evolution of CO WDs that accrete CO-rich material by employing the stellar evolution code MESA on the basis of the thick-disc assumption. We found that the final outcomes of CO WDs strongly depend on Macc (Msun/yr) based on the thick-disc assumption, which can be divided into four regions: (1) explosive carbon ignition in the centre, then SNe Ia (Macc < 2.45*10^-6); (2) OSi cores, then neutron stars (2.45*10^-6 < Macc < 4.5*10^-6); (3) ONe cores, then e-capture SNe (4.5*10^-6 < Macc < 1.05*10^-5); (4) off-centre oxygen and neon ignition, then off-centre explosion or Si-Fe cores (Macc > 1.05*10^-5). Our results indicate that the final fates of double CO WD mergers are strongly dependent on the merging processes (e.g. slow merger, fast merger, composite merger, violent merger, etc.).

It has been suggested that SNe Ia could be produced in the condition of the violent merger scenario of the double-degenerate model, in which a thermonuclear explosion could be produced when the merging of double carbon-oxygen white dwarfs (CO WDs) is still ongoing. It has been recently found that the nucleus of the bipolar planetary nebula Henize 2-428 consists of double CO WDs that have a total mass of ~1.76Msun, a mass ratio of ~1 and an orbital period of ~4.2 hours, which is the first and only discovered progenitor candidate of SNe Ia predicted by the violent merger scenario. In this work, we aim to reproduce the evolutionary history of the central double CO WDs of Henize 2-428. We find that the planetary nebula Henize 2-428 may originate from a primordial binary that have a ~5.4Msun primary and a ~2.7Msun secondary with an initial orbital period of ~15.9 days. The double CO WDs are formed after the primordial binary experiencing two Roche-lobe overflows and two common-envelope ejection processes. According to our calculations, it takes about ~840 Myr for the double CO WDs to merge and form an SN Ia driven by the gravitational wave radiation after their birth. To produce the current status of Henize 2-428, a large common-envelope parameter is needed. We also estimate that the rate of SNe Ia from the violent merger scenario is at most 2.9*10-4 yr-1, and that the delay time is in the range of ~90 Myr to the Hubble time.

We devise a fully self-consistent simulation pipeline for the first time to study the interaction between dark matter and dark energy. We perform convergence tests and show that our code is accurate on different scales. Using the parameters constrained by Planck, Type Ia Supernovae, Baryon Acoustic Oscillations (BAO) and Hubble constant observations, we perform cosmological N-body simulations. We calculate the resulting matter power spectra and halo mass functions for four different interacting dark energy models. In addition to the dark matter density distribution, we also show the inhomogeneous density distribution of dark energy. With this new simulation pipeline, we can further refine and constrain interacting dark energy models.

The IceCube Collaboration has published four years of data and the observed neutrino flux is significantly in excess of the expected atmospheric background. Due to the steeply falling atmospheric background spectrum, events at the highest energies are most likely extraterrestrial. In our previous approach we have studied blazars as the possible origin of the High-Energy Starting Events (HESE) neutrino events at PeV energies. In this work we extend our study to include all HESE neutrinos (which does not include IC 170922A) at or above a reconstructed energy of 100 TeV, but below 1 PeV. We study the X-ray and $\gamma$-ray data of all ($\sim200$) 3LAC blazars that are positionally consistent with the neutrino events above 100 TeV to determine the maximum neutrino flux from these sources. This larger sample allows us to better constrain the scaling factor between the observed and maximum number of neutrino events. We find that when we consider a realistic neutrino spectrum and other factors, the number of neutrinos is in good agreement with the detected number of IceCube HESE events. We also show that there is no direct correlation between \Fermi-LAT $\gamma$-ray flux and the IceCube neutrino flux and that the expected number of neutrinos is consistent with the non-detection of individual bright blazars.

A recent observation points to an excess in the expected 21-cm brightness temperature from cosmic dawn. In this paper, we present an alternative explanation of this phenomenon, an interaction in the dark sector. Interacting dark energy models have been extensively studied recently and there is a whole variety of such in the literature. Here we particularize to a specific model in order to make explicit the effect of an interaction.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

We study scalar condensations around asymptotically Anti-de Sitter(AdS) regular reflecting shells. We show that the charged scalar field can condense around charged reflecting shells with both Dirichlet and Neumann boundary conditions. In particular, the radii of the asymptotically AdS hairy shells are discrete, which is similar to cases in asymptotically flat spacetimes. We also provide upper bounds for the radii of the hairy Dirichlet reflecting shells and above the bound, the scalar field cannot condense around the shell.

It has been suggested that the accretion-induced collapse (AIC) of an oxygen-neon white dwarf (ONe WD) to a neutron star is a theoretically predicted outcome in stellar evolution, likely relating to the formation of some neutron star systems. However, the progenitor models of AIC events are still not well studied, and recent studies indicated that CO WD+He star systems may also contribute to the formation of neutron star systems through AIC process when off-centre carbon ignition happens on the surface of the CO WD. In this work, I studied the single-degenerate (SD) model of AIC events in a systematic way, including the contribution of the CO WD+He star channel and the ONe WD+MS/RG/He star channels. Firstly, I gave the initial parameter space of these SD channels for producing AIC events in the orbital period--secondary mass plane based on detailed binary evolution computations. Then, according to a binary population synthesis approach, I gave the rates and delay times of AIC events for these SD channels based on their initial parameter space. I found that the rates of AIC events in our galaxy are in the range of $\sim0.3-0.9\times10^{-3}$\,yr$^{-1}$,and that their delay times are $>$30\,Myr. I also found that the ONe WD+He star channel is the main way to produce AIC events, and that the CO WD+He star channel cannot be ignored when studying the progenitors of AIC events.

We analyze simulated maps of the Cosmology Large Angular Scale Surveyor (CLASS) experiment and recover a nearly cosmic-variance limited estimate of the reionization optical depth $\tau$. We use a power spectrum-based likelihood to simultaneously clean foregrounds and estimate cosmological parameters in multipole space. Using software specifically designed to constrain $\tau$, the amplitude of scalar fluctuations $A_s$, and the tensor-to-scalar ratio $r$, we demonstrate that the CLASS experiment will be able to estimate $\tau$ within a factor of two of the full-sky cosmic variance limit allowed by cosmic microwave background polarization measurements. Additionally, we discuss the role of CLASS's $\tau$ constraint in conjunction with gravitational lensing of the CMB on obtaining a $\gtrsim4\sigma$ measurement of the sum of the neutrino masses.

The Cosmology Large Angular Scale Surveyor consists of four instruments performing a CMB polarization survey. Currently, the 40 GHz and first 90 GHz instruments are deployed and observing, with the second 90 GHz and a multichroic 150/220 GHz instrument to follow. The receiver is a central component of each instrument's design and functionality. This paper describes the CLASS receiver design, using the first 90 GHz receiver as a primary reference. Cryogenic cooling and filters maintain a cold, low-noise environment for the detectors. We have achieved receiver detector temperatures below 50 mK in the 40 GHz instrument for 85% of the initial 1.5 years of operation, and observed in-band efficiency that is consistent with pre-deployment estimates. At 90 GHz, less than 26% of in-band power is lost to the filters and lenses in the receiver, allowing for high optical efficiency. We discuss the mounting scheme for the filters and lenses, the alignment of the cold optics and detectors, stray light control, and magnetic shielding.

The search for inflationary primordial gravitational waves and the measurement of the optical depth to reionization, both through their imprint on the large angular scale correlations in the polarization of the cosmic microwave background (CMB), has created the need for high sensitivity measurements of polarization across large fractions of the sky at millimeter wavelengths. These measurements are subject to instrumental and atmospheric $1/f$ noise, which has motivated the development of polarization modulators to facilitate the rejection of these large systematic effects. Variable-delay polarization modulators (VPMs) are used in the Cosmology Large Angular Scale Surveyor (CLASS) telescopes as the first element in the optical chain to rapidly modulate the incoming polarization. VPMs consist of a linearly polarizing wire grid in front of a movable flat mirror. Varying the distance between the grid and the mirror produces a changing phase shift between polarization states parallel and perpendicular to the grid which modulates Stokes U (linear polarization at $45^\circ$) and Stokes V (circular polarization). The CLASS telescopes have VPMs as the first optical element from the sky; this simultaneously allows a lock-in style polarization measurement and the separation of sky polarization from any instrumental polarization further along in the optical path. The Q-band CLASS VPM was the first VPM to begin observing the CMB full time, starting in the Spring of 2016. The first W-band CLASS VPM was installed in the Spring of 2018.

The Cosmology Large Angular Scale Surveyor (CLASS) aims to detect and characterize the primordial B-mode signal and make a sample-variance-limited measurement of the optical depth to reionization. CLASS is a ground-based, multi-frequency microwave polarimeter that surveys 70% of the microwave sky every day from the Atacama Desert. The focal plane detector arrays of all CLASS telescopes contain smooth-walled feedhorns that couple to transition-edge sensor (TES) bolometers through symmetric planar orthomode transducer (OMT) antennas. These low noise polarization-sensitive detector arrays are fabricated on mono-crystalline silicon wafers to maintain TES uniformity and optimize optical efficiency throughout the wafer. In this paper, we discuss the design and characterization of the first CLASS 93 GHz detector array. We measure the dark parameters, bandpass, and noise spectra of the detectors and report that the detectors are photon-noise limited. With current array yield of 82%, we estimate the total array noise-equivalent power (NEP) to be 2.1 aW$\sqrt[]{\mathrm{s}}$.

Coronal jets are always produced by magnetic reconnection between emerging flux and pre-existing overlying magnetic fields. When the overlying field is vertical/obilique or horizontal, the coronal jet will appear as anemone type or two-sided-loop type. Most of observational jets are of the anemone type, and only a few of two-sided-loop jets have been reported. Using the high-quality data from New Vacuum Solar Telescope, Interface Region Imaging Spectrograph, and Solar Dynamics Observatory, we present an example of two-sided-loop jets simultaneously observed in the chromosphere, transition region, and corona. The continuous emergence of magnetic flux brought in successively emerging of coronal loops and the slowly rising of an overlying horizontal filament threads. Sequentially, there appeared the deformation of the loops, the plasmoids ejection from the loop top, and pairs of loop brightenings and jet moving along the untwisting filament threads. All the observational results indicate there exist magnetic reconnection between the emerging loops and overlying horizontal filament threads, and it is the first example of two-sided-loop jets associated with ejected plasmoids and twisted overlying fields.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

Inhomogeneous cosmological models such as the Stephani universes could, in principle, provide an explanation for the observed accelerated expansion of the Universe. Working with a concrete, popular model of the Stephani cosmology -- the Stephani-Dabrowski model, we found that it is entropically viable. We also comment on the energy conditions and the two-sheeted geometry of the spacetime. However, similar to the LTB models, despite satisfying the holographic principle, Stephani cosmology has difficulty satisfying all the constraints from observations.

Coronal loops interconnecting two active regions, called as interconnecting loops (ILs), are prominent large-scale structures in the solar atmosphere. They carry a significant amount of magnetic flux, therefore are considered to be an important element of the solar dynamo process. Earlier observations show that eruptions of ILs are an important source of CMEs. It is generally believed that ILs are formed through magnetic reconnection in the high corona (>150-200"), and several scenarios have been proposed to explain their brightening in soft X-rays (SXRs). Yet, the detailed IL formation process has not been fully explored and the associated energy release in the corona still remains unresolved. Here we report the complete formation process of a set of ILs connecting two nearby active regions, with successive observations by STEREO-A on the far side of the Sun and SDO and Hinode on the Earth side. We conclude that ILs are formed by gradual reconnection high in the corona, in line with earlier postulations. In addition, we show evidence supporting that ILs become brightened in SXRs and EUVs through heating at or close to the reconnection site in the corona (i.e., through direct heating process of reconnection), a process that has been largely overlooked in earlier studies on ILs

Extreme ultraviolet (EUV) waves, spectacular horizontally propagating disturbances in the low solar corona, always trigger horizontal secondary waves (SWs) when they encounter ambient coronal structure. We present a first example of upward SWs in a streamer-like structure after the passing of an EUV wave. The event occurred on 2017 June 1. The EUV wave happened during a typical solar eruption including a filament eruption, a CME, a C6.6 flare. The EUV wave was associated with quasi-periodic fast propagating (QFP) wave trains and a type II radio burst that represented the existence of a coronal shock. The EUV wave had a fast initial velocity of $\sim$1000 km s$^{-1}$, comparable to high speeds of the shock and the QFP wave trains. Intriguingly, upward SWs rose slowly ($\sim$80 km s$^{-1}$) in the streamer-like structure after the sweeping of the EUV wave. The upward SWs seemed to originate from limb brightenings that were caused by the EUV wave. All the results show the EUV wave is a fast-mode magnetohydrodynamic shock wave, likely triggered by the flare impulses. We suggest that part of the EUV wave was probably trapped in the closed magnetic fields of streamer-like structure, and upward SWs possibly resulted from the release of trapped waves in the form of slow-mode. It is believed that an interplay of the strong compression of the coronal shock and the configuration of the streamer-like structure is crucial for the formation of upward SWs.

Context. The companions of the exploding carbon-oxygen white dwarfs (CO WDs) for producing type Ia supernovae (SNe Ia) are still not conclusively confirmed. A red-giant (RG) star has been suggested to be the mass donor of the exploding WD, named as the symbiotic channel. However, previous studies on the this channel gave a relatively low rate of SNe Ia. Aims. We aim to systematically investigate the parameter space, Galactic rates and delay time distributions of SNe Ia from the symbiotic channel by employing a revised mass-transfer prescription. Methods. We adopted an integrated mass-transfer prescription to calculate the mass-transfer process from a RG star onto the WD. In this prescription, the mass-transfer rate varies with the local material states. Results. We evolved a large number of WD+RG systems, and found that the parameter space of WD+RG systems for producing SNe Ia is significantly enlarged. This channel could produce SNe Ia with intermediate and old ages, contributing to at most 5% of all SNe Ia in the Galaxy. Our model increases the SN Ia rate from this channel by a factor of 5. We suggest that the symbiotic systems RS Oph and T CrB are strong candidates for the progenitors of SNe Ia.

It has been suggested that accretion-induced collapse (AIC) is a non-negligible path for the formation of the observed neutron stars (NSs). An ONe white dwarf (WD) that accretes material from a He star may experience AIC process and eventually produce intermediate-mass binary pulsars (IMBPs), named as the ONe WD+He star scenario. Note that previous studies can only account for part of the observed IMBPs with short orbital periods. In this work, we investigate the evolution of about 900 ONe WD+He star binaries to explore the distribution of IMBPs. We found that the ONe WD+He star scenario could form IMBPs including pulsars with 5-340 ms spin periods and 0.75-1.38 Msun WD companions, in which the orbital periods range from 0.04 to 900 d. Compared with the 20 observed IMBPs, this scenario can cover the parameters of 13 sources in the final orbital period-WD mass plane and the Corbet diagram, most of which has short orbital periods. We found that the ONe WD+He star scenario can explain almost all the observed IMBPs with short orbital periods. This work can well match the observed parameters of PSR J1802-2124 (one of the two precisely observed IMBPs), providing a possible evolutional path for its formation. We also speculate that the compact companion of HD 49798 (a hydrogen depleted sdO6 star) may be not a NS based on the present work.

Type Ia supernovae (SNe Ia) play a prominent role in understanding the evolution of Universe. They are thought to be thermonuclear explosions of mass-accreting carbon-oxygen white dwarfs (CO WDs) in binaries, although the mass donors of the accreting WDs are still not well determined. In this article, I review recent studies on mass-accreting WDs, including H- and He-accreting WDs. I also review currently most studied progenitor models of SNe Ia, i.e., the single-degenerate model (including the WD+MS channel, the WD+RG channel and the WD+He star channel), the double-degenerate model (including the violent merger scenario) and the sub-Chandrasekhar mass model. Recent progress on these progenitor models is discussed, including the initial parameter space for producing SNe Ia, the binary evolutionary paths to SNe Ia, the progenitor candidates of SNe Ia, the possible surviving companion stars of SNe Ia, and some observational constraints, etc. Some other potential progenitor models of SNe Ia are also summarized, including the hybrid CONe WD model, the core-degenerate model, the double WD collision model, the spin-up/spin-down model, and the model of WDs near black holes.To date, it seems that two or more progenitor models are needed to explain the observed diversity among SNe Ia.

This paper presents the latest observations from the newly-built solar radio spectrograph at the \emph{Chashan Solar Observatory}. On July 18 2016, the spectrograph records a solar spike burst event, which has several episodes showing harmonic structures, with the second, third, and fourth harmonics. The lower harmonic radio spike emissions are observed later than the higher harmonic bands, and the temporal delay of the second (third) harmonic relative to the fourth harmonic is about 30\ --\ 40 (10) ms. Based on the electron cyclotron maser emission mechanism, we analyze possible causes of the temporal delay and further infer relevant coronal parameters, such as the magnetic field strength and the electron density at the radio source.

We constrain interacting dark matter and dark energy (IDMDE) models using a 450-degree-square cosmic shear data from the Kilo Degree Survey (KiDS) and the angular power spectra from Planck's latest cosmic microwave background measurements. We revisit the discordance problem in the standard Lambda cold dark matter ($\Lambda$CDM) model between weak lensing and Planck datasets and extend the discussion by introducing interacting dark sectors. The IDMDE models are found to be able to alleviate the discordance between KiDS and Planck as previously inferred from the $\Lambda$CDM model, and moderately favored by a combination of the two datasets.

The final outcomes of accreting ONe white dwarfs (ONe WDs) have been studied for several decades, but there are still some issues not resolved. Recently, some studies suggested that the deflagration of oxygen would occur for accreting ONe WDs with Chandrasekhar masses. In this paper, we aim to investigate whether ONe WDs can experience accretion-induced collapse (AIC) or explosions when their masses approach the Chandrasekhar limit. Employing the stellar evolution code modules for experiments in stellar astrophysics (MESA), we simulate the long-term evolution of ONe WDs by accreting CO material. The ONe WDs undergo weak multicycle carbon flashes during the mass-accretion process, leading to the mass increase of the WDs. We found that different initial WD masses and mass-accretion rates have influence on the evolution of central density and temperature. However, the central temperature cannot reach the explosive oxygen ignition temperature due to the neutrino cooling. This work implies that the final outcome of accreting ONe WDs is electron-capture induced collapse rather than thermonuclear explosion.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

This paper discusses three matter-of-principle methods for measuring the general relativity correction to the Newtonian values of the position of collinear Lagrangian points L1 and L2 of the Sun-Earth-satellite system. All approaches are based on time measurements. The first approach exploits a pulsar emitting signals and two receiving antennas located at L1 and L2, respectively. The second method is based on a relativistic positioning system based on the Lagrangian points themselves. These first two methods depend crucially on the synchronization of clocks at L1 and L2. The third method combines a pulsar and an artificial emitter at the stable points L4 or L5 forming a basis for the positioning of the collinear points L1 and L2. Further possibilities are mentioned and the feasibility of the measurements is considered.

Many studies have shown that either the nearby astrophysical source or dark matter (DM) annihilation/decay is required to explain the origin of high energy cosmic ray (CR) $e^\pm$, which are measured by many experiments, such as PAMELA and AMS-02. Recently, the Dark Matter Particle Explorer (DAMPE) collaboration has reported its first result of the total CR $e^\pm$ spectrum from $25 \,\mathrm{GeV}$ to $4.6 \,\mathrm{TeV}$ with high precision. In this work, we study the DM annihilation and pulsar interpretations of the DAMPE high energy $e^\pm$ spectrum. In the DM scenario, the leptonic annihilation channels to $\tau^+\tau^-$, $4\mu$, $4\tau$, and mixed charged lepton final states can well fit the DAMPE result, while the $\mu^+\mu^-$ channel has been excluded. In addition, we find that the mixed charged leptons channel would lead to a sharp drop at $\sim$ $\mathrm{TeV}$. However, these DM explanations are almost excluded by the observations of gamma-ray and CMB, unless some complicated DM models are introduced. In the pulsar scenario, we analyze 21 nearby known pulsars and assume that one of them is the primary source of high energy CR $e^\pm$.Considering the constraint from the Fermi-LAT observation of the $e^\pm$ anisotropy, we find that two pulsars are possible to explain the DAMPE data. Our results show that it is difficult to distinguish between the DM annihilation and single pulsar explanations of high energy $e^\pm$ with the current DAMPE result.

We estimate uncertainties of cosmological parameters for phenomenological interacting dark energy models using weak lensing convergence power spectrum and bispectrum. We focus on the bispectrum tomography and examine how well the weak lensing bispectrum with tomography can constrain the interactions between dark sectors, as well as other cosmological parameters. Employing the Fisher matrix analysis, we forecast parameter uncertainties derived from weak lensing bispectra with a two-bin tomography and place upper bounds on strength of the interactions between the dark sectors. The cosmic shear will be measured from upcoming weak lensing surveys with high sensitivity, thus it enables us to use the higher order correlation functions of weak lensing to constrain the interaction between dark sectors and will potentially provide more stringent results with other observations combined.

The double-degenerate (DD) model, involving the merging of massive double carbon-oxygen white dwarfs (CO WDs) driven by gravitational wave radiation, is one of the classical pathways for the formation of type Ia supernovae (SNe Ia). Recently, it has been proposed that the WD+He subgiant channel has a significant contribution to the production of massive double WDs, in which the primary WD accumulates mass by accreting He-rich matter from a He subgiant. We evolved about 1800 CO WD+He star systems and obtained a large and dense grid for producing SNe Ia through the DD model. We then performed a series of binary population synthesis simulations for the DD model, in which the WD+He subgiant channel is calculated by interpolations in this grid. According to our standard model, the Galactic birthrate of SNe Ia is about 2.4*10^{-3} yr^{-1} for the WD+He subgiant channel of the DD model; the total birthrate is about 3.7*10^{-3} yr^{-1} for all channels, reproducing that of observations. Previous theoretical models still have deficit with the observed SNe Ia with delay times <1 Gyr and >8 Gyr. After considering the WD+He subgiant channel, we found that the delay time distributions is comparable with the observed results. Additionally, some recent studies proposed that the violent WD mergers are more likely to produce SNe Ia based on the DD model. We estimated that the violent mergers through the DD model may only contribute to about 16% of all SNe Ia.

We report the detection of widespread CH$_2$OHCHO and HOCH$_2$CH$_2$OH emission in Galactic center giant molecular cloud Sagittarius B2 using the Shanghai Tianma 65m Radio Telescope. Our observations show for the first time that the spatial distribution of these two important prebiotic molecules extends over 15 arc-minutes, corresponding to a linear size of approximately 36 pc. These two molecules are not just distributed in or near the hot cores. The abundance of these two molecules seems to decrease from the cold outer region to the central region associated with star-formation activity. Results present here suggest that these two molecules are likely to form through a low temperature process. Recent theoretical and experimental studies demonstrated that prebiotic molecules can be efficiently formed in icy grain mantles through several pathways. However, these complex ice features cannot be directly observed, and most constraints on the ice compositions come from millimeter observations of desorbed ice chemistry products. These results, combined with laboratory studies, strongly support the existence of abundant prebiotic molecules in ices.

We report the discovery of two new unbound hypervelocity stars (HVSs) from the LAMOST spectroscopic surveys. They are respectively a B2V type star of ~ 7 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 502 km/s at a Galactocentric radius of ~ 21 kpc and a B7V type star of ~ 4 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 408 km/s at a Galactocentric radius of ~ 30 kpc. The origins of the two HVSs are not clear given their currently poorly measured proper motions. However, the future data releases of Gaia should provide proper motion measurements accurate enough to solve this problem. The ongoing LAMOST spectroscopic surveys are expected to yield more HVSs to form a statistical sample, providing vital constraint on understanding the nature of HVSs and their ejection mechanisms.

Type Ia supernovae (SNe Ia) play a crucial role in studying cosmology and galactic chemical evolution. They are thought to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs) when their masses reach the Chandrasekar mass limit in binaries. Previous studies have suggested that He novae may be progenitor candidates of SNe Ia. However, the mass retention efficiencies during He nova outbursts are still uncertain. In this article, we aim to study the mass retention efficiencies of He nova outbursts and to investigate whether SNe Ia can be produced through He nova outbursts. Using the stellar evolution code Modules for Experiments in Stellar Astrophysics, we simulated a series of multicycle He-layer flashes, in which the initial WD masses range from 0.7 to 1.35 Msun with various accretion rates. We obtained the mass retention efficiencies of He nova outbursts for various initial WD masses, which can be used in the binary population synthesis studies. In our simulations, He nova outbursts can increase the mass of the WD to the Chandrasekar mass limit and the explosive carbon burning can be triggered in the center of the WD; this suggests that He nova outbursts can produce SNe Ia. Meanwhile, the mass retention efficiencies in the present work are lower than those of previous studies, which leads to a lower birthrates of SNe Ia through the WD + He star channel. Furthermore, we obtained the elemental abundances distribution at the moment of explosive carbon burning, which can be used as the initial input parameters in studying explosion models of SNe Ia.

He accretion onto carbon-oxygen white dwarfs (CO WDs) plays a fundamental role when studying the formation of type Ia supernovae (SNe Ia). Employing the MESA stellar evolution code, we calculated the long-term evolution of He-accreting CO WDs. Previous studies usually supposed that a WD can grow in mass to the Chandrasekhar limit in the stable He burning region and finally produce a SN Ia. However, in this study we find that off-centre carbon ignition occurs in the stable He burning region if the accretion rate is above a critical value (~2.05*10^{-6}{M$_\odot$}{yr}^-1), resulting in accretion-induced collapse rather than a SN Ia. If the accretion rate is below the critical value, explosive carbon ignition will eventually happen in the centre producing a SN Ia. Taking into account the possibility of off-centre carbon ignition, we have re-determined the initial parameter space that produces SNe Ia in the He star donor channel, one of the promising channels to produce SNe Ia in young populations. Since this parameter space is smaller than was found in the previous study of Wang et al. (2009), the SN Ia rates are also correspondingly smaller. We also determined the chemical abundance profile of the He-accreting WDs at the moment of explosive carbon ignition, which can be used as initial input for SN Ia explosion models.

Context: Although type Ia supernovae (SNe Ia) play a key role in astrophysics, the companions of the exploding carbon-oxygen white dwarfs (CO WDs) are still not completely identified. It has been suggested recently that a He-rich WD (a He WD or a hybrid HeCO WD) merges with a CO WD may produce an SN Ia. This theory was based on the double-detonation model, in which the shock compression in the CO core caused by the surface explosion of the He-rich shell might lead to the explosion of the whole CO WD. However, so far, very few binary population synthesis (BPS) studies have been made on the merger scenario of a CO WD and a He-rich WD in the context of SNe Ia. Aims: We aim to systematically study the Galactic birthrates and delay-time distributions of SNe Ia based on the merger scenario of a CO WD and a He-rich WD. Methods: We performed a series of Monte Carlo BPS simulations to investigate the properties of SNe Ia from the merging of a CO WD and a He-rich WD based on the Hurley rapid binary evolution code. We also considered the influence of different metallicities on the final results. Results: From our simulations, we found that no more than 15% of all SNe Ia stem from the merger scenario of a CO WD and a He-rich WD, and their delay times range from ~110 Myr to the Hubble time. This scenario mainly contributes to SN Ia explosions with intermediate and long delay times. The present work indicates that the merger scenario of a CO WD and a He-rich WD can roughly reproduce the birthrates of SN 1991bg-like events, and cover the range of their delay times. We also found that SN Ia birthrates from this scenario would be higher for the cases with low metallicities.

Double coronal hard X-ray (HXR) sources are believed to be critical observational evidence of bi-directional energy release through magnetic reconnection in a large-scale current sheet in solar ares. Here we present a study on double coronal sources observed in both HXR and microwave regimes, revealing new characteristics distinct from earlier reports. This event is associated with a footpoint-occulted X1.3-class flare (25 April 2014, starting at 00:17 UT) and a coronal mass ejection that are likely triggered by the magnetic breakout process, with the lower source extending upward from the top of the partially-occulted flare loops and the upper source co-incident with rapidly squeezing-in side lobes (at a speed of ~250 km/s on both sides). The upper source can be identified at energies as high as 70-100 keV. The X-ray upper source is characterized by flux curves different from the lower source, a weak energy dependence of projected centroid altitude above 20 keV, a shorter duration and a HXR photon spectrum slightly-harder than those of the lower source. In addition, the microwave emission at 34 GHz also exhibits a similar double source structure and the microwave spectra at both sources are in line with gyro-synchrotron emission given by non- thermal energetic electrons. These observations, especially the co-incidence of the very-fast squeezing-in motion of side lobes and the upper source, indicate that the upper source is associated with (possibly caused by) this fast motion of arcades. This sheds new lights on the origin of the corona double-source structure observed in both HXRs and microwaves.

Type-I bursts (i.e. noise storms) are the earliest-known type of solar radio emission at the metre wavelength. They are believed to be excited by non-thermal energetic electrons accelerated in the corona. The underlying dynamic process and exact emission mechanism still remain unresolved. Here, with a combined analysis of extreme ultraviolet (EUV), radio and photospheric magnetic field data of unprecedented quality recorded during a type-I storm on 30 July 2011, we identify a good correlation between the radio bursts and the co-spatial EUV and magnetic activities. The EUV activities manifest themselves as three major brightening stripes above a region adjacent to a compact sunspot, while the magnetic field there presents multiple moving magnetic features (MMFs) with persistent coalescence or cancelation and a morphologically similar three-part distribution. We find that the type-I intensities are correlated with those of the EUV emissions at various wavelengths with a correlation coefficient of 0.7-0.8. In addition, in the region between the brightening EUV stripes and the radio sources there appear consistent dynamic motions with a series of bi-directional flows, suggesting ongoing small-scale reconnection there. Mainly based on the induced connection between the magnetic motion at the photosphere and the EUV and radio activities in the corona, we suggest that the observed type-I noise storms and the EUV brightening activities are the consequence of small-scale magnetic reconnection driven by MMFs. This is in support of the original proposal made by Bentely et al. (Solar Phys. 193, 227, 2000).

We estimated the ortho-{\rm{H$_2$O}} abundances of G267.9--1.1, G268.4--0.9, G333.1--0.4 and G336.5--1.5, four of the brightest ortho-{\rm{H$_2$O}} sources in the southern sky observed by the Submillimeter Wave Astronomy Satellite (ortho-{\rm{H$_2$O}} 1$_{10}$ -- 1$_{01}$ line, 556.936~GHz). The typical molecular clumps in our sample have H$_2$ column densities of $10 ^{22}$ to $10 ^{23}${\,}cm$^{-2}$ and ortho-{\rm{H$_2$O}} abundances of 10$^{-10}$. Compared with previous studies, the ortho-{\rm{H$_2$O}} abundances are at a low level, which can be caused by the low temperatures of these clumps. To estimate the ortho-{\rm{H$_2$O}} abundances, we used the CS $J = 2 \to 1$ line (97.98095~GHz) and CS $J = 5 \to 4$ (244.93556~GHz) line observed by{ the} Swedish-ESO 15\,m Submillimeter Telescope (SEST) to calculate the temperatures of the clumps and the 350~$\upmu$m dust continuum observed by{ the} Caltech Submillimeter Observatory (CSO) telescope to estimate the H$_2$ column densities. The observations of {\rm{N$_2$H$^+$}} ($J = 1 \to 0$) for these clumps were also acquired by SEST and the corresponding abundances were estimated. The {\rm{N$_2$H$^+$}} abundance in each clump shows a common decreasing trend toward the center and {a} typical abundance range from 10$^{-11}$ to 10$^{-9}$.

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

The result presented by the BOSS-SDSS Collaboration measuring the baryon acoustic oscillations of the Lyman-$\alpha$ forest from high-redshift quasars indicates a $2.5\sigma$ departure from the standard $\Lambda$-cold-dark-matter model. This is the first time that the evolution of dark energy at high redshifts has been measured, and the current results cannot be explained by simple generalizations of the cosmological constant. We show here that a simple phenomenological interaction in the dark sector provides a good explanation for this deviation, naturally accommodating the Hubble parameter obtained by BOSS, $H(z=2.34)=222 \pm 7 ~\mathrm{km~s^{-1}~Mpc^{-1}}$. By performing a global fit of the parameters with the inclusion of this new data set together with the Planck data for the interacting model, we are able to show that some interacting models have constraints for $H(2.34)$ and $D_\mathrm{A}(2.34)$ that are compatible with the ones obtained by the BOSS Collaboration, showing a better concordance than $\Lambda$CDM. We also show that the interacting models that have a small positive coupling constant, which helps alleviate the coincidence problem, are compatible with the cosmological observations. Adding the likelihood of these new baryon acoustic oscillations data shows an improvement in the global fit, although it is not statistically significant. The coupling constant could not be fully constrained by the data sets used, but the dark energy equation of state shows a slight preference for a value different from a cosmological constant.

Using the high-quality observations of the Solar Dynamics Observatory, we present the interaction of two filaments (F1 and F2) in a long filament channel associated with twin coronal mass ejections (CMEs) on 2016 January 26. Before the eruption, a sequence of rapid cancellation and emergence of the magnetic flux has been observed, which likely triggered the ascending of the west filament (F1). The east footpoints of rising F1 moved toward the east far end of the filament channel, accompanying with post-eruption loops and flare ribbons. It likely indicated a large-scale eruption involving the long filament channel, resulted from the interaction between F1 and the east filament (F2). Some bright plasma flew over F2, and F2 stayed at rest during the eruption, likely due to the confinement of its overlying lower magnetic field. Interestingly, the impulsive F1 pushed its overlying magnetic arcades to form the first CME, and F1 finally evolved into the second CME after the collision with the nearby coronal hole. We suggest that the interaction of F1 and the overlying magnetic field of F2 led to the merging reconnection that form a longer eruptive filament loop. Our results also provide a possible picture of the origin of twin CMEs, and show the large-scale magnetic topology of the coronal hole is important for the eventual propagation direction of CMEs.

We investigate phenomenological interactions between dark matter and dark energy and constrain these models by employing the most recent cosmological data including the cosmic microwave background radiation anisotropies from Planck 2015, Type Ia supernovae, baryon acoustic oscillations, the Hubble constant and redshift-space distortions. We find that the interaction in the dark sector parameterized as an energy transfer from dark matter to dark energy is strongly suppressed by the whole updated cosmological data. On the other hand, an interaction between dark sectors with the energy flow from dark energy to dark matter is proved in better agreement with the available cosmological observations. This coupling between dark sectors is needed to alleviate the coincidence problem.

The AMS-02 has measured the cosmic ray electron (plus positron) spectrum up to ~TeV with an unprecedent precision. The spectrum can be well described by a power law without any obvious features above 10 GeV. The satellite instrument Dark Matter Particle Explorer (DAMPE), which was launched a year ago, will measure the electron spectrum up to 10 TeV with a high energy resolution. The cosmic electrons beyond TeV may be attributed to few local cosmic ray sources, such as supernova remnants. Therefore, spectral features, such as cutoff and bumps, can be expected at high energies. In this work we give a careful study on the perspective of the electron spectrum beyond TeV. We first examine our astrophysical source models on the latest leptonic data of AMS-02 to give a self-consistent picture. Then we focus on the discussion about the candidate sources which could be electron contributors above TeV. Depending on the properties of the local sources (especially on the nature of Vela), DAMPE may detect interesting features in the electron spectrum above TeV in the future.

Growing evidence reveals universal hardening on various cosmic ray spectra, e.g. proton, positron, as well as antiproton fraction. Such universality may indicate they have a common origin. In this paper, we argue that these widespread excesses can be accounted for by a nearby supernova remnant surrounded by a giant molecular cloud. Secondary cosmic rays ($\rm p$, $\rm e^+$) are produced through the collisions between the primary cosmic ray nuclei from this supernova remnant and the molecular gas. Different from the background, which is produced by the ensemble of large amount of sources in the Milky Way, the local injected spectrum can be harder. The time-dependent transport of particles would make the propagated spectrum even harder. Under this scenario, the anomalies of both primary ($\rm p$, $\rm e^-$) and secondary ($\rm e^+$, $\rm \bar{p}/p$) cosmic rays can be properly interpreted. We further show that the TeV to sub-PeV anisotropy of proton is consistent with the observations if the local source is relatively young and lying at the anti-Galactic center direction.

Large-scale structure has been shown as a promising cosmic probe for distinguishing and constraining dark energy models. Using the growth index parametrization, we obtain an analytic formula for the growth rate of structures in a coupled dark energy model in which the exchange of energy-momentum is proportional to the dark energy density. We find that the evolution of $f \sigma_8$ can be determined analytically once we know the coupling, the dark energy equation of state, the present value of the dark energy density parameter and the current mean amplitude of dark matter fluctuations. After correcting the growth function for the correspondence with the velocity field through the continuity equation in the interacting model, we use our analytic result to compare the model's predictions with large-scale structure observations.

The origin of the progenitors of type Ia supernovae (SNe Ia) is still uncertain. The core-degenerate (CD) scenario has been proposed as an alternative way for the production of SNe Ia. In this scenario, SNe Ia are formed at the final stage of common-envelope evolution from a merger of a carbon-oxygen white dwarf (CO WD) with the CO core of an asymptotic giant branch companion. However, the birthrates of SNe Ia from this scenario are still not well determined. In this work, we performed a detailed investigation on the CD scenario based on a binary population synthesis approach. The SN Ia delay times from this scenario are basically in the range of 90Myr-2500Myr, mainly contributing to the observed SNe Ia with short and intermediate delay times although this scenario can also produce some old SNe Ia. Meanwhile, our work indicates that the Galactic birthrates of SNe Ia from this scenario are no more than 20% of total SNe Ia due to more careful treatment of mass transfer. Although the SN Ia birthrates in the present work are lower than those in Ilkov & Soker, the CD scenario cannot be ruled out as a viable mechanism for the formation of SNe Ia. Especially, SNe Ia with circumstellar material from this scenario contribute to 0.7-10% of total SNe Ia, which means that the CD scenario can reproduce the observed birthrates of SNe Ia like PTF 11kx. We also found that SNe Ia happen systemically earlier for a high value of metallicity and their birthrates increase with metallicity.

Hot channel (HC) structure, observed in the high-temperature passbands of the AIA/SDO, is regarded as one candidate of coronal flux rope which is an essential element of solar eruptions. Here we present the first radio imaging study of an HC structure in the metric wavelength. The associated radio emission manifests as a moving type-IV (t-IVm) burst. We show that the radio sources co-move outwards with the HC, indicating that the t-IV emitting energetic electrons are efficiently trapped within the structure. The t-IV sources at different frequencies present no considerable spatial dispersion during the early stage of the event, while the sources spread gradually along the eruptive HC structure at later stage with significant spatial dispersion. The t-IV bursts are characterized by a relatively-high brightness temperature ($\sim$ 10$^{7}$ $-$ 10$^{9}$ K), a moderate polarization, and a spectral shape that evolves considerably with time. This study demonstrates the possibility of imaging the eruptive HC structure at the metric wavelength and provides strong constraints on the t-IV emision mechanism, which, if understood, can be used to diagnose the essential parameters of the eruptive structure.

The physical connection among and formation mechanisms of various components of the prominence-horn cavity system remain elusive. Here we present observations of such a system, focusing on a section of the prominence that rises and separates gradually from the main body. This forms a configuration sufficiently simple to yield clues to the above issues. It is characterized by embedding horns, oscillations, and a gradual disappearance of the separated material. The prominence-horn structure exhibits a large amplitude longitudinal oscillation with a period of ~150 minutes and an amplitude of ~30 Mm along the trajectory defined by the concave horn structure. The horns also experience a simultaneous transverse oscillation with a much smaller amplitude (~3 Mm) and shorter period (~10-15 minutes), likely representative of a global mode of the large-scale magnetic structure. The gradual disappearance of the structure indicates that the horn, an observational manifestation of the field-aligned transition region separating the cool and dense prominence from the hot and tenuous corona, is formed due to the heating and diluting process of the central prominence mass, while most previous studies suggest that it is the opposite process, i.e., the cooling and condensation of coronal plasmas, to form the horn. This study also demonstrates how the prominence transports magnetic flux to the upper corona, a process essential for the gradual build-up of pre-eruption magnetic energy.

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

Solar coronal radio bursts are enhanced radio emission excited by energetic electrons accelerated during solar eruptions, studies on which are important for investigating the origin and physical mechanism of energetic particles and further diagnosing coronal parameters. Earlier studies suffered from a lack of simultaneous high-quality imaging data of the radio burst and the eruptive structure in the inner corona. Here we present a study on a complex solar radio eruption consisting of a type II and three reversely-drifting type III bursts, using simultaneous EUV and radio imaging data. It is found that the type II burst is closely associated with a propagating and evolving CME-driven EUV shock structure, originated initially at the northern shock flank and later transferred to the top part of the shock. This source transfer is co-incident with the presence of shock decay and enhancing signatures observed at the corresponding side of the EUV front. The electron energy accelerated by the shock at the flank is estimated to be $\sim$ 0.3 c by examining the imaging data of the fast-drifting herringbone structure of the type II burst. The reversely-drifting type III sources are found to be within the ejecta and correlated with a likely reconnection event therein. Implications on further observational studies and relevant space-weather forecasting techniques are discussed.

Recent observations suggest that some type Ia supernovae (SNe Ia) originate from the merging of two carbon-oxygen white dwarfs (CO WDs). Meanwhile, recent hydrodynamical simulations have indicated that the accretion-induced collapse may be avoided under certain conditions when double WDs merge violently. However, the properties of SNe Ia from this violent merger scenario are highly dependent on a particular mass-accretion stage, the so-called WD + He subgiant channel, during which the primary WD is able to increase its mass by accreting He-rich material from a He subgiant before the systems evolves into a double WD system. In this article, we aim to study this particular evolutionary stage systematically and give the properties of violent WD mergers. By employing the Eggleton stellar evolution code, we followed a large number of binary calculations and obtained the regions in parameter space for producing violent mergers based on the WD + He subgiant channel. According to these simulations, we found that the primary WDs can increase their mass by ~0.10-0.45 Msun during the mass-accretion stage. We then conducted a series of binary population synthesis calculations and found that the Galactic SN Ia birthrate from this channel is about 0.01-0.4*10-3 yr-1. This suggests that the violent WD mergers from this channel may only contribute to ~0.3%-10% of all SNe Ia in our Galaxy. The delay times of violent WD mergers from this channel are >1.7 Gyr, contributing to the SNe Ia in old populations. We also found that the WD + He subgiant channel is the dominent way for producing violent WD mergers that may be able to eventually explode as SNe Ia.

We perform a phase space analysis of a generalized modified gravity theory with nonminimally coupling between geometry and matter. We apply the dynamical system approach to this generalized model and find that in the cosmological context, different choices of Lagrangian density will apparently result in different phases of the Universe. By carefully choosing the variables, we prove that there is an attractor solution to describe the late time accelerating universe when the modified gravity is chosen in a simple power-law form of the curvature scalar. We further examine the temperature evolution based on the thermodynamic understanding of the model. Confronting the model with supernova type Ia data sets, we find that the nonminimally coupled theory of gravity is a viable model to describe the late time Universe acceleration.

With the observations of the Solar Dynamics Observatory, we present the slipping magnetic reconnections with multiple flare ribbons (FRs) during an X1.2 eruptive flare on 2014 January 7. A center negative polarity was surrounded by several positive ones, and there appeared three FRs. The three FRs showed apparent slipping motions, and hook structures formed at their ends. Due to the moving footpoints of the erupting structures, one tight semi-circular hook disappeared after the slippage along its inner and outer edge, and coronal dimmings formed within the hook. The east hook also faded as a result of the magnetic reconnection between the arcades of a remote filament and a hot loop that was impulsively heated by the under flare loops. Our results are accordant with the slipping magnetic reconnection regime in 3D standard model for eruptive flares. We suggest that complex structures of the flare is likely a consequence of the more complex flux distribution in the photosphere, and the eruption involves at least two magnetic reconnections.

Models where Dark Matter and Dark Energy interact with each other have been proposed to solve the coincidence problem. We review the motivations underlying the need to introduce such interaction, its influence on the background dynamics and how it modifies the evolution of linear perturbations. We test models using the most recent observational data and we find that the interaction is compatible with the current astronomical and cosmological data. Finally, we describe the forthcoming data sets from current and future facilities that are being constructed or designed that will allow a clearer understanding of the physics of the dark sector.

Photometric and spectroscopic observations of a slowly declining, luminous Type Ia supernova (SN Ia) SN 2011hr in the starburst galaxy NGC 2691 are presented. SN 2011hr is found to peak at $M_{B}=-19.84 \pm 0.40\,\rm{mag}$, with a post-maximum decline rate $\Delta$m$_{15}$(B) = 0.92 $\pm$ 0.03\,$\rm{mag}$. From the maximum-light bolometric luminosity, $L=(2.30 \pm 0.90) \times 10^{43}\,\rm{erg\,s^{-1}}$, we estimate the mass of synthesized \Nifs\ in SN 2011hr to be $M(\rm{^{56}Ni})=1.11 \pm 0.43\,M_{\sun}$. SN 2011hr appears more luminous than SN 1991T at around maximum light, and the absorption features from its intermediate-mass elements (IMEs) are noticeably weaker than the latter at similar phases. Spectral modeling suggests that SN 2011hr has the IMEs of $\sim$\,0.07 M$_{\sun}$ in the outer ejecta, which is much lower than the typical value of normal SNe Ia (i.e., 0.3 -- 0.4 M$_{\sun}$) and is also lower than the value of SN 1991T (i.e., $\sim$\,0.18 M$_{\sun}$). These results indicate that SN 2011hr may arise from a Chandrasekhar-mass white dwarf progenitor that experienced a more efficient burning process in the explosion. Nevertheless, it is still possible that SN 2011hr may serve as a transitional object connecting the SN 1991T-like SNe Ia with the superluminous subclass like SN 2007if given that the latter also shows very weak IMEs at all phases.

We study a class of early dark energy models which has substantial amount of dark energy in the early epoch of the universe. We examine the impact of the early dark energy fluctuations on the growth of structure and the CMB power spectrum in the linear approximation. Furthermore we investigate the influence of the interaction between the early dark energy and the dark matter and its effect on the structure growth and CMB. We finally constrain the early dark energy model parameters and the coupling between dark sectors by confronting to different observations.

Supernovae of type Ia (SNe Ia) are believed to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs). However, the mass accretion process onto CO WDs is still not completely understood. In this paper, we study the accretion of He-rich matter onto CO WDs and explore a scenario in which a strong wind forms on the surface of the WD if the total luminosity exceeds the Eddington limit. Using a stellar evolution code called modules for experiments in stellar astrophysics (MESA), we simulated the He accretion process onto CO WDs for WDs with masses of 0.6-1.35Msun and various accretion rates of 10^{-8}-10^{-5}Msun/yr. If the contribution of the total luminosity is included when determining the Eddington accretion rate, then a super-Eddington wind could be triggered at relatively lower accretion rates than those of previous studies based on steady-state models. The super-Eddington wind can prevent the WDs with high accretion rates from evolving into red-giant-like He stars. We found that the contributions from thermal energy of the WD are non-negligible, judging by our simulations, even though the nuclear burning energy is the dominating source of luminosity. We also provide the limits of the steady He-burning regime in which the WDs do not lose any accreted matter and increase their mass steadily, and calculated the mass retention efficiency during He layer flashes for various WD masses and accretion rates. These obtained results can be used in future binary population synthesis computations.

The tensor-vector-scalar (TeVeS) model is considered a viable theory of gravity. It produces the Milgrom's modified Newtonian dynamics in the nonrelativistic weak field limit and is free from ghosts. This model has been tested against various cosmological observations. Here we investigate whether new observations such as the galaxy velocity power spectrum measured by 6dF and the kinetic Sunyaev Zel'dovich effect power spectrum measured by ACT/SPT can put further constraints on the TeVeS model. Furthermore, we perform the test of TeVeS cosmology with a sterile neutrino by confronting to Planck data, and find that it is ruled out by cosmic microwave background measurements from the Planck mission.

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

Type Iax supernovae (SNe Iax) are proposed as one new sub-class of SNe Ia since they present observational properties that are sufficiently distinct from the bulk of SNe Ia. SNe Iax are the most common of all types of peculiar SNe by both number and rate, with an estimated rate of occurrence of about 5-30% of the total SN Ia rate. However, the progenitor systems of SNe Iax are still uncertain. Analyzing pre-explosion images at SN Iax positions provides a direct way to place strong constraints on the nature of progenitor systems of SNe Iax. In this work, we predict pre-explosion properties of binary companion stars in a variety of potential progenitor systems by performing detailed binary evolution calculations with the one-dimensional stellar evolution code STARS. This will be helpful for constraining progenitor systems of SNe Iax from their pre-explosion observations. With our binary evolution calculations, it is found that the non-degenerate helium (He) companion star to both a massive C/O WD (> 1.1 solar mass) and a hybrid C/O/Ne WD can provide an explanation for the observations of SN~2012Z-S1, but the hybrid WD+He star scenario is more favorable.

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

Full-sky maps of the cosmic microwave background temperature reveal a 7% asymmetry of fluctuation power between two halves of the sky. A common phenomenological model for this asymmetry is an overall dipole modulation of statistically isotropic fluctuations, which produces particular off-diagonal correlations between multipole coefficients. We compute these correlations and construct corresponding estimators for the amplitude and direction of the dipole modulation. Applying these estimators to various cut-sky temperature maps from Planck and WMAP data shows consistency with a dipole modulation, differing from a null signal at 2.5$\sigma$, with an amplitude and direction consistent with previous fits based on the temperature fluctuation power. The signal is scale dependent, with a statistically significant amplitude at angular scales larger than 2 degrees. Future measurements of microwave background polarization and gravitational lensing can increase the significance of the signal. If the signal is not a statistical fluke in an isotropic Universe, it requires new physics beyond the standard model of cosmology.

Hyper-velocity stars are believed to be ejected out from the Galactic center through dynamical interactions between (binary) stars and the central massive black hole(s). In this paper, we report 19 low mass F/G/K type hyper-velocity star candidates from over one mil- lion stars of the first data release of the LAMOST general survey. We determine the unbound probability for each candidate using a Monte-Carlo simulation by assuming a non-Gaussian proper-motion error distribution, Gaussian heliocentric distance and radial velocity error dis- tributions. The simulation results show that all the candidates have unbound possibilities over 50% as expected, and one of them may even exceed escape velocity with over 90% probabili- ty. In addition, we compare the metallicities of our candidates with the metallicity distribution functions of the Galactic bulge, disk, halo and globular cluster, and conclude that the Galactic bulge or disk is likely the birth place for our candidates.

IceCube has reported the detection of a diffuse TeV-PeV neutrino emission, for which the flat spectrum radio quasars (FSRQs) have been proposed to be the candidate sources. Here we assume that the neutrino flux from FSRQs is proportional to their gamma-ray ones, and obtain the gamma-ray/neutrino flux ratio by the diffuse gamma-ray flux from Fermi-LAT measurement of FSRQs and the diffuse neutrino flux detected by IceCube. We apply this ratio to individual FSRQs and hence predict their neutrino flux. We find that a large fraction of candidate FSRQs from the northern sky in the IceCube point source search has predicted neutrino flux above the IceCube upper limit; and for the sample of stacking search for neutrinos by IceCube, the predicted stacked flux is even larger than the upper limit of stacked flux by orders of magnitude. Therefore the IceCube limit from stacking searches, combined with the Fermi-LAT observations, already rejects FSRQs as the main sources of IceCube-detected diffuse neutrinos: FSRQs can only account for <10% (4%) of the IceCube-detected diffuse neutrino flux, according to the stacking searches from the whole (northern) sky. The derived neutrino/gamma-ray flux ratio is <<1, also implying that the gamma-ray emission from FSRQs cannot be produced by the secondary leptons and photons from the pion production processes. The caveat in the assumptions is discussed.

The Square Kilometer Array (SKA) is an international effort to build the world's largest radio telescope, with one square kilometer collecting area. Besides its ambitious scientific objectives, such as probing the cosmic dawn and cradle of life, SKA also demands several revolutionary technological breakthroughs, with ultra-high precision synchronisation of the frequency references for thousands of antennas being one of them. In this report, aimed at applications to SKA, we demonstrate a frequency reference synchronization and dissemination scheme with the phase noise compensation function placed at the client site. Hence, one central hub can be linked to a large number of client sites, forming a star-shaped topology. As a performance test, the 100 MHz reference signal from a Hydrogen maser clock is disseminated and recovered at two remote sites. Phase noise characteristics of the recovered reference frequency signal coincides with that of the hydrogen-maser source and satisfies SKA requirement.

HD 49798 (a hydrogen depleted subdwarf O6 star) with its massive white dwarf (WD) companion has been suggested to be a progenitor candidate of type Ia supernovae (SNe Ia). However, it is still uncertain whether the companion of HD 49798 is a carbon-oxygen (CO) WD or an oxygen-neon (ONe) WD. A CO WD will explode as an SN Ia when its mass grows approach to Chandrasekhar mass, while the outcome of an accreting ONe WD is likely to be a neutron star. We followed a series of Monte Carlo binary population synthesis approach to simulate the formation of ONe WD + He star systems. We found that there is almost no orbital period as large as HD 49798 with its WD companion in these ONe WD + He star systems based on our simulations, which means that the companion of HD 49798 might not be an ONe WD. We suggest that the companion of HD 49798 is most likely a CO WD, which can be expected to increase its mass to the Chandrasekhar mass limit by accreting He-rich material from HD 49798. Thus, HD 49798 with its companion may produce an SN Ia in its future evolution.

The core-degenerate (CD) scenario has been suggested to be a possible progenitor model of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (CO WD) merges with the hot CO core of a massive asymptotic giant branch (AGB) star during their common-envelope phase. However, the SN Ia birthrates for this scenario are still uncertain. We conducted a detailed investigation into the CD scenario and then gave the birthrates for this scenario using a detailed Monte Carlo binary population synthesis approach. We found that the delay times of SNe Ia from this scenario are ~70Myrs-1400Myrs, which means that the CD scenario contributes to young SN Ia populations. The Galactic SN Ia birthrates for this scenario are in the range of ~7.4*10^{-5} yr^{-1}-3.7*10^{-4}yr^{-1}, which roughly accounts for ~2-10% of all SNe Ia. This indicates that, under the assumptions made here, the CD scenario only contributes a small portion of all SNe Ia, which is not consistent with the results of Ilkov & Soker (2013).

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

The super-Eddington wind scenario has been proposed as an alternative way for producing type Ia supernovae (SNe Ia). The super-Eddington wind can naturally prevent the carbon--oxygen white dwarfs (CO WDs) with high mass-accretion rates from becoming red-giant-like stars. Furthermore, it works in low-metallicity environments, which may explain SNe Ia observed at high redshifts. In this article, we systematically investigated the most prominent single-degenerate WD+MS channel based on the super-Eddington wind scenario. We combined the Eggleton stellar evolution code with a rapid binary population synthesis (BPS) approach to predict SN Ia birthrates for the WD+MS channel by adopting the super-Eddington wind scenario and detailed mass-accumulation efficiencies of H-shell flashes on the WDs. Our BPS calculations found that the estimated SN Ia birthrates for the WD+MS channel are ~0.009-0.315*10^{-3}{yr}^{-1} if we adopt the Eddington accretion rate as the critical accretion rate, which are much lower than that of the observations (<10% of the observed SN Ia birthrates). This indicates that the WD+MS channel only contributes a small proportion of all SNe Ia. The birthrates in this simulation are lower than previous studies, the main reason of which is that new mass-accumulation efficiencies of H-shell flashes are adopted. We also found that the critical mass-accretion rate has a significant influence on the birthrates of SNe Ia. Meanwhile, the results of our BPS calculations are sensitive to the values of the common-envelope ejection efficiency.

Since the discovery of a neutrino flux in excess of the atmospheric background by the IceCube Collaboration, searches for the astrophysical sources have been ongoing. Due to the steeply falling background towards higher energies, the PeV events detected in three years of IceCube data are the most likely ones to be of extraterrestrial origin. Even excluding the PeV events detected so far, the neutrino flux is well above the atmospheric background, so it is likely that a number of sub-PeV events originate from the same astrophysical sources that produce the PeV events. We study the high-energy properties of AGN that are positionally coincident with the neutrino events from three years of IceCube data and show the results for event number 4. IC 4 is a event with a low angular error (7.1$^\circ$) and a large deposited energy of 165 TeV. We use multiwavelength data, including Fermi/LAT and X-ray data, to construct broadband spectra and present parametrizations of the broadband spectral energy distributions with logarithmic parabolas. Assuming the X-ray to {\gamma}-ray emission in blazars originates in the photoproduction of pions by accelerated protons, their predicted neutrino luminosity can be estimated. The measurements of the diffuse extragalactic background by Fermi/LAT gives us an estimate of the flux contributions from faint unresolved blazars. Their contribution increases the number of expected events by a factor of $\sim$2. We conclude that the detection of the IceCube neutrinos IC4, IC14, and IC20 can be explained by the integral emission of blazars, even though no individual source yields a sufficient energy output.

SNe Iax are proposed as one new sub-class of SNe Ia. SNe Iax have been estimated to account for ~5-30% of the total SN Ia rate, and most SNe Iax have been discovered in late-type galaxies. In addition, observations constrain the progenitor systems of some SN Iax progenitors have ages of <80Myr. Recently, the weak deflagration explosions of Ch-mass CO WDs seem to provide a viable physical scenario for SNe Iax. Here, comparing theoretical predictions from BPS calculations with SN Iax observations, we put constraints on the SD Ch-mass model as a possible Iax progenitor. A long delay times of >3Gyr and low SN rates of ~3e-5/yr are found in our RG donor channel, indicating that this channel is unlikely to produce SNe Iax. We predict that the Galactic rate from the MS (He) donor channel is ~1.5e-3/yr (3e-4/yr), which is consistent with the observed Iax rate. The short delay times in the He channel support the young host environments of SNe Iax. However, the relative long delay times in the MS donor channel are less favourable for the observational constraints on the Iax progenitor ages. Finally, we set an upper limit on the pre-SN mass-loss rate at < 10e-4 M_sun/yr (for V_wind=100 km/s). The delay times in the SD Ch-mass model do not account for that most SNe Iax are located in late-type galaxies. However, at least one Iax event (SN 2008ge) is hosted by a S0 galaxy with no signs of star formation. Moreover, current X-ray observations cannot rule out the SD Ch-mass model. Taking all these into account and considering the uncertainty of the observed rate for SNe Iax, we suggest that some SNe Iax may be produced from weak deflagrations of Ch-mass CO WDs in SD progenitors, especially in the He-star channel. This is consistent with recent analysis of HST observations for the SN Iax SN2012Z. However, this SD deflagration model is still unlikely to be the most common progenitor scenario for SNe Iax.

Many current and future dark matter and neutrino detectors are designed to measure scintillation light with a large array of photomultiplier tubes (PMTs). The energy resolution and particle identification capabilities of these detectors depend in part on the ability to accurately identify individual photoelectrons in PMT waveforms despite large variability in pulse amplitudes and pulse pileup. We describe a Bayesian technique that can identify the times of individual photoelectrons in a sampled PMT waveform without deconvolution, even when pileup is present. To demonstrate the technique, we apply it to the general problem of particle identification in single-phase liquid argon dark matter detectors. Using the output of the Bayesian photoelectron counting algorithm described in this paper, we construct several test statistics for rejection of backgrounds for dark matter searches in argon. Compared to simpler methods based on either observed charge or peak finding, the photoelectron counting technique improves both energy resolution and particle identification of low energy events in calibration data from the DEAP-1 detector and simulation of the larger MiniCLEAN dark matter detector.

The IceCube (IC) collaboration recently reported the detection of TeV-PeV extraterrestrial neutrinos whose origin is yet unknown. By the photon-neutrino connection in $pp$ and $p\gamma$ interactions, we use the \fermi-LAT observations to constrain the origin of the IC detected neutrinos. We find that Galactic origins, i.e., the diffuse Galactic neutrinos due to cosmic ray (CR) propagation in the Milky Way, and the neutrinos from the Galactic point sources, may not produce the IC neutrino flux, thus these neutrinos should be of extragalactic origin. Moreover, the extragalactic gamma-ray bursts (GRBs) may not account for the IC neutrino flux, the jets of active galactic nuclei may not produce the IC neutrino spectrum, but the starburst galaxies (SBGs) may be promising sources. As suggested by the consistency between the IC detected neutrino flux and the Waxman-Bahcall bound, GRBs in SBGs may be the sources of both the ultrahigh energy, $>10^{19}$eV, CRs and the $1-100$~PeV CRs that produce the IC detected TeV-PeV neutrinos.

We consider scenario in which our $(3+1)$-dim universe is actually a dense stack of multiple parallel $(2+1)$-dim branes. For this purpose, we generalize the DGP model to a multi-brane case. We solve for the propagation of the scalar field and gravity in this setup. At short distances (high momenta) along the branes interactions follow the $(2+1)$-dim laws, while at large distances (low momenta) interactions follow the usual $(3+1)$-dim laws. This feature is inherited from the original DGP model. In the direction perpendicular to the brane, we show that interactions become $(3+1)$-dim at low momenta which are unable to resolve the interbrane separation. Thus, this is one of the explicit constructs of the "vanishing dimensions" scenario where high energy physics appears to be lower dimensional rather than higher dimensional.

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

The accelerating expansion of the universe at recent epochs is encoded in the cosmic microwave background: a few percent of the total temperature fluctuations are generated by evolving gravitational potentials which trace the large-scale structures in the universe. This signature of dark energy, the Integrated Sachs-Wolfe Effect, has been detected by averaging temperatures in the WMAP sky maps corresponding to the directions of superstructures in the Sloan Digital Sky Survey data release 6. We model the maximum average peak signal expected in the standard $\Lambda$CDM cosmological model, using Gaussian random realizations of the microwave sky, including correlations between different physical contributions to the temperature fluctuations and between different redshift ranges of the evolving gravitational potentials. We find good agreement with the mean temperature peak amplitude from previous theoretical estimates based on large-scale structure simulations, but with larger statistical uncertainties. We apply our simulation pipeline to four different foreground-cleaned microwave temperature maps from Planck and WMAP data, finding a mean temperature peak signal at previously identified sky locations which exceeds our theoretical mean signal at a statistical significance of about $2.5\sigma$ and which differs from a null signal at $3.5\sigma$.

The hybrid CONe white dwarfs (WDs) have been suggested to be possible progenitors of type Ia supernovae (SNe Ia). In this article, we systematically studied the hybrid CONe WD + He star scenario for the progenitors of SNe Ia, in which a hybrid CONe WD increases its mass to the Chandrasekhar mass limit by accreting He-rich material from a non-degenerate He star. According to a series of detailed binary population synthesis simulations, we obtained the SN Ia birthrates and delay times for this scenario. The SN Ia birthrates for this scenario are ~0.033-0.539*10^(-3)yr^(-1), which roughly accounts for 1-18% of all SNe Ia. The estimated delay times are ~28Myr-178Myr, which are the youngest SNe Ia predicted by any progenitor model so far. We suggest that SNe Ia from this scenario may provide an alternative explanation of type Iax SNe. We also presented some properties of the donors at the point when the WDs reach the Chandrasekhar mass. These properties may be a good starting point for investigating the surviving companions of SNe Ia, and for constraining the progenitor scenario studied in this work.

The potential importance of the angular momentum which is gained by accreting white dwarfs (WDs) has been increasingly recognized in the context of type Ia supernova (SN Ia) single-degenerate model. The expectation that the spin of the WD can delay the explosion should help the single-degenerate model to be consistent with the observed properties of most SNe Ia, in particular by avoiding hydrogen contamination. In this article, we attempt to study the most prominent single-degenerate supersoft (WD + MS) channel when the rotation of accreting WDs is considered. We present a detailed binary population synthesis study to examine the predicted population of SNe Ia for this channel. For our standard model, we find that 77% of these SNe Ia explode with WD masses which are low enough to be supported by solid-body rotation (1.378-1.5Msun); this is a substantially higher proportion than found by previous work. Only 2% have WD explosion masses >2.0Msun; these require the initial WD mass to be larger than 1.0 Msun. We further discuss the possible origin of the diversity of SNe Ia from the pre- and post- accretion properties of the WDs in this population. We also suggest that some SN Ia progenitors with substantial circumstellar hydrogen, including some apparent type IIn SNe, might be related to WDs which required support from differential rotation to avoid explosion, since these can still be accreting from hydrogen-rich donors with a relatively high mass-transfer rate at the time of the SN explosion.

Among various factors affecting the space weather effects of a coronal mass ejection (CME), its propagation trajectory in the interplanetary space is an important one determining whether and when the CME will hit the Earth. Many direct observations have revealed that a CME may not propagate along a straight trajectory in the corona, but whether or not a CME also experiences a deflected propagation in the interplanetary space is a question, which has never been fully answered. Here by investigating the propagation process of an isolated CME from the corona to interplanetary space during 2008 September 12 -- 19, we present solid evidence that the CME was deflected not only in the corona but also in the interplanetary space. The deflection angle in the interplanetary space is more than 20 degrees toward the west, resulting a significant change in the probability the CME encounters the Earth. A further modeling and simulation-based analysis suggest that the cause of the deflection in the interplanetary space is the interaction between the CME and the solar wind, which is different from that happening in the corona.

Interacting Dark Energy and Dark Matter is used to go beyond the standard cosmology. We base our arguments on Planck data and conclude that an interaction is compatible with the observations and can provide a strong argument towards consistency of different values of cosmological parameters.

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

Extensive optical and ultraviolet (UV) observations of the type Ia supernova (SN Ia) 2012fr are presented in this paper. It has a relatively high luminosity, with an absolute $B$-band peak magnitude of about $-19.5$ mag and a smaller post-maximum decline rate than normal SNe Ia [e.g., $\Delta m _{15}$($B$) $= 0.85 \pm 0.05$ mag]. Based on the UV and optical light curves, we derived that a $^{56}$Ni mass of about 0.88 solar masses was synthesized in the explosion. The earlier spectra are characterized by noticeable high-velocity features of \ion{Si}{2} $\lambda$6355 and \ion{Ca}{2} with velocities in the range of $\sim22,000$--$25,000$ km s$^{-1}$. At around the maximum light, these spectral features are dominated by the photospheric components which are noticeably narrower than normal SNe Ia. The post-maximum velocity of the photosphere remains almost constant at $\sim$12,000 km s$^{-1}$ for about one month, reminiscent of the behavior of some luminous SNe Ia like SN 1991T. We propose that SN 2012fr may represent a subset of the SN 1991T-like SNe Ia viewed in a direction with a clumpy or shell-like structure of ejecta, in terms of a significant level of polarization reported in Maund et al. (2013).

A scenario is proposed to explain the preferential heating of minor ions and differential streaming velocity between minor ions and protons observed in the solar corona and in the solar wind. It is demonstrated by test particle simulations that minor ions can be nearly fully picked up by intrinsic Alfv\'en-cyclotron waves observed in the solar wind based on the observed wave energy density. Both high frequency ion-cyclotron waves and low frequency Alfv\'en waves play crucial roles in the pickup process. A minor ion can first gain a high magnetic moment through the resonant wave-particle interaction with ion-cyclotron waves, and then this ion with a large magnetic moment can be trapped by magnetic mirror-like field structures in the presence of the lower-frequency Alfv\'en waves. As a result, the ion is picked up by these Alfv\'en-cyclotron waves. However, minor ions can only be partially picked up in the corona due to low wave energy density and low plasma beta. During the pickup process, minor ions are stochastically heated and accelerated by Alfv\'en-cyclotron waves so that they are hotter and flow faster than protons. The compound effect of Alfv\'en waves and ion-cyclotron waves is important on the heating and acceleration of minor ions. The kinetic properties of minor ions from simulation results are generally consistent with in situ and remote features observed in the solar wind and solar corona.

The direct search for dark matter is entering a period of increased sensitivity to the hypothetical Weakly Interacting Massive Particle (WIMP). One such technology that is being examined is a scintillation only noble liquid experiment, MiniCLEAN. MiniCLEAN utilizes over 500 kg of liquid cryogen to detect nuclear recoils from WIMP dark matter and serves as a demonstration for a future detector of order 50 to 100 tonnes. The liquid cryogen is interchangeable between argon and neon to study the A$^{2}$ dependence of the potential signal and examine backgrounds. MiniCLEAN utilizes a unique modular design with spherical geometry to maximize the light yield using cold photomultiplier tubes in a single-phase detector. Pulse shape discrimination techniques are used to separate nuclear recoil signals from electron recoil backgrounds. MiniCLEAN will be spiked with additional $^{39}$Ar to demonstrate the effective reach of the pulse shape discrimination capability. Assembly of the experiment is underway at SNOLAB and an update on the project is given.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

The interaction between Dark Matter and Dark Energy has been proposed as a mechanism to alleviate the coincidence problem. We analyze the effect of the interaction on the evolution of the gravitational field and propose two new observables based on its effect on the matter peculiar velocity field. We find that for different model parameters the matter peculiar velocity could be factor 2 times larger or 5 times smaller than the amplitude of velocity perturbations in the concordance LCDM cosmological model at the same scales. We compare the effect on the peculiar velocities with those on the Integrated Sachs-Wolfe effect and we show that velocities can potentially provide constraints on the strength of the interaction stronger than those currently available. We show that the current upper limits on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum provide constraints on the coupling within the dark sectors that are consistent with those obtained previously from the CMB and galaxy clusters. In particular, we show that current upper limits from the Atacama Cosmology Telescope and the South Pole Telescope favor Dark Energy decaying into Dark Matter, as required to solve the coincidence problem.

The most favored progenitor scenarios for Type Ia supernovae (SNe Ia) involve the single-degenerate (SD) scenario and the double-degenerate scenario. The absence of stripped hydrogen (H) in the nebular spectra of SNe Ia challenges the SD progenitor models. Recently, it was shown that pure deflagration explosion models of Chandrasekhar-mass white dwarfs ignited off-center reproduce the characteristic observational features of 2002cx-like SNe Ia very well. In this work we predict, for the first time, the amount of stripped H for the off-center pure deflagration explosions. We find that their low kinetic energies lead to inefficient H mass stripping (less than 0.01 M_sun), indicating that the stripped H may be hidden in (observed) late-time spectra of SN 2002cx-like SNe Ia.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. However, their progenitors are still uncertain. We study whether the population properties of this class might be understood if the events originate from a subset of sub-Chandrasekhar mass explosions. In this potential progenitor population, a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star; ignition of that helium layer then leads to ignition of the CO WD. We incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain rates and delay times for such events. The predicted Galactic event rate of these explosions is ~1.5\times10^{-3}{yr}^{-1} according to our standard model, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are ~70Myr-800Myr, consistent with the fact that most of SNe Iax have been discovered in late-type galaxies. If the explosions are assumed to be double-detonations -- following current model expectations -- then based on the CO WD masses at explosion we also estimate the distribution of resulting SN brightness (-13 \gtrsim M_{bol} \gtrsim -19mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not. We suggest that the high mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might be necessary to produce SN 2002cx-like characteristics, perhaps even by changing the nature of the CO ignition.

In the version of the SD scenario of SNe Ia studied here, a CO WD explodes close to the Chandrasekhar limit after accreting material from a non-degenerate He companion. In the present study, we employ the Stellar GADGET code to perform 3D hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion taking into account its orbital motion and spin. It is found that only 2%--5% of the initial companion mass are stripped off from the outer layers of He companions due to the SN impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted in good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about > 10e-3 M_sun for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power law relation in good approximation.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. We show that the population properties of this class can be understood if the events originate from helium double-detonation sub-Chandrasekhar mass explosions, in which a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star. We have incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain birthrates and delay times for such events. The predicted Galactic event rate is $\sim$$0.6-1.8\times 10^{-3}\,{\rm yr}^{-1}$, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are $\sim$70\,Myr$-$710\,Myr, consistent with the fact that SNe Iax have so far only been discovered in late-type galaxies. Based on the CO WD mass at explosion and previous detonation models, we also estimate the distribution of resulting SN brightness ($-13 \gtrsim M_{\rm bol} \gtrsim -19$ mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not and suggest that the higher mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might provide the explanation.

Type Ia supernovae (SN Ia) are the most important standard candles for measuring the expansion history of the universe. The thermonuclear explosion of a white dwarf can explain their observed properties, but neither the progenitor systems nor any stellar remnants have been conclusively identified. Underluminous SN Ia have been proposed to originate from a so-called double-detonation of a white dwarf. After a critical amount of helium is deposited on the surface through accretion from a close companion, the helium is ignited causing a detonation wave that triggers the explosion of the white dwarf itself. We have discovered both shallow transits and eclipses in the tight binary system CD-30 11223 composed of a carbon/oxygen white dwarf and a hot helium star, allowing us to determine its component masses and fundamental parameters. In the future the system will transfer mass from the helium star to the white dwarf. Modelling this process we find that the detonation in the accreted helium layer is sufficiently strong to trigger the explosion of the core. The helium star will then be ejected at so large a velocity that it will escape the Galaxy. The predicted properties of this remnant are an excellent match to the so-called hypervelocity star US 708, a hot, helium-rich star moving at more than 750 km/s, sufficient to leave the Galaxy. The identification of both progenitor and remnant provides a consistent picture of the formation and evolution of underluminous type Ia supernovae.

In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

We collected 19.4 days of data from four 750 g TeO2 bolometers, and in three of them we were able to set the energy threshold around 3 keV using a new analysis technique. We found a background rate ranging from 25 cpd/keV/kg at 3 keV to 2 cpd/keV/kg at 25 keV, and a peak at 4.7 keV. The origin of this peak is presently unknown, but its presence is confirmed by a reanalysis of 62.7 kg.days of data from the finished CUORICINO experiment. Finally, we report the expected sensitivities of the CUORE0 (52 bolometers) and CUORE (988 bolometers) experiments to a WIMP annual modulation signal.

We reconstruct an $f(R)$ gravity model that gives rise to the particular $\Lambda$CDM background evolution of the universe. We find well-defined, real-valued analytical forms for the $f(R)$ model to describe the universe both in the early epoch from the radiation to matter dominated eras and the late time acceleration period. We further examine the viability of the derived $f(R)$ model and find that it is viable to describe the evolution of the universe in the past and there does not exist the future singularity in the Lagrangian.

The identity of the progenitor systems of SNe Ia is still uncertain. In the single-degenerate (SD) scenario, the interaction between the SN blast wave and the outer layers of a main sequence (MS) companion star strips off H-rich material which is then mixed into the ejecta. Strong contamination of the SN ejecta with stripped material could lead to a conflict with observations of SNe Ia. This constrains the SD progenitor model. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. We use Eggleton's stellar evolution code including the optically thick accretion wind model and the possibility of the effects of accretion disk instabilities to obtain realistic models of companions for different progenitor systems. The impact of the SN blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations employing the SPH code GADGET3. We find that the stripped masses range from 0.11 to 0.18 M_sun. The kick velocity is between 51 and 105 km/s. We find that the stripped mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion. They can be fitted by a power law for a given companion model. However, the structure of the companion star is also important for the amount of stripped material. With more realistic companion star models than in previous studies, our simulations show that the H masses stripped from companions are inconsistent with the best observational limits (< 0.01 M_sun) derived from nebular spectra. However, a rigorous forward modeling based on impact simulations with radiation transfer is required to reliably predict observable signatures of the stripped H and to conclusively assess the viability of the considered SN Ia progenitor scenario.

Waves play a crucial role in diagnosing the plasma properties of various structures in the solar corona and coronal heating. Slow magneto-acoustic (MA) waves are one of the important magnetohydrodynamic waves. In past decades, numerous slow MA waves were detected above the active regions and coronal holes, but rarely found elsewhere. Here, we investigate a `tornado'-like structure consisting of quasi-periodic streaks within a dark cavity at about 40--110 Mm above the quiet-Sun region on 2011 September 25. Our analysis reveals that these streaks are actually slow MA wave trains. The properties of these wave trains, including the phase speed, compression ratio, kinetic energy density, etc., are similar to those of the reported slow MA waves, except that the period of these waves is about 50 s, much shorter than the typical reported values (3--5 minutes).

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

We present a new measurement of the $\alpha$-spectroscopic factor ($S_\alpha$) and the asymptotic normalization coefficient (ANC) for the 6.356 MeV 1/2$^+$ subthreshold state of $^{17}$O through the $^{13}$C($^{11}$B, $^{7}$Li)$^{17}$O transfer reaction and we determine the $\alpha$-width of this state. This is believed to have a strong effect on the rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction, the main neutron source for {\it slow} neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. Based on the new width we derive the astrophysical S-factor and the stellar rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction. At a temperature of 100 MK our rate is roughly two times larger than that by \citet{cau88} and two times smaller than that recommended by the NACRE compilation. We use the new rate and different rates available in the literature as input in simulations of AGB stars to study their influence on the abundances of selected $s$-process elements and isotopic ratios. There are no changes in the final results using the different rates for the $^{13}$C($\alpha$, $n$)$^{16}$O reaction when the $^{13}$C burns completely in radiative conditions. When the $^{13}$C burns in convective conditions, as in stars of initial mass lower than $\sim$2 $M_\sun$ and in post-AGB stars, some changes are to be expected, e.g., of up to 25% for Pb in our models. These variations will have to be carefully analyzed when more accurate stellar mixing models and more precise observational constraints are available.

Type Ia supernovae (SNe Ia) play an important role in astrophysics and are crucial for the studies of stellar evolution, galaxy evolution and cosmology. They are generally thought to be thermonuclear explosions of accreting carbon-oxygen white dwarfs (CO WDs) in close binaries, however, the nature of the mass donor star is still unclear. In this article, we review various progenitor models proposed in the past years and summarize many observational results that can be used to put constraints on the nature of their progenitors. We also discuss the origin of SN Ia diversity and the impacts of SN Ia progenitors on some fields. The currently favourable progenitor model is the single-degenerate (SD) model, in which the WD accretes material from a non-degenerate companion star. This model may explain the similarities of most SNe Ia. It has long been argued that the double-degenerate (DD) model, which involves the merger of two CO WDs, may lead to an accretion-induced collapse rather than a thermonuclear explosion. However, recent observations of a few SNe Ia seem to support the DD model, and this model can produce normal SN Ia explosion under certain conditions. Additionally, the sub-luminous SNe Ia may be explained by the sub-Chandrasekhar mass model. At present, it seems likely that more than one progenitor model, including some variants of the SD and DD models, may be required to explain the observed diversity of SNe Ia.

We study the interaction between dark sectors by considering the momentum transfer caused by the dark matter scattering elastically within the dark energy fluid. Describing the dark scattering analogy to the Thomson scattering which couples baryons and photons, we examine the impact of the dark scattering in CMB observations. Performing global fitting with the latest observational data, we find that for a dark energy equation of state $w<-1$, the CMB gives tight constraints on dark matter-dark energy elastic scattering. Assuming a dark matter particle of proton mass, we derive an elastic scattering cross section of $\sigma_D < 3.295 \times 10^{-10} \sigma_T$ where $\sigma_T$ is the cross section of Thomson scattering. For $w>-1$, however, the constraints are poor. For $w=-1$, $\sigma_D$ can formally take any value.

We review the conformal equivalence in describing the background expansion of the universe by $f(R)$ gravity both in the Jordan frame and the Einstein frame. In the Jordan frame, we present the general analytic expression for $f(R)$ models that have the same expansion history as the $\Lambda$CDM model. This analytic form can provide further insights on how cosmology can be used to test the $f(R)$ gravity at the largest scales. Moreover we present a systematic and self-consistent way to construct the viable $f(R)$ model in Jordan frame using the mass dilation rate function from the Einstein frame through the conformal transformation. In addition, we extend our study to the linear perturbation theories and we further exhibit the equivalence of the $f(R)$ gravity presented in the Jordan frame and Einstein frame in the perturbed space-time. We argue that this equivalence has solid physics root.

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $\alpha_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We examine the conformal equivalence between the $f(R)$ gravity and the interacting dark sector model. We review the well-known result that the conformal transformation physically corresponds to the mass dilation which marks the strength of interaction between dark sectors. Instead of modeling f(R) gravity in the Jordan frame, we construct the $f(R)$ gravity in terms of mass dilation function in the Einstein frame. We find that the condition to keep $f(R)$ gravity consistent with CMB observations ensures the energy flow from dark energy to dark matter in the corresponding interacting model, which meets the requirement to alleviate the coincidence problem in the Einstein framework.

Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (WD) accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. Based on a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (\sim100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

We present a comprehensive study of the massless scalar field perturbation in the Reissner-Nordstrom anti-de Sitter (RNAdS) spacetime and compute its quasinormal modes (QNM). For the lowest lying mode, we confirm and extend the dependence of the QNM frequencies on the black hole charge got in previous works. In near extreme limit, under the scalar perturbation we find that the imaginary part of the frequency tends to zero, which is consistent with the previous conjecture based on electromagnetic and axial gravitational perturbations. For the extreme value of the charge, the asymptotic field decay is dominated by a power-law tail, which shows that the extreme black hole can still be stable to scalar perturbations. We also study the higher overtones for the RNAdS black hole and find large variations of QNM frequencies with the overtone number and black hole charge. The nontrivial dependence of frequencies on the angular index is also discussed.

We study the interacting dark energy model with time varying dark energy equation of state. We examine the stability in the perturbation formalism and the degeneracy among the coupling between dark sectors, the time-dependent dark energy equation of state and dark matter abundance in the cosmic microwave background radiation. Further we discuss the possible ways to break such degeneracy by doing global fitting using the latest observational data and we get a tight constraint on the interaction between dark sectors.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

Cosmological analysis based on currently available observations are unable to rule out a sizeable coupling between dark energy and dark matter. However, the signature of the coupling is not easy to grasp, since the coupling is degenerate with other cosmological parameters, such as the dark energy equation of state and the dark matter abundance. We discuss the possible ways to break such degeneracy. Based on the perturbation formalism, we carry out the global fitting by using latest observational data and get a tight constraint on the interaction between dark sectors. We find that the appropriate interaction can alleviate the coincidence problem.

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

Recent discovery of several overluminous type Ia supernovae (SNe Ia) indicates that the explosive masses of white dwarfs may significantly exceed the canonical Chandrasekhar mass limit. Rapid differential rotation may support these massive white dwarfs. Based on the single-degenerate scenario, and assuming that the white dwarfs would differentially rotate when the accretion rate $\dot{M}>3\times 10^{-7}M_{\odot}\rm yr^{-1}$, employing Eggleton's stellar evolution code we have performed the numerical calculations for $\sim$ 1000 binary systems consisting of a He star and a CO white dwarf (WD). We present the initial parameters in the orbital period - helium star mass plane (for WD masses of $1.0 M_{\odot}$ and $1.2 M_{\odot}$, respectively), which lead to super-Chandrasekhar mass SNe Ia. Our results indicate that, for an initial massive WD of $1.2 M_{\odot}$, a large number of SNe Ia may result from super-Chandrasekhar mass WDs, and the highest mass of the WD at the moment of SNe Ia explosion is 1.81 $M_\odot$, but very massive ($>1.85M_{\odot}$) WDs cannot be formed. However, when the initial mass of WDs is $1.0 M_{\odot}$, the explosive masses of SNe Ia are nearly uniform, which is consistent with the rareness of super-Chandrasekhar mass SNe Ia in observations.

We study the effects of electron screening on nuclear reaction rates occurring during the Big Bang nucleosynthesis epoch. The sensitivity of the predicted elemental abundances on electron screening is studied in details. It is shown that electron screening does not produce noticeable results in the abundances unless the traditional Debye-H\"uckel model for the treatment of electron screening in stellar environments is enhanced by several orders of magnitude. The present work rules out electron screening as a relevant ingredient to Big Bang nucleosynthesis and ruling out exotic possibilities for the treatment of screening, beyond the mean-field theoretical approach.

Based on perturbation theory, we study the dynamics of how dark matter and dark energy in the collapsing system approach dynamical equilibrium while interacting. We find that the interaction between dark sectors cannot ensure the dark energy to fully cluster along with dark, leading to the energy non-conservation problem in the collapsing system We examine the cluster number counts dependence on the interaction between dark sectors. Furthermore, we analyze how dark energy inhomogeneities affect cluster abundances. It is shown that cluster number counts can provide specific signature of dark sectors interaction and dark energy inhomogeneities.

We examine the effectiveness of the weak gravity conjecture in constraining the dark energy by comparing with observations. For general dark energy models with plausible phenomenological interactions between dark sectors, we find that although the weak gravity conjecture can constrain the dark energy, the constraint is looser than that from the observations.

The nature of type Ia supernovae (SNe Ia) is still unclear. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of SNe Ia, in which a carbon-oxygen white dwarf accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. According to a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (~100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

The universe, with large-scale homogeneity, is locally inhomogeneous, clustering into stars, galaxies and larger structures. Such property is described by the smoothness parameter $\alpha$ which is defined as the proportion of matter in the form of intergalactic medium. If we take consideration of the inhomogeneities in small scale, there should be modifications of the cosmological distances compared to a homogenous model. Dyer and Roeder developed a second-order ordinary differential equation (D-R equation) that describes the angular diameter distance-redshift relation for inhomogeneous cosmological models. Furthermore, we may obtain the D-R equation for observational $H(z)$ data (OHD). The density-parameter $\Omega_{\rm M}$, the state of dark energy $\omega$, and the smoothness-parameter $\alpha$ are constrained by a set of OHD in a spatially flat $\Lambda$CDM universe as well as a spatially flat XCDM universe. By using of $\chi^2$ minimization method we get $\alpha=0.81^{+0.19}_{-0.20}$ and $\Omega_{\rm M}=0.32^{+0.12}_{-0.06}$ at $1\sigma$ confidence level. If we assume a Gaussian prior of $\Omega_{\rm M}=0.26\pm0.1$, we get $\alpha=0.93^{+0.07}_{-0.19}$ and $\Omega_{\rm M}=0.31^{+0.06}_{-0.05}$. For XCDM model, $\alpha$ is constrained to $\alpha\geq0.80$ but $\omega$ is weakly constrained around -1, where $\omega$ describes the equation of the state of the dark energy ($p_{\rm X}=\omega\rho_{\rm X}$). We conclude that OHD constrains the smoothness parameter more effectively than the data of SNe Ia and compact radio sources.

Recent observations show that our universe is accelerating by dark energy, so it is important to investigate the thermodynamical properties of it. The Undulant Universe is a model with equation of state $\omega(a)=-\cos(b\ln a)$ for dark energy, where we show that there neither the event horizon nor the particle horizon exists. However, as a boundary of keeping thermodynamical properties, the apparent horizon is a good holographic screen. The Universe has a thermal equilibrium inside the apparent horizon, so the Unified First Law and the Generalized Second Law of thermodynamics are satisfied. As a thermodynamical whole, the evolution of the Undulant Universe behaves very well in the current phase. However, when considering the unification theory, the failure of conversation law at the epoch of the matter dominated or near singularity need some more consideration for the form of the Undulant Universe.

The paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

HD 49798 is a hydrogen depleted subdwarf O6 star and has an X-ray pulsating companion (RX J0648.0-4418). The X-ray pulsating companion is a massive white dwarf. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we find that the hot subdwarf HD 49798 and its X-ray pulsating companion could produce a type Ia supernova (SN Ia) in future evolution. This implies that the binary system is a likely candidate of SN Ia progenitors. We also discussed the possibilities of some other WD + He star systems (e.g. V445 Pup and KPD 1930+2752) for producing SNe Ia.

We examine in some detail the influence of the systematics in different data sets including type Ia supernova sample, baryon acoustic oscillation data and the cosmic microwave background information on the fitting results of the Chevallier-Polarski-Linder parametrization. We find that the systematics in the data sets does influence the fitting results and leads to different evolutional behavior of dark energy. To check the versatility of Chevallier-Polarski-Linder parametrization, we also perform the analysis on the Wetterich parametrization of dark energy. The results show that both the parametrization of dark energy and the systematics in data sets influence the evolutional behavior of dark energy.

Type Ia supernovae (SNe Ia) play an important role in the study of cosmic evolution, especially in cosmology. There are several progenitor models for SNe Ia proposed in the past years. In this paper, by considering the effect of accretion disk instability on the evolution of white dwarf (WD) binaries, we performed detailed binary evolution calculations for the WD + red-giant (RG) channel of SNe Ia, in which a carbon-oxygen WD accretes material from a RG star to increase its mass to the Chandrasekhar mass limit. According to these calculations, we mapped out the initial and final parameters for SNe Ia in the orbital period--secondary mass plane for various WD masses for this channel. We discussed the influence of the variation of the duty cycle value on the regions for producing SNe Ia. Similar to previous studies, this work also indicates that the long-period dwarf novae offer a possible ways for producing SNe Ia. Meanwhile, we find that the surviving companion stars from this channel have a low mass after SN explosion, which may provide a means for the formation of the population of single low-mass WDs (<0.45Msun).

Context: The nature of type Ia supernovae (SNe Ia) is still unclear. Metallicities may have an important effect on their properties. Aims: In this paper, we study the He star donor channel towards SNe Ia comprehensively and systematically at various metallicities. Methods: Employing Eggleton's stellar evolution code with the optically thick wind assumption, we calculated about 10 000 WD + He star systems and obtained SN Ia production regions of the He star donor channel with metallicities Z=0.03, 0.02, 0.004 and 0.0001. According to a detailed binary population synthesis approach, we also obtained SN Ia birthrates at various metallicities. Results: Our study shows that both the initial mass of the He donor star and the initial orbital period for SNe Ia increase with metallicity, while the minimum initial mass of the carbon--oxygen white dwarfs producing SNe Ia decreases with metallicity. For a constant star-formation galaxy, SN Ia birthrates increase with metallicity. If a single starburst is assumed, SNe Ia occur systemically earlier and the peak value of the birthrate is larger for a high Z, and the He star donor channel with different metallicities can produce the young SNe Ia with delay times $\sim$45-220Myr.

Recent investigations of the WD + MS channel of Type Ia supernovae (SNe Ia) imply that this channel may be the main contribution to the old population (>1Gyr) of SNe Ia. In the WD + MS channel, the WD could accrete material from a main-sequence or a slightly evolved star until it reaches the Chandrasekhar mass limit. The companions in this channel would survive after SN explosion and show distinguishing properties. In this Letter, based on SN Ia production regions of the WD + MS channel and three formation channels of WD + MS systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. The properties can be verified by future observations. We find that the surviving companions of the old SNe Ia have a low mass, which provides a possible way to explain the formation of the population of single low-mass WDs (<0.45Msun).

Type Ia supernovae (SNe Ia) play an important role in diverse areas of astrophysics, from the chemical evolution of galaxies to observational cosmology. However, the nature of the progenitors of SNe Ia is still unclear. In this paper, according to a detailed binary population synthesis study, we obtained SN Ia birthrates and delay times from different progenitor models, and compared them with observations. We find that the Galactic SN Ia birthrate from the double-degenerate (DD) model is close to those inferred from observations, while the birthrate from the single-degenerate (SD) model accounts for only about 1/2-2/3 of the observations. If a single starburst is assumed, the distribution of the delay times of SNe Ia from the SD model is a weak bimodality, where the WD + He channel contributes to the SNe Ia with delay times shorter than 100Myr, and the WD + MS and WD + RG channels to those with age longer than 1Gyr.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Cosmic rays are hampered by the Moon and a deficit in its direction is expected (the so-called Moon shadow). The Moon shadow is an important method to determine the performance of an air shower array. In fact, the westward displacement of the shadow centre, due to the propagation of cosmic rays in the geomagnetic field, allows to calibrate the energy scale of the primary particles observed by the detector. In addition, the shape of the shadow allows a measurement of the angular resolution and the position of the deficit at high energy allows the evaluation of the pointing accuracy of the detector. In this paper we present the observation of the galactic cosmic rays Moon shadowing effect performed by the ARGO-YBJ experiment in the multi-TeV energy region. The measured angular resolution as a function of the shower size is compared with the expectations from a MC simulation.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

Employing Eggleton's stellar evolution code with an optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon-oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is $\sim$$0.3\times 10^{-3} {\rm yr}^{-1}$, and that this channel can produce SNe Ia with short delay times ($\sim45-$140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity ($>$400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

{Context} Recent investigations of the white dwarf (WD) + He star channel of Type Ia supernovae (SNe Ia) imply that this channel can produce SNe Ia with short delay times. The companion stars in this channel would survive and be potentially identifiable. {Aims} In this Letter, we study the properties of the companion stars of this channel at the moment of SN explosion, which can be verified by future observations. {Methods} According to SN Ia production regions of the WD + He star channel and three formation channels of WD + He star systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. {Results} We obtained the distributions of many properties of the companion stars of this channel at the moment of SN explosion. We find that the surviving companion stars have a high spatial velocity (>400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

The nature of the progenitors of Type Ia supernovae (SNe Ia) is still unclear. In this paper, by considering the effect of the instability of accretion disk on the evolution of white dwarf (WD) binaries, we performed binary evolution calculations for about 2400 close WD binaries, in which a carbon--oxygen WD accretes material from a main-sequence star or a slightly evolved subgiant star (WD + MS channel), or a red-giant star (WD + RG channel) to increase its mass to the Chandrasekhar (Ch) mass limit. According to these calculations, we mapped out the initial parameters for SNe Ia in the orbital period--secondary mass ($\log P^{\rm i}-M^{\rm i}_2$) plane for various WD masses for these two channels, respectively. We confirm that WDs in the WD + MS channel with a mass as low as $0.61 M_\odot$ can accrete efficiently and reach the Ch limit, while the lowest WD mass for the WD + RG channel is $1.0 \rm M_\odot$. We have implemented these results in a binary population synthesis study to obtain the SN Ia birthrates and the evolution of SN Ia birthrates with time for both a constant star formation rate and a single starburst. We find that the Galactic SN Ia birthrate from the WD + MS channel is $\sim$$1.8\times 10^{-3} {\rm yr}^{-1}$ according to our standard model, which is higher than previous results. However, similar to previous studies, the birthrate from the WD + RG channel is still low ($\sim$$3\times 10^{-5} {\rm yr}^{-1}$). We also find that about one third of SNe Ia from the WD + MS channel and all SNe Ia from the WD + RG channel can contribute to the old populations ($\ga$1 Gyr) of SN Ia progenitors.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Dark energy interacting with dark matter is a promising model to solve the cosmic coincidence problem. We study the signature of such interaction on large scale cosmic microwave background (CMB) temperature anisotropies. Based on the detail analysis in perturbation equations of dark energy and dark matter when they are in interaction, we find that the large scale CMB, especially the late Integrated Sachs Wolfe effect, is a useful tool to measure the coupling between dark sectors. We also discuss the possibility to detect the coupling by cross-correlating CMB maps with tracers of the large scale structure. We finally perform the global fitting to constrain the coupling by using the CMB power spectrum data together with other observational data. We find that in the $1\sigma$ range, the constrained coupling between dark sectors can solve the coincidence problem.

We present a general formalism to study the growth of dark matter perturbations when dark energy perturbations and interactions between dark sectors are present. We show that dynamical stability of the growth of structure depends on the type of coupling between dark sectors. By taking the appropriate coupling to ensure the stable growth of structure, we observe that the effect of the dark sectors' interaction overwhelms that of dark energy perturbation on the growth function of dark matter perturbation. Due to the influence of the interaction, the growth index can differ from the value without interaction by an amount within the observational sensibility, which provides a possibility to disclose the interaction between dark sectors through future observations on the growth of large structure.

The single-degenerate model is currently a favourable progenitor model for Type Ia supernovae (SNe Ia). Recent investigations on the WD + He star channel of the single-degenerate model imply that this channel is noteworthy for producing SNe Ia. In this paper we studied SN Ia birthrates and delay times of this channel via a detailed binary population synthesis approach. We found that the Galactic SN Ia birthrate from the WD + He star channel is $\sim 0.3\times 10^{-3} {\rm yr}^{-1}$ according to our standard model, and that this channel can explain SNe Ia with short delay times ($\sim4.5\times10^7 - 1.4\times10^8$ yr). Meanwhile, these WD + He star systems may be related to the young supersoft X-ray sources prior to SN Ia explosions.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.

We examine the effect of the thermal vacuum on the power spectrum of inflation by using the thermal field dynamics. We find that the thermal effect influences the CMB anisotropy at large length scale. After removing the divergence by using the holographic cutoff, we observe that the thermal vacuum explains well the observational CMB result at low multipoles. This shows that the temperature dependent factor should be considered in the study of power spectrum in inflation, especially at large length scale.

We investigate the influence of an interaction between dark energy and dark matter upon the dynamics of galaxy clusters. We obtain the general Layser-Irvine equation in the presence of interactions, and find how, in that case, the virial theorem stands corrected. Using optical, X-ray and weak lensing data from 33 relaxed galaxy clusters, we put constraints on the strength of the coupling in the dark sector. We find that this coupling is small but positive, indicating that dark energy can be decaying into dark matter.

We study perturbations of black holes absorbing dark energy. Due to the accretion of dark energy, the black hole mass changes. We observe distinct perturbation behaviors for absorption of different forms of dark energy into the black holes. This provides the possibility of extracting information whether dark energy lies above or below the cosmological constant boundary $w=-1$. In particular, we find in the late time tail analysis that, differently from the other dark energy models, the accretion of phantom energy exhibits a growing mode in the perturbation tail. The instability behavior found in this work is consistent with the Big Rip scenario, in which all of the bound objects are torn apart with the presence of the phantom dark energy.

Type Ia supernovae (SNe Ia) play an important role in astrophysics, especially in the study of cosmic evolution. There are several progenitor models for SNe Ia proposed in the past years. In this paper, we have carried out a detailed study of the He star donor channel, in which a carbon-oxygen white dwarf (CO WD) accretes material from a He main sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. Employing Eggleton's stellar evolution code with an optically thick wind assumption, and adopting the prescription of Kato & Hachisu (2004) for the mass accumulation efficiency of the He-shell flashes onto the WDs, we performed binary evolution calculations for about 2600 close WD binary systems. According to these calculations, we mapped out the initial parameters for SNe Ia in the orbital period--secondary mass ($\log P^{\rm i}-M^{\rm i}_2$) plane for various WD masses from this channel. The study shows that the He star donor channel is noteworthy for producing SNe Ia (i.e. $\sim 1.2\times10^{-3} {\rm yr}^{-1}$ in the Galaxy), and that the progenitors from this channel may appear as supersoft X-ray sources. Importantly, this channel can explain SNe Ia with short delay times ($\la 10^{8}$ yr), which is consistent with recent observational implications of young populations of SN Ia progenitors.

We use observations of Sunyaev-Zel'dovich effect and X-ray surface brightness to reconstruct the radial profiles of gas temperature and density under the assumption of a spherically symmetric distribution of the gas. The method of reconstruction, first raised by Silk & White, depends directly on the observations of the Sunyaev-Zel'dovich effect and the X-ray surface brightness, without involving additional assumptions such as the equation of state of the gas or the conditions of hydrostatic equilibrium. We applied this method to the cluster RX J1347.5-1145, which has both the Sunyaev-Zel'dovich effect and X-ray observations with relative high precision. It is shown that it will be an effective method to obtain the gas distribution in galaxy clusters. Statistical errors of the derived temperature and density profiles of gas were estimated according to the observational uncertainties.

In this paper, we use a set of observational $H(z)$ data (OHD) to constrain the $\Lambda$CDM cosmology. This data set can be derived from the differential ages of the passively evolving galaxies. Meanwhile, the $\mathcal {A}$-parameter, which describes the Baryonic Acoustic Oscillation (BAO) peak, and the newly measured value of the Cosmic Microwave Background (CMB) shift parameter $\mathcal {R}$ are used to present combinational constraints on the same cosmology. The combinational constraints favor an accelerating flat universe while the flat $\Lambda$CDM cosmology is also analyzed in the same way. We obtain a result compatible with that by many other independent cosmological observations. We find that the observational $H(z)$ data set is a complementarity to other cosmological probes.

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations; however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density.

We examine different phenomenological interaction models for Dark Energy and Dark Matter by performing statistical joint analysis with observational data arising from the 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and age estimates of 35 galaxies. Including the time-dependent observable, we add sensitivity of measurement and give complementary results for the fitting. The compatibility among three different data sets seem to imply that the coupling between dark energy and dark matter is a small positive value, which satisfies the requirement to solve the coincidence problem and the second law of thermodynamics, being compatible with previous estimates.

We examine the effects of possible phenomenological interactions between dark energy and dark matter on cosmological parameters and their efficiency in solving the coincidence problem. We work with two simple parameterizations of the dynamical dark energy equation of state and the constant dark energy equation of state. Using observational data coming from the new 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, we perform a statistical joint analysis of different forms of phenomenological interactions between dark energy and dark matter.

Assuming that dark energy may be treated as a fluid with a well defined temperature, close to equilibrium, we argue that if nowadays there is a transfer of energy between dark energy and dark matter, it must be such that the latter gains energy from the former and not the other way around.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

The amount of $^{56}$Ni produced in type Ia supernova (SN Ia) explosion is probably the most important physical parameter underlying the observed correlation of SN Ia luminosities with their light curves. Based on an empirical relation between the $^{56}$Ni mass and the light curve parameter $\triangle m_{15}$, we obtained rough estimates of the $^{56}$Ni mass for a large sample of nearby SNe Ia with the aim of exploring the diversity in SN Ia. We found that the derived $^{56}$Ni masses for different SNe Ia could vary by a factor of ten (e.g., $M_{\rm Ni}=0.1 - 1.3$ $M_{\odot}$), which cannot be explained in terms of the standard Chandrasekhar-mass model (with a $^{56}$Ni mass production of 0.4 -- 0.8 $M_{\odot}$). Different explosion and/or progenitor models are clearly required for various SNe Ia, in particular, for those extremely nickel-poor and nickel-rich producers. The nickel-rich (with $M_{\rm Ni}$ $>$ 0.8 $M_{\odot}$) SNe Ia are very luminous and may have massive progenitors exceeding the Chandrasekhar-mass limit since extra progenitor fuel is required to produce more $^{56}$Ni to power the light curve. This is also consistent with the finding that the intrinsically bright SNe Ia prefer to occur in stellar environments of young and massive stars. For example, 75% SNe Ia in spirals have $\Delta m_{15} < 1.2$ while this ratio is only 18% in E/S0 galaxies. On the other hand, the nickel-poor SNe Ia (with $M_{\rm Ni}$ $<$ 0.2 $M_{\odot}$) may invoke the sub-Chandrasekhar model, as most of them were found in early-type E/S0 galaxies dominated by the older and low-mass stellar populations. This indicates that SNe Ia in spiral and E/S0 galaxies have progenitors of different properties.

We present a thermodynamical description of the interaction between holographic dark energy and dark matter. If holographic dark energy and dark matter evolve separately, each of them remains in thermodynamic equilibrium. A small interaction between them may be viewed as a stable thermal fluctuation that brings a logarithmic correction to the equilibrium entropy. From this correction we obtain a physical expression for the interaction which is consistent with phenomenological descriptions and passes reasonably well the observational tests.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

We investigate the effect of the bulk contents in the DGP braneworld on the evolution of the universe. We find that although the pure DGP model cannot accommodate the transition of the effective equation of state of dark energy, once the bulk matter T^5_5 is considered, the modified model can realize the w_{eff} crossing -1. However this transition of the equation of state cannot be realized by just considering bulk-brane energy exchange or the GB effect while the bulk matter contribution is not included. T^5_5 plays the major role in the modified DGP model to have the w crossing -1 behavior. We show that our model can describe the super-acceleration of our universe with the equation of state of the effective dark energy and the Hubble parameter in agreement with observations.

Using the data coming from the new 182 Gold type Ia supernova samples, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and the H(z) data, we have performed a statistical joint analysis of the DGP brane-world model with a high curvature Gauss-Bonnet term in the bulk. Consistent parameters estimations show that the Gauss-Bonnet-Induced Gravity model is a viable candidate to explain the observed acceleration of our universe.

Using the data coming from the new 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, $H(z)$ and lookback time measurements, we have performed a statistical joint analysis of the interacting holographic dark energy model. Consistent parameter estimations show us that the interacting holographic dark energy model is a viable candidate to explain the observed acceleration of our universe.

We investigate the effect of the bulk content in the general Gauss-Bonnet braneworld on the evolution of the universe. We find that the Gauss-Bonnet term and the combination of the dark radiation and the matter content of the bulk play a crucial role in the universe evolution. We show that our model can describe the super-acceleration of our universe with the equation of state of the effective dark energy in agreement with observations.

We have investigated the thermodynamical properties of dark energy. Assuming that the dark energy temperature $T\sim a^{-n}$ and considering that the volume of the Universe enveloped by the apparent horizon relates to the temperature, we have derived the dark energy entropy. For dark energy with constant equation of state $w>-1$ and the generalized Chaplygin gas, the derived entropy can be positive and satisfy the entropy bound. The total entropy, including those of dark energy, the thermal radiation and the apparent horizon, satisfies the generalized second law of thermodynamics. However, for the phantom with constant equation of state, the positivity of entropy, the entropy bound, and the generalized second law cannot be satisfied simultaneously.

Several observational constraints are imposed on the interacting holographic model of Dark Energy and Dark Matter. First we use the age parameter today, as given by the WMAP results. Subsequently, we explained the reason why it is possible, as recently observed, for an old quasar to be observed in early stages of the universe. We discuss this question in terms of the evolution of the age parameter as well as in terms of the structure formation. Finally, we give a detailed discussion of the constraints implied by the observed CMB low $\ell$ suppression. As a result, the interacting holographic model has been proved to be robust and with reasonable bounds predicts a non vanishing interaction of Dark Energy and Dark Matter.

In this paper we first obtain Friedmann equations for the $(n-1)$-dimensional brane embedded in the $(n+1)$-dimensional bulk, with intrinsic curvature term of the brane included in the action (DGP model). Then, we show that one can always rewrite the Friedmann equations in the form of the first law of thermodynamics, $dE=TdS+WdV$, at apparent horizon on the brane, regardless of whether there is the intrinsic curvature term on the brane or a cosmological constant in the bulk. Using the first law, we extract the entropy expression of the apparent horizon on the brane. We also show that in the case without the intrinsic curvature term, the entropy expressions are the same by using the apparent horizon on the brane and by using the bulk geometry. When the intrinsic curvature appears, the entropy of apparent horizon on the brane has two parts, one part follows the $n$-dimensional area formula on the brane, and the other part is the same as the entropy in the case without the intrinsic curvature term. As an interesting result, in the warped DGP model, the entropy expression in the bulk and on the brane are not the same. This is reasonable, since in this model gravity on the brane has two parts, one induced from the $(n+1)$-dimensional bulk gravity and the other due to the intrinsic curvature term on the brane.

We have investigated the thermodynamical properties of the Universe with dark energy. Adopting the usual assumption in deriving the constant co-moving entropy density that the physical volume and the temperature are independent, we observed some strange thermodynamical behaviors. However, these strange behaviors disappeared if we consider the realistic situation that the physical volume and the temperature of the Universe are related. Based on the well known correspondence between the Friedmann equation and the first law of thermodynamics of the apparent horizon, we argued that the apparent horizon is the physical horizon in dealing with thermodynamics problems. We have concentrated on the volume of the Universe within the apparent horizon and considered that the Universe is in thermal equilibrium with the Hawking temperature on the apparent horizon. For dark energy with $w\ge -1$, the holographic principle and the generalized second law are always respected.

We have studied the wave dynamics and the Hawking radiation for the scalar field as well as the brane-localized gravitational field in the background of the braneworld black hole with tidal charge containing information of the extra dimension. Comparing with the four-dimensional black holes, we have observed the signature of the tidal charge which presents the signals of the extra dimension both in the wave dynamics and the Hawking radiation.

We have studied the uncertainty on the determination of the dark energy equation of state due to a non-vanishing spatial curvature by considering some fundamental observables. We discussed the sensitivity of these observables to the value and redshift history of the equation of state and the spatial curvature and investigated whether these different observables are complementary and can help to reduce the cosmic confusion.

We investigate the question of the suppression of the CMB power spectrum for the lowest multipoles in closed Universes. The intrinsic reason for a lowest cutoff in closed Universes, connected with the discrete spectrum of the wavelength, is shown not to be enough to explain observations. We thus extend the holographic cosmic duality to closed universes by relating the dark energy equation of state and the power spectrum, showing a suppression behavior which describes the low l features extremely well. We also explore the possibility to disclose the nature of the dark energy from the observed small l CMB spectrum by employing the holographic idea.

We consider a brane cosmological model with energy exchange between brane and bulk. Parameterizing the energy exchange term by the scale factor and Hubble parameter, we are able to exactly solve the modified Friedmann equation on the brane. In this model, the equation of state for the effective dark energy has a transition behavior changing from $w_{de}^{eff}>-1$ to $w_{de}^{eff}<-1$, while the equation of state for the dark energy on the brane has $w>-1$. Fitting data from type Ia supernova, Sloan Digital Sky Survey and Wilkinson Microwave Anisotropy Probe, our universe is predicted now in the state of super-acceleration with $w_{de0}^{eff}=-1.21$.

Using the evolution history of the universe, one can make constraint on the parameter space of dynamic dark energy models. We discuss two different parameterized dark energy models. Our results further restrict the combined constraints obtained from supernova and WMAP observations. From the allowed parameter space, it is found that our universe will experience an eternal acceleration. We also estimate the bound on the physically relevant regions both in the re-inflationary and inflationary phases.

Recently the 1/R gravity has been proposed in order to explain the accelerated expansion of the universe. However, it was argued that the 1/R gravity conflicts with solar system tests. While this statement is true if one views the 1/R gravity as an effective theory, we find that this difficulty might be avoided if one treats the 1/R term as a correction to the scalar curvature term in the high curvature limit $R \gg \mu^2$.

We have reexamined the holographic dark energy model by considering the spatial curvature. We have refined the model parameter and observed that the holographic dark energy model does not behave as phantom model. Comparing the holographic dark energy model to the supernova observation alone, we found that the closed universe is favored. Combining with the Wilkinson Microwave Anisotropy Probe (WMAP) data, we obtained the reasonable value of the spatial curvature of our universe.

A model of holographic dark energy with an interaction with matter fields has been investigated. Choosing the future event horizon as an IR cutoff, we have shown that the ratio of energy densities can vary with time. With the interaction between the two different constituents of the universe, we observed the evolution of the universe, from early deceleration to late time acceleration. In addition, we have found that such an interacting dark energy model can accommodate a transition of the dark energy from a normal state where $w_D>-1$ to $w_D<-1$ phantom regimes. Implications of interacting dark energy model for the observation of dark energy transition has been discussed.

We present the first study on the mass functions of Jenkins et al (J01) and an estimate of their corresponding largest virialized dark halos in the Universe for a variety of dark-energy cosmological models with a running spectral index. Compared with the PL-CDM model, the RSI-CDM model can raise the mass abundance of dark halos for small mass halos at lower redshifts, but it is not apparent on scales of massive mass halos. Particularly, this discrepancy increases largely with the decrease of redshift, and the RSI-CDM model can suppress the mass abundance on any scale of halo masses at higher redshift. As for the largest mass of virialized halos, the spatially flat $\Lambda$CDM models give more massive mass of virialized objects than other models for both of PL-CDM and RSI-CDM power spectral indexs, and the RSI-CDM model can enhance the mass of largest virialized halos for all of models considered in this paper. So we probably distinguish the PL-CDM and RSI-CDM models by the largest virialized halos in the future survey of cluster of galaxies.

Employing the Quark Mass Denisity- and temperature- dependent model and the Hartle's method, We have studied the slowly rotating strange star with uniform angular velocity. The mass-radius relation, the moment of inertia and the frame dragging for different frequencies are given. We found that we cannot use the strange star to solve the challenges of Stella and Vietri for the horizontal branch oscillations and the moment of inertia $I_{45}/(M/M_ s)>2.3$. Furthermore, we extended the Hartle's method to study the differential rotating strange star and found that the differential rotation is an effective way to get massive strange star.

Using the holographic principle we constrained the Friedmann equation, modified by brane-cosmology inspired terms which accommodate dark energy contributions in the context of extra dimensions.

Using the recently obtained holographic cosmic duality, we reached a reasonable quantitative agreement between predictions of the Cosmic Microwave Background Radiation at small l and the WMAP observations, showing the power of the holographic idea. We also got constraints on the dark energy and its behaviour as a function of the redshift upon relating it to the small l CMB spectrum. For a redshift independent dark energy, our constraint is consistent with the supernova results, which again shows the correctness of the cosmic duality prescription. We have also extended our study to the redshift dependence of the dark energy.

We have studied the evolution of the massless scalar field propagating in time-dependent charged Vaidya black hole background. A generalized tortoise coordinate transformation were used to study the evolution of the massless scalar field. It is shown that, for the slowest damped quasinormal modes, the approximate formulae in stationary Reissner-Nordstr\"{o}m black hole turn out to be a reasonable prescription, showing that results from quasinormal mode analysis are rather robust.

We study the possibility to probe the spatial geometry of the Universe by supernova measurement of the cubic correction to the luminosity distance. We illustrate with an accelerating universe model with infinite-volume extra dimensions, for which the 1$\sigma$ level supernova results indicate that the Universe is closed.

Based solely on the arguments relating Friedmann equation and the Cardy formula we derive a bound for the number of e-folds during inflation for a standard Friedmann-Robertson-Walker as well as non-standard four dimensional cosmology induced by a Randall-Sundrum-type model.

In light of WMAP results, we examine the observational constraint on the P-term inflation. With the tunable parameter $f$, P-term inflation contains richer physics than D-term and F-term inflationary models. We find the logarithmic derivative spectral index with $n>1$ on large scales and $n<1$ on small scales in agreement to observation. We obtained a reasonable range for the choice of the gauge coupling constant $g$ in order to meet the requirements of WMAP observation and the expected number of the e-foldings. Although tuning $f$ and $g$ we can have larger values for the logarithmic derivative of the spectral index, it is not possible to satisfy all observational requirements for both, the spectral index and its logarithmic derivative at the same time.

Recent observations of Type Ia supernovae provide evidence for the acceleration of our universe, which leads to the possibility that the universe is entering an inflationary epoch. We simulate it under a ``big bounce'' model, which contains a time variable cosmological ``constant'' that is derived from a higher dimension and manifests itself in 4D spacetime as dark energy. By properly choosing the two arbitrary functions contained in the model, we obtain a simple exact solution in which the evolution of the universe is divided into several stages. Before the big bounce, the universe contracts from a $\Lambda $-dominated vacuum, and after the bounce, the universe expands. In the early time after the bounce, the expansion of the universe is decelerating. In the late time after the bounce, dark energy (i.e., the variable cosmological ``constant'') overtakes dark matter and baryons, and the expansion enters an accelerating stage. When time tends to infinity, the contribution of dark energy tends to two third of the total energy density of the universe, qualitatively in agreement with observations.

The big bounce singularity of a simple 5D cosmological model is studied. Contrary to the standard big bang space-time singularity, this big bounce singularity is found to be an event horizon at which the scale factor and the mass density of the universe are finite, while the pressure undergoes a sudden transition from negative infinity to positive infinity. By using coordinate transformation it is also shown that before the bounce the universe contracts deflationary, and the universe has been existed, according to the proper-time, for an infinitely long time.

We study the massless scalar wave propagation in the time-dependent Schwarzschild black hole background. We find that the Kruskal coordinate is an appropriate framework to investigate the time-dependent spacetime. A time-dependent scattering potential is derived by considering dynamical black hole with parameters changing with time. It is shown that in the quasinormal ringing both the decay time-scale and oscillation are modified in the time-dependent background.

The Sunyaev Zel'dovich effect (SZ effect) is a complete probe of ionized baryons, the majority of which are likely hiding in the intergalactic medium. We ran a $512^3$ $\Lambda$CDM simulation using a moving mesh hydro code to compute the statistics of the thermal and kinetic SZ effect such as the power spectra and measures of non-Gaussianity. The thermal SZ power spectrum has a very broad peak at multipole $l\sim 2000-10^4$ with temperature fluctuations $\Delta T \sim 15\mu$K. The power spectrum is consistent with available observations and suggests a high $\sigma_8\simeq 1.0$ and a possible role of non-gravitational heating. The non-Gaussianity is significant and increases the cosmic variance of the power spectrum by a factor of $\sim 5$ for $l<6000$. We explore optimal driftscan survey strategies for the AMIBA CMB interferometer and their dependence on cosmology. For SZ power spectrum estimation, we find that the optimal sky coverage for a 1000 hours of integration time is several hundred square degrees. One achieves an accuracy better than 40% in the SZ measurement of power spectrum and an accuracy better than 20% in the cross correlation with Sloan galaxies for $2000<l<5000$. For cluster searches, the optimal scan rate is around 280 hours per square degree with a cluster detection rate 1 every 7 hours, allowing for a false positive rate of 20% and better than 30% accuracy in the cluster SZ distribution function measurement.

We calculated the entropy of a class of inhomogeneous dust universes. Allowing spherical symmetry, we proposed a holographic principle by reflecting all physical freedoms on the surface of the apparent horizon. In contrast to flat homogeneous counterparts, the principle may break down in some models, though these models are not quite realistic. We refined fractal parabolic solutions to have a reasonable entropy value for the present observable universe and found that the holographic principle always holds in the realistic cases.

We propose that cooling outflows from star-forming dwarf galaxies that are spatially correlated with but unbound to bright galaxies may account for the low column QSO Ly$\alpha$ forests, providing an explanation for both the heavy elements recently detected in them and the large sizes associated with them. Thus, we suggest that structure formations in the early universe on scales where matters collapse gravitationally and the gas forms stars are responsible for the Ly$\alpha$ forests (cf. the mini-halo model). Supernovae shock the interstellar gas to high temperature, and the resulting gas outflows from the dwarfs owing to their shallow potential wells. The gas cools radiatively as the radiative cooling time shortens because of adiabatic expansion. Subsequent thermal instability results in condensation of clouds. The clouds inherit the kinetic energy of the flow and coast to large distances, providing regions of the absorbing gas with sizes of up to hundreds of kpc. We calculate the Ly$\alpha$ absorption line profiles and the neutral column distribution, and show that they are consistent with observations. We suggest that the Ly$\alpha$ forests are caused by the faint blue galaxies found in the deep extragalactic surveys.

We study steady, radial gas outflows from galaxies in an effort to understand the way tenuous and hot gas is transported to large distances away from galaxies. In particular, we obtain solutions for outflow problems, and study the outflow topology, effects of the galaxy mass, the size of outflow regions, the efficiency of radiative cooling, and the fate of the cooled gas. Under general power-law forms for the cooling function and the gravitational field of galaxies, we show that the outflow solutions are determined by the two-parameter initial conditions. In an analogy with stellar wind or accretion problems, we demonstrate that there exists no transonic flow, but either subsonic or supersonic flows are obtainable. Solutions of the supersonic outflows are studied in detail as they are most likely to carry gas to large distances away from galaxies. We find that if gravity is weak, the outflow is characterized by the ratio of the radiative cooling time to the flow time, $t_c/t_f$. The importance of the galactic gravitational field is characterized by the fractional energy lost radiatively within the flow time in outflows with velocity equal to the circular velocity of the galaxy; if the fraction is small, gravity stops the outflow before the gas has a chance to cool radiatively. In the case the gas does cool radiatively, the cooled gas is most likely to form clouds via various instabilities. The clouds coast farther away from the galaxy because of the finite kinetic energy they inherit. We find that the hot gas in dwarf galaxies can either flow out as galactic winds, or cools radiatively to form clouds. In the latter case, the clouds escape the galaxies. In contrast, massive galaxies like our own tend to confine the gas. We present the surface brightness in various x-ray energy bands. We also estimate the mean

We present the calculations of steady, spherical gas outflows from galaxies in an effort to understand the effects of the galaxy mass on the flow properties such as the size of the outflow regions, the efficiency of radiative cooling, and the fate of the cooled gas. We show that there exist no transonic flows but either subsonic or supersonic flows are obtainable, in an analogy with the stellar wind or accretion problems. Solutions of the supersonic outflows are studied in detail as they are most likely to carry gas to large distances away from galaxies. In case the gas does cool radiatively, the cooled gas is most likely to form clouds via various instabilities. The clouds coast farther away from the galaxy because of the finite kinetic energy they inherit. Depending on the initial energies, the clouds can either leave the galaxy or fall back ballistically. Applications of the calculations to dwarf and normal galaxies are made, and we argue that normal galaxies like our own should have relatively small gaseous halos formed through this mechanism. We discuss the implications of cold clouds from dwarf galaxies on recent observations of the QSO absorption line systems. For massive galaxies, satellites provide alternative ways to form large gaseous halos. We show how dynamical friction can increase tidal stripping of the gas as the satellite spirals in toward the primary, and that the gas ejected from satellites may account for the QSO heavy element absorption lines.