We investigate observational constraints on the running vacuum model (RVM) of $\Lambda=3\nu (H^{2}+K/a^2)+c_0$ in the spatially curved universe, where $\nu$ is the model parameter, $K$ corresponds to the spatial curvature constant, and $c_{0}$ is a constant defined by the boundary conditions. As $\dot{\Lambda} \ne 0$, there are energy exchanges between vacuum, matter and radiation in RVM. We study the "geometrical degeneracy" of RVM on the CMB power spectra. By fitting the cosmological data, we find that the values of $\chi^2$ in RVM and $\Lambda$CDM are similar to each other for the non-flat universe. Explicitly, we obtain the constraints of $\nu\leq O(10^{-4})$ (68 $\%$ C.L.) and $|\Omega_K|\leq O(10^{-2})$ (95 $\%$ C.L.) in our study. In addition, we show that the cosmological constraints of $\Sigma m_{\nu}=0.416^{+0.311}_{-0.407}$ (RVM) and $\Sigma m_{\nu}=0.497^{+0.335}_{-0.387}$ ($\Lambda$CDM) at 95$\%$ C.L. for the neutrino mass sum are relaxed in both models in the spatially curved universe.

In the light of the latest measurements on the total $e^+ + e^-$ flux by CALET and DAMPE experiments, we revisit the multicomponent leptonically decaying dark matter (DM) explanations to the cosmic-ray electron/positron excesses observed previously. Especially, we use the single and double-component DM models to explore the compatibility of the AMS-02 positron fraction with the new CALET or DAMPE data. It turns out that neither single nor double-component DM models are able to fit the AMS-02 positron fraction and DAMPE total $e^+ + e^-$ flux data simultaneously. On the other hand, for the combined AMS-02 and CALET dataset, both the single and double-component DM models can provide reasonable fits. If we further take into the diffuse $\gamma$-ray constraints from Fermi-LAT, only the double-component DM models are allowed.

We investigate the cosmological viability of the generalized proca theory. We first implement the background and linear perturbation equations of motion in the Boltzmann code and then study the constraints on the parameters of the generalized proca theory after running MCMC against the cosmological data set. With Planck + HST data, we obtain the constraint $h=0.7334_{-0.0269}^{+0.0246}$, which indicates that the tension between early universe and late time universe within this theory is removed. By adding other late-time data sets (BAO, RSD, etc.) we show that the tension is reduced, as the 2$\sigma$ allowed region for $h$ in Proca, $h=0.7041_{-0.0087}^{+0.0094}$, overlaps with the 2$\sigma$ region of the HST data.

We study a special running vacuum model (RVM) with $\Lambda = 3 \alpha H^2+3\beta H_0^4 H^{-2}+\Lambda_0$, where $\alpha$, $\beta$ and $\Lambda_0$ are the model parameters and $H$ is the Hubble one. This RVM has non-analytic background solutions for the energy densities of matter and radiation, which can only be evaluated numerically. From the analysis of the CMB power spectrum and baryon acoustic oscillation along with the prior of $\alpha>0$ to avoid having a negative dark energy density, we find that $\alpha<2.83\times 10^{-4}$ and $\beta=(-0.2^{+3.9}_{-4.5})\times 10^{-4}$ (95$\%$ C.L.). We show that the RVM fits the cosmological data comparably to the $\Lambda$CDM. In addition, we relate the fluctuation amplitude $\sigma_8$ to the neutrino mass sum $\Sigma m_\nu$.

We apply the first law of thermodynamics to the apparent horizon of the universe with the power-law corrected and non-extensive Tsallis entropies rather than the Bekenstein-Hawking one. We examine the cosmological properties in the two entropy models by using the CosmoMC package. In particular, the first numerical study for the cosmological observables with the power-law corrected entropy is performed. We also show that the neutrino mass sum has a non-zero central value with a relaxed upper bound in the Tsallis entropy model comparing with that in the $\Lambda$CDM one.

Due to the universality of gravitational interactions, it is generally expected that a stochastic gravitational wave (GW) background could form during the reheating period when the inflaton perturbatively decays with the emission of gravitons. Previously, only models in which the inflaton dominantly decays into a pair of light scalar and/or fermion particles were considered in the literature. In the present paper, we focus on the cases with a vector particle pair in the final decay product. The differential decay rates for the three-body gravitational inflaton decays are presented for two typical couplings between the inflaton and vector fields, from which we predict their respective GW frequency spectra. It turns out that, similar to the scalar and fermion cases, the obtained GW spectra is too high in frequency to be observed by the current and near-future GW detection experiments and calls for a new design of high-frequency GW detectors.

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.

Large High Altitude Air Shower Observatory(LHAASO) is a composite cosmic ray observatory consisting of three detector arrays: kilometer square array (KM2A) which includes the electromagnetic detector array and muon detector array, water Cherenkov detector array (WCDA) and wide field of view Cherenkov telescope array (WFCTA). One of the main scientific objectives of LHAASO is to precisely measure the cosmic rays energy spectrum of individual components from 1014 eV to 1018 eV. The hybrid observation will be employed by LHAASO experiment, in which the lateral and longitudinal distributions of the extensive air shower can be observed simultaneously. Thus many kinds of parameters can be used for primary nuclei identification. In this paper, high purity cosmic ray simulation samples of light nuclei component are obtained through Multi-Variable Analysis. The apertures of 1/4 LHAASO array for pure proton and mixed proton and helium (H&He) samples are 900 m2Sr and 1800 m2Sr respectively. A prospect of proton and H&He spectra from 100 TeV to 4 PeV is discussed.

We study the $\Lambda(\alpha)$CDM models with $\Lambda(\alpha)$ being a function of the time-varying fine structure constant $\alpha$. We give a close look at the constraints on two specific $\Lambda(\alpha)$CDM models with one and two model parameters, respectively, based on the cosmological observational measurements along with 313 data points for the time-varying $\alpha$. We find that the model parameters are constrained to be around $10^{-4}$, which are similar to the results discussed previously but more accurately.

The cosmic-ray (CR) knee and the compositions contain abundant information for probing the CR's origin, acceleration and propagation mechanisms, as well as the frontier of the fundamental physics. Realizing that major proposals toward the knee's shape can be divided into two categories: the rigidity-dependent (also Z-dependent) knee and the mass-dependent (also A-dependent) knee, where the former one relates to the acceleration or the propagation mechanisms, and the other one is often associated with the new physics, it is essential to precisely measure the individual compositions. Benefit from the high altitude and hybrid detection methods, the LHAASO experiment has the ability in determining the individual component and brings us an opportunity in discriminating these two models. We test this expected ability of LHAASO from 100 TeV to 10 PeV with 3-year observation. And find the dominant component at the knee is essential to this issue, while much heavier nuclei occupying the knee leads to higher significance. In the analysis, the He-dominant knee under the A-dependent case can be recognized at the significance about 6.6 $\sigma$, while the P-dominant knee under the Z-dependent case will be classified with 2 $\sigma$ significance.

We examine the running vacuum model with $\Lambda (H) = 3 \nu H^2 + \Lambda_0$, where $\nu$ is the model parameter and $\Lambda_0$ is the cosmological constant. From the data of the cosmic microwave background radiation, weak lensing and baryon acoustic oscillation along with the time dependent Hubble parameter $H(z)$ and weighted linear growth $f (z)\sigma_8(z)$ measurements, we find that $\nu=(1.37^{+0.72}_{-0.95})\times 10^{-4}$ with the best fitted $\chi^2$ value slightly smaller than that in the $\Lambda$CDM model.

We investigate the interacting dark energy models by using the diagnostics of statefinder hierarchy and growth rate of structure. We wish to explore the deviations from $\Lambda$CDM and to differentiate possible degeneracies in the interacting dark energy models with the geometrical and structure growth diagnostics. We consider two interacting forms for the models, i.e., $Q_1=\beta H\rho_c$ and $Q_2=\beta H\rho_{de}$, with $\beta$ being the dimensionless coupling parameter. Our focus is the I$\Lambda$CDM model that is a one-parameter extension to $\Lambda$CDM by considering a direct coupling between the vacuum energy ($\Lambda$) and cold dark matter (CDM), with the only additional parameter $\beta$. But we begin with a more general case by considering the I$w$CDM model in which dark energy has a constant $w$ (equation-of-state parameter). For calculating the growth rate of structure, we employ the "parametrized post-Friedmann" theoretical framework for interacting dark energy to numerically obtain the $\epsilon(z)$ values for the models. We show that in both geometrical and structural diagnostics the impact of $w$ is much stronger than that of $\beta$ in the I$w$CDM model. We thus wish to have a closer look at the I$\Lambda$CDM model by combining the geometrical and structural diagnostics. We find that the evolutionary trajectories in the $S^{(1)}_3$--$\epsilon$ plane exhibit distinctive features and the departures from $\Lambda$CDM could be well evaluated, theoretically, indicating that the composite null diagnostic $\{S^{(1)}_3, \epsilon\}$ is a promising tool for investigating the interacting dark energy models.

We present large scale 3D particle-in-cell (PIC) simulations to examine particle energization in magnetic reconnection of relativistic electron-positron (pair) plasmas. The initial configuration is set up as a relativistic Harris equilibrium without a guide field. These simulations are large enough to accommodate a sufficient number of tearing and kink modes. Contrary to the non-relativistic limit, the linear tearing instability is faster than the linear kink instability, at least in our specific parameters. We find that the magnetic energy dissipation is first facilitated by the tearing instability and followed by the secondary kink instability. Particles are mostly energized inside the magnetic islands during the tearing stage due to the spatially varying electric fields produced by the outflows from reconnection. Secondary kink instability leads to additional particle acceleration. Accelerated particles are, however, observed to be thermalized quickly. The large amplitude of the vertical magnetic field resulting from the tearing modes by the secondary kink modes further help thermalizing the non-thermal particles generated from the secondary kink instability. Implications of these results for astrophysics are briefly discussed.