By the end of 2014, a cosmic ray muon telescope was installed at Zhongshan Station in Antarctic and has been continuously collecting data since then. It is the first surface muon telescope to be built in Antarctic. In June 2015, five CMEs were ejected towards the Earth initiating a big large Forbush decrease (FD) event. We conduct a comprehensive study of the galactic cosmic ray intensity fluctuations during the FD using the data from cosmic ray detectors of multiple stations (Zhongshan, McMurdo, South Polar and Nagoya) and he solar wind measurements from ACE and WIND. A pre-increase before the shock arrival was observed. Distinct differences exist in the timelines of the galactic cosmic ray recorded by the neutron monitors and the muon telescopes. FD onset for Zhongshan muon telescope is delayed (2.5h) with respect to SSC onset. This FD had a profile of four-step decrease. The traditional one- or two-step classification of FDs was inadequate to explain this FD.

We study the classical dynamics of the universe undergoing a transition from contraction to expansion through a big bang singularity. The dynamics is described by a system of differential equations for a set of physical quantities, such as the scale factor $a$, the Hubble parameter $H$, the equation of state parameter $w$, and the density parameter $\Omega$. The solutions of the dynamical system have a singularity at the big bang. We study if the solutions can be regularized at the singularity in the sense of whether they have unique branch extensions through the singularity. In particular, we consider the model in which the contracting universe is dominated by a scalar field with a time varying equation of state $w$, which approaches a constant value $w_c$ near the singularity. We prove that, for $w_c > 1$, the solutions are regularizable only for a discrete set of $w_c$ values that satisfy a coprime number condition. Our result implies that the evolution of a bouncing universe through the big bang singularity does not have a continuous classical limit unless the equation of state is extremely fine-tuned.

In bouncing cosmology, the primordial fluctuations are generated in a cosmic contraction phase before the bounce into the current expansion phase. For a nonsingular bounce, curvature and anisotropy grow rapidly during the bouncing phase, raising questions about the reliability of perturbative analysis. In this paper, we study the evolution of adiabatic perturbations in a nonsingular bounce by nonperturbative methods including numerical simulations of the nonsingular bounce and the covariant formalism for calculating nonlinear perturbations. We show that the bounce is disrupted in regions of the universe with significant inhomogeneity and anisotropy over the background energy density, but is achieved in regions that are relatively homogeneous and isotropic. Sufficiently small perturbations, consistent with observational constraints, can pass through the nonsingular bounce with negligible alteration from nonlinearity. We follow scale invariant perturbations generated in a matter-like contraction phase through the bounce. Their amplitude in the expansion phase is determined by the growing mode in the contraction phase, and the scale invariance is well preserved across the bounce.

We introduce a new set of flow parameters to describe the time dependence of the equation of state and the speed of sound in single field cosmological models. A scale invariant power spectrum is produced if these flow parameters satisfy specific dynamical equations. We analyze the flow of these parameters and find four types of fixed points that encompass all known single field models. Moreover, near each fixed point we uncover new models where the scale invariance of the power spectrum relies on having simultaneously time varying speed of sound and equation of state. We describe several distinctive new models and discuss constraints from strong coupling and superluminality.

We consider bouncing cosmologies in which an ekpyrotic contraction phase with w >> 1 is followed by a bouncing phase with w < -1 that violates the null energy condition. The bouncing phase, induced by ghost condensation, is designed to produce a classically nonsingular bounce at a finite value of the scale factor. We show that the initial curvature and anisotropy, though diluted during the ekpyrotic phase, grow back exponentially during the bouncing phase. Moreover, curvature perturbations and anisotropy are generated by quantum fluctuations during the ekpyrotic phase. In the bouncing phase, however, an adiabatic curvature perturbation grows to dominate and contributes a blue spectrum that spoils the scale-invariance. Meanwhile, a scalar shear perturbation grows nonlinear and creates an overwhelming anisotropy that disrupts the nonsingular bounce altogether.

We consider non-singular bouncing cosmologies, such as the new ekpyrotic model, in which the universe undergoes a slow contraction phase with equation of state $w \gg 1$, followed by a bounce that occurs at a finite scale factor when quantum gravity corrections are still negligible. Such a non-singular bounce requires a violation of the null energy condition in which $w$ falls below -1 at some time before the bounce. In this paper, we show that a component of the adiabatic curvature perturbations, though decaying and negligible during the ekpyrotic phase, is exponentially amplified just before $w$ approaches -1, enough to spoil the scale-invariant perturbation spectrum. We discuss how the problem may be avoided, for example, in singular bounces.

We propose to extract the charge information of high energy muons in very inclined extensive air showers by analyzing their relative lateral positions in the shower transverse plane. We calculate the muon lateral deviation under the geomagnetic field and compare it to dispersive deviations from other causes. By our criterion of resolvability, positive and negative muons with energies above $10^4$ GeV will be clearly separated into two lobes if the shower zenith angle is larger than $70^\circ$. Thus we suggest a possible approach to measure the $\mu^+ / \mu^-$ ratio for high energy muons.