By using recent data, we directly determine the dark matter (DM) induced $e^\pm$ spectrum at the source from experimental measurements at the earth, without reference to specific particle physics models. The DM induced gamma rays emitted via inverse Compton scattering are then obtained in a model independent way. However the results depend on the choice of the astrophysical $e^\pm$ background, which is not reliably known. Nevertheless, we calculate, as an illustration, the fluxes of gamma rays from the Fornax cluster in the decaying DM scenario with various astrophysical $e^\pm$ backgrounds. Without any assumptions on details of the DM model, the predictions turn out to be either in disagreement with or only marginally below the upper limits measured recently by the Fermi-LAT Collaboration. In addition, these DM induced ICS gamma rays in the GeV range are shown to be almost independent of choices of cosmic ray propagation model and of DM density profile, when a given astrophysical $e^\pm$ background is assumed. This provides a strong constraint on decaying DM scenario as the gamma rays may be produced in other processes besides inverse Compton scattering, such as the bremsstrahlung and neutral pion decays.

Recently, the Fermi-LAT collaboration reported upper limits on the GeV gamma-ray flux from nearby clusters of galaxies. Motivated by these limits, we study corresponding constraints on gamma-ray emissions from two specific decaying dark matter models, one via grand unification scale suppressed operators and the other via R-parity violating operators. Both can account for the PAMELA and Fermi-LAT excesses of e^\pm. For GUT decaying dark matter, the gamma-rays from the M49 and Fornax clusters, with energy in the range of 1 to 10 GeV, lead to the most stringent constraints to date. As a result, this dark matter is disfavored with conventional model of e^\pm background. In addition, it is likely that some tension exists between the Fermi-LAT e^\pm excess and the gamma ray constraints for any decaying dark matter model, provided conventional model of e^\pm background is adopted. Nevertheless, the GUT decaying dark matter can still solely account for the PAMELA positron fraction excess without violating the gamma-ray constraints. For the gravitino dark matter model with R-parity violation, cluster observations do not give tight constraints. This is because a different e^\pm background has been adopted which leads to relatively light dark matter mass around 200 GeV.

Motivated by recent observations from Pamela, Fermi and H.E.S.S., we consider dark matter decays in the framework of supersymmetric SU(5) grand unification theories. An SU(5) singlet S is assumed to be the main component of dark matters, which decays into visible particles through dimension six operators suppressed by the grand unification scale. Under certain conditions, S decays dominantly into a pair of sleptons with universal coupling for all generations. Subsequently, electrons and positrons are produced from cascade decays of these sleptons. These cascade decay chains smooth the electron/positron spectrum, which permit naturally a good fit to the Fermi LAT data. The observed positron fraction upturn by PAMELA can be reproduced simultaneously. We have also calculated diffuse gamma-ray spectra due to the electron/positron excesses and compared them with the preliminary Fermi LAT data from 0.1 GeV to 10 GeV in the region 0<l <360, 10<|b|<20. The photon spectrum of energy above 100 GeV, mainly from final state radiations, may be checked in the near future.

Quintessence field is a widely-studied candidate of dark energy. There is "tracker solution" in quintessence models, in which evolution of the field $\phi$ at present times is not sensitive to its initial conditions. When the energy density of dark energy is neglectable ($\Omega_\phi\ll1$), evolution of the tracker solution can be well analysed from "tracker equation". In this paper, we try to study evolution of the quintessence field from "full tracker equation", which is valid for all spans of $\Omega_\phi$. We get stable fixed points of $w_\phi$ and $\Omega_\phi$ (noted as $\hat w_\phi$ and $\hat\Omega_\phi$) from the "full tracker equation", i.e., $w_\phi$ and $\Omega_\phi$ will always approach $\hat w_\phi$ and $\hat\Omega_\phi$ respectively. Since $\hat w_\phi$ and $\hat\Omega_\phi$ are analytic functions of $\phi$, analytic relation of $\hat w_\phi\sim\hat\Omega_\phi$ can be obtained, which is a good approximation for the $w_\phi\sim\Omega_\phi$ relation and can be obtained for the most type of quintessence potentials. By using this approximation, we find that inequalities $\hat w_\phi<w_\phi$ and $\hat\Omega_\phi<\Omega_\phi$ are statisfied if the $w_\phi$ (or $\hat w_\phi$) is decreasing with time. In this way, the potential $U(\phi)$ can be constrained directly from observations, by no need of solving the equations of motion numerically.

We analyze the existing solar neutrino experiment data and show the allowed regions. The result from SNO's salt phase itself restricts quite a lot the allowed region's area. Reactor neutrinos play an important role in determining oscillation parameters. KamLAND gives decisive conclusion on the solution to the solar neutrino puzzle, in particular, the spectral distortion in the 766.3 Ty KamLAND data gives another new improvement in the constraint of solar MSW-LMA solutions. We confirm that at 99.73% C.L. the high-LMA solution is excluded.

Non-minimally coupled theories of special potentials are analyzed numerically. Such theories yield equations of state $\omega<-1$ and oscillations of the cosmological expansion, which are favored by the recent analysis of observations. Fitting these theories to the Gold SnIa dataset, we obtain results comparable with other models. A potential of the form $V(\phi)=V_0e^{a_1\phi^2}$ yields $\chi^2_{min}=170.127$. Similar results are obtained for potentials of the form $V(\phi)=V_0+a_1\phi^n$.

We report results on the broad scale anisotropy of cosmic ray arrival directions in the energy rage from $2 \times 10^{17} eV$ to $10^{20} eV$. The data was taken by the Fly's Eye detector in both monocular and stereo modes of operation. We look for dependence on galactic latitude or supergalactic latitude by fitting the data to a Wdowczyk and Wolfendale plane enhancement function and a N-S gradient functional form. We report a small but statistically significant galactic plane enhancement in the energy range between $2\times 10^{17}eV$ and $3.2 \times 10^{18} eV$. The probability that this anisotropy is due to fluctuations of an isotropic distribution is less than 0.06%. The most significant galactic plane enhancement factor $f_{E} = 0.104 \pm 0.036$ is in the energy range $0.4-1.0 \times 10^{18} eV$. No statistically significant evidence for a N-S gradient is found. There is no sign of significant deviation from isotropic background when the data is analyzed in terms of supergalactic latitude distributions