Using Stergioulas's RNS code for investigating fast pulsars with Equation of States (EOSs) on the causality surface (where the speed of sound equals to that of light) of the high-density EOS parameter space satisfying all known constraints from both nuclear physics and astrophysics, we show that the GW190814's secondary component of mass $(2.50-2.67)$ M$_{\odot}$ can be a super-fast pulsar spinning faster than 971 Hz about 42\% below its Kepler frequency. There is a large and physically allowed EOS parameter space below the causality surface where pulsars heavier than 2.50 M$_{\odot}$ are supported if they can rotate even faster with critical frequencies depending strongly on the high-density behavior of nuclear symmetry energy.

Based on the dynamics of single scalar field slow-roll inflation and the theory of reheating, we investigate the generalized natural inflationary (GNI) model. Concretely, we give constraints on the scalar spectral index $n_{s}$ and tensor-to scalar ratio $r$ for $\Lambda$CDM $+r$ model according to the latest data from Plack 2018 TT,TE,EE+lowE+lensing (P18) and BICEP2/Keck 2015 season (BK15), i.e., $n_{s}=0.9659\pm0.0044$ at $68\%$ confidence level (CL) and $r<0.0623$ at $95\%$CL. We find that the GNI model is favored by P18 plus BK15 in the ranges of $\log_{10}(f/M_{p})=0.62^{+0.17}_{-0.18}$ and $m=0.35^{+0.13}_{-0.23}$ at $68\%$CL. In addition, the corresponding predictions of the general and two-phase reheating are respectively discussed. It follows that the parameter $m$ has the significant effect on the model behaviors.

Orbital resonances play an important role in the dynamics of planetary systems. Classical theoretical analyses found in textbooks report that libration widths of first order mean motion resonances diverge for nearly circular orbits. Here we examine the nature of this divergence with a non-perturbative analysis of a few first order resonances interior to a Jupiter-mass planet. We show that a first order resonance has two branches, the pericentric and the apocentric resonance zone. As the eccentricity approaches zero, the centers of these zones diverge away from the nominal resonance location but their widths shrink. We also report a novel finding of "bridges" between adjacent first order resonances: at low eccentricities, the apocentric libration zone of a first order resonance smoothly connects with the pericentric libration zone of the neighboring first order resonance. These bridges may facilitate resonant migration across large radial distances in planetary systems, entirely in the low eccentricity regime.

Muons in neutron stars (NSs) play especially important roles in addressing several interesting new physics questions associated with detecting as well as understanding interactions and astrophysical effects of muonphilic dark matter particles. The key model inputs for studying the latter are the total muon mass $M_{\mu}$, the muon mass fraction $M_{\mu}/M_{\rm NS}$ over the NS mass $M_{\rm NS}$ and the muon radial density profile $\rho_{\mu}(r)$ in NSs of varying masses. We investigate these quantities within a minimum model for the core of NSs consisting of neutrons, protons, electrons, and muons using an explicitly isospin-dependent parametric Equation of State (EOS) constrained by available nuclear laboratory experiments and the latest astrophysical observations of NS masses, radii and tidal deformabilities. We found that the absolutely maximum muon mass $M_{\mu}$ and its mass fraction $M_{\mu}/M_{\rm NS}$ in the most massive NSs allowed by causality are about 0.025 $M_\odot$ and 1.1\%, respectively. For the most massive NS of mass 2.14 $M_\odot$ observed so far, they reduce to about 0.020 $M_\odot$ and 0.9\%, respectively. We also study respective effects of individual parameters describing the EOS of high-density neutron-rich nucleonic matter on the muon contents in NSs with varying masses. We found that the most important but uncertain nuclear physics ingredient for determining the muon contents in NSs is the high-density nuclear symmetry energy.

The constraints on tidal deformability $\Lambda$ of neutron stars are first extracted from GW170817 by LIGO and Virgo Collaborations but the relation between radius $R$ and tidal deformability $\Lambda$ is still nuder debate. Using an isospin-dependent parameterized equation of state (EOS), we study the relation between $R$ and $\Lambda$ of neutron stars and its dependence on parameters of symmetry energy $E_{\rm sym}$ and EOS of symmetric nuclear matter $E_0$ when the mass is fixed as $1.4$ $M_\odot$, $1.0$ $M_\odot$, and $1.8$ $M_\odot$, respectively. We find that, though the changes of high order parameters of $E_{\rm sym}$ and $E_0$ can shift the individual values of $R_{1.4}$ and $\Lambda_{1.4}$ to different values, the $R_{1.4}\sim\Lambda_{1.4}$ relation approximately locates at the same fitted curve. The slope of symmetry energy $L$ plays the dominated role in determining the $R_{1.4}\sim\Lambda_{1.4}$ relation. By checking the mass dependence of $R\sim\Lambda$ relation, the well fitted $R\sim\Lambda$ relation for 1.4 $M_\odot$ is broken for massive neutron stars.

When photons from distant galaxies and stars pass through our neighboring environment, the wavelengths of the photons would be shifted by our local gravitational potential. This local gravitational redshift effect can potentially have an impact on the measurement of cosmological distance-redshift relation. Using available supernovae data, Wojtak et al [1] found seemingly large biases of cosmological parameters for some extended models (non-flat $\Lambda$CDM, $w$CDM, etc.). Huang [2] pointed out that, however, the biases can be reduced to a negligible level if cosmic microwave background (CMB) data are added to break the strong degeneracy between parameters in the extended models. In this article we forecast the cosmological bias due to local gravitational redshifts for a future WFIRST-like supernovae survey. We find that the local gravitational redshift effect remains negligible, provided that CMB data or some future redshift survey data are added to break the degeneracy between parameters.

Determining the Equation of State (EOS) of dense neutron-rich nuclear matter is a shared goal of both nuclear physics and astrophysics. Except possible phase transitions, the density dependence of nuclear symmetry \esym is the most uncertain part of the EOS of neutron-rich nucleonic matter especially at supra-saturation densities. Much progresses have been made in recent years in predicting the symmetry energy and understanding why it is still very uncertain using various microscopic nuclear many-body theories and phenomenological models. Simultaneously, significant progresses have also been made in probing the symmetry energy in both terrestrial nuclear laboratories and astrophysical observatories. In light of the GW170817 event as well as ongoing or planned nuclear experiments and astrophysical observations probing the EOS of dense neutron-rich matter, we review recent progresses and identify new challenges to the best knowledge we have on several selected topics critical for understanding astrophysical effects of the nuclear symmetry energy.

We study implications of the very recently reported mass $M=2.17^{+0.11}_{-0.10}$M$_\odot$ of PSR~J0740+6620 on the Equation of State (EOS) of super-dense neutron-rich nuclear matter with respect to existing constraints on the EOS based on the mass $M=2.01\pm 0.04$M$_\odot$ of PSR~J0348+0432, the maximum tidal deformability of GW170817 and earlier results of various terrestrial nuclear laboratory experiments. The lower limit of the skewness $J_0$ measuring the stiffness of super-dense isospin-symmetric nuclear matter is raised raised from about -220 MeV to -150 MeV, reducing significantly its current uncertainty range. The lower bound of the high-density symmetry energy also increases appreciably leading to a rise of the minimum proton fraction in neutron stars at $\beta$-equilibrium from about 0 to 5\% around three times the saturation density of nuclear matter. The difficulties for some of the most widely used and previously well tested model EOSs to predict simultaneously both a maximum mass higher than 2.17 M$_\odot$ and a pressure consistent with that extracted from GW170817 present some interesting new challenges for nuclear theories.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

Within a minimum model for neutron stars consisting of nucleons, electrons and muons at $\beta$-equilibrium using about a dozen Equation of States (EOSs) from microscopic nuclear many-body theories and 40,000 EOSs randomly generated using an explicitly isospin-dependent parametric EOS model for high-density neutron-rich nucleonic matter within its currently known uncertainty range, we study correlations among the f-mode frequency, its damping time and the tidal deformability as well as the compactness of neutron stars. Except for quark stars, both the f-mode frequency and damping time of canonical neutron stars are found to scale with the tidal deformability independent of the EOSs used. Applying the constraint on the tidal deformability of canonical neutron stars $\Lambda_{1.4}=190^{+390}_{-120}$ extracted by the LIGO+VIRGO Collaborations from their improved analyses of the GW170817 event, the f-mode frequency and its damping time of canonical neutron stars are limited to 1.67 kHz - 2.18 kHz and 0.155 s - 0.255 s, respectively, providing a useful guidance for the ongoing search for gravitational waves from the f-mode oscillations of isolated neutron stars. Moreover, assuming either or both the f-mode frequency and its damping time will be measured precisely in future observations with advanced gravitational wave detectors, we discuss how information about the mass and/or radius as well as the still rather elusive nuclear symmetry energies at supra-saturation densities may be extracted.

By numerically inverting the Tolman-Oppenheimer-Volkov (TOV) equation using an explicitly isospin-dependent parametric Equation of State (EOS) of dense neutron-rich nucleonic matter, a restricted EOS parameter space is established using observational constraints on the radius, maximum mass, tidal polarizability and causality condition of neutron stars (NSs). The constraining band obtained for the pressure as a function of energy (baryon) density is in good agreement with that extracted recently by the LIGO+Virgo Collaborations from their improved analyses of the NS tidal polarizability in GW170817. Rather robust upper and lower boundaries on nuclear symmetry energies are extracted from the observational constraints up to about twice the saturation density $\rho_0$ of nuclear matter. More quantitatively, the symmetry energy at $2\rho_0$ is constrained to $E_{\rm{sym}}(2\rho_0)=46.9\pm10.1$ MeV excluding many existing theoretical predictions scattered between $E_{\rm{sym}}(2\rho_0)=15$ and 100 MeV. Moreover, by studying variations of the causality surface where the speed of sound equals that of light at central densities of the most massive neutron stars within the restricted EOS parameter space, the absolutely maximum mass of neutron stars is found to be 2.40 M$_{\odot}$ approximately independent of the EOSs used. This limiting mass is consistent with findings of several recent analyses and numerical general relativity simulations about the maximum mass of the possible super-massive remanent produced in the immediate aftermath of GW170817.

We derive in detail the equations of motion for the tensorial modes of primordial metric perturbations in the Coupled-Scalar-Tachyon Bounce Universe. We solve for the gravitational wave equations in the pre-bounce contraction and the post-bounce expansion epochs. To match the solutions of the tensor perturbations, we idealise the bounce process yet retaining the essential physical properties of the bounce universe. We put forward two matching conditions: one ensures the continuity of the gravitational wave functions and the other respects the symmetric nature of the bounce dynamics. The matching conditions connect the two independent modes of gravitational waves solutions before and after the bounce. We further analyse the scale dependence and time dependence of the gravitational waves spectra in the bounce universe and compare them with the primordial spectrum in the single field inflation scenario. We discuss the implications to early universe physics and present model independent observational signatures extracted from the bounce universe.

What can we learn about the density dependence of nuclear symmetry energy $E_{\rm{sym}}(\rho )$ from precise measurements of the radius ($R_{\rm{1.4}}$) and/or tidal polarizability ($\Lambda_{1.4}$) of canonical neutron stars (NSs) with a mass of 1.4 M$_\odot$? With the $E_{\rm{sym}}(\rho )$ parameterized using three parameters $L$, $K_{\rm{sym}}$, and $J_{\rm{sym}}$ which have the asymptotic meaning of being respectively the slope, curvature, and skewness of symmetry energy at saturation density, we found that, while both the $R_{\rm{1.4}}$ and $\Lambda_{1.4}$ depend strongly on the slope $L$, the $K_{\rm{sym}}$ and $J_{\rm{sym}}$ parameters characterizing the high-density behavior of $E_{\rm{sym}}(\rho )$ also play appreciable roles. Thus, there is not a simple relation between the $\Lambda_{\rm{1.4}}$/$R_{\rm{1.4}}$ and $L$ alone. Precise measurements of just the $\Lambda_{\rm{1.4}}$ and $R_{\rm{1.4}}$ can not completely determine the $E_{\rm{sym}}(\rho )$ but limit combinations of its parameters. In particular, stringent constraints approximately independent of the $J_{\rm{sym}}$ on the $L$-$K_{\rm{sym}}$ correlations can be obtained. However, infinite combinations of the larger (smaller) $L$ and smaller (larger) $K_{\rm{sym}}$ can lead to the same $\Lambda_{\rm{1.4}}$ and $R_{\rm{1.4}}$. Additional observables including those from terrestrial nuclear experiments are thus necessary to break this degeneracy in order to completely determine the density dependence of nuclear symmetry energy $E_{\rm{sym}}(\rho )$.

By assuming that only gravitation acts between dark matter (DM) and normal matter (NM), we studied DM admixed neutron stars (DANSs) using the two-fluid TOV equations. The NM and DM of compact stars are simulated by the relativistic mean field (RMF) theory and non-self-annihilating self-interacting fermionic model, respectively. The effects of the particle mass of fermionic DM $m_f$ and the interaction strength parameter $y$ on the properties of DANSs are investigated in detail. $m_f$ and $y$ are considered as the free parameters due to the lack of information about the particle nature of DM so far. For a DANS, we suggest a simple universal relationship $M_D^{\max}=(0.267 y +0.627-3.21\frac{M_N}{\M_{\odot}})( \frac{1\GeV}{{m_f}})^2 \M_{\odot}$ for $y>100$, where $M_D^{\max}$ is the maximum mass of DM existing in DANSs and $M_N$ is the mass of the neutron star without DM. For free fermion DM model ($y$=0), the relationship becomes $ M_D^{\max}=(0.627-0.027\frac{M_N^2}{\M_{\odot}^2}) ( \frac{1\GeV}{{m_f}})^2 \M_{\odot}$. The radius of DM $R_D$ shows a linear relationship with $M_D^{\max}$ in DANSs, namely $R_D=(7.02 \frac{M_D^{\max}}{ \M_{\odot}}+1.36)$~km. These conclusions are independent of the different NM EOSs from RMF theory. Such a kind of universal relationship connecting the nature of DM particle and mass of stars might shed light on the constraining the nature of the DM by indirect method.

Within the parameter space of equation of state (EOS) of dense neutron-rich matter limited by existing constraints mainly from terrestrial nuclear experiments, we investigate how the neutron star maximum mass $M_{\rm{max}}>2.01\pm0.04$ M$_\odot$, radius $10.62<R_{\rm{1.4}}< 12.83$ km and tidal deformability $\Lambda_{1.4}\leq800$ of canonical neutron stars all together constrain the EOS of dense neutron-rich nucleonic matter. While the 3-D parameter space of $K_{\rm{sym}}$ (curvature of nuclear symmetry energy), $J_{\rm{sym}}$ and $J_0$ (skewness of the symmetry energy and EOS of symmetric nuclear matter, respectively) are narrowed down significantly by the observational constraints, more data are needed to pin down the individual values of $K_{\rm{sym}}$, $J_{\rm{sym}}$ and $J_0$ with quantified uncertainties. The $J_0$ largely controls the maximum mass of neutron stars. While the EOS with $J_0=0$ is sufficiently stiff to support neutron stars as massive as 2.37 M$_{\odot}$, to support the hyperthetical ones as massive as 2.74 M$_{\odot}$ (composite mass of GW170817) requires $J_0$ to be larger than its currently known maximum value of about 400 MeV and beyond the causality limit. The upper limit on the tidal deformability of $\Lambda_{1.4}=800$ from the recent observation of GW170817 is found to provide upper limits on some EOS parameters consistent with but far less restrictive than the existing constraints of other observables studied.

For the propose of confronting updated pulsar observations with developed neutron star equation of states (EoSs), we employ four unified EoSs for both the core and the crust, namely BCPM, BSk20, BSk21, Shen-TM1, as well as two non-unified EoSs widely used in the literature, i.e. APR and GM1 EoS, which are commonly matched with the Negele-Vautherin and the Baym-Pethick-Sutherland crust EoS. All the core EoSs satisfy the recent observational constraints of the two massive pulsars whose masses are precisely measured. We show that the NS mass-equatorial radius relations are slightly affected by the smoothness at the core-crust matching interface. Moreover, the uncertainties in the crust EoS and the matching interface bring insignificant changes, even at maximally rotating (Keplerian) configurations. We also find that for all four unified EoSs, rotation can increase the star's gravitational mass up to $18\%-19\%$ and the equatorial radius by $29\%-36\%$, which are consistent with the previous calculations using non-unified EoSs. For stars as heavy as 1.4 M$_{\odot}$, the radius increase is more pronounced, reaching $41\%-43\%$, i.e. $5-6$ km. Moreover, by confronting the results using unified EoSs, which give the correct empirical values at saturation density, with two controversial determinations of the radius for the fast rotator 4U 1608-52, we address that a small radius may be better justified for this source.

We present $Chandra$ observations of seven broad absorption line (BAL) quasars at $z=0.863$-2.516 with redshifted BAL troughs (RSBALs). Five of our seven targets were detected by $Chandra$ in 4-13 ks exposures with ACIS-S. The $\alpha_{\rm ox}$ values, $\Delta\alpha_{\rm ox}$ values, and spectral energy distributions of our targets demonstrate they are all X-ray weak relative to expectations for non-BAL quasars, and the degree of X-ray weakness is consistent with that of appropriately-matched BAL quasars generally. Furthermore, our five detected targets show evidence for hard X-ray spectral shapes with a stacked effective power-law photon index of $\Gamma_{\rm eff}=0.5^{+0.5}_{-0.4}$. These findings support the presence of heavy X-ray absorption ($N_{\rm H}\approx 2 \times 10^{23}$ cm$^{-2}$) in RSBAL quasars, likely by the shielding gas found to be common in BAL quasars more generally. We use these X-ray measurements to assess models for the nature of RSBAL quasars, finding that a rotationally-dominated outflow model is favored while an infall model also remains plausible with some stipulations. The X-ray data disfavor a binary quasar model for RSBAL quasars in general.

A supramassive, strongly-magnetized millisecond neutron star (NS) has been proposed to be the candidate central engine of at least some short gamma-ray bursts (SGRBs), based on the "internal plateau" commonly observed in the early X-ray afterglow. While a previous analysis shows a qualitative consistency between this suggestion and the Swift SGRB data, the distribution of observed break time $t_b$ is much narrower than the distribution of the collapse time of supramassive NSs for the several NS equations-of-state (EoSs) investigated. In this paper, we study four recently-constructed "unified" NS EoSs, as well as three developed strange quark star (QS) EoSs within the new confinement density-dependent mass model. All the EoSs chosen here satisfy the recent observational constraints of the two massive pulsars whose masses are precisely measured. We construct sequences of rigidly rotating NS/QS configurations with increasing spinning frequency $f$, from non-rotating ($f = 0$) to the Keplerian frequency ($f = f_{\rm K}$), and provide convenient analytical parametrizations of the results. Assuming that the cosmological NS-NS merger systems have the same mass distribution as the Galactic NS-NS systems, we demonstrate that all except the BCPM NS EoS can reproduce the current $22\%$ supramassive NS/QS fraction constraint as derived from the SGRB data. We simultaneously simulate the observed quantities (the break time $t_b$, the break time luminosity $L_b$ and the total energy in the electromagnetic channel $E_{\rm total}$) of SGRBs, and find that while equally well reproducing other observational constraints, QS EoSs predict a much narrower $t_b$ distribution than that of the NS EoSs, better matching the data. We therefore suggest that the post-merger product of NS-NS mergers might be fast-rotating supramassive QSs rather than NSs.

N70 and N185 are two large ($\ge$100 pc in diameter) shell nebulae in the Large Magellanic Cloud (LMC). Their high [\ion{S}{2}]/H$\alpha$ ratios rival those of supernova remnants (SNRs), but they are not confirmed as SNRs. To study their physical nature, we have obtained \emph{XMM-Newton} X-ray observations and high-dispersion long-slit echelle spectroscopic observations of these two nebulae. The X-ray spectra of both nebulae can be well interpreted with an optically thin thermal ($\sim$0.2 keV) plasma with the average LMC abundance in a collisional ionization equilibrium. N70 encompasses the OB association LH114. Although N70 has a modest expansion velocity and essentially thermal radio emission, its diffuse X-ray luminosity ($\sim6.1\times10^{35}$ erg s$^{-1}$) is higher than that from a quiescent superbubble with N70's density, size, and expansion velocity; thus, N70 is most likely a superbubble that is recently energized by an interior SNR. N185 does not contain any known OB association, and its X-ray luminosity is an order of magnitude lower than expected if it is a quiescent superbubble. N185 has nonthermal radio emission and has high-velocity material expanding at nearly 200 km s$^{-1}$, similar to many known SNRs in the LMC. Its X-ray luminosity ($\sim1.9\times10^{35}$ erg s$^{-1}$) is also consistent with that of an evolved SNR. We therefore suggest that N185 is energized by a recent supernova.

We study the key factor to determine the dimensionless spin parameter $j\equiv cJ/(GM^2)$ of different kinds of uniformly rotating compact stars, including the traditional neutron stars, hyperonic neutron stars, and hybrid stars, and check the reliability of the results on various types of equations of state of dense matter. The equations of state from the relativistic mean field theory and the MIT bag model are adopted to simulate compact stars. Numerical calculations of rigidly rotating neutron stars are performed using the RNS code in the framework of general relativity by solving the Einstein equations for stationary axis-symmetric spacetime. The crust structure of compact stars is found to be a key factor to determine the maximum value of the spin parameter $j_{\rm max}$. For the stars with inclusion of the crust, $j_{\rm max}\sim 0.7$ is sustained for various kinds of compact stars with $M>0.5 M_{\odot}$, and is found to be insensitive to the mass of star and selected equations of state. For the traditional neutron stars, hyperonic neutron stars and hybrid stars without crust structure, the value $j_{\rm max}$ lies in the range of $[0.7, 1.0]$. Thus, not $j>0.7$ but $j>1$ could be treated as the candidate criterion to distinguish the strange quark stars from the other kinds of compact stars. Furthermore, a universal formula $j=0.48(f/f_K)^3-0.42(f/f_K)^2+0.63(f/f_K)$ is suggested to calculate the spin parameter at any rotational frequency for all kinds of compact stars with crust structure and $M>0.5M_{\odot}$.

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.

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.

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.

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 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 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.

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.

A new BeppoSAX broad-band (0.6-150 keV) spectrum of the Seyfert 2 galaxy Mkn 3 is presented. The spectrum provides a direct measurement of a large, neutral column of gas with Nh~10**24 cm**-2 in the source direction. The source, as bright as 3C 273 above 10 keV, has a steep (Gamma~1.8) spectrum without any evidence of a high-energy cutoff up to at least 150 keV. At lower energies, the data are best modeled with the addition of an unabsorbed reflection component. Combining these data with previous Ginga and ASCA observations, the Fe K$_{\alpha}$ and reflection continuum indicate that the reprocessed emission is responding slower than the intrinsic continuum variations suggesting a size of the reprocessor >= 2 pc. Identifying such a reprocessor with a (close to edge-on) obscuring torus, the overall result fits well into unified models since, presumably, one can interpret the strong absorption as due to transmission through the rim of the torus and the unabsorbed (directly viewed) reflection component as due to reprocessing from the torus inner surface.

We report broad band (0.1-200 keV) X-ray observations, made by BeppoSAX, of a sample of bright Seyfert 2 galaxies: NGC7172, NGC2110, NGC4507, Mkn 3 and NGC7674. These spectra provide a better understanding of the effects of X-ray reprocessing by cold material in the source and allow to put tighter constraints on the various spectral parameters involved. In particular, the data are used to determine, with less ambiguities than in the past, the shape of the intrinsic continuum emission by means of the high energy data. Within the small sample both Compton thin and Compton thick sources are found according to the expectations of the unified theory.