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.

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.

We use the FIRE-2 cosmological simulations to study the formation of a virial temperature, quasi-static gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that virialization of the inner CGM coincides with abrupt changes in the properties of the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results hence suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed.

We study the escape fraction of ionizing photons (f_esc) in two cosmological zoom-in simulations of galaxies in the reionization era with halo mass M_halo~10^10 and 10^11 M_sun (stellar mass M*~10^7 and 10^9 M_sun) at z=5 from the Feedback in Realistic Environments project. These simulations explicitly resolve the formation of proto-globular clusters (GCs) self-consistently, where 17-39% of stars form in bound clusters during starbursts. Using post-processing Monte Carlo radiative transfer calculations of ionizing radiation, we compute f_esc from cluster stars and non-cluster stars formed during a starburst over ~100 Myr in each galaxy. We find that the averaged f_esc over the lifetime of a star particle follows a similar distribution for cluster stars and non-cluster stars. Clusters tend to have low f_esc in the first few Myrs, presumably because they form preferentially in more extreme environments with high optical depths; the f_esc increases later as feedback starts to disrupt the natal cloud. On the other hand, non-cluster stars formed between cluster complexes or in the compressed shell at the front of a superbubble can also have high f_esc. We find that cluster stars on average have comparable f_esc to non-cluster stars. This result is robust across several star formation models in our simulations. Our results suggest that the fraction of ionizing photons from proto-GCs to cosmic reionization is comparable to the cluster formation efficiency in high-redshift galaxies and hence proto-GCs likely contribute an appreciable fraction of photons but are not the dominant sources for reionization.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

In this work, we achieve the determination of the cosmic curvature $\Omega_K$ in a cosmological model-independent way, by using the Hubble parameter measurements $H(z)$ and type Ia supernovae (SNe Ia). In our analysis, two nonlinear interpolating tools are used to reconstruct the Hubble parameter, one is the Artificial Neural Network (ANN) method, and the other is the Gaussian process (GP) method. We find that $\Omega_K$ based on the GP method can be greatly influenced by the prior of $H_0$, while the ANN method can overcome this. Therefore, the ANN method may have more advantages than GP in the measurement of the cosmic curvature. Based on the ANN method, we find a spatially open universe is preferred by the current $H(z)$ and SNe Ia data, and the difference between our result and the value inferred from Planck CMB is $1.6\sigma$. In order to test the reliability of the ANN method, and the potentiality of the future gravitational waves (GW) standard sirens in the measurement of the cosmic curvature, we constrain $\Omega_K$ using the simulated Hubble parameter and GW standard sirens in a model-independent way. We find that the ANN method is reliable and unbiased, and the error of $\Omega_K$ is $\sim0.186$ when 100 GW events with electromagnetic counterparts are detected, which is $\sim56\%$ smaller than that constrained from the Pantheon SNe Ia. Therefore, the data-driven method based on ANN has potential in the measurement of the cosmic curvature.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

In this work, we propose a new nonparametric approach for reconstructing a function from observational data using an Artificial Neural Network (ANN), which has no assumptions about the data and is a completely data-driven approach. We test the ANN method by reconstructing functions of the Hubble parameter measurements $H(z)$ and the distance-redshift relation $D_L(z)$ of Type Ia supernovae. We find that both $H(z)$ and $D_L(z)$ can be reconstructed with high accuracy. Furthermore, we estimate cosmological parameters using the reconstructed functions of $H(z)$ and $D_L(z)$ and find the results are consistent with those obtained using the observational data directly. Therefore, we propose that the function reconstructed by ANN can represent the actual distribution of observational data and can be used for parameter estimation in further cosmological research. In addition, we present a new strategy for training and evaluating the neural network, and a code for reconstructing functions using ANN has been developed and is available for download.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We present a systematic analysis of the phase lags associated with the type-C QPOs in GRS 1915+105 using RXTE data. Our sample comprises of 620 RXTE observations with type-C QPOs ranging from ~0.4 Hz to ~6.3 Hz. Based on our analysis, we confirm that the QPO phase lags decrease with QPO frequency, and change sign from positive to negative at a QPO frequency of ~2 Hz. In addition, we find that the slope of this relation is significantly different between QPOs below and above 2 Hz. The relation between the QPO lags and QPO rms can be well fitted with a broken line: as the QPO lags go from negative to positive, the QPO rms first increases, reaching its maximum at around zero lag, and then decreases. The phase-lag behaviour of the subharmonic of the QPO is similar to that of the QPO fundamental, where the subharmonic lags decrease with subharmonic frequency and change sign from positive to negative at a subharmonic frequency of ~1 Hz; on the contrary, the second harmonic of the QPO shows a quite different phase-lag behaviour, where all the second harmonics show hard lags that remain more or less constant. For both the QPO and its (sub)harmonics, the slope of the lag-energy spectra shows a similar evolution with frequency as the average phase lags. This suggests that the lag-energy spectra drives the average phase lags. We discuss the possibility for the change in lag sign, and the physical origin of the QPO lags.

We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

We present a concrete realization of the sound speed resonance mechanism for primordial black hole formation within a specific model of Dirac-Born-Infeld (DBI) inflation. We perform a perturbative approach to phenomenologically construct such a viable DBI inflation model that involves the non-oscillating stage and the oscillating stage, with a type of specific forms of the warp factor and the potential. We show that the continuous but non-smooth conjunction of sound speed between two stages does not yield manifest effects on the phenomenology of SSR, and thus, our model gives rise to the same PBH mass spectrum as the original predictions of SSR. Making use of observational data, we derive various cosmological constraints on the parameter space. Our analyses show that, the predicted tensor-to-scalar ratio is typically small, while the amplitude of primordial non-Gaussianity can meet with CMB bounds, and additionally, the consistency relation for single-field slow-roll inflation is softly violated in our case due to the small sound speed variations.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

We use Hubble Space Telescope (HST) imaging and near-infrared spectroscopy from Keck/MOSFIRE to study the sub-structure around the progenitor of a Milky Way-mass galaxy in the Hubble Frontier Fields (HFF). Specifically, we study an $r_e = 40^{+70}_{-30}$pc, $M_{\star} \sim 10^{8.2} M_{\odot}$ rest-frame ultra-violet luminous "clump" at a projected distance of $\sim$100~pc from a $M_{\star} \sim 10^{9.8}$M$_{\odot}$ galaxy at $z = 2.36$ with a magnification $\mu = 5.21$. We measure the star formation history of the clump and galaxy by jointly modeling the broadband spectral energy distribution from HST photometry and H$\alpha$ from MOSFIRE spectroscopy. Given our inferred properties (e.g., mass, metallicity, dust) of the clump and galaxy, we explore scenarios in which the clump formed \emph{in-situ} (e.g., a star forming complex) or \emph{ex-situ} (e.g., a dwarf galaxy being accreted). If it formed \emph{in-situ}, we conclude that the clump is likely a single entity as opposed to a aggregation of smaller star clusters, making it one of the most dense star clusters cataloged. If it formed \emph{ex-situ}, then we are witnessing an accretion event with a 1:40 stellar mass ratio. However, our data alone are not informative enough to distinguish between \emph{in-situ} and \emph{ex-situ} scenarios to a high level of significance. We posit that the addition of high-fidelity metallicity information, such as [OIII]4363\AA, which can be detected at modest S/N with only a few hours of JWST/NIRSpec time, may be a powerful discriminant. We suggest that studying larger samples of moderately lensed sub-structures across cosmic time can provide unique insight into the hierarchical formation of galaxies like the Milky Way.

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically-thin) and moments-based M1 (exact when optically-thick). The most important RFB channels on galaxy scales are photo-ionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (~1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyper-dense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyper-clustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, ~10% of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with A_v~1.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

Several of the many proposed Dark Matter candidate particles, already investigated with lower exposure and a higher software energy threshold, are further analyzed including the first DAMA/LIBRA--phase2 data release, with an exposure of 1.13 ton $\times$ yr and a lower software energy threshold (1 keV). The cumulative exposure above 2 keV considering also DAMA/NaI and DAMA/LIBRA--phase1 results is now 2.46 ton $\times$ yr. The analysis permits to constraint the parameters' space of the considered candidates restricting their values -- with respect to previous analyses -- thanks to the increase of the exposure and to the lower energy threshold.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

We use cosmological hydrodynamical simulations of Milky-Way-mass galaxies from the FIRE project to evaluate various strategies for estimating the mass of a galaxy's stellar halo from deep, integrated-light images. We find good agreement with integrated-light observations if we mimic observational methods to measure the mass of the stellar halo by selecting regions of an image via projected radius relative to the disk scale length or by their surface density in stellar mass . However, these observational methods systematically underestimate the accreted stellar component, defined in our (and most) simulations as the mass of stars formed outside of the host galaxy, by up to a factor of ten, since the accreted component is centrally concentrated and therefore substantially obscured by the galactic disk. Furthermore, these observational methods introduce spurious dependencies of the estimated accreted stellar component on the stellar mass and size of galaxies that can obscure the trends in accreted stellar mass predicted by cosmological simulations, since we find that in our simulations the size and shape of the central galaxy is not strongly correlated with the assembly history of the accreted stellar halo. This effect persists whether galaxies are viewed edge-on or face-on. We show that metallicity or color information may provide a way to more cleanly delineate in observations the regions dominated by accreted stars. Absent additional data, we caution that estimates of the mass of the accreted stellar component from single-band images alone should be taken as lower limits.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The Hard X-ray Modulation Telescope ($Insight$-HXMT) was successfully launched on June 15th, 2017. It performs broad band X-ray scan survey of the Galactic Plane to detect new black holes and other objects in active states. It also observes X-ray binaries to study their X-ray variabilities. Here we will introduce the Science Operations of $Insight$-HXMT, which is responsible for collecting and evaluating observation proposals, scheduling observations, and monitoring the working status of the payloads.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning and gamma ray burst (GRB) observations and , based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present the observational results from a detailed timing analysis of the black hole candidate Swift J1658.2-4242 during its 2018 outburst with the observations of Hard X-ray Modulation Telescope (Insight-HXMT), Neutron Star Interior Composition Explorer (NICER) and AstroSat in 0.1-250keV. The evolution of intensity, hardness and integrated fractional root mean square (rms) observed by Insight-HXMT and NICER are presented in this paper. Type-C quasi-periodic oscillations (QPOs) observed by NICER (0.8-3.5Hz) and Insight-HXMT (1-1.6Hz) are also reported in this work. The features of the QPOs are analysed with an energy range of 0.5-50keV. The relations between QPO frequency and other characteristics such as intensity, hardness and QPO rms are carefully studied. The timing and spectral properties indicate that Swift J1658.2-4242 is a black hole binary system. Besides, the rms spectra of the source calculated from the simultaneous observation of Insight-HXMT, NICER and AstroSat support the Lense-Thirring origin of the QPOs. The relation between QPO phase lag and the centroid frequency of Swift J1658.2-4242 reveals a near zero constant when < 4Hz and a soft phase lag at 6.68 Hz. This independence follows the same trend as the high inclination galactic black hole binaries such as MAXI J1659-152.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formation rate variations sufficiently large to reproduce the observed Balmer break of MACS1149-JD1. We find that several independent galaxy simulations are unable to account for Balmer breaks of the inferred size, suggesting that MACS1149-JD1 either must be a very rare type of object or that our simulations are missing some key ingredient. We present predictions of spectroscopic Balmer break strength distributions for $z\approx 7-9$ galaxies that may be tested through observations with the upcoming James Webb Space Telescope and also discuss the impact that various assumptions on dust reddening, Lyman continuum leakage and deviations from a standard stellar initial mass function would have on the results.

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimeter broadband flux to total infrared luminosity (L_IR), and hence star formation rate, in high-z galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g., redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (M* > 10^10 Msun) z=2-6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. At z~2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature does not strongly evolve with redshift over z=2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. The mass-weighted temperature is important for accurately estimating the total dust mass. We also analyse an 'equivalent' dust temperature for converting (sub)millimeter flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature ('warmer dust') do not necessarily have higher mass-weighted temperatures. A 'two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

We present the first detailed study of the spatially-resolved dust continuum emission of simulated galaxies at 1<z<5. We run the radiative transfer code SKIRT on a sample of submillimeter-bright galaxies drawn from the Feedback in Realistic Environments (FIRE) project. These simulated galaxies reach Milky Way masses by z=2. Our modelling provides predictions for the full rest-frame far-ultraviolet-to-far-infrared spectral energy distributions of these simulated galaxies, as well as 25-pc-resolution maps of their emission across the wavelength spectrum. The derived morphologies are notably different in different wavebands, with the same galaxy often appearing clumpy and extended in the far-ultraviolet yet an ordered spiral at far-infrared wavelengths. The observed-frame 870-$\mu$m half-light radii of our FIRE-2 galaxies are ~0.5-4kpc, consistent with existing ALMA observations of galaxies with similarly high redshifts and stellar masses. In both simulated and observed galaxies, the dust continuum emission is generally more compact than the cold gas and the dust mass, but more extended than the stellar component. The most extreme cases of compact dust emission seem to be driven by particularly compact recent star-formation, which generates steep dust temperature gradients. Our results confirm that the spatial extent of the dust continuum emission is sensitive to both the dust mass and SFR distributions.

We studied the long-term evolution of the spectral-temporal correlated properties of the black hole candidate Swift J1753.5-0127 from the onset of its outburst until 2011 with the Rossi X-ray Timing Explorer (RXTE). The source stayed most of its lifetime during hard state, with occasionally transitioned to the hard intermediate state. Similar to typical black hole transients, Swift J1753.5-0127 traces a clear hard line in absolute rms - intensity diagram during the low hard state, with expected highest absolute rms, while shows a clear turn during the hard intermediate state, accompanied by lower absolute rms. Different from Cyg X-1, we found that frequency-dependent time lag increased significantly in the 0.02 -- 3.2 Hz band during state transition in this source. The X-ray time lags in 0.02 -- 3.2 Hz can therefore be used as indicators of state transition in this source. Type-C QPO frequency is positively related with its fractional rms and X-ray photon index, suggesting a moving inwards disk/corona scenario. We discussed the physical interpretation of our results in our paper.

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies.

We investigate the bursty star formation histories (SFHs) of dwarf galaxies using the distribution of log($L_{H\alpha}/L_{UV}$) of 185 local galaxies. We expand on the work of Weisz et al. 2012 to consider a wider range of SFHs and stellar metallicities, and show that there are large degeneracies in a periodic, top-hat burst model. We argue that all galaxies of a given mass have similar SFHs and we can therefore include the $L_{H\alpha}$ distributions (subtracting the median trend with stellar mass, referred to as $\Delta \text{log}(L_{H\alpha})$) in our analyses. $\Delta \text{log}(L_{H\alpha})$ traces the amplitude of the bursts, and log($L_{H\alpha}/L_{UV}$) is a function of timescale, amplitude, and shape of the bursts. We examine the 2-dimensional distribution of these two indicators constrain the SFHs. We use exponentially rising/falling bursts to determine timescales ($e$-folding time, $\tau$). We find that galaxies below $10^{7.5}$ M$_{\odot}$ undergo large (amplitudes of $\sim 100$) and rapid ($\tau < 30$ Myr) bursts, while galaxies above $10^{8.5}$ M$_{\odot}$ experience smaller (maximum amplitudes $\sim 10$), slower ($\tau \gtrsim 300$ Myr) bursts. We compare to the FIRE-2 hydrodynamical simulations and find that the burst amplitudes agree with observations, but they are too rapid in more massive galaxies ($M_* > 10^8$ M$_{\odot}$). Finally, we confirm that stochastic sampling of the stellar mass function can not reproduce the observed distributions unless the standard assumptions of cluster and stellar mass functions are changed. With the next generation of telescopes, measurements of $L_{UV}$ and $L_{H\alpha}$ will become available for dwarf galaxies at high-redshift, enabling similar analyses of galaxies in the early universe.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

The first model independent results obtained by the DAMA/LIBRA-phase2 experiment are presented. The data have been collected over 6 annual cycles corresponding to a total exposure of 1.13 ton $\times$ yr, deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. The DAMA/LIBRA-phase2 apparatus, $\simeq$ 250 kg highly radio-pure NaI(Tl), profits from a second generation high quantum efficiency photomultipliers and of new electronics with respect to DAMA/LIBRA-phase1. The improved experimental configuration has also allowed to lower the software energy threshold. New data analysis strategies are presented. The DAMA/LIBRA-phase2 data confirm the evidence of a signal that meets all the requirements of the model independent Dark Matter (DM) annual modulation signature, at 9.5 $\sigma$ C.L. in the energy region (1-6) keV. In the energy region between 2 and 6 keV, where data are also available from DAMA/NaI and DAMA/LIBRA-phase1 (exposure $1.33$ ton $\times$ yr, collected over 14 annual cycles), the achieved C.L. for the full exposure (2.46 ton $\times$ yr) is 12.9 $\sigma$; the modulation amplitude of the single-hit scintillation events is: $(0.0103 \pm 0.0008)$ cpd/kg/keV, the measured phase is $(145 \pm 5)$ days and the measured period is $(0.999 \pm 0.001)$ yr, all these values are well in agreement with those expected for DM particles. No systematics or side reaction able to mimic the exploited DM signature (i.e. to account for the whole measured modulation amplitude and to simultaneously satisfy all the requirements of the signature), has been found or suggested by anyone throughout some decades thus far.

The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

Lyman-alpha (Ly{\alpha}) photons from ionizing sources and cooling radiation undergo a complex resonant scattering process that generates unique spectral signatures in high-redshift galaxies. We present a detailed Ly{\alpha} radiative transfer study of a cosmological zoom-in simulation from the Feedback In Realistic Environments (FIRE) project. We focus on the time, spatial, and angular properties of the Ly{\alpha} emission over a redshift range of z = 5-7, after escaping the galaxy and being transmitted through the intergalactic medium (IGM). Over this epoch, our target galaxy has an average stellar mass of $M_{\rm star} \approx 5 \times 10^8 {\rm M}_\odot$. We find that many of the interesting features of the Ly{\alpha} line can be understood in terms of the galaxy's star formation history. The time variability, spatial morphology, and anisotropy of Ly{\alpha} properties are consistent with current observations. For example, the rest frame equivalent width has a ${\rm EW}_{{\rm Ly}\alpha,0} > 20 {\rm \AA}$ duty cycle of 62% with a non-negligible number of sightlines with $> 100 {\rm \AA}$, associated with outflowing regions of a starburst with greater coincident UV continuum absorption, as these conditions generate redder, narrower (or single peaked) line profiles. The lowest equivalent widths correspond to cosmological filaments, which have little impact on UV continuum photons but efficiently trap Ly{\alpha} and produce bluer, broader lines with less transmission through the IGM. We also show that in dense self-shielding, low-metallicity filaments and satellites Ly{\alpha} radiation pressure can be dynamically important. Finally, despite a significant reduction in surface brightness with increasing redshift, Ly{\alpha} detections and spectroscopy of high-$z$ galaxies with the upcoming James Webb Space Telescope is feasible.

The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without AGN, and show that various proposed "non-AGN" solution mechanisms in the literature, including Type Ia supernovae, shocked AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields, Spitzer-Braginskii conduction, or "morphological quenching" do not halt or substantially reduce cooling flows nor maintain "quenched" galaxies in this mass range. We show that stellar feedback (including cosmic rays from SNe) alters the balance of cold/warm gas and the rate at which the cooled gas within the galaxy turns into stars, but not the net baryonic inflow. If anything, outflowing metals and dense gas promote additional cooling. Conduction is important only in the most massive halos, as expected, but even at $\sim 10^{14}\,{\rm M}_{\odot}$ reduces inflow only by a factor $\sim 2$ (owing to saturation effects and anisotropic suppression). Changing the morphology of the galaxies only slightly alters their Toomre-$Q$ parameter, and has no effect on cooling (as expected), so has essentially no effect on cooling flows or maintaining quenching. This all supports the idea that additional physics, e.g., AGN feedback, must be important in massive galaxies.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

We present a suite of cosmological zoom-in simulations at z>5 from the Feedback In Realistic Environments project, spanning a halo mass range M_halo~10^8-10^12 M_sun at z=5. We predict the stellar mass-halo mass relation, stellar mass function, and luminosity function in several bands from z=5-12. The median stellar mass-halo mass relation does not evolve strongly at z=5-12. The faint-end slope of the luminosity function steepens with increasing redshift, as inherited from the halo mass function at these redshifts. Below z~6, the stellar mass function and ultraviolet (UV) luminosity function slightly flatten below M_star~10^4.5 M_sun (fainter than M_1500~-12), owing to the fact that star formation in low-mass halos is suppressed by the ionizing background by the end of reionization. Such flattening does not appear at higher redshifts. We provide redshift-dependent fitting functions for the SFR-M_halo, SFR-M_star, and broad-band magnitude-stellar mass relations. We derive the star formation rate density and stellar mass density at z=5-12 and show that the contribution from very faint galaxies becomes more important at z>8. Furthermore, we find that the decline in the z~6 UV luminosity function brighter than M_1500~-20 is largely due to dust attenuation. Approximately 37% (54%) of the UV luminosity from galaxies brighter than M_1500=-13 (-17) is obscured by dust at z~6. Our results broadly agree with current data and can be tested by future observations.

Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypothesis of present-day metal-poor globular clusters. We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ~> 1e5-6 Msun collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. Our simulation then allows us to trace how clusters could become virialized and tightly-bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly-bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

We study the morphologies and sizes of galaxies at z>5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frame UV light; sizes measured from stellar mass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellar mass. The B-band light traces stellar mass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z>5, the size of galaxies at fixed stellar mass/luminosity evolves as (1+z)^{-m}, with m~1-2. For galaxies less massive than M_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% for more massive galaxies. A galaxy's "observed" size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z>5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z=6.

We use hydrodynamic cosmological zoom-in simulations from the FIRE project to explore the morphologies and kinematics of fifteen Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90% of their stellar mass within 2.5 kpc to well-ordered disks that reach $\gtrsim15$ kpc. The gas in our galaxies always forms a thin, rotation-supported disk at $z=0$, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disk-like orbits and with the radial extent of the stellar disk. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disk after $z\sim1$ develop well-ordered stellar disks. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-$z$ mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact disks in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in $z=0$ MW-mass galaxies (even $z=0$ bulge stars) form in a disk: $\gtrsim$60-90% of stars begin their lives rotationally supported.

Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these processes in cosmological simulations of dwarf galaxies with $z=0$ stellar masses $M_{\ast}\sim 10^{4}-3\times10^{6}\,M_\odot$ (halo masses $\sim 10^{9}-10^{10}\,M_\odot$). We show that the discrete nature of individual SNe (as opposed to a model in which their energy/momentum deposition is continuous over time, similar to stellar winds) is crucial in generating a reasonable ISM structure and galactic winds and in regulating dwarf stellar masses. However, once SNe are discretized, accounting for the effects of IMF sampling on continuous mechanisms such as radiative feedback and stellar mass-loss (as opposed to adopting IMF-averaged rates) has weak effects on galaxy-scale properties. We also consider the effects of rare HNe events with energies $\sim 10^{53}\,{\rm erg}$. The effects of HNe are similar to the effects of clustered explosions of SNe -- which are already captured in our default simulation setup -- and do not quench star formation (provided that the HNe do not dominate the total SNe energy budget), which suggests that HNe yield products should be observable in ultra-faint dwarfs today.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon excess appeared at a time coincident with the solar flare of 8 of November 2000. This muon excess corresponds to above 40 GeV primary protons which came from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation is estimated to be about 0.1%. It has been convinced that solar protons could be accelerated to tens of GeV in those Class X solar flares which usually arose solar cosmic ray ground level enhancement (GLE) events. However, whether a Class M solar flare like the non-GLE event of 8 November 2000 may also accelerate solar protons to such high energies? It is interesting and noteworthy.

One of the key observations regarding the stellar initial mass function (IMF) is its near-universality in the Milky Way (MW), which provides a powerful way to constrain different star formation models that predict the IMF. However, those models are almost universally "cloud-scale" or smaller -- they take as input or simulate single molecular clouds (GMCs), clumps, or cores, and predict the resulting IMF as a function of the cloud properties. Without a model for the progenitor properties of all clouds which formed the stars at different locations in the MW (including ancient stellar populations formed in high-redshift, likely gas-rich dwarf progenitor galaxies that looked little like the Galaxy today), the predictions cannot be fully explored, nor safely applied to "live" cosmological calculations of the IMF in different galaxies at different cosmological times. We therefore combine a suite of high-resolution cosmological simulations (from the Feedback In Realistic Environments project), which form MW-like galaxies with reasonable star formation properties and explicitly resolve massive GMCs, with various proposed cloud-scale IMF models. We apply the models independently to every star particle formed in the simulations to synthesize the predicted IMF in the present-day galaxy. We explore models where the IMF depends on Jeans mass, sonic or turbulent Bonner-Ebert mass, fragmentation with a polytropic equation-of-state, or where it is self-regulated by protostellar feedback. We show that all of these models, except the feedback-regulated ones, predict far more variation ($\sim 0.6-1$ dex $1\,\sigma$ scatter in the IMF turnover mass) in the simulations than is observed in the MW.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback in Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer-Braginskii conduction and viscosity, and sub-grid metal diffusion from unresolved turbulence. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological "zoom-in" simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modeled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs) . In contrast, removing stellar feedback leads to SFRs being over-predicted by factors of $\sim 10 -100$. Without feedback, neither galactic winds nor volume filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scale-heights with unphysically dense clumps congregating at the galactic center. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions and for the investigated halo mass range $10^{10}-10^{13} M_{\odot}$, the additional physics investigated here (MHD, conduction, viscosity, metal diffusion) have only weak ($\sim10\%$-level) effects on regulating SFR and altering the balance of phases, outflows, or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities, and feedback. We add the caveat that AGN feedback is not included in the present work.

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height increases with stellar age at all radii, because (1) stars that formed earlier were thicker "at birth", and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well-described by two components, with scale heights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disk in cosmological simulations even with strong stellar feedback.

Passive scalar mixing (metals, molecules, etc.) in the turbulent interstellar medium (ISM) is critical for abundance patterns of stars and clusters, galaxy and star formation, and cooling from the circumgalactic medium. However, the fundamental scaling laws remain poorly understood in the highly supersonic, magnetized, shearing regime relevant for the ISM. We therefore study the full scaling laws governing passive-scalar transport in idealized simulations of supersonic turbulence. Using simple phenomenological arguments for the variation of diffusivity with scale based on Richardson diffusion, we propose a simple fractional diffusion equation to describe the turbulent advection of an initial passive scalar distribution. These predictions agree well with the measurements from simulations, and vary with turbulent Mach number in the expected manner, remaining valid even in the presence of a large-scale shear flow (e.g. rotation in a galactic disk). The evolution of the scalar distribution is not the same as obtained using simple, constant "effective diffusivity" as in Smagorinsky models, because the scale-dependence of turbulent transport means an initially Gaussian distribution quickly develops highly non-Gaussian tails. We also emphasize that these are mean scalings that only apply to ensemble behaviors (assuming many different, random scalar injection sites): individual Lagrangian "patches" remain coherent (poorly-mixed) and simply advect for a large number of turbulent flow-crossing times.

We performed a global study on the timing and spectral properties of type-B quasi-periodic oscillations (QPOs) in the outbursts of black hole X-ray binaries. The sample is built based on the observations of {\it Rossi X-ray Timing Explorer}, via searching in the literature in RXTE era for all the identified type-B QPOs. To enlarge the sample, we also investigated some type-B QPOs that are reported but not yet fully identified. Regarding to the time lag and hard/soft flux ratio, we found that the sources with type-B QPOs behave in two subgroups. In one subgroup, type-B QPO shows a hard time lag that firstly decreases and then reverse into a soft time lag along with softening of the energy spectrum. In the other subgroup, type-B QPOs distribute only in a small region with hard time lag and relatively soft hardness. These findings may be understood with a diversity of the homogeneity showing up for the hot inner flow of different sources. We confirm the universality of a positive relation between the type-B QPO frequency and the hard component luminosity in different sources. We explain the results by considering that the type-B QPO photons are produced in the inner accretion flow around the central black hole, under a local Eddington limit. Using this relationship, we derived an mass estimation of $9.3-27.1$ $\rm{M_{\odot}}$ for the black hole in H 1743-322.

Recent spatially resolved observations of galaxies at z=0.6-3 reveal that high-redshift galaxies show complex kinematics and a broad distribution of gas-phase metallicity gradients. To understand these results, we use a suite of high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, which include physically motivated models of the multi-phase ISM, star formation, and stellar feedback. Our simulations reproduce the observed diversity of kinematic properties and metallicity gradients, broadly consistent with observations at z=0-3. Strong negative metallicity gradients only appear in galaxies with a rotating disk, but not all rotationally supported galaxies have significant gradients. Strongly perturbed galaxies with little rotation always have flat gradients. The kinematic properties and metallicity gradient of a high-redshift galaxy can vary significantly on short time-scales, associated with starburst episodes. Feedback from a starburst can destroy the gas disk, drive strong outflows, and flatten a pre-existing negative metallicity gradient. The time variability of a single galaxy is statistically similar to the entire simulated sample, indicating that the observed metallicity gradients in high-redshift galaxies reflect the instantaneous state of the galaxy rather than the accretion and growth history on cosmological time-scales. We find weak dependence of metallicity gradient on stellar mass and specific star formation rate (sSFR). Low-mass galaxies and galaxies with high sSFR tend to have flat gradients, likely due to the fact that feedback is more efficient in these galaxies. We argue that it is important to resolve feedback on small scales in order to produce the diverse metallicity gradients observed.

Supernova remnants (SNRs) are believed to be the most important acceleration sites for cosmic rays (CRs) below $\sim10^{15}$ eV in the Galaxy. High energy photons, either directly from the shocks of the SNRs or indirectly from the interaction between SNRs and the nearby clouds, are crucial probes for the CR acceleration. Big progresses on observations of SNRs have been achieved by space- and ground-based $\gamma$-ray facilities. However, whether $\gamma$-rays come from accelerated hadrons or not, as well as their connection with the CRs observed at Earth, remains in debate. Large High Altitude Air Shower Observatory (LHAASO), the next generation experiment, is designed to survey the northern part of the very high energy $\gamma$-ray sky from $\sim 0.3$ TeV to PeV with the sensitivity of $\lesssim1\%$ of the Crab nebula flux. In this paper, we indicate that LHAASO will be dedicated to enlarging the $\gamma$-ray SNR samples and improving the spectral and morphological measurements. These measurements, especially at energies above 30 TeV, will be important for us to finally understand the CR acceleration in SNRs.

Air radon concentration measurement is useful for research on geophysical effects, but it is strongly sensitive to site geology and many geophysical and microclimatic processes such as wind, ventilation, air humidity and so on that induce very big fluctuations on the concentration of radon in air. On the contrary, monitoring the radon concentration in soil by measuring the thermal neutron flux reduces environmental effects. In this paper we report some experimental results on the natural thermal neutron flux as well as the concentration of air radon and its variations at 4300 m a.s.l. These results were obtained with unshielded thermal neutron scintillation detectors (en-detectors) and radon monitors located inside the ARGO-YBJ experimental hall. The correlation of these variations with the lunar month and 1-year period is undoubtedly confirmed. A method for earthquakes prediction provided by a global net of the en-detectors is currently under study.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

In recent decades, inorganic crystals have been widely used in dark matter direct search experiments. To contribute to the understanding of the capabilities of $CsI(Na)$ and $CaF_{2}(Eu)$ crystals, a mono-energetic neutron beam is utilized to study the properties of nuclear recoils, which are expected to be similar to signals of dark matter direct detection. The quenching factor of nuclear recoils in $CsI(Na)$ and $CaF_{2}(Eu)$, as well as an improved discrimination factor between nuclear recoils and ${\gamma}$ backgrounds in $CsI(Na)$, are reported.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.

Supernovae of type Ia (SNe Ia) are believed to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs). However, the mass accretion process onto CO WDs is still not completely understood. In this paper, we study the accretion of He-rich matter onto CO WDs and explore a scenario in which a strong wind forms on the surface of the WD if the total luminosity exceeds the Eddington limit. Using a stellar evolution code called modules for experiments in stellar astrophysics (MESA), we simulated the He accretion process onto CO WDs for WDs with masses of 0.6-1.35Msun and various accretion rates of 10^{-8}-10^{-5}Msun/yr. If the contribution of the total luminosity is included when determining the Eddington accretion rate, then a super-Eddington wind could be triggered at relatively lower accretion rates than those of previous studies based on steady-state models. The super-Eddington wind can prevent the WDs with high accretion rates from evolving into red-giant-like He stars. We found that the contributions from thermal energy of the WD are non-negligible, judging by our simulations, even though the nuclear burning energy is the dominating source of luminosity. We also provide the limits of the steady He-burning regime in which the WDs do not lose any accreted matter and increase their mass steadily, and calculated the mass retention efficiency during He layer flashes for various WD masses and accretion rates. These obtained results can be used in future binary population synthesis computations.

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

In the present paper the results obtained in the investigation of possible diurnal effects for low-energy single-hit scintillation events of DAMA/LIBRA-phase1 (1.04 ton $\times$ yr exposure) have been analysed in terms of an effect expected in case of Dark Matter (DM) candidates inducing nuclear recoils and having high cross-section with ordinary matter, which implies low DM local density in order to fulfill the DAMA/LIBRA DM annual modulation results. This effect is due to the different Earth depths crossed by those DM candidates during the sidereal day.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

The super-Eddington wind scenario has been proposed as an alternative way for producing type Ia supernovae (SNe Ia). The super-Eddington wind can naturally prevent the carbon--oxygen white dwarfs (CO WDs) with high mass-accretion rates from becoming red-giant-like stars. Furthermore, it works in low-metallicity environments, which may explain SNe Ia observed at high redshifts. In this article, we systematically investigated the most prominent single-degenerate WD+MS channel based on the super-Eddington wind scenario. We combined the Eggleton stellar evolution code with a rapid binary population synthesis (BPS) approach to predict SN Ia birthrates for the WD+MS channel by adopting the super-Eddington wind scenario and detailed mass-accumulation efficiencies of H-shell flashes on the WDs. Our BPS calculations found that the estimated SN Ia birthrates for the WD+MS channel are ~0.009-0.315*10^{-3}{yr}^{-1} if we adopt the Eddington accretion rate as the critical accretion rate, which are much lower than that of the observations (<10% of the observed SN Ia birthrates). This indicates that the WD+MS channel only contributes a small proportion of all SNe Ia. The birthrates in this simulation are lower than previous studies, the main reason of which is that new mass-accumulation efficiencies of H-shell flashes are adopted. We also found that the critical mass-accretion rate has a significant influence on the birthrates of SNe Ia. Meanwhile, the results of our BPS calculations are sensitive to the values of the common-envelope ejection efficiency.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

This paper summarizes in a simple and intuitive way why the neutrons, the muons and the solar neutrinos cannot give any significant contribution to the DAMA annual modulation results. A number of these elements have already been presented in individual papers; they are recalled here. Afterwards, few simple considerations are summarized which already demonstrate the incorrectness of the claim reported in PRL 113 (2014) 081302.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

The results obtained in the search for possible diurnal effect in the single-hit low energy data collected by DAMA/LIBRA-phase1 (total exposure: 1.04 ton x yr) deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. are presented. At the present level of sensitivity the presence of any significant diurnal variation and of diurnal time structures in the data can be excluded for both the cases of solar and sidereal time. In particular, the diurnal modulation amplitude expected, because of the Earth diurnal motion, on the basis of the DAMA Dark Matter annual modulation results is below the present sensitivity.

Using all the observations from Rossi X-ray Timing Explorer for accreting neutron star XTE J1701-462, we carry out a systematic study on the cross-correlation between its soft and hard light curves. The anti-correlations appear at the hard vertex and on the upper normal branch (NB) in the Cyg-like interval, occur on the horizontal branch (HB) and/or upper NB in the first and second Sco-like intervals, and display at the upper flaring branch (FB) in the third Sco-like interval. We suggest that the anti-correlation might evolve with luminosity, because with decreasing of luminosity, the source evolved from a Cyg-like Z source, via a Sco-like Z source, to an atoll source. In the Cyg-like interval, the positive correlations are mostly distributed on the HB, which is not consistent with that of the Cyg-like Z source GX 5-1 and Cyg X-2 whose HBs host ambiguous correlations and anti-correlations. From our spectral analyses, it is found that the ratio of the hard emission to the soft emission basically keeps unvaried from the HB to the hard vertex, which might result in the positive correlation. Additionally, we compare the spectra in the third Sco-like interval with those of atoll source 4U 1735-44 and find that their spectral evolution along the tracks in the color-color diagrams are similar, indicating that in this interval, XTE J1701-462 approaches an atoll source. Both hard time lags and soft time lags are detected and the possible origins of them are discussed.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

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.

The results obtained with the total exposure of 1.04 ton x yr collected by DAMA/LIBRA-phase1 deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. during 7 annual cycles (i.e. adding a further 0.17 ton x yr exposure) are presented. The DAMA/LIBRA-phase1 data give evidence for the presence of Dark Matter (DM) particles in the galactic halo, on the basis of the exploited model independent DM annual modulation signature by using highly radio-pure NaI(Tl) target, at 7.5 sigma C.L.. Including also the first generation DAMA/NaI experiment (cumulative exposure 1.33 ton x yr, corresponding to 14 annual cycles), the C.L. is 9.3 sigma and the modulation amplitude of the single-hit events in the (2-6) keV energy interval is: (0.0112 \pm 0.0012) cpd/kg/keV; the measured phase is (144 \pm 7) days and the measured period is (0.998 \pm 0.002) yr, values well in agreement with those expected for DM particles. No systematic or side reaction able to mimic the exploited DM signature has been found or suggested by anyone over more than a decade.

The accretion of hydrogen-rich material onto carbon-oxygen white dwarfs (CO WDs) is crucial for understanding type Ia supernova (SN Ia) from the single-degenerate model, but this process has not been well understood due to the numerical difficulties in treating H and He flashes during the accretion. For the CO WD masses from 0.5 to $1.378\,{M}_\odot$ and accretion rates in the range from $10^{-8}$ to $10^{-5}\,{M}_\odot\,\mbox{yr}^{-1}$, we simulated the accretion of solar-composition material onto CO WDs using the state-of-the-art stellar evolution code of {\sc MESA}. For comparison with the steady-state models (e.g \citet{nskh07}), we firstly ignored the contribution from nuclear burning to the luminosity when determining the Eddington accretion rate and found that the properties of H burning in our accreting CO WD models are similar to those from the steady-state models, except that the critical accretion rates at which the WDs turn into red giants or H-shell flashes occur on their surfaces are slightly higher than those from the steady-state models. However, the super-Eddington wind is triggered at much lower accretion rates, than previously thought, when the contribution of nuclear burning to the total luminosity is included. This super-Eddington wind naturally prevents the CO WDs with high accretion rates from becoming red giants, thus presenting an alternative to the optically thick wind proposed by \cite{hkn96}. Furthermore, the super-Eddington wind works in low-metallicity environments, which may explain SNe Ia observed at high redshifts.

Some experiments have found multi-core events in cosmic rays extensive air showers which should be interpreted by hadronic interaction theory. In this paper, the multi-core events are reproduced by Monte Carlo simulation with CORSIKA. The origin of each sub-cores is tracked back from the observation level. The interaction mechanism and original particles of sub-core are studied in this paper.

The rich statistics of galaxy rotations as captured by the velocity function (VF) provides invaluable constraints on galactic baryon physics and the nature of dark matter (DM). However, the comparison of observed galaxy rotations against cosmological models is prone to subtle caveats that can easily lead to misinterpretations. Our analysis reveals full statistical consistency between ~5000 galaxy rotations, observed in line-of-sight projection, and predictions based on the standard cosmological model (LCDM) at the mass-resolution of the Millennium simulation (HI line-based circular velocities above ~50 km/s). Explicitly, the HI linewidths in the HI Parkes All Sky Survey (HIPASS) are found consistent with those in S3-SAX, a post-processed semi-analytic model for the Millennium simulation. Previously found anomalies in the VF can be plausibly attributed to (1) the mass-limit of the Millennium simulation, (2) confused sources in HIPASS, (3) inaccurate inclination measurements for optically faint sources, and (4) the non-detectability of gas-poor early-type galaxies. These issues can be bypassed by comparing observations and models using linewidth source counts rather than VFs. We investigate if and how well such source counts can constrain the temperature of DM.

Experimental observations and theoretical arguments at Galaxy and larger scales have suggested that a large fraction of the Universe is composed by Dark Matter particles. This has motivated the DAMA experimental efforts to investigate the presence of such particles in the galactic halo by exploiting a model independent signature and very highly radiopure set-ups deep underground. Few introductory arguments are summarized before presenting a review of the present model independent positive results obtained by the DAMA/NaI and DAMA/LIBRA set-ups at the Gran Sasso National Laboratory of the INFN. Implications and model dependent comparisons with other different kinds of results will be shortly addressed. Some arguments put forward in literature will be confuted.

Bound electron contribution to the Doppler-shift of gamma-ray spectra in the positron-electron annihilation process of molecular methane has been studied in gas phase. Two accurate ab initio quantum mechanical schemes, i.e. the delocalized molecular orbital (MO) and the localized natural bond orbital (NBO) schemes, are applied to study the multi-centred methane molecule. The present ab initio calculations of methane indicate that the C-H bonds are polarized with the partial negative charge of -0.36 a.u. on the carbon atom and the partial positive charge of +0.09 a.u. on each of the hydrogen atoms. The positively charged hydrogen atoms produce repulsive Coulomb potentials to a positron. Both the MO and NBO schemes further reveal that the 2a1 electrons of methane, that is, the 2a1 electron component of the C-H bonds rather than the whole C-H bonds of methane, predominates the positron-electron annihilation gamma-ray spectra of the molecule. Electrons of a molecule which are dominant the positron-electron annihilation processes are called positrophilic electrons in the present study. It is further shown that the negative electrostatic potential (ESP) of methane facilitates with the density of the positrophilic 2a1 electrons of methane. Other valence electrons (e.g. 1t2) in the C-H bonds play a minor spectator role in the annihilation process of methane.

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.

The DAMA/LIBRA experiment, running at the Gran Sasso National Laboratory of the I.N.F.N. in Italy, has a sensitive mass of about 250 kg highly radiopure NaI(Tl). It is mainly devoted to the investigation of Dark Matter (DM) particles in the Galactic halo by exploiting the model independent DM annual modulation signature. The present DAMA/LIBRA experiment and the former DAMA/NaI one (the first generation experiment having an exposed mass of about 100 kg) have released so far results corresponding to a total exposure of 1.17 ton yr over 13 annual cycles. They provide a model independent evidence of the presence of DM particles in the galactic halo at 8.9 sigma C.L.. A short summary of the obtained results is presented and future perspectives of the experiment mentioned.

By combining the newly infrared photometric data from the All-Sky Data Release of the Wide-field Infrared Survey Explorer with the spectroscopic data from the Seventh Data Release of the Sloan Digital Sky Survey, we study the covering factor of warm dust ($\CF$) for a large quasar sample, as well as the relations between $\CF$ and other physical parameters of quasars. We find a strong correlation between the flux ratio of mid-infrared to near-ultraviolet and the slope of near-ultraviolet spectra, which is interpreted as the dust extinction effect. After correcting for the dust extinction utilizing the above correlation, we examine the relations between $\CF$ and AGN properties: bolometric luminosity ($\Lbol$), black hole mass ($\MBH$) and Eddington ratio ($L/L_{\rm Edd}$). We confirm the anti-correlation between $\CF$ and $\Lbol$. Further we find that $\CF$ is anti-correlated with $\MBH$, but is independent of $L/L_{\rm Edd}$. Radio-loud quasars are found to follow the same correlations as for radio-quiet quasars. Monte Carlo simulations show that the anisotropy of UV-optical continuum of the accretion disc can significantly affect, but is less likely to dominate the $\CF$--$\Lbol$ correlation.

The abstract is additional with repect to the paper published in JGCD. Ordinary Earth satellites are usually placed into five categories of special orbits: sun-synchronous orbits, orbits at the critical inclination, frozen orbits, repeating ground track orbits, and geostationary orbits. This paper investigates their counterparts around Mars and examines the basic nature of these orbits, which are of special interest for missions conducted around Mars, including Mars reconnaissance. Mars' gravity field is much more complicated, with relatively smaller J2, compared to Earth's, which makes the behaviors of these Martian orbits different from those of Earth. Analytical formulations and numerical simulations are used to analyze these Martian orbits and compare them with their Earth counterparts. First, mean element theory is employed to describe variations of orbital elements and give the constraint conditions for achieving these special orbits. Then, numerical verifications based on the PSODE algorithm (particle swarm optimization combined with differential evolution) are adopted to provide more accurate conditions for achieving these orbits when considering an Mars gravity field. Using the numerical method can significantly improve the design in the full gravity field, and it is therefore possible to select these usable orbits for Mars that can reduce or eliminate the need for stationkeeping.

The triple asteroids and triple Kuiper belt objects (collectively called the triple minor planets) in the Solar system are of particular interest to the scientific community since the discovery of the first triple asteroid system in 2004. In this paper, the Hill stability of the nine known triple minor planets in the Solar system is investigated. Seven of the systems are of large size ratio, i.e. they consist of a larger primary and two moonlets, while the other two systems have components of comparable size. Each case is treated separately. For the triple minor planets that have large size ratio, the sufficient condition for Hill stability is expressed in closed form. This is not possible for the systems with comparable size components, for which the Hill stability is assessed by a combination of analytical and numerical means. It is found that all the known triple minor planets are Hill stable, except 3749 Balam, for which the incomplete orbital parameters make the Hill stability of the system uncertain. This suggests that there might be more such stable triple minor planets in the Solar system yet to be observed. It is also shown that the Hill stability regions increase as the mutual inclination between the inner orbit and outer orbit decreases, the semimajor axis ratio of the inner orbit with respect to the outer orbit decreases, and the mass ratio of the outer satellite with respect to the inner satellite increases. This study therefore provides useful information about dynamical properties of the triple minor planets in the Solar system.

A few final comments on arXiv:1210.7548 are given to confute incorrect arguments claimed there.

Theoretical gamma-ray spectra of molecule hexane have been calculated and compared with the experimental results in both gas (Surko et al, 1997) and liquid (Kerr et al, 1965) phases. The present study reveals that in gas phase not all valence electrons of hexane exhibit the same probability to annihilate a positron. Only the positrophilic electrons in the valence space dominate the gamma-ray spectra, which are in good agreement with the gas phase measurement. When hexane is confined in liquid phase, however, the intermolecular interactions ultimately eliminate the free molecular orientation and selectivity for the positrophilic electrons in the gas phase. As a result, the gamma-ray spectra of hexane become an averaged contribution from all valence electrons, which is again in agreement with liquid phase measurement. The roles of the positrophilic electrons in annihilation process for gas and liquid phases of hexane have been recognized for the first time in the present study.

We briefly remind references and arguments, already discussed in the past, which confute erroneous claims in arXiv:1210.5501.

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.

As is well known, a varying effective gravitational "constant" is one of the common features of most modified gravity theories. Of course, as a modified gravity theory, $f(T)$ theory is not an exception. Noting that the observational constraint on the varying gravitational "constant" is very tight, in the present work we try to constrain $f(T)$ theories with the varying gravitational "constant". We find that the allowed model parameter $n$ or $\beta$ has been significantly shrunk to a very narrow range around zero. In fact, the results improve the previous constraints by an order of magnitude.

This paper gathers arguments and reasons why muons surviving the Gran Sasso mountain cannot mimic the Dark Matter annual modulation signature exploited by the DAMA/NaI and DAMA/LIBRA experiments. A number of these items have already been presented in individual papers. Further arguments have been addressed here in order to present a comprehensive collection and to enable a wider community to correctly approach this point.

This study investigates the problem of areostationary orbits around Mars in the three-dimensional space. Areostationary orbits are expected to be used to establish a future telecommunication network for the exploration of Mars. However, no artificial satellites have been placed in these orbits thus far. In this paper, the characteristics of the Martian gravity field are presented, and areostationary points and their linear stability are calculated. By taking linearized solutions in the planar case as the initial guesses and utilizing the Levenberg-Marquardt method, families of periodic orbits around areostationary points are shown to exist. Short-period orbits and long-period orbits are found around linearly stable areostationary points, and only short-period orbits are found around unstable areostationary points. Vertical periodic orbits around both linearly stable and unstable areostationary points are also examined. Satellites in these periodic orbits could depart from areostationary points by a few degrees in longitude, which would facilitate observation of the Martian topography. Based on the eigenvalues of the monodromy matrix, the evolution of the stability index of periodic orbits is determined. Finally, heteroclinic orbits connecting the two unstable areostationary points are found, providing the possibility for orbital transfer with minimal energy consumption.

The results of Monte-Carlo simulations of Extensive Air Shower are presented to show the difference of hadronic component content at various altitudes with the aim to choose an optimal altitude for the PRISMA-like experiment. CORSIKA program for EAS simulations with QGSJET and GHEISHA models was used to calculate the number of hadrons reaching the observational level inside a ring of 50 m radius around the EAS axis. Then the number of neutrons produced by the hadronic component was calculated using an empirical relationship between the two components. We have tested the results with the ProtoPRISMA array at sea level, and recorded neutrons are close to the simulation results.

In the current study, a double-averaged analytical model including the action of the perturbing body's inclination is developed to study third-body perturbations. The disturbing function is expanded in the form of Legendre polynomials truncated up to the second-order term, and then is averaged over the periods of the spacecraft and the perturbing body. The efficiency of the double-averaged algorithm is verified with the full elliptic restricted three-body model. Comparisons with the previous study for a lunar satellite perturbed by Earth are presented to measure the effect of the perturbing body's inclination, and illustrate that the lunar obliquity with the value 6.68\degree is important for the mean motion of a lunar satellite. The application to the Mars-Sun system is shown to prove the validity of the double-averaged model. It can be seen that the algorithm is effective to predict the long-term behavior of a high-altitude Martian spacecraft perturbed by Sun. The double-averaged model presented in this paper is also applicable to other celestial systems.

Some theoretical and experimental aspects regarding the direct dark matter field are mentioned. In particular some arguments, which play a relevant role in the evaluation of model dependent interpretations of experimental results and in comparisons, are shortly addressed.

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.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

In analogy to $f(R)$ theory, recently a new modified gravity theory, namely the so-called $f(T)$ theory, has been proposed to drive the current accelerated expansion without invoking dark energy. In the present work, by extending Bisabr's idea, we try to constrain $f(T)$ theories with the varying fine structure "constant", $\alpha\equiv e^2/\hbar c$. We find that the constraints on $f(T)$ theories from the observational $\Delta\alpha/\alpha$ data are very severe. In fact, they make $f(T)$ theories almost indistinguishable from $\Lambda$CDM model.

This paper investigates the dynamics of a particle orbiting around a rotating homogeneous cube, and shows fruitful results that have implications for examining the dynamics of orbits around non-spherical celestial bodies. This study can be considered as an extension of previous research work on the dynamics of orbits around simple shaped bodies, including a straight segment, a circular ring, an annulus disk, and simple planar plates with backgrounds in celestial mechanics. In the synodic reference frame, the model of a rotating cube is established, the equilibria are calculated, and their linear stabilities are determined. Periodic orbits around the equilibria are computed using the traditional differential correction method, and their stabilities are determined by the eigenvalues of the monodromy matrix. The existence of homoclinic and heteroclinic orbits connecting periodic orbits around the equilibria is examined and proved numerically in order to understand the global orbit structure of the system. This study contributes to the investigation of irregular shaped celestial bodies that can be divided into a set of cubes.

Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80\times80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre are kept at 270 degrees because J3 has the same sign with J2; while for a minority of ones with low altitude and low inclination, argument of pericentre are possible to be kept at 90 degrees because of the effect of the higher degree odd zonals; for the critical inclinations cases, argument of pericentre can also be kept at 90 degrees. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.

In the current study, the existence of periodic orbits around a fixed homogeneous cube is investigated, and the results have powerful implications for examining periodic orbits around non-spherical celestial bodies. In the two different types of symmetry planes of the fixed cube, periodic orbits are obtained using the method of the Poincar\'e surface of section. While in general positions, periodic orbits are found by the homotopy method. The results show that periodic orbits exist extensively in symmetry planes of the fixed cube, and also exist near asymmetry planes that contain the regular Hex cross section. The stability of these periodic orbits is determined on the basis of the eigenvalues of the monodromy matrix. This paper proves that the homotopy method is effective to find periodic orbits in the gravity field of the cube, which provides a new thought of searching for periodic orbits around non-spherical celestial bodies. The investigation of orbits around the cube could be considered as the first step of the complicated cases, and helps to understand the dynamics of orbits around bodies with complicated shapes. The work is an extension of the previous research work about the dynamics of orbits around some simple shaped bodies, including a straight segment, a circular ring, an annulus disk, and simple planar plates.

The critical inclination is of special interest in artificial satellite theory. The critical inclination can maintain minimal deviations of eccentricity and argument of pericentre from the initial values, and orbits at this inclination have been applied to some space missions. Most previous researches about the critical inclination were made under the assumption that the oblateness term J2 is dominant among the harmonic coefficients. This paper investigates the extension of the critical inclination where the concept of the critical inclination is different from that of the traditional sense. First, the study takes the case of Venus for instance, and provides some preliminary results. Then for general cases, given the values of argument of pericentre and eccentricity, the relationship between the multiplicity of the solutions for the critical inclination and the values of J2 and J4 is analyzed. Besides, when given certain values of J2 and J4, the relationship between the multiplicity of the solutions for the critical inclination and the values of semimajor axis and eccentricity is studied. The results show that for some cases, the value of the critical inclination is far away from that of the traditional sense or even has multiple solutions. The analysis in this paper could be used as starters of correction methods in the full gravity field of celestial bodies

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.