We present the observational results from the 43-GHz Very Long Baseline Array (VLBA) observations of 124 compact radio-loud active galactic nuclei (AGNs) that were conducted between 2014 November and 2016 May. The typical dimensions of the restoring beam in each image are about 0.5 mas $\times$ 0.2 mas. The highest resolution of 0.2 mas corresponds to a physical size of 0.02 pc for the lowest redshift source in the sample. The 43-GHz very long baseline interferometry (VLBI) images of 97 AGNs are presented for the first time. We study the source compactness on milli-arcsec (mas) and sub-mas scales, and suggest that 95 sources in our sample are suitable for future space VLBI observations. By analyzing our data supplemented with other VLBA AGN surveys from literature, we find that the core brightness temperature increases with increasing frequency below a break frequency ~ 7 GHz, and decreases between ~7--240~GHz but increases again above~240 GHz in the rest frame of the sources. This indicates that the synchrotron opacity changes from optically thick to thin. We also find a strong statistical correlation between radio and $\gamma$-ray flux densities. Our correlation is tighter than those in literature derived from lower-frequency VLBI data, suggesting that the $\gamma$-ray emission is produced more co-spatially with the 43-GHz VLBA core emission. This correlation can also be extrapolated to the un-beamed AGN population, implying that a universal $\gamma$-ray production mechanism might be at work for all types of AGNs.

We analyzed 199 GRBs detected by Fermi-LAT during the years 2008-2019. We found 67 photons at energies >=10 GeV, which come from 34 GRBs. Out of these 34 GRBs, Fermi-LAT detects significant (TS>=4) afterglow 0.1-200 GeV photons from 25 GRBs. We present time-integrated 0.1-200 GeV spectra of these 25 GRBs. The spectra of a significant fraction (9/25) of these GRBs revealed a harder component above an energy break of 0.3-2 GeV. While shock synchrotron may account for the photons at the lower energy end, high energy photons above the break is naturally explained by synchrotron self-Compton (SSC) emission. We perform broadband model fit to the X-ray-LAT emission of GRB 131231A. Comparing the afterglow spectra of these 25 GRBs with other Fermi-LAT detected GRBs, we found that the power-law index distribution is similar for the two populations. This may indicate that the additional high-energy component may also exist in Fermi-LAT GRBs in general.

We investigate the nature of the unidentified very-high-energy (VHE) gamma-ray object, HESS J1832-093, in a multi-wavelength context. Based on X-ray variability and spectral index ($\Gamma_X\sim\,1.5$), and its broad-band spectrum (which was remarkably similar to HESS J0632+057, a confirmed gamma-ray binary), HESS J1832-093 has been considered to be a strong gamma-ray binary candidate in previous works. In this work, we provide further evidence for this scenario. We obtained a spectrum of its IR counterpart using Gemini/Flamingo, finding absorption lines that are usually seen in massive stars, in particular O stars. We also obtained a rather steep ATCA spectrum ($\alpha=-1.18^{+1.04}_{-0.88}$) which prefers a gamma-ray binary over an AGN scenario. Based on spatial-spectral analysis and variability search, we found that 4FGL J1832.9-0913 is possible to be associated with SNR G22.7-0.2 rather than with HESS J1832-093 only.

The supernova remnant (SNR) HESS J1731-347 is a young SNR which displays a non-thermal X-ray and TeV shell structure. A molecular cloud at a distance of 3.2 kpc is spatially coincident with the western part of the SNR, and it is likely hit by the SNR. The X-ray emission from this part of the shell is much lower than from the rest of the SNR. Moreover, a compact GeV emission region coincident with the cloud has been detected with a soft spectrum. These observations seem to imply a shock-cloud collision scenario at this area, where the stalled shock can no longer accelerate super-TeV electrons or maintain strong magnetic turbulence downstream, while the GeV cosmic rays (CRs) are released through this stalled shock. To test this hypothesis, we have performed a detailed Fermi-LAT reanalysis of the HESS J1731-347 region with over 9 years of data. We find that the compact GeV emission region displays a spectral power-law index of -2.4, whereas the GeV emission from the rest of the SNR (excluding the cloud region) has an index of -1.8. A hadronic model involving a shock-cloud collision scenario is built to explain the -ray emission from this area. It consists of three CR sources: run-away super-TeV CRs that have escaped from the fast shock, leaked GeV CRs from the stalled shock, and the local CR sea. The X-ray and -ray emission of the SNR excluding the shock-cloud interaction region is explained in a one-zone leptonic model. Our shock-cloud collision model explains well the GeV-TeV observations from both cloud regions around HESS J1731-347, i.e. from the cloud in contact with the SNR and from the more distant cloud which is coincident with the nearby TeV source HESS J1729-345. We find however that the leaked GeV CRs from the shock-cloud collision do not necessarily dominate the GeV emission from the clouds, due to a comparable contribution from the local CR sea.

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.

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.

GRB 190114C is a famous Gamma-Ray Burst (GRB) due to its detection at sub-TeV energies by MAGIC, seen at redshift z = 0.42. This burst is one of the brightest GRB detected by \fermi. We present a detailed analysis of GRB 190114C prompt emission, using the two \fermi detectors: GBM and LAT. The \emph{LAT low energy events} (LLE) data is also considered. A joint GBM-LAT analysis reveals a sub-$GeV$ spectral cutoff. A similar high energy cutoff was reported in GRB 160509A and GRB 100724B earlier, and a handful of other sources. The cutoff can be explained by the intrinsic opacity due to pair production within the emitting region. Such morphology in these GRBs suggests that they belong to one specific class having a similar source of the radiation mechanism. GRB 190114C shows a transition from non-thermal to a quasi-thermal-like spectrum along with radiation due to external shock. From spectrum analysis and $Lorentz$ factor evolution from the trigger time to late emission, considering the fact that sub-TeV photons are detected in MAGIC, we are able to draw an emission mechanism picture, where the prompt emission spectrum is more consistent with spectrum via photospheric dissipation with presence of external shock emission simultaneously.

The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.

Galactic cosmic rays are the high-energy particles that stream into our solar system from distant corners of our Galaxy and some low energy particles are from the Sun which are associated with solar flares. The Earth atmosphere serves as an ideal detector for the high energy cosmic rays which interact with the air molecule nuclei causing propagation of extensive air showers. In recent years, there are growing interests in the applications of the cosmic ray measurements which range from the space/earth weather monitoring, homeland security based on the cosmic ray muon tomography, radiation effects on health via air travel, etc. A simulation program (based on the GEANT4 software package developed at CERN) has been developed at Georgia State University for studying the cosmic ray showers in atmosphere. The results of this simulation study will provide unprecedented knowledge of the geo-position-dependent cosmic ray shower profiles and significantly enhance the applicability of the cosmic ray applications. In the paper, we present the computational challenges and the opportunities for carrying out the cosmic ray shower simulations at the global scale using various computing resources including XSEDE.

GeV flares from PSR B1259-63/LS 2883 were seen starting around 30 days after the two periastron passages in 2010 and 2014. The flares are clearly delayed compared to the occurrence of the X-ray and TeV flux peaks during the post-periastron disk crossing. While several attempts have been put forward to explain this phenomenon, the origin of these GeV flares remains a puzzle. Here we present a detailed analysis of the observational data taken by the Fermi and Swift observatories over the 2017 September periastron passage. For the first time, we find short-lived but powerful GeV flares on time scales of down to three hours. The onset of the GeV flaring period in 2017 is also delayed compared to those seen in 2011 and 2014. Supplemented by a re-analysis of previous data, we compare the Fermi/LAT, Swift/XRT and Swift/UVOT light curves in 2017 with those taken over the 2010 and 2014 periastrons, and difference in UVOT light curves are noted.

The angular distribution of GRB jets is not yet clear. The observed luminosity of GRB 170817A is the lowest among all known short GRBs, which is best explained by the fact that our line of sight is outside of the jet opening angle, $\theta_{obs}$ $>$ $\theta_j$, where $\theta_{obs}$ is the angle between our line of sight and jet axis. Inferred by gravitational wave observations, as well as radio and X-ray afterglow modeling of GRB 170817A, it is likely that $\theta_{obs}$ $\sim$ 20$^{\circ}$ -- 28$^{\circ}$. In this work, we quantitatively consider two scenarios of angular energy distribution of GRB ejecta: top-hat jet, and structured jet with a power law index $s$. For top-hat jet model, we get a large $\theta_j$ (e.g., $\theta_{j}$$>$$10^{\circ}$), a rather high local (i.e., z$<$0.01) short GRB rate $\sim$8--15$\times10^3$ $Gpc^{-3}$$yr^{-1}$ (estimated 90$\sim$1850 $ Gpc^{-3}$$yr^{-1}$ in Fong et al. 2015) and an extremely high $E_{peak,0}$(on-axis,rest-frame)$>$7.5$\times10^4$keV($\sim$500keV for typical short GRB). For structured jet model, we use $\theta_{obs}$ to give limits on $s$ and $\theta_{j}$ for a typical on-axis luminosity short GRB (e.g., $10^{49}$ erg s$^{-1}$$\sim$$10^{51}$ erg s$^{-1}$), and a low on-axis luminosity case (e.g., $10^{49}$ erg s$^{-1}$) gives more reasonable values of $s$. A structured jet model is more feasible for GRB 170817A than a top-hat jet model, due to the rather high local short GRB rate and the extremely high on-axis $E_{peak,0}$ almost rules out the top-hat jet model. GRB 170817A is likely a low on-axis luminosity GRB ($10^{49}$ erg s$^{-1}$) with a structured jet.

Main sequence turn-off (MSTO) stars have advantages as indicators of Galactic evolution since their ages could be robustly estimated from atmospheric parameters. Hundreds of thousands of MSTO stars have been selected from the LAMOST Galactic sur- vey to study the evolution of the Galaxy, and it is vital to derive accurate stellar parameters. In this work, we select 150 MSTO star candidates from the MSTO stars sample of Xiang that have asteroseismic parameters and determine accurate stellar parameters for these stars combing the asteroseismic parameters deduced from the Kepler photometry and atmospheric parameters deduced from the LAMOST spectra.With this sample, we examine the age deter- mination as well as the contamination rate of the MSTO stars sample. A comparison of age between this work and Xiang shows a mean difference of 0.53 Gyr (7%) and a dispersion of 2.71 Gyr (28%). The results show that 79 of the candidates are MSTO stars, while the others are contaminations from either main sequence or sub-giant stars. The contamination rate for the oldest stars is much higher than that for the younger stars. The main cause for the high contamination rate is found to be the relatively large systematic bias in the LAMOST surface gravity estimates.

We present the high-energy emission properties of GRB 160509A, from its prompt mission to late afterglow phase. GRB 160509A contains two emission episodes: 0-40s and 280-420s after the burst onset (t0). The relatively high fluence of GRB 160509A allows us to establish an evolving spectrum above 100 MeV. During the first emission episode, the >100 MeV spectrum is soft with \Gamma=>3.0, which can be smoothly connected to keV energies with a Band function with a high-energy cutoff. The >100 MeV spectrum rapidly changes to a hard spectrum with \Gamma<=1.5 after t0+40s. The existence of very energetic photons, e.g., a 52 GeV that arrives t0+77 seconds, and a 29 GeV that arrives t0+70 ks, is hard to reconcile by the synchrotron emission from forward-shock electrons, but likely due to inverse Compton mechanism (e.g., synchrotron self-Compton emission). A soft spectrum (\Gamma~2) between 300s and 1000s after the burst onset is also found at a significance of about 2 standard deviations, which suggests a different emission mechanism at work for this short period of time. GRB 160509A represents the latest example where inverse Compton emission has to be taken into account in explaining the afterglow GeV emission, which had been suggested long before the launch of Fermi LAT.

We explore simple Higgs-portal models of dark matter (DM) with spin 1/2, 3/2, and 1, respectively, applying to them constraints from the LUX and PandaX-II direct detection experiments and from LHC measurements on the 125-GeV Higgs boson. With only one Higgs doublet, we find that the spin-1/2 DM having a purely scalar effective coupling to the doublet is viable only in a narrow range of mass near the Higgs pole, whereas the vector DM is still allowed if its mass is also close to the Higgs pole or exceeds 1.4 TeV, both in line with earlier analyses. Moreover, the spin-3/2 DM is in a roughly similar situation to the spin-1/2 DM, but has surviving parameter space which is even more restricted. We also consider the two-Higgs-doublet extension of each of the preceding models, assuming that the expanded Yukawa sector is that of the two-Higgs-doublet model of type II. We show that in these two-Higgs-doublet-portal models significant portions of the DM mass regions excluded in the simplest scenarios by direct search bounds can be reclaimed due to suppression of the effective DM interactions with nucleons at some ratios of the $CP$-even Higgs bosons' couplings to the up and down quarks. The regained parameter space contains areas which can yield a DM-nucleon scattering cross-section that is far less than its current experimental limit or even goes below the neutrino-background floor.

Direct searches for dark matter (DM) by the LUX and PandaX-II Collaborations employing xenon-based detectors have recently come up with the most stringent limits to date on the spin-independent elastic scattering of DM off nucleons. For Higgs-portal scalar DM models, the new results have precluded any possibility of accommodating low-mass DM as suggested by the DAMA and CDMS II Si experiments utilizing other target materials, even after invoking isospin-violating DM interactions with nucleons. In the simplest model, SM+D, which is the standard model plus a real singlet scalar named darkon acting as the DM candidate, the LUX and PandaX-II limits rule out DM masses roughly from 4 to 450 GeV, except a small range around the resonance point at half of the Higgs mass where the interaction cross-section is near the neutrino-background floor. In the THDM II+D, which is the type-II two-Higgs-doublet model combined with a darkon, the region excluded in the SM+D by the direct searches can be recovered due to suppression of the DM effective interactions with nucleons at some values of the ratios of Higgs couplings to the up and down quarks, making the interactions significantly isospin-violating. However, in either model, if the 125-GeV Higgs boson is the portal between the dark and SM sectors, DM masses less than 50 GeV or so are already ruled out by the LHC constraint on the Higgs invisible decay. In the THDM II+D, if the heavier $CP$-even Higgs boson is the portal, theoretical restrictions from perturbativity, vacuum stability, and unitarity requirements turn out to be important instead and exclude much of the region below 100 GeV. For larger DM masses, the THDM II+D has plentiful parameter space that corresponds to interaction cross-sections under the neutrino-background floor and therefore is likely to be beyond the reach of future direct searches without directional sensitivity.

A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy.

There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy.

The frequency ratios and of HD 49933 exhibit an increase at high frequencies. This behavior also exists in the ratios of other stars, which is considered to result from the low signal-to-noise ratio and the larger line width at the high-frequency end and could not be predicted by stellar models in previous work. Our calculations show that the behavior not only can be reproduced by stellar models, but can be predicted by asymptotic formulas of the ratios. The frequency ratios of the Sun, too, can be reproduced well by the asymptotic formulas. The increased behavior derives from the fact that the gradient of mean molecular weight at the bottom of the radiative region hinders the propagation of p-modes, while the hindrance does not exist in the convective core. This behavior should exist in the ratios of stars with a large convective core. The characteristic of the ratios at high frequencies provides a strict constraint on stellar models and aids in determining the size of the convective core and the extent of overshooting. Observational constraints point to a star with $M=1.28\pm0.01 M_{\odot}$, $R=1.458\pm0.005 R_{\odot}$, $t=1.83\pm0.1$ Gyr, $r_{cc}=0.16\pm0.02 R_{\odot}$, $\alpha=1.85\pm0.05$, and $\delta_{ov}=0.6\pm0.2$ for HD 49933.

We report in this letter for the first time the numerical simulations of muon and neutron flux variations at the surface of the earth with varying air densities in the troposphere and stratosphere. The simulated neutron and muon flux variations are in very good agreement with the measured neutron flux variation in Oulu and the muon flux variation in Atlanta. We conclude from this study that the stratosphere air density variation dominates the effects on the muon flux changes while the density variation in troposphere mainly influences the neutron count variation. These results pave a new path for systematically studying the global temperature evolution using worldwide cosmic ray data.

Seesaw mechanism provides a natural explanation of light neutrino masses through suppression of heavy seesaw scale. In inverse seesaw models the seesaw scale can be much lower than that in the usual seesaw models. If terms inducing seesaw masses are further induced by loop corrections, the seesaw scale can be lowered to be in the range probed by experiments at the LHC without fine tuning. In this paper we construct models in which inverse seesaw neutrino masses are generated at two loop level. These models also naturally have dark matter candidates. Although the recent data from Xenon100 put stringent constraint on the models, they can be consistent with data on neutrino masses, mixing, dark matter relic density and direct detection. These models also have some interesting experimental signatures for collider and flavor physics.

The most recent results of searches at the LHC for the Higgs boson h have turned up possible hints of such a particle with mass m_h about 125 GeV consistent with standard model (SM) expectations. This has many potential implications for the SM and beyond. We consider some of them in the contexts of a simple Higgs-portal dark matter (DM) model, the SM plus a real gauge-singlet scalar field D as the DM candidate, and a couple of its variations. In the simplest model with one Higgs doublet and three or four generations of fermions, for D mass m_D < m_h/2 the invisible decay h -> DD tends to have a substantial branching ratio. If future LHC data confirm the preliminary Higgs indications, m_D will have to exceed m_h/2. To keep the DM lighter than m_h/2, one will need to extend the model and also satisfy constraints from DM direct searches. The latter can be accommodated if the model provides sizable isospin violation in the DM-nucleon interactions. We explore this in a two-Higgs-doublet model combined with the scalar field D. This model can offer a 125-GeV SM-like Higgs and a light DM candidate having isospin-violating interactions with nucleons at roughly the required level, albeit with some degree of fine-tuning.

The model-independent piecewise parametrizations (0-spline, linear-spline and cubic-spline) are used to estimate constraints of equation of state of dark energy ($w_{de}$) from current observational data (including SNIa, BAO and Hubble parameter) and the simulated future data. A combination of fitting results of $w_{de}$ from these three spline methods reveal essential properties of real equation of state $w_{de}$. It is shown that $w_{de}$ beyond redshift $z\sim0.5$ is poorly constrained from current data, and the mock future $\sim2300$ supernovae data give poor constraints of $w_{de}$ beyond $z\sim1$. The fitting results also indicate that there might exist a rapid transition of $w_{de}$ around $z\sim0.5$. The difference between three spline methods in reconstructing and constraining $w_{de}$ has also been discussed.

Reconnection of the self-generated magnetic fields in laser-plasma interaction was first investigated experimentally by Nilson {\it et al.} [Phys. Rev. Lett. 97, 255001 (2006)] by shining two laser pulses a distance apart on a solid target layer. An elongated current sheet (CS) was observed in the plasma between the two laser spots. In order to more closely model magnetotail reconnection, here two side-by-side thin target layers, instead of a single one, are used. It is found that at one end of the elongated CS a fan-like electron outflow region including three well-collimated electron jets appears. The ($>1$ MeV) tail of the jet energy distribution exhibits a power-law scaling. The enhanced electron acceleration is attributed to the intense inductive electric field in the narrow electron dominated reconnection region, as well as additional acceleration as they are trapped inside the rapidly moving plasmoid formed in and ejected from the CS. The ejection also induces a secondary CS.

We present the results of our optical monitoring of the BL Lac object S5 0716+714 on seven nights in 2006 December. The monitoring was carried out simultaneously at three optical wavelengths with a novel photometric system. The object did not show large-amplitude internight variations during this period. Intranight variations were observed on four nights and probably on one more. Strong bluer-when-brighter chromatism was detected on both intranight and internight timescales. The intranight variation amplitude decreases in the wavelength sequence of B', R', and V'. Cross correlation analyses revealed that the variability at the $B'$ and $V'$ bands lead that at the $R'$ band by about 30 minutes on one night.

Consistent descriptions of the equation of states, and information about transport coefficients of deuterium-tritium mixture are demonstrated through quantum molecular dynamic (QMD) simulations (up to a density of 600 g/cm$^{3}$ and a temperature of $10^{4}$ eV). Diffusion coefficients and viscosity are compared with one component plasma model in different regimes from the strong coupled to the kinetic one. Electronic and radiative transport coefficients, which are compared with models currently used in hydrodynamic simulations of inertial confinement fusion, are evaluated up to 800 eV. The Lorentz number is also discussed from the highly degenerate to the intermediate region.

X-ray loud M dwarfs are a major source of by-product (contamination) in the X-ray band of the multiwavelength quasar survey (MWQS). As a by-product, the low dispersion spectra of 22 M dwarfs are obtained in which the spectra of 16 sources are taken for the first time. The spectral types and distance of the sample are given based on spectral indices CaH2, CaH3, and TiO5. The parameter {\zeta}TiO/CaH is calculated to make the metallicity class separation among dwarfs, subdwarfs and extreme subdwarfs. We also discuss the distributions in the diagrams of Log(Lx/Lbol) versus spectral type and infrared colors.

A model relating radiative seesaw and minimal dark matter mass scales without beyond the standard model (SM) gauge symmetry (R$\nu$MDM) is constructed. In addition to the SM particles, the R$\nu$MDM contains, a Majorana fermion multiplet $N_R$ and scalar multiplet $\chi$ that transform respectively as $(1,5,0)$ and $(1,6,-1/2)$ under the SM gauge group $SU(3)_C\times SU(2)_L\times U(1)_Y$. The neutral component $N_R^0$ plays the role of dark matter with a mass in the range of 9 to 10 TeV. This scale also sets the lower limit for the scale for the heavy degrees of freedom in $N_R$ and $\chi$ which generate light neutrino masses through the radiative seesaw mechanism. The model predicts an $N_R^0$-nucleus scattering cross section that would be accessible with future dark matter direct detection searches as well as observable effects in present and searches for charged lepton flavor violating processes, such as $l_i\to l_j \gamma$ and $\mu - e$ conversion.

We study the thermophysical properties of dense helium plasmas by using quantum molecular dynamics and orbital-free molecular dynamics simulations, where densities are considered from 400 to 800 g/cm$^{3}$ and temperatures up to 800 eV. Results are presented for the equation of state. From the Kubo-Greenwood formula, we derive the electrical conductivity and electronic thermal conductivity. In particular, with the increase in temperature, we discuss the change in the Lorenz number, which indicates a transition from strong coupling and degenerate state to moderate coupling and partial degeneracy regime for dense helium.

An NUV-optical diagram made for sources from the secend Galaxy Evolution Explorer (GALEX) Ultraviolet Variability (GUVV-2) Catalog provide us a method to tentatively classify the unknown GUVV2 sources by their NUV-optical magnitudes. On the purpose of testing the correctness and generality of the method, we carry out a program on the spectroscopic observations of the unidentified GUVV2 sources. The spectroscopic identification for these 37 sources are 19 type -A to -F stars, 10 type -G to -K stars and 7 M dwarf stars together with an AGN. We also present the light curves in R-band for two RR Lyrae star candidates selected from the NUV-optical diagram, both of which perform cyclic variations. Combining there light curves and colors, we classify them as RR Lyrae stars. To confirm the results, we shows a color-color diagram for the 37 newly spectroscopically identified objects compared with the previously identified ones, which manifests good consistence with our former results, indicating that the ultroviolet variable sources can be initially classified by their NUV/optical color-color diagram.

The Standard Model (SM) plus a real gauge-singlet scalar field dubbed darkon (SM+D) is the simplest model possessing a weakly interacting massive particle (WIMP) dark-matter candidate. In this model, the parameters are constrained from dark matter relic density and direct searches. The fact that interaction between darkon and SM particles is only mediated by Higgs boson exchange may lead to significant modifications to the Higgs boson properties. If the dark matter mass is smaller than a half of the Higgs boson mass, the Higgs boson can decay into a pair of darkons resulting in a large invisible branching ratio. The Higgs boson will be searched for at the LHC and may well be discovered in the near future. If a Higgs boson with a small invisible decay width will be found, the SM+D model with small dark matter mass will be in trouble. We find that by extending the SM+D to a two-Higgs-doublet model plus a darkon (THDM+D) it is possible to have a Higgs boson with a small invisible branching ratio and at the same time the dark matter can have a low mass. We also comment on other implications of this model.

The direct-search experiment for dark matter performed by the CDMS II Collaboration has observed two candidate events. Although these events cannot be interpreted as significant evidence for the presence of weakly interacting massive particle (WIMP) dark matter (DM), the total CDMS II data have led to an improved upper-limit on the WIMP-nucleon spin-independent cross-section. We study some implications of these results for the simplest WIMP DM model, the SM+D, which extends the standard model (SM) by the addition of a real SM-singlet scalar field dubbed darkon to play the role of the DM. We find that, although the CDMS II data rule out a sizable portion of parameter space of the model, a large part of the parameter space is still allowed. We obtain strong correlations among the darkon mass, darkon-nucleon cross-section, mass of the Higgs boson, and branching ratio of its invisible decay. We point out that measurements of the Higgs invisible branching-ratio at the LHC can lift some possible ambiguities in determining the darkon mass from direct DM searches.

The standard model (SM) plus a real gauge-singlet scalar field dubbed darkon (SM+D) is the simplest model possessing a weakly interacting massive particle (WIMP) dark-matter candidate. The upper limits for the WIMP-nucleon elastic cross-section as a function of WIMP mass from the recent XENON10 and CDMS-II experiments rule out darkon mass ranges from 10 to (50,70,75) GeV for Higgs-boson masses of (120,200,350) GeV, respectively. This may exclude the possibility of the darkon providing an explanation for the gamma-ray excess observed in the EGRET data. We show that by extending the SM+D to a two-Higgs-doublet model plus a darkon the experimental constraints on the WIMP-nucleon interactions can be circumvented due to suppression occurring at some values of the product tan(alpha)tan(beta), with alpha being the neutral-Higgs mixing angle and tan(beta) the ratio of vacuum expectation values of the Higgs doublets. We also comment on the implication of the darkon model for Higgs searches at the LHC.

We show that black holes can be quantized in an intuitive and elegant way with results in agreement with conventional knowledge of black holes by using Bohr's idea of quantizing the motion of an electron inside the atom in quantum mechanics. We find that properties of black holes can be also derived from an Ansatz of quantized entropy $\Delta S=4\pi k {\Delta R / \lambdabar}$, which was suggested in a previous work to unify the black hole entropy formula and Verlinde's conjecture to explain gravity as an entropic force. Such an Ansatz also explains gravity as an entropic force from quantum effect. This suggests a way to unify gravity with quantum theory. Several interesting and surprising results of black holes are given from which we predict the existence of primordial black holes ranging from Planck scale both in size and energy to big ones in size but with low energy behaviors.

We study entropic force effects on black holes and photons. We find that application of an entropic analysis restricts the radial change $\Delta R$ of a black hole of radius $R_{\mathrm{H}}$, due to a test particle of a Schwartzchild radius $R_{h}$ moving towards the black hole by $\Delta x$ near black body surface, to be given by a relation $R_{\mathrm{H}} \Delta R= R_h \Delta x/2$, or ${\Delta R}/{\lambdabar_M} = {\Delta x}/{2 \lambdabar_m}$. We suggest a new rule regarding entropy changes in different dimensions, $\Delta S= 2\pi k D \Delta l /\lambdabar,$ which unifies Verlinde's conjecture and the black hole entropy formula. We also propose to extend the entropic force idea to massless particles such as a photon. We find that there is an entropic force on a photon of energy $E_\gamma$, with $F=G M m_{\gamma}/R^2$, and therefore the photon has an effective gravitational mass $m_\gamma = E_\gamma/c^2$.

Dark matter annihilation is one of the leading explanations for the recently observed $e^\pm$ excesses in cosmic rays by PAMELA, ATIC, FERMI-LAT and HESS. Any dark matter annihilation model proposed to explain these data must also explain the fact that PAMELA data show excesses only in $e^\pm$ spectrum but not in anti-proton. It is interesting to ask whether the annihilation mode into anti-proton is completely disallowed or only suppressed at low energies. Most models proposed have negligible anti-protons in all energy ranges. We show that the leptocentric $U(1)_{B-3L_i}$ dark matter model can explain the $e^\pm$ excesses with suppressed anti-proton mode at low energies, but at higher energies there are sizable anti-proton excesses. Near future data from PAMELA and AMS can provide crucial test for this type of models. Cosmic $\gamma$ ray data can further rule out some of the models. We also show that this model has interesting cosmic neutrino signatures.

Recently data from PAMELA, ATIC, FERMI-LAT and HESS show that there are $e^{\pm}$ excesses in the cosmic ray energy spectrum. PAMELA observed excesses only in $e^+$, but not in anti-proton spectrum. ATIC, FERMI-LAT and HESS observed excesses in $e^++e^-$ spectrum, but the detailed shapes are different which requires future experimental observations to pin down the correct data set. Nevertheless a lot of efforts have been made to explain the observed $e^\pm$ excesses, and also why PAMELA only observed excesses in $e^+$ but not in anti-proton. In this brief review we discuss one of the most popular mechanisms to explain the data, the dark matter annihilation. It has long been known that about 23% of our universe is made of relic dark matter. If the relic dark matter was thermally produced, the annihilation rate is constrained resulting in the need of a large boost factor to explain the data. We will discuss in detail how a large boost factor can be obtained by the Sommerfeld and Briet-Wigner enhancement mechanisms. Some implications for particle physics model buildings will also be discussed.

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.

Unparticle $\U$ with scaling dimension $d_\U$ has peculiar thermal properties due to its unique phase space structure. We find that the equation of state parameter $\omega_\U$, the ratio of pressure to energy density, is given by $1/(2d_\U +1)$ providing a new form of energy in our universe. In an expanding universe, the unparticle energy density $\rho_\U(T)$ evolves dramatically differently from that for photons. For $d_\U >1$, even if $\rho_\U(T_D)$ at a high decoupling temperature $T_D$ is very small, it is possible to have a large relic density $\rho_\U(T^0_\gamma)$ at present photon temperature $T^0_\gamma$, large enough to play the role of dark matter. We calculate $T_D$ and $\rho_\U(T^0_\gamma)$ using photon-unparticle interactions for illustration.

We study baryon number violating nucleon decays induced by unparticle interactions with the standard model particles. We find that the lowest dimension operators which cause nucleon decays can arise at dimension 6 + (d_s-3/2) with the unparticles being a spinor of dimension d_s=d_\U +1/2. For scalar and vector unparticles of dimension d_\U, the lowest order operatoers arise at 6+d_\U and 7+d_\U dimensions,respectively. Comparing the spinor unparticle induced n \to O^s_\U and experimental bound on invisible decay of a neutron from KamLAND, we find that the scale for unparticle physics is required to be larger than 10^{10} GeV for d_\U < 2 if the couplings are set to be of order one. For scalar and vector unparticles, the dominant baryon number violating decay modes are n\to \bar \nu + O_\U (O^\mu_\U) and p \to e^+ + O_\U (O^\mu_\U). The same experimental bound puts the scales for scalar and vector unparticle to be larger than 10^{7} and 10^{5} GeV for d_\U <2 with couplings set to be of order one. Data on, p \to e^+ invisible, puts similar constraints on unparticle interactions.

The ROSAT X-ray source 1RXS J114003.0+124112 was identified as a starburst galaxy at redshift 0.177 by He et al. (2001). Meanwhile, the authors also noted that the source is almost two orders of magnitude brighter in X-ray than the X-ray-brightest starburst galaxy and it seems to be in a merging system, making this source an enigmatic system for further observations. This paper reports the re-identification of 1RXS J114003.0+124112 with the observations on the 2.6 m telescope at Byurakan Astrophysical Observatory, Armenia and with the SDSS data. The results indicate that the starburst activity is associated with the brighter object in the system, while the fainter object is a typical Seyfert 1 galaxy at redshift 0.282. Therefore, the two objects are not in a merging system, and the Seyfert 1 galaxy naturally accounts for the high X-ray flux. Three more objects reside in the vicinity, but they are all too faint to be responsible for the high X-ray flux.

The well-documented E1 relationships are first extended to infrared color $\alpha(60,25)$ and flux ratio [OIII]/H$\beta_{\rm{n}}$ by comparing emission line properties to continuum properties in infrared wavelength. Both direct correlations and a principle component analysis are used in a sample of 50 IRAS IR-selected Seyfert 1.5 galaxies. In addition to confirm the correlations of E1 in Boroson & Green (1992), our Eigenvector 1 turns out to be dominated by mid-infrared color $\alpha$(60,25), and most strongly effected by RFe, [OIII]/H$\beta_{\rm{n}}$, and EW(H$\beta_{\rm{b}}$). Our analysis indicate that the objects with large E1 tend to co-existent with relatively young nuclear stellar populations, which implies that the E1 is related with nuclear starformation history. The IR-dominated Eigenvector 1 can be, therefore, inferred to be interpreted as ``age'' of AGN. In confirmation of Xu et al. (2003), it is clear that the extreme Seyfert galaxies with both large RFe and large [OIII]/H$\beta_{\rm{n}}$ are rare in our Universe.

We study the Berry phase in neutrino oscillations for both Dirac and Majorana neutrinos. In order to have a Berry phase, the neutrino oscillations must occur in a varying medium, the neutrino-background interactions must depend on at least two independent densities, and also there must be CP violation if the neutrino interactions with matter are mediated only by the standard model W and Z boson exchanges which implies that there must be at least three generations of neutrinos. The CP violating Majorana phases do not play a role in generating a Berry phase. We show that a natural way to satisfy the conditions for the generation of a Berry phase is to have sterile neutrinos with active-sterile neutrino mixing, in which case at least two active and one sterile neutrinos are required. If there are additional new CP violating flavor changing interactions, it is also possible to have a non-zero Berry phase with just two generations.

The combination of the first-year Wilkinson Microwave Anisotropy Probe (WMAP) data with other finer scale cosmic microwave background (CMB) experiments (CBI and ACBAR) and two structure formation measurements (2dFGRS and Lyman $\alpha$ forest) suggest a $\Lambda$CDM cosmological model with a running spectral power index of primordial density fluctuations. Motivated by this new result on the index of primordial power spectrum, we present the first study on the predicted lensing probabilities of image separation in a spatially flat $\Lambda$CDM model with a running spectral index (RSI-$\Lambda$CDM model). It is shown that the RSI-$\Lambda$CDM model suppress the predicted lensing probabilities on small splitting angles of less than about 4$^{''}$ compared with that of standard power-law $\Lambda$CDM (PL-$\Lambda$CDM) model.

The variability of optical FeII blends in NGC 4051 is examined from spectra extracted from the AGN Watch program.In our analysis, the FeII complex are subtracted and measured with the following results. Firstly, the FeII variations were detected in NGC 4051 during a 3-year period. The identified FeII variations followed the variations in continuum closely. Secondly, the EW of FeII is reported to increase with the rising continuum flux, which is consistent with previous claims that there is no convincing Baldwin Effect in optical FeII. Thirdly, by comparing the variations of H$\beta$ and FeII, we find that RFe scales with continuum flux as $\rm{R_{Fe}\propto (5.0\pm0.8) \log(L/M)}$, which is significantly different from the theoretical expectations. Finally, in six selected Seyferts, four out of five Narrow-Line Seyfert 1 galaxies present positive correlation between RFe and the continuum flux. The negative correlations are identified in the remaindng two objects that have relatively broad profiles of H$\beta$ ($\rm{FWHM>1500 km\ s^{-1}}$).We argue that the difference of electron density of broad line clouds and/or variability behavior of incident high-energy radiation can explain the dichotomy in variability behavior of RFe.

Although the existence of asymmetrical profile of [OIII]$\lambda$5007 has been discovered for ages, its filiation and physics are poorly understood. Two new spectra of I ZWI taken on Nov 16, 2001 and on Dec 3, 2002 were compared with the spectra taken by BG92. Following results are obtained. 1)The certain variations of broad [OIII] during about 10 years separating the observations are identified. The inferred length scale of broad [OIII] emitting region ranges from 0.3pc to 3pc. By assuming a Keplerian motion in emitting region, the material emitting broad [OIII] is likely to be located at transient emission line region, between BLR and NLR. 2)We find a positive relation between the FeII emission and flux of H$\beta$(or continuum). On the other hand, the parameter RFe decreases with ionizing continuum marginally. 3)We detect a low ionized NLR in I ZWI, because of the low flux ratios $\rm{[OIII]_{n}/H\beta_{n}}$($\sim1.7$).

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

Recent measurements of type Ia supernovae as well as other concordant observations suggest that the expansion of our universe is accelerating. A dark energy component has usually been invoked as the most feasible mechanism for the acceleration. However, the effects arising from possible extra dimensions can mimic well the role of a dark energy through a modified Friedmann equation. In this work, we investigate some observational constraints on a scenario in which this modification is given by $H^2 = {8\pi G \over 3} (\rho +C\rho^n)$. We mainly focus our attention on the constraints from recent measurements of the dimensionless coordinate distances to type Ia supernovae and Fanaroff-Riley type IIb radio galaxies compiled by Daly and Djorgovski (2003) and the X-ray gas mass fractions in clusters of galaxies published by Allen et al. (2002,2003). We obtain the confidence region on the power index $n$ of the modificative term and the density parameter $\Omega_m$ of the universe from a combined analysis of these databases. It is found that $n=0.06^{+0.22}_{-0.18}$ and $\Omega_m=0.30^{+0.02}_{-0.02}$, at the 95.4% confidence level, which is consistent within the errors with the standard $\Lambda$CDM model. These parameter ranges give a universe whose expansion swithes from deceleration to acceleration at a redshift between 0.52 to 0.73.

A sample of close major-merger pairs (projected separation ${\rm 5 \leq r \leq 20 h^{-1}}$ kpc, ${\rm K_s}$ band magnitude difference $\delta {\rm K_s} \leq 1$ mag) is selected from the matched 2MASS-2dFGRS catalog of Cole et al. (2001). The pair primaries are brighter than ${\rm K_s} = 12.5$ mag. After corrections for various biases, the comparison between counts in the paired galaxy sample and counts in the parent sample shows that for the local `M* galaxies' sampled by flux limited surveys, the fraction of galaxies in the close major-merger pairs is 1.70$\pm 0.32%$. Using 38 paired galaxies in the sample, a ${\rm K_s}$ band luminosity function (LF) is calculated. This is the first unbiased LF for a sample of objectively defined interacting/merging galaxies in the local universe, while all previously determined LFs of paired galaxies are biased by mistreating paired galaxies as singles. A stellar mass function (MF) is translated from the LF. Compared to the LF/MF of 2MASS galaxies, a differential pair fraction function is derived. The results suggest a trend in the sense that less massive galaxies may have lower chance to be involved in close major-merger pairs than more massive galaxies. The algorithm presented in this paper can be easily applied to much larger samples of 2MASS galaxies with redshifts in near future.

Using recent measurements of angular size of high-z milliarcsecond compact radio sources compiled by Gurvits, Kellermann and Frey (1999) and X-ray gas mass fraction of galaxy clusters published by Allen et al. (2002,2003), we explore their bounds on the equation of state, omega_x = p_x/rho_x, of the dark energy, whose existence has been congruously suggested by various cosmological observations. We relaxe the usual constraint omega_x >= -1, and find that combining the two databases yields a nontrivial lower bound on omega_x. Under the assumption of a flat universe, we obtain a bound -2.22 < omega_x < -0.62 at 95.4% confidence level. The 95.4% confidence bound goes to -1 =< omega_x < -0.60 when the constraint omega_x >= -1 is imposed.

Casimir vacuum energy is divergent. It needs to be regularized. The regularization introduces a renormalization scale which may lead to a scale dependent cosmological constant. We show that the requirement of physical cosmological constant is renormalization scale independent provides important constraints on possible particle contents and their masses in particle physics models. In the Standard Model of strong and electroweak interactions, besides the Casimir vacuum energy there is also vacuum energy induced from spontaneous symmetry breaking. The requirement that the total vacuum energy to be scale independent dictates the Higgs mass to be $m_H^4 = 8\sum_i N_i m^4_i - 12m^4_W - 6 m^4_Z$, where the summation is over fermions and $N_i$ equals to 3 and 1 for quarks and leptons, respectively. The Higgs mass is predicted to be approximately 382 GeV.

We study bound object formation in a background universe accelerated by cosmological constant and by quintessence. If the acceleration lasts forever, due to the existence of event horizon, one would have naively expected the universe to approach a state of cold death. However, we find that many local regions in the universe can in fact be protected by their own gravity to form mini-universes, provided that their present matter densities exceed some critical value. In the case with cosmological constant ($\Lambda$CDM cosmology) the condition of forming mini-universe is that the ratio of present density parameters $\Omega_{i}^0/\Omega_\Lambda$ ought to be larger than a critical value 3.63. Such mini-universes typically weigh less than $2\times 10^{14}$ solar masses, with the lighter ones having tight and compact configurations. In the case with quintessence, the final ratio $\Omega_i/\Omega_q$ of a mini-universe is found to be always larger than $w_q^2-w_q$, where $w_q$ is the equation of state parameter.

Using semi-analytic approach, we present an estimate of the properties of the largest virialized dark halos in the present universe for three different scenarios of structure formation: SCDM, LCDM and OCDM models. The resulting virial mass and temperature increase from the lowest values of $1.6 \times 10^{15}h^{-1}M_{\odot}$ and 9.8 keV in OCDM, the mid-range values of $9.0 \times 10^{15}h^{-1}M_{\odot}$ and 31 keV in LCDM, to the highest values of $20.9 \times 10^{15}h^{-1}M_{\odot}$, 65 keV in SCDM. As compared with the largest virialized object seen in the universe, the richest clusters of galaxies, we can safely rule out the OCDM model. In addition, the SCDM model is very unlikely because of the unreasonably high virial mass and temperature. Our computation favors the prevailing LCDM model in which superclusters may be marginally regarded as dynamically-virialized systems.

It has been argued in the literature that if a universe is expanding with an accelerating rate indefinitely, it presents a challenge to string theories due to the existence of event horizons. We study the fate of a currently accelerating universe. We show that the universe will continue to accelerate indefinitely if the parameter $\omega = p/\rho$ of the equation of state is a constant, no matter how many different types of energy (matter, radiation, quintessence, cosmological constant and etc) are contained in the universe. This type of universe will always exhibit an event horizon indicating that such a universe may not be derived from string theories. We also comment on some related issues.

We study gravitational lensing and the bending of light in low energy scale (M_S) gravity theories with extra space-time dimensions n. We find that due to the presence of spin-2 Kaluza-Klein states from compactification, a correction to the deflection angle with a strong quadratic dependence on the photon energy is introduced. No deviation from the Einstein General Relativity prediction for the deflection angle for photons grazing the Sun in the visible band with 15% accuracy (90% c.l.) implies that the scale M_S has to be larger than 1.4(2/(n-2))^{1/4} TeV and approximately 4 TeV for n=2. This lower bound is comparable with that from collider physics constraints. Gravitational lensing experiments with higher energy photons can provide stronger constraints.