Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

As a next-generation complex extensive air shower array with a large field of view, the large high altitude air shower observatory (LHAASO) is very sensitive to the very high energy gamma-rays from $\sim$ 300 GeV to 1 PeV, and may thus serve as an important probe for the heavy dark matter (DM) particles. In this study, we make a forecast for the LHAASO sensitivities to the gamma-ray signatures resulting from DM decay in dwarf spheroidal satellite galaxies (dSphs) within the LHAASO field of view. Both individual and combined limits for 19 dSphs incorporating the uncertainties of the DM density profile are explored. Owing to the large effective area and strong capability of the photon-proton discrimination, we find that LHASSSO is sensitive to the signatures from decaying DM particles above $\mathcal{O}(1)$ TeV. The LHAASO sensitivity to the DM decay lifetime reaches $\mathcal{O} (10^{26}) \sim \mathcal{O} (10^{28})$ s for several decay channels at the DM mass scale from 1 TeV to 100 TeV.

Conventional cosmic-ray propagation models usually assume an isotropic diffusion coefficient to account for the random deflection of cosmic rays by the turbulent interstellar magnetic field. Such a picture is very successful in explaining many observational phenomena related to the propagation of Galactic cosmic rays, such as broken power-law energy spectra, secondary-to-primary ratios, etc. However, the isotropic diffusion presupposition is facing severe challenges from recent observations. In particular, such observations on the large-scale anisotropy of TeV cosmic rays show that the dipole direction differs from the prediction of the conventional model. One possible reason is that the large-scale regular magnetic field, which leads to an anisotropic diffusion of cosmic rays, has not been included in the model provided by the public numerical packages. In this work, we propose two numerical schemes to solve the $3$-dimensional anisotropic transport equation: the pseudo source method and Hundsdorfer-Verwer scheme. Both methods are verified by reproducing the measured B/C and proton spectrum and the radial variation of spectral index expected by former 2D simulation. As a demonstration of the prediction capability, dipole anisotropy is also calculated by a toy simulation with a rough magnetic field.

We investigate the axion like particle (ALP)-photon oscillation effect in the high energy $\gamma$-ray spectra of PG 1553+113 and PKS 2155-304 measured by Fermi-LAT and H.E.S.S.. The choice of extragalactic background light (EBL) model, which induces the attenuate effect in observed $\gamma$-ray spectra, would affect the ALP implication. For the ordinary EBL model that prefers a null hypothesis, we set constraint on the ALP-photon coupling constant at 95% C.L. as $g_{a\gamma}\lesssim 5\times 10^{-11} ~\rm{GeV}^{-1}$ for the ALP mass $\sim 10$ neV. We also consider the CIBER observation of the cosmic infrared radiation, which shows an excess at the wave wavelength of $\sim 1~\mu$m after the substraction of foregrounds. The high energy gamma-rays from extragalactic sources at high redshifts would suffer from a more significant attenuate effect caused by this excess. In this case, we find that the ALP-photon oscillation would improve the fit to the observed spectra of PKS 2155-304 and PG 1553+113 and find a favored parameter region at 95% C.L..

The deuterium fractionation in starless cores gives us a clue to estimate their lifetime scales, thus allowing us to distinguish between different dynamical theories of core formation. Cores also seem to be subject to a differential N2 and CO depletion which was not expected from models. We aim to make a survey of 10 cores to estimate their lifetime scales and depletion profiles in detail. After L183, in Serpens, we present the second cloud of the series, L1512 in Auriga. To constrain the lifetime scale, we perform chemical modeling of the deuteration profiles across L1512 based on dust extinction measurements from near-infrared observations and non-local thermal equilibrium radiative transfer with multiple line observations of N2H+, N2D+, DCO+, C18O, and 13CO, plus H2D+ (1$_{10}$--1$_{11}$). We find a peak density of 1.1$\times$10$^5$ cm$^{-3}$ and a central temperature of 7.5$\pm$1 K, which are respectively higher and lower compared with previous dust emission studies. The depletion factors of N2H+ and N2D+ are 27$^{+17}_{-13}$ and 4$^{+2}_{-1}$ in L1512, intermediate between the two other more advanced and denser starless core cases, L183 and L1544. These factors also indicate a similar freeze-out of N2 in L1512, compared to the two others despite a peak density one to two orders of magnitude lower. Retrieving CO and N2 abundance profiles with the chemical model, we find that CO has a depletion factor of $\sim$430-870 and the N2 profile is similar to that of CO unlike towards L183. Therefore, L1512 has probably been living long enough so that N2 chemistry has reached steady state. N2H+ modeling remains compulsory to assess the precise physical conditions in the center of cold starless cores, rather than dust emission. L1512 is presumably older than 1.4 Myr. Therefore, the dominating core formation mechanism should be ambipolar diffusion for this source.

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

The Large High Altitude Air Shower Observatory (LHAASO) is a next-generation observatory for high energy gamma rays and cosmic rays with wide field of view. It will detect gamma rays with high sensitivity in the energy range from 300 GeV to 1 PeV. Therefore, it is promising for LHAASO to search for the high-energy gamma rays induced by dark matter (DM) self-annihilation in dwarf spheroidal satellite galaxies (dSphs), which are ideal objects for the DM indirect detection. In this work, we investigate the LHAASO sensitivity to DM self-annihilation signatures for 19 dSphs and take the uncertainties on the spatial DM distribution of dSphs into account. We perform a joint likelihood analysis for the 19 dSphs and find that the LHAASO sensitivity to the DM annihilation cross section will reach $\mathcal{O}(10^{-24})\sim \mathcal{O}(10^{-25})$ cm$^3$ s$^{-1}$ at the mass scale above TeV for several annihilation modes, which is larger than the canonical thermal relic cross section by a factor of 10 to 100.

We present the results of mapping observations toward an outflow-shocked region, OMC-2 FIR 4 using the Nobeyama 45-m telescope. We observed the area in $^{13}$CO ($J=1-0$), C$^{18}$O ($J=1-0$), N$_2$H$^+$ ($J=1-0$), CCS ($J_N=8_7-7_6$), HCO$^+$ ($J=1-0$), H$^{13}$CO$^+$ ($J=1-0$), HN$^{13}$C ($J=1-0$), H$^{13}$CN ($J=1-0$), DNC ($J=1-0$), N$_2$D$^+$ ($J=1-0$), and DC$_3$N ($J=9-8$). We detected a dense molecular clump that contains FIR 4/5. We also detected in $^{13}$CO blueshifted and redshifted components driven presumably by protostellar outflows in this region. The axes of the FIR 3 and VLA 13 outflows, projected on the plane of the sky, appear to point toward the FIR 4 clump, suggesting that the clump may be compressed by protostellar outflows from Class I sources, FIR 3 and VLA 13. Applying the hyperfine fit of N$_2$H$^+$ lines, we estimated the excitation temperature to be $\sim$ 20 K. The high excitation temperature is consistent with the fact that the clump contains protostars. The CCS emission was detected in this region for the first time. Its abundance is estimated to be a few $\times 10^{-12}$, indicating that the region is chemically evolved at $\sim 10^5$ years, which is comparable to the typical lifetime of the Class I protostars. This timescale is consistent with the scenario that star formation in FIR 4 is triggered by dynamical compression of the protostellar outflows. The [HNC]/[HCN] ratio is evaluated to be $\sim 0.5$ in the dense clump and the outflow lobes, whereas it is somewhat larger in the envelope of the dense clump. The small [HNC]/[HCN] ratio indicates that the HNC formation was prevented due to high temperatures. Such high temperatures seem to be consistent with the scenario that either protostellar radiation or outflow compression, or both, affected the thermal properties of this region.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

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

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

We investigate the implications on the dark matter (DM) signal from the AMS-02 cosmic antiproton flux. Global fits to the data are performed under different propagation and hadronic interaction models. The uncertainties from the injection spectrum, propagation effects and solar modulation of the cosmic rays are taken into account comprehensively. Since we need to investigate extended parameter regions with multiple free parameters in the fit, the machine learning method is adopted to maintain a realistic time cost. We find all the effects considered in the fitting process interplay with each other, among which the hadronic interaction model is the most important factor affecting the result. In most hadronic interaction and CR propagation models no DM signal is found with significance larger than $2\sigma$ except that the EPOS-LHC interaction model requires a more than $3\sigma$ DM signal with DM mass around $1\,\mathrm{TeV}$. For the diffusive reacceleration propagation model there is a highly significant DM signal with mass around $100\,\mathrm{GeV}$. However, the signal becomes less than $1\sigma$ if we take a charge dependent solar modulation potential in the analysis.

A study on cosmogenic activation in germanium was carried out to evaluate the cosmogenic background level of natural and $^{70}$Ge depleted germanium detectors. The production rates of long-lived radionuclides were calculated with Geant4 and CRY. Results were validated by comparing the simulated and experimental spectra of CDEX-1B detector. Based on the validated codes, the cosmogenic background level was predicted for further tonne-scale CDEX experiment. The suppression of cosmogenic background level could be achieved by underground germanium crystal growth and high-purity germanium detector fabrication to reach the sensitivity requirement for direct detection of dark matter. With the low cosmogenic background, new physics channels, such as solar neutrino research and neutrinoless double-beta decay experiments, were opened and the corresponding simulations and evaluations were carried out.

In this work, we propose a new source for gravitational wave (GW) radiation associated with the quantum chromodynamics (QCD) phase transition in the inner cores of neutron stars. The mechanism is based on the bubble dynamics during the first-order phase transition from nuclear matter to quark matter. We identify the characteristic frequency to be of order $\omega_c\sim 10^6~{\rm rad/s}$ for this kind of sources and the strain magnitude ($h\sim 10^{-24}$ for a neutron star at a distance of $0.1~{\rm Mpc}$) reachable by future GW detectors. The GW spectra are shown to be useful to check the transition nature at high baryon chemical potential as well as to constrain the radius and density of the inner cores, which are still indistinct up to now.

Wave dark matter ($\psi$DM), which satisfies the Schr\"odinger-Poisson equation, has recently attracted substantial attention as a possible dark matter candidate. Numerical simulations have in the past provided a powerful tool to explore this new territory of possibility. Despite their successes to reveal several key features of $\psi$DM, further progress in simulations is limited, in that cosmological simulations so far can only address formation of halos below $\sim 2\times 10^{11} M_\odot$ and substantially more massive halos have become computationally very challenging to obtain. For this reason, the present work adopts a different approach in assessing massive halos by constructing wave-halo solutions directly from the wave distribution function. This approach bears certain similarity with the analytical construction of particle-halo (cold dark matter model). Instead of many collisionless particles, one deals with one single wave that has many non-interacting eigenstates. The key ingredient in the wave-halo construction is the distribution function of the wave power, and we use several halos produced by structure formation simulations as templates to determine the wave distribution function. Among different models, we find the fermionic King model presents the best fits and we use it for our wave-halo construction. We have devised an iteration method for constructing the nonlinear halo, and demonstrate its stability by three-dimensional simulations. A Milky-Way-sized halo has also been constructed, and the inner halo is found flatter than the NFW profile. These wave-halos have small-scale interferences both in space and time producing time-dependent granules. While the spatial scale of granules varies little, the correlation time is found to increase with radius by one order of magnitude across the halo.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

The cosmic antideuteron is a promising probe for the dark matter annihilation signature. In order to determine the DM signature, the background astrophysical antideuteron flux should be carefully studied. In this work we provide a new calculation of the secondary antideuteron flux, and pay special attention to the uncertainties from hadronic interaction models by using several Monte Carlo generators. The uncertainties from propagation effects are also carefully investigated for both the astrophysical background and DM annihilation signature in several scenarios, which are constrained by the latest B/C ratio measured by AMS-02. Considering these uncertainties, we find that the secondary antideuteron flux is hard to detect in the near future detectors. However, the antideuteron signature from dark matter annihilation will be detectable even considering the constraint from the AMS-02 observation of the $\bar{p}/p$ ratio.

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

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

As one of the most popular applications of lidar systems, the atmospheric visibility is defined to be inversely proportional to the atmospheric extinction coefficient at 0.55 um. Since the laser at 1.5 um shows the highest maximum permissible exposure in the wavelength ranging from 0.3 um to 10 um, the eye-safe 1.5 um lidar can be deployed in urban areas. In such a case, the measured extinction coefficient at 1.5 um should be converted to that at 0.55 um for visibility retrieval. Although several models have been established since 1962, the accurate wavelength conversion remains a challenge. An adaptive inversion algorithm for 1.5 um visibility lidar is proposed and demonstrated by using the in situ Angstrom wavelength exponent, which is derived from either a sun photometer or an aerosol spectrometer. The impact of the particle size distribution of atmospheric aerosols and the Rayleigh backscattering of atmospheric molecules are taken into account. In comparison, the Angstrom wavelength exponent derived from the sun photometer is 7.7% higher than that derived from the aerosol spectrometer. Then, using 1.5 um visibility lidar, the visibility with a temporal resolution of 5 min is detected over 48 hours in Hefei (31.83 N, 117.25 E). The average visibility error between the new method and a visibility sensor (Vaisala, PWD52) is 5.2% with the R-square value of 0.96, while the relative error between another reference visibility lidar at 532 nm and the visibility sensor is 6.7% with the R-square value of 0.91. All results agree with each other well, demonstrating the accuracy and stability of the algorithm.

The unprecedented depth and area surveyed by the Subaru Strategic Program with the Hyper Suprime-Cam (HSC-SSP) have enabled us to construct and publish the largest distant cluster sample out to $z\sim 1$ to date. In this exploratory study of cluster galaxy evolution from $z=1$ to $z=0.3$, we investigate the stellar mass assembly history of brightest cluster galaxies (BCGs), evolution of stellar mass and luminosity distributions, stellar mass surface density profile, as well as the population of radio galaxies. Our analysis is the first high redshift application of the top N richest cluster selection, which is shown to allow us to trace the cluster galaxy evolution faithfully. Over the 230 deg$^2$ area of the current HSC-SSP footprint, selecting the top 100 clusters in each of the 4 redshift bins allows us to observe the buildup of galaxy population in descendants of clusters whose $z\approx 1$ mass is about $2\times 10^{14}\,M_\odot$. Our stellar mass is derived from a machine-learning algorithm, which is found to be unbiased and accurate with respect to the COSMOS data. We find very mild stellar mass growth in BCGs (about 35% between $z=1$ and 0.3), and no evidence for evolution in both the total stellar mass-cluster mass correlation and the shape of the stellar mass surface density profile. We also present the first measurement of the radio luminosity distribution in clusters out to $z\sim 1$, and show hints of changes in the dominant accretion mode powering the cluster radio galaxies at $z\sim 0.8$.

We present clustering properties from 579,492 Lyman break galaxies (LBGs) at z~4-6 over the 100 deg^2 sky (corresponding to a 1.4 Gpc^3 volume) identified in early data of the Hyper Suprime-Cam (HSC) Subaru strategic program survey. We derive angular correlation functions (ACFs) of the HSC LBGs with unprecedentedly high statistical accuracies at z~4-6, and compare them with the halo occupation distribution (HOD) models. We clearly identify significant ACF excesses in 10"<$\theta$<90", the transition scale between 1- and 2-halo terms, suggestive of the existence of the non-linear halo bias effect. Combining the HOD models and previous clustering measurements of faint LBGs at z~4-7, we investigate dark-matter halo mass (Mh) of the z~4-7 LBGs and its correlation with various physical properties including the star-formation rate (SFR), the stellar-to-halo mass ratio (SHMR), and the dark matter accretion rate (dotMh) over a wide-mass range of Mh/M$_\odot$=4x10^10-4x10^12. We find that the SHMR increases from z~4 to 7 by a factor of ~4 at Mh~1x10^11 M$_\odot$, while the SHMR shows no strong evolution in the similar redshift range at Mh~1x10^12 M$_\odot$. Interestingly, we identify a tight relation of SFR/dotMh-Mh showing no significant evolution beyond 0.15 dex in this wide-mass range over z~4-7. This weak evolution suggests that the SFR/dotMh-Mh relation is a fundamental relation in high-redshift galaxy formation whose star formation activities are regulated by the dark matter mass assembly. Assuming this fundamental relation, we calculate the cosmic SFR densities (SFRDs) over z=0-10 (a.k.a. Madau-Lilly plot). The cosmic SFRD evolution based on the fundamental relation agrees with the one obtained by observations, suggesting that the cosmic SFRD increase from z~10 to 4-2 (decrease from z~4-2 to 0) is mainly driven by the increase of the halo abundance (the decrease of the accretion rate).

The cosmic-ray (CR) $e^\pm$ excess observed by AMS-02 can be explained by dark matter (DM) annihilation. However, the DM explanation requires a large annihilation cross section which is strongly disfavored by other observations, such as the Fermi-LAT gamma-ray observation of dwarf galaxies and the Planck observation of the cosmic microwave background (CMB). Moreover, the DM annihilation cross section required by the CR $e^\pm$ excess is also too large to generate the correct DM relic density with thermal production. In this work we use the Breit-Wigner mechanism with a velocity dependent DM annihilation cross section to reconcile these tensions. If DM particles accounting for the CR $e^\pm$ excess with $v\sim \mathcal{O}(10^{-3})$ are very close to a resonance in the physical pole case, their annihilation cross section in the Galaxy reaches a maximal value. On the other hand, the annihilation cross section would be suppressed for DM particles with smaller relative velocities in dwarf galaxies and at recombination, which may affect the gamma-ray and CMB observations, respectively. We find a proper parameter region that can simultaneously explain the AMS-02 results and the thermal relic density, while satisfying the Fermi-LAT and Planck constraints.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

In this work we use the newly reported Boron-to-Carbon ratio (B/C) from AMS-02 and the time-dependent proton fluxes from PAMELA and AMS-02 to constrain the source and propagation parameters of cosmic rays in the Milky Way. A linear correlation of the solar modulation parameter with solar activities is assumed to account for the time-varying cosmic ray fluxes. A comprehensive set of propagation models, with/without reacceleration or convection, have been discussed and compared. We find that only the models with reacceleration can self-consistently fit both the proton and B/C data. The rigidity dependence slope of the diffusion coefficient, $\delta$, is found to be about $0.38-0.50$ for the diffusion-reacceleration models. The plain diffusion and diffusion-convection models fit the data poorly. We compare different model predictions of the positron and antiproton fluxes with the data. We find that the diffusion-reacceleration models over-produce low energy positrons, while non-reacceleration models give better fit to the data. As for antiprotons, reacceleration models tend to under-predict low energy antiproton fluxes, unless a phenomenological modification of the velocity-dependence of the diffusion coefficient is applied. Our results suggest that there could be important differences of the propagation for nuclei and leptons, in either the Milky Way or the solar heliosphere.

In the holographic description of cosmic preheating proposed in an accompanied Letter [1], the energy transfer between the inflaton and matter field at strong coupling is suggested to be mimicked by superfluid and normal components of a superconductor on Friedmann-Robertson-Walker (FRW) boundary in an asymptotically Anti-de Sitter (AdS) spacetime. In this paper we investigated two aspects of the scenario of holographic preheating that are not included in the accompanied work. Firstly, we study in detail the evolution of the quasi-normal modes (QNMs) surrounding a metastable hairy black hole. This analysis can quantitatively describe the preheating process of the matter field that is produced continuously in the case of strong coupling. Secondly, we present a detailed analysis of the holographic renormalization for the AdS-FRW background, which allows us to extract operator expectation values for studying cosmological implications.

We propose a holographic description of cosmic preheating at strong coupling. In this scenario the energy transfer between the inflaton and matter field is mimicked by a model of holographic superconductor. An exponential amplification of the matter field during preheating can be described by the quasi-normal modes of a metastable "black hole" in the bulk spacetime with an expanding boundary. Our results reveal that the matter field can be produced continuously at strong coupling in contrast to the case of weak coupling with a discontinuous matter growth as inflaton oscillates. Furthermore, the amplification of matter field has an enhanced dependence on the vacuum expectation value of the inflaton at strong coupling. By virtue of the proposed mechanism, physics of the very early universe at an extremely high temperature right after inflation may become accessible.

Recently the AMS-02 collaboration has published the measurement of the cosmic antiproton to proton ratio $\bar{p}/p$ and the $\bar{p}$ flux with a high precision up to $\sim 450\,\mathrm{GeV}$. In this work, we perform a systematic analysis of the secondary antiproton flux generated by the cosmic ray interaction with the interstellar gas. The uncertainty of the prediction originates from the cosmic ray propagation process and the hadronic interaction models. Although the cosmic ray propagation parameters have been well controlled by the AMS-02 $B/C$ ratio data for a specified model, different propagation models can not be discriminated by the $B/C$ data. The $\bar{p}$ flux is also calculated for several hadronic interaction models, which are generally adopted by the cosmic ray community. However, the results for different hadronic models do not converge. We find the EPOS LHC model, which seems to fit the collider data very well, predicts a slightly lower $\bar{p}/p$ ratio than the AMS-02 data at the high energy end. Finally we derive the constraints on the dark matter annihilation cross section from the AMS-02 $\bar{p}/p$ ratio for different propagation and hadronic interaction models.

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

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

We present an alternative Corrector-ADC design for GMT. The design consists of just 3 silica lenses, of maximum size 1.51m, and includes only a single low-precision asphere for 20' field-of-view, and none for 10'. The polychromatic (360nm-1300nm) image quality is d80<0.043" at zenith and d80<0.20" for ZD<60 degrees. The monochromatic image quality is d80<0.1" everywhere, and typically ~0.05". The ADC action is achieved by tilt and translation of all three lenses; L1 and L2 via simple slide mechanisms each using a single encoded actuator, and L3 via a novel 'tracker-ball' support and three actuators. There is also a small motion of M2 via the hexapod, automatically generated by the AGWS system. The ADC action causes a small non-telecentricity, but this is much less than the unavoidable chromatic effects shared with the baseline design. The ADC action also changes the distortion pattern of the telescope, but this can be used positively, to reduce the maximum image motion due to differential refraction by a factor of three. The transmission is superb at all wavelengths, because of the reduced number of air/glass surfaces, and the use only of fused silica.

L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

In this work we study the constraints on dark matter (DM) annihilation/decay from the Fermi-LAT Isotropic Gamma-Ray Background (IGRB) observation. We consider the contributions from both extragalactic and galactic DM components. For DM annihilation, the evolutions of extragalactic DM halos are taken into account. We find that the IGRB constraints under some DM subhalo models can be comparable to those derived from the observations of dwarf spheroidal galaxies. We also use the IGRB results to constrain the parameter regions accounting for the latest AMS-02 electron-positron anomaly. We find that the majority of DM annihilation/decay channels are strongly disfavored by the latest Fermi-LAT IGRB observation; only DM annihilation/decay to $\mu^+\mu^-$ may be valid.

We investigate the observational constraints on the modified gravity, which combines the $R^{2-q}$ inflation with the power-law (exponential) type of the viable $f(R)$ dark energy models. We discuss the difference between the combined model and $R^{2-q}$ gravity in the inflationary epoch and obtain the constraints on the deviation power $q$ as well as the parameters in $f(R)$ by using the CosmoMC package. The allowed ranges of the spectral index and tensor-to-scalar ratio from the Planck data are highly restricted, resulting in $q < 2.66 \times 10^{-2}$ and $2.17 \times 10^{-2}$ for the power-law and exponential types of $f(R)$ gravity, respectively.

The AMS-02 collaboration has just released the cosmic antiproton to proton ratio $\bar{p}/p$ with a high precision up to $\sim 450$ GeV. In this work, we calculate the secondary antiprotons generated by cosmic ray interactions with the interstellar medium taking into account the uncertainties from the cosmic ray propagation. The $\bar{p}/p$ ratio predicted by these processes shows some tension with the AMS-02 data in some regions of propagation parameters, but the excess is not significant. We then try to derive upper bounds on the dark matter annihilation cross section from the $\bar{p}/p$ data or signal regions favored by the data. It is shown that the constraint derived by the AMS-02 data is similar to that from Fermi-LAT observations of dwarf galaxies. The signal region for dark matter is usually required $m_\chi \sim O(10)$ TeV and $\left<\sigma v\right>\sim\mathcal{O}(10^{-23})~\cm^3~\sec^{-1}$.

The AMS-02 has just published the unprecedentedly precise measurement of the cosmic electron and positron spectra. In this paper we try to give a quantitative study on the AMS-02 results by a global fitting to the electron and positron spectra, together with the updated positron fraction data. The Markov Chain Monte Carlo algorithm is adopted to do the fitting. The primary electron spectrum and the parameters for pulsars or dark matter which contribute extra positrons are determined simultaneously. We find that there is a hardening of the primary electron spectrum at $\sim 60$ GeV. With such a new feature at the background spectrum both the pulsars and dark matter can explain the AMS-02 results very well. The dark matter scenario shows a drop at positron fraction at $\sim 300$ GeV, however, suffers very strong constraints from Fermi $\gamma$-ray observations. The fitting results also suggest that the propagation model with convection may be more favored by the lepton data than the reacceleration model.

The inorganic CsI(Tl) crystal scintillator is a candidate served as an anti-compton detector for China Dark matter EXperiment (CDEX). Studying the intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The timing, energy and spatial correlations, as well as the capability of pulse shape discrimination provide powerful methods for the measurement of intrinsic radiopurities. The experimental design, detector performance and event-selection algorithms are described. A total of 359*3 kg-days data from three prototypes of CsI(Tl) crystals were taken at China Jinping underground laboratory where offers a good shielding environment. The contamination levels of internal isotopes from Cs-137, Th-232 and U-238 series, as well as the upper bounds of U-235 series are reported. Identification of the whole alpha8 peaks from U/Th decay chains and derivation of those corresponding quenching factors are achieved.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

The AMS-02 collaboration has just released its first result of the cosmic positron fraction $e^+/(e^-+e^+)$ with high precision up to $\sim 350$ GeV. The AMS-02 result shows the same trend with the previous PAMELA result, which requires extra electron/positron sources on top of the conventional cosmic ray background, either from astrophysical sources or from dark matter annihilation/decay. In this paper we try to figure out the nature of the extra sources by fitting to the AMS-02 $e^+/(e^-+e^+)$ data, as well as the electron and proton spectra by PAMELA and the $(e^-+e^+)$ spectrum by Fermi and HESS. We adopt the GALPROP package to calculate the propagation of the Galactic cosmic rays and the Markov Chain Monte Carlo sampler to do the fit. We find that the AMS-02 data have implied essential difference from the PAMELA data. There is {\rm tension} between the AMS-02 $e^+/(e^-+e^+)$ data and the Fermi/HESS $(e^-+e^+)$ spectrum, that the AMS-02 data requires less contribution from the extra sources than Fermi/HESS. Then we redo the fit without including the Fermi/HESS data. In this case both the pulsars and dark matter annihilation/decay can explain the AMS-02 data. The pulsar scenario has a soft inject spectrum with the power-law index $\sim 2$, while the dark matter scenario needs $\tau^+\tau^-$ final state with mass $\sim 600$ GeV and a boost factor $\sim 200$.

We have built and commissioned gas absorption cells for precision spectroscopic radial velocity measurements in the near-infrared in the H and K bands. We describe the construction and installation of three such cells filled with 13CH4, 12CH3D, and 14NH3 for the CSHELL spectrograph at the NASA Infrared Telescope Facility (IRTF). We have obtained their high-resolution laboratory Fourier Transform spectra, which can have other practical uses. We summarize the practical details involved in the construction of the three cells, and the thermal and mechanical control. In all cases, the construction of the cells is very affordable. We are carrying out a pilot survey with the 13CH4 methane gas cell on the CSHELL spectrograph at the IRTF to detect exoplanets around low mass and young stars. We discuss the current status of our survey, with the aim of photon-noise limited radial velocity precision. For adequately bright targets, we are able to probe a noise floor of 7 m/s with the gas cell with CSHELL at cassegrain focus. Our results demonstrate the feasibility of using a gas cell on the next generation of near-infrared spectrographs such as iSHELL on IRTF, iGRINS, and an upgraded NIRSPEC at Keck.

We have built and commissioned a prototype agitated non-circular core fiber scrambler for precision spectroscopic radial velocity measurements in the near-infrared H band. We have collected the first on-sky performance and modal noise tests of these novel fibers in the near-infrared at H and K bands using the CSHELL spectrograph at the NASA InfraRed Telescope Facility (IRTF). We discuss the design behind our novel reverse injection of a red laser for co-alignment of star-light with the fiber tip via a corner cube and visible camera. We summarize the practical details involved in the construction of the fiber scrambler, and the mechanical agitation of the fiber at the telescope. We present radial velocity measurements of a bright standard star taken with and without the fiber scrambler to quantify the relative improvement in the obtainable blaze function stability, the line spread function stability, and the resulting radial velocity precision. We assess the feasibility of applying this illumination stabilization technique to the next generation of near-infrared spectrographs such as iSHELL on IRTF and an upgraded NIRSPEC at Keck. Our results may also be applied in the visible for smaller core diameter fibers where fiber modal noise is a significant factor, such as behind an adaptive optics system or on a small < 1 meter class telescope such as is being pursued by the MINERVA and LCOGT collaborations.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.

We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

An energy threshold of (220$\pm$10) eV was achieved at an efficiency of 50% with a four-channel ultra-low-energy germanium detector each with an active mass of 5 g\cite{wimppaper}. This provides a unique probe to WIMP dark matter with mass below 10 GeV. With low background data taken at the Kuo-Sheng Laboratory, constraints on WIMPs in the galactic halo were derived. Both spin-independent WIMP-nucleon and spin-dependent WIMP-neutron bounds improve over previous results for WIMP mass between 3$-$6 GeV. These results, together with those on spin-dependent couplings, will be presented. Sensitivities for full-scale experiments were projected. This detector technique makes the unexplored sub-keV energy window accessible for new neutrino and dark matter experiments.

The circumstellar dust disk of the Herbig Ae star AB Aur has been found to exhibit complex spiral-like structures in the near-IR image obtained with the Subaru Telescope. We present maps of the disk in both 12CO (3-2) and dust continuum at 345 GHz with the SMA at an angular resolution of 1.0"x0.7" (144AU x 100AU). The continuum emission traces a dust disk with a central depression and a maximum overall dimension of 450AU (FWHM). This dust disk exhibits several distinct peaks that appear to coincide with bright features in the near-IR image, in particular the brightest inner spiral arm. The CO emission traces a rotating gas disk of size 530AU x 330AU with a deprojected maximum velocity of 2.8km/s at 450AU. In contrast to the dust disk, the gas disk exhibits an intensity peak at the stellar position. Furthermore, the CO emission in several velocity channels traces the innermost spiral arm seen in the near-IR. We compare the observed spatial-kinematic structure of the CO emission to a simple model of a disk in Keplerian rotation, and find that only the emission tracing the main spiral arm clearly lies outside the confines of our model. This emission has a net outward radial motion compared with the radial velocity predicted by the model at the location of the main spiral arms. The disk of AB Aur is therefore quite different from the Keplerian disks seen around many Herbig Ae stars. The spiral-like structures of the disk with non-Keplerian motions we revealed in 12CO (3-2), together with the central depression of the dust disk, may be explained to be driven by the possible existence of a giant planet forming in the disk.

A new scalar-tensor theory of gravity induced by dynamically broken scale invariance is proposed, and its cosmological implications are discussed. It is found that the model admits an inflation via the Hawking-Moss bubbling, but the inflation rate remains undetermined due to the strong gravity limit. In light of this, scale-invariant metric perturbations having a dominant tensor component can be generated without slow-rollover. In addition, the deviation from the standard hot big-bang is vanishingly small after inflation.