We have searched for exotic neutrino-electron interactions that could be produced by a neutrino millicharge, by a neutrino magnetic moment, or by dark photons using solar neutrinos in the XMASS-I liquid xenon detector. We observed no significant signals in 711 days of data. We obtain an upper limit for neutrino millicharge of 5.4$\times$10$^{-12} e$ at 90\% confidence level assuming all three species of neutrino have common millicharge. We also set flavor dependent limits assuming the respective neutrino flavor is the only one carrying a millicharge, $7.3 \times 10^{-12} e$ for $\nu_e$, $1.1 \times 10^{-11} e$ for $\nu_{\mu}$, and $1.1 \times 10^{-11} e$ for $\nu_{\tau}$. These limits are the most stringent yet obtained from direct measurements. We also obtain an upper limit for the neutrino magnetic moment of 1.8$\times$10$^{-10}$ Bohr magnetons. In addition, we obtain upper limits for the coupling constant of dark photons in the $U(1)_{B-L}$ model of 1.3$\times$10$^{-6}$ if the dark photon mass is 1$\times 10^{-3}$ MeV$/c^{2}$, and 8.8$\times$10$^{-5}$ if it is 10 MeV$/c^{2}$.

A search for Weakly Interacting Massive Particles (WIMPs) was conducted with the single-phase liquid-xenon detector XMASS through inelastic scattering in which $^{129}$Xe nuclei were excited, using an exposure ($\rm 327\; kg \times 800.0 \; days$) 48 times larger than that of our previous study. The inelastic excitation sensitivity was improved by detailed evaluation of background, event classification based on scintillation timing that distinguished $\gamma$-rays and $\beta$-rays, and simultaneous fitting of the energy spectra of $\gamma$-like and $\beta$-like samples. No evidence of a WIMP signal was found. Thus, we set the upper limits of the inelastic channel cross section at 90\% confidence level, for example, $4.1\times 10^{-39} \;{\rm cm^2}$ for a $200\; {\rm GeV}/c^2$ WIMP. This result provides the most stringent limits on the SD WIMP-neutron interaction and is better by a factor of 7.7 at $200\;{\rm GeV}/c^2$ than the existing experimental limit.

XMASS-I is a single-phase liquid xenon detector whose purpose is direct detection of dark matter. To achieve the low background requirements necessary in the detector, a new model of photomultiplier tubes (PMTs), R10789, with a hexagonal window was developed based on the R8778 PMT used in the XMASS prototype detector. We screened the numerous component materials for their radioactivity. During development, the largest contributions to the reduction of radioactivity came from the stem and the dynode support. The glass stem was exchanged to the Kovar alloy one and the ceramic support were changed to the quartz one. R10789 is the first model of Hamamatsu Photonics K. K. that adopted these materials for low background purposes and provided a groundbreaking step for further reductions of radioactivity in PMTs. Measurements with germanium detectors showed 1.2$\pm$0.3 mBq/PMT of $^{226}$Ra, less than 0.78 mBq/PMT of $^{228}$Ra, 9.1$\pm$2.2 mBq/PMT of $^{40}$K, and 2.8$\pm$0.2 mBq/PMT of $^{60}$Co. In this paper, the radioactive details of the developed R10789 are described together with our screening methods and the components of the PMT.

A search for dark matter using an underground single-phase liquid xenon detector was conducted at the Kamioka Observatory in Japan, particularly for Weakly Interacting Massive Particles (WIMPs). We have used 705.9 live days of data in a fiducial volume containing 97 kg of liquid xenon at the center of the detector. The event rate in the fiducial volume after the data reduction was ${\rm (4.2 \pm 0.2) \times 10^{-3} \, day^{-1}kg^{-1} keV_{ee}^{-1}}$ at ${\rm 5 \, keV_{ee}}$, with a signal efficiency of ${\rm 20\%}$. All the remaining events are consistent with our background evaluation, mostly of the "mis-reconstructed events" originated from $^{210}$Pb in the copper plates lining the detector's inner surface. The obtained upper limit on a spin-independent WIMP-nucleon cross section was ${\rm 2.2 \times 10^{-44} \, cm^{2}}$ for a WIMP mass of ${\rm 60 \, GeV/c^{2}}$ at the $90\%$ confidence level, which was the most stringent limit among results from single-phase liquid xenon detectors.

A search for dark matter (DM) with mass in the sub-GeV region (0.32-1 GeV) was conducted by looking for an annual modulation signal in XMASS, a single-phase liquid xenon detector. Inelastic nuclear scattering accompanied by bremsstrahlung emission was used to search down to an electron equivalent energy of 1 keV. The data used had a live time of 2.8 years (3.5 years in calendar time), resulting in a total exposure of 2.38 ton-years. No significant modulation signal was observed and 90% confidence level upper limits of $1.6 \times 10^{-33}$ cm$^2$ at 0.5 GeV was set for the DM-nucleon cross section. This is the first experimental result of a search for DM mediated by the bremsstrahlung effect. In addition, a search for DM with mass in the multi-GeV region (4-20 GeV) was conducted with a lower energy threshold than previous analysis of XMASS. Elastic nuclear scattering was used to search down to a nuclear recoil equivalent energy of 2.3 keV, and upper limits of 2.9 $\times$10$^{-42}$ cm$^2$ at 8 GeV was obtained.

Hidden photons and axion-like particles are candidates for cold dark matter if they were produced non-thermally in the early universe. We conducted a search for both of these bosons using 800 live-days of data from the XMASS detector with 327 kg of liquid xenon in the fiducial volume. No significant signal was observed, and thus we set constraints on the $\alpha' / \alpha$ parameter related to kinetic mixing of hidden photons and the coupling constant $g_{Ae}$ of axion-like particles in the mass range from 40 to 120 keV/$c^2$, resulting in $\alpha' / \alpha < 6 \times 10^{-26}$ and $g_{Ae} < 4 \times 10^{-13}$. These limits are the most stringent derived from both direct and indirect searches to date.

An annual modulation signal due to the Earth orbiting around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for this annual modulation with our single-phase liquid xenon XMASS-I detector. That analysis resulted in a slightly negative modulation amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with the same 832 kg target mass. When we assume weakly interacting massive particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of 2 to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ at 90\% confidence level with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1$-$20 keV energy region. This search probed this region with an exposure that was larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days in the energy region between 1 to 6 keV.

We analyze in this paper the existence of the "out-of-plane" equilibrium points in the restricted three-body problem with oblateness. From the series expansion of the potential function of an oblate asteroid, we show analytically all equilibrium points locate on the orbital plane of primaries and how artificial equilibrium points may arise due to an inappropriate application of the potential function. Using the closed form of the potential of a triaxial ellipsoid, we analytically demonstrate that the gravitational acceleration in $z$-direction is always pointing toward the equatorial plane, thus it could not be balanced out at any value of $z\neq 0$ and the out-of-plane equilibrium points cannot exist. The out-of-plane equilibrium points appear only when additional acceleration other than the gravitation from primaries is taken into account. We suggest that special attention must be paid to the application of the spherical harmonics expansion of potential to find the equilibrium points, especially when these points may be very close to the celestial body.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

We numerically determine the cosmological branch of the free function in a nonlocal metric-based modification of gravity which provides a relativistic generalization of Milgrom's Modified Newtonian Dynamics. We are able to reproduce the $\Lambda$CDM expansion history over virtually all of cosmic history, including the era of radiation domination during Big Bang Nucleosynthesis, the era of matter domination during Recombination, and most of the era of vacuum energy domination. The very late period of $0 \leq z < 0.0880$, during which the model deviates from the $\Lambda$CDM expansion history, is interesting because it causes the current value of the Hubble parameter to be about 4.5\% larger than it would be for the $\Lambda$CDM model. This may resolve the tension between inferences of $H_0$ which are based on data from large redshift and inferences based on Hubble plots.

A search for dark matter was conducted by looking for an annual modulation signal due to the Earth's rotation around the Sun using XMASS, a single phase liquid xenon detector. The data used for this analysis was 359.2 live days times 832 kg of exposure accumulated between November 2013 and March 2015. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the target nuclei, the exclusion upper limit of the WIMP-nucleon cross section 4.3$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ was obtained and we exclude almost all the DAMA/LIBRA allowed region in the 6 to 16 GeV/c$^2$ range at $\sim$10$^{-40}$cm$^2$. The result of a simple modulation analysis, without assuming any specific dark matter model but including electron/$\gamma$ events, showed a slight negative amplitude. The $p$-values obtained with two independent analyses are 0.014 and 0.068 for null hypothesis, respectively. we obtained 90\% C.L. upper bounds that can be used to test various models. This is the first extensive annual modulation search probing this region with an exposure comparable to DAMA/LIBRA.

In the late stages of nuclear burning for massive stars ($M>8~M_{\sun}$), the production of neutrino-antineutrino pairs through various processes becomes the dominant stellar cooling mechanism. As the star evolves, the energy of these neutrinos increases and in the days preceding the supernova a significant fraction of emitted electron anti-neutrinos exceeds the energy threshold for inverse beta decay on free hydrogen. This is the golden channel for liquid scintillator detectors because the coincidence signature allows for significant reductions in background signals. We find that the kiloton-scale liquid scintillator detector KamLAND can detect these pre-supernova neutrinos from a star with a mass of $25~M_{\sun}$ at a distance less than 690~pc with 3$\sigma$ significance before the supernova. This limit is dependent on the neutrino mass ordering and background levels. KamLAND takes data continuously and can provide a supernova alert to the community.

We developed an algorithm to find and characterize gravitationally lensed galaxies (arcs) to perform a comparison of the observed and simulated arc abundance. Observations are from the Cluster Lensing And Supernova survey with Hubble (CLASH). Simulated CLASH images are created using the MOKA package and also clusters selected from the high resolution, hydrodynamical simulations, MUSIC, over the same mass and redshift range as the CLASH sample. The algorithm' s arc elongation accuracy, completeness and false positive rate are determined and used to compute an estimate of the true arc abundance. We derive a lensing efficiency of $4 \pm 1$ arcs (with length $\ge 6"$ and length-to-width ratio $\ge 7$) per cluster for the X-ray selected CLASH sample, $4 \pm 1$ arcs per cluster for the MOKA simulated sample and $3 \pm 1$ arcs per cluster for the MUSIC simulated sample. The observed and simulated arc statistics are in full agreement. We measure the photometric redshifts of all detected arcs and find a median redshift $z_s = 1.9$ with 33% of the detected arcs having $z_s > 3$. We find that the arc abundance does not depend strongly on the source redshift distribution but is sensitive to the mass distribution of the dark matter halos (e.g. the $c-M$ relation). Our results show that consistency between the observed and simulated distributions of lensed arc sizes and axial ratios can be achieved by using cluster-lensing simulations that are carefully matched to the selection criteria used in the observations.

We report a measurement of the neutrino-electron elastic scattering rate of 862 keV 7Be solar neutrinos based on a 165.4 kton-day exposure of KamLAND. The observed rate is 582 +/- 90 (kton-day)^-1, which corresponds to a 862 keV 7Be solar neutrino flux of (3.26 +/- 0.50) x 10^9 cm^-2s^-1, assuming a pure electron flavor flux. Comparing this flux with the standard solar model prediction and further assuming three flavor mixing, a nu_e survival probability of 0.66 +/- 0.14 is determined from the KamLAND data. Utilizing a global three flavor oscillation analysis, we obtain a total 7Be solar neutrino flux of (5.82 +/- 0.98) x 10^9 cm^-2s^-1, which is consistent with the standard solar model predictions.

We search for electron anti-neutrinos ($\overline{\nu}_e$) from long and short-duration gamma-ray bursts~(GRBs) using data taken by the KamLAND detector from August 2002 to June 2013. No statistically significant excess over the background level is found. We place the tightest upper limits on $\overline{\nu}_e$ fluence from GRBs below 7 MeV and place first constraints on the relation between $\overline{\nu}_e$ luminosity and effective temperature.

We study the origin of non-thermal emissions from the Galactic black hole X-ray binary Cygnus X--1 which is a confirmed high mass microquasar. By analogy with those methods used in studies of active galactic nuclei, we propose a two-dimensional, time-dependent radiation model from the microquasar Cygnus X--1. In this model, the evolution equation for relativistic electrons in a conical jet are numerically solved by including escape, adiabatic and various radiative losses. The radiative processes involved are synchrotron emission, its self-Compton scattering, and inverse Compton scatterings of an accretion disk and its surrounding stellar companion. This model also includes an electromagnetic cascade process of an anisotropic $\gamma$-$\gamma$ interaction. We study the spectral properties of electron evolution and its emission spectral characteristic at different heights of the emission region located in the jet. We find that radio data from Cygnus X--1 are reproduced by the synchrotron emission, the \emph{Fermi} LAT measurements by the synchrotron emission and Comptonization of photons of the stellar companion, the TeV band emission fluxes by the Comptonization of the stellar photons. Our results show that: (1) Radio emission region extends from the binary system scales to the termination of the jet. (2) The GeV band emissions should originate from the distance close to the binary system scales. (3) The TeV band emissions could be inside the binary system, and these emissions could be probed by the upcoming CTA telescope. (4) The MeV tail emissions, which produces a strongly linearly polarized signal, are emitted inside the binary system. The location of the emissions is very close to the inner region of the jet.

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

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

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

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

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

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

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

The tight correlation between galaxy bulges and their central black hole masses likely emerges in a phase of rapid collapse and starburst at high redshift, due to the balance of gravity on gas with the feedback force from starbursts and the wind from the black hole; the average gravity on per unit mass of gas is ~ 2 x 10^-10 m/sec^2 during the star burst phase. This level of gravity could come from the real r^{-1} cusps of Cold Dark Matter (CDM) halos, but the predicted gravity would have a large scatter due to dependence on cosmological parameters and formation histories. Better agreement is found with the gravity from the scalar field in some co-variant versions of MOND, which can create the mirage of a Newtonian effective dark halo of density Pi r^{-1} near the center, where the characteristic surface density Pi=130alpha^{-1} Msun pc^{-2} and alpha is a fundamental constant of order unity fixed by the Lagrangian of the co-variant theory if neglecting environmental effects. We show with a toy analytical model and a hydrodynamical simulation that a constant background gravity due to MOND/TeVeS scalar field implies a critical pressure synchronizing starbursts and the formation of galaxy bulges and ellipticals. A universal threshold for the formation of the brightest regions of galaxies in a MONDian universe suggests that the central BHs, bulges and ellipticals would respect tight correlations like the M_{bulge}-M_{BH}-sigma relations. In general MOND tends to produce tight correlations in galaxy properties because its effective halo has less freedom and scatter than CDM halos.

Gravity is nearly a universal constant in the cusp of an NFW galaxy halo. Inside this external field an isothermal gas sphere will collapse and trigger a starburst if above a critical central pressure. Thus formed spheroidal stellar systems have Sersic-profile and satisfy the Faber-Jackson relation. The process is consistent with observed starbursts. We also recover the M_BH vs. velocity dispersion relation, if the gas collapse is regulated or resisted by the feedback from radiation from the central BH.

We present a "feedback compression" model to describe the galactic spheroid formation and its relation with the central nuclear activity. We suggest that the star formation itself can serve as the "positive feedback" in some extremely dense region to trigger the starburst. The star formation rate as well as the related stellar feedback-induced turbulence will be maximized under the regulation of the background dark halo's gravity. There is also stellar feedback acting inward to confine and obscure the central black hole (BH) till the BH grows sufficiently large to satisfy a balance condition between the accretion disk wind and the inward stellar feedback. The extremely vigorous star formation activity, the BH - bulge relations, the maximum velocity dispersion as well as the maximum BH mass are investigated based on such scenario, and are found to be consistent with observations.

We present a simple scenario where the formation of galactic bulges was regulated by the dark halo gravity and regulated the growth of the central supermassive black hole. Assuming the angular momentum is low, we suggest that bulges form in a runaway collapse due to the "gravothermal" instability once the central gas density or pressure exceeds certain threshold (Xu & Zhao 2007). We emphasize that the threshold is nearly universal, set by the background NFW dark matter gravity $g_{DM} \sim 1.2 \times 10^{-8}{\rm cm} {\rm sec}^{-2}$ in the central cusps of halos. Unlike known thresholds for gradual formation of galaxy disks, we show that the universal "halo-regulated" star formation threshold for spheroids matches the very high star formation rate and star formation efficiency shown in high-redshift observations of central starburst regions. The starburst feedback also builds up a pressure shortly after the collapse. This large pressure could both act outward to halt further infall of gas from larger scale, and act inward to counter the Compton-thick wind launched from the central black hole in an Eddington accretion. Assuming the feedback balancing inward and outward forces, our scenario naturally gives rise to the black hole-bulge relationships observed in the local universe.