The phantom dark energy remarkably boosts our prehension of the accelerating Universe. Various models are widely discussed in the phantom Universe without bulk viscosity. From the hydrodynamics' point of view, it is natural to introduce the nonperfect fluid in the study of the Universe, as an ideal fluid just an approximation to the real world after all and using the generalized equation of state (EoS) with bulk viscosity, the early inflationary universe and the accelerated expansion of the late-time universe are described by many authors. In this paper, in order to investigate how the viscosity will influence the evolution of the Universe, we study a class of phantom dark energy models with bulk viscosity by the method of dynamical analysis technique. We show that there are different cosmic late-time behaviors and the stability also brings some constraints on the models. We also plot the evolutionary trajectories of this model in the statefinder parameter-planes to see the different behaviors of the models from the statefinder viewpoint.

Recently, a new kind of $f(z)$ theory is proposed to provide a different perspective for the development of reliable alternative models of gravity in which the $f(R)$ Lagrangian terms are reformulated as polynomial parameterizations $f(z)$. In the previous study, the parameters in the $f(z)$ models have been constrained by using cosmological data. In this paper, these models will be tested by the observations in the solar system. After solving the Ricci scalar as a function of the redshift, one could obtain $f(R)$ that could be used to calculate the standard Parameterized-Post-Newtonian (PPN) parameters. First, we fit the parametric models with the latest cosmological observational data. Then the tests are performed by solar system observations. And last we combine the constraints of solar system and cosmology together and reconstruct the $f(R)$ actions of the $f(z)$ parametric models.

The spectral energy distribution (SED) sequence for type Ia supernovae (SN Ia) is modeled by an artificial neural network. The SN Ia luminosity is characterized as a function of phase, wavelength, a color parameter and a decline rate parameter. After training and testing the neural network, the SED sequence could give both the spectrum with wavelength range from 3000\AA~to 8000\AA~ and the light curve with phase from 20 days before to 50 days after the maximum luminosity for the supernovae with different colors and decline rates. Therefore, we call this the Artificial Neural Network Spectral Light Curve Template (ANNSLCT) model. We retrain the Joint Light-curve Analysis (JLA) supernova sample by using the ANNSLCT model and obtain the parameters for each supernova to make a constraint on the cosmological $\Lambda$CDM model. We find that the best fitting values of these parameters are almost the same as those from the JLA sample trained with the Spectral Adaptive Lightcurve Template 2 (SALT2) model. So we believe that the ANNSLCT model could be used to analyze a large number of SN Ia multi-color light curves measured in the current and future observational projects.

There is undetermined potential function $V(\phi)$ in the action of mimetic gravity which should be resolved through physical means. In general relativity(GR), the static spherically symmetric(SSS) solution to the Einstein equation is a benchmark and its deformation also plays a crucial role in mimetic gravity. The equation of motion is provided with high nonlinearity, but we can reduce primal nonlinearity to a frequent Riccati form in the SSS case of mimetic gravity. In other words, we obtain an expression of solution to the functional differential equation of motion with any potential function. Remarkably, we proved rigorously that there is a zero point of first order for the metric function $\beta(r)$ if another metric function $\alpha(r)$ possesses a pole of first order within mimetic gravity. The zero point theorem may be regarded as the generalization of Birkhoff theorem $\alpha\beta=1$ in GR. As a corollary, we show that there is a modified black hole solution for any given $V(\phi)$, which can pass the test of solar system. As another corollary, the zero point theorem provides a dynamical mechanism for the maximum size of galaxies. Especially, there are two analytic solutions which provide good fits to the rotation curves of galaxies without the demand for particle dark matter.

The currently accelerated expansion of our Universe is unarguably one of the most intriguing problems in today's physics research. Two realistic non-minimal torsion-matter coupling $f(T)$ models have been established and studied in our previous papers [Phys. Rev. D92, 104038(2015) and Eur. Phys. J. C77, 504(2017)] aiming to explain this "dark energy" problem. In this paper, we study the generalized power-law torsion-matter coupling $f(T)$ model. Dynamical system analysis shows that the three expansion phases of the Universe, i.e. the radiation dominated era, the matter dominated era and the dark energy dominated era, can all be reproduced in this generalized model. By using the statefinder and $Om$ diagnostics, we find that the different cases of the model can be distinguished from each other and from other dark energy models such as the two models in our previous papers, $\Lambda$CDM, quintessence and Chaplygin gas. Furthermore, the analyses also show that all kinds of generalized power-law torsion-matter coupling model are able to cross the $w=-1$ divide from below to above, thus the decrease of the energy density resulting from the crossing of $w$ will make the catastrophic fate of the Universe avoided and a de Sitter expansion fate in the future will be approached.

We construct and train an artificial neural network called the back-propagation neural network to describe the evolution of the type Ia supernova spectrum by using the data from the CfA Supernova Program. This network method has many attractive features, and one of them is that the constructed model is differentiable. Benefitting from this, we calculate the absorption velocity and its variation. The model we constructed can well describe not only the spectrum of SNe Ia with wavelength range from $3500\AA$ to $8000\AA$, but also the light-curve evolution with phase time from $-15$ to $50$ with different colors. Moreover, the number of parameters needed during the training process is much less than the usual methods.

It is known that the Cardassian universe is successful in describing the accelerated expansion of the universe, but its dynamical equations are hard to get from the action principle. In this paper, we establish the connection between the Cardassian universe and $f(T, \mathcal{T})$ gravity, where $T$ is the torsion scalar and $\mathcal{T}$ is the trace of the matter energy-momentum tensor. For dust matter, we find that the modified Friedmann equations from $f(T, \mathcal{T})$ gravity can correspond to those of Cardassian models, and thus, a possible origin of Cardassian universe is given. We obtain the original Cardassian model, the modified polytropic Cardassian model, and the exponential Cardassian model from the Lagrangians of $f(T,\mathcal{T})$ theory. Furthermore, by adding an additional term to the corresponding Lagrangians, we give three generalized Cardassian models from $f(T,\mathcal{T})$ theory. Using the observation data of type Ia supernovae, cosmic microwave background radiation, and baryon acoustic oscillations, we get the fitting results of the cosmological parameters and give constraints of model parameters for all of these models.

Using the classical top-hat profile, we study the non-linear growth of spherically symmetric density perturbation and structure formation in $f(T)$ gravities. In particular, three concrete models, which have been tested against the observation of large-scale evolution and linear perturbation of the universe in the cosmological scenario, are investigated in this framework, covering both minimal and nonminimal coupling cases of $f(T)$ gravities. Moreover, we consider the virialization of the overdense region in the models after they detach from the background expanding universe and turn around to collapse. We find that there are constraints in the magnitude and occurring epoch of the initial perturbation. The existence of these constraints indicates that a perturbation that is too weak or occurs too late will not be able to stop the expanding of the overdense region. The illustration of the evolution of the perturbation shows that in $f(T)$ gravities, the initial perturbation within the constraints can eventually lead to clustering and form structure. The evolution also shows that nonminimal coupling models collapse slower than the minimal coupling one.

We present a detailed study of the Vaidya solution and its generalization in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored with the St\"{u}ckelberg fields $\phi^a$ introduced, there is a new invariant $I^{ab}=g^{\mu \nu}\partial_\mu \phi^a\partial_\nu \phi^b$ in the massive gravity, which adds to the ones usually encountered in general relativity. There is no conventional Vaidya solution if we choose unitary gauge. In this paper, we obtain three types of self-consistent ansatz with some nonunitary gauge, and find accordingly the Vaidya, generalized Vaidya and furry Vaidya solution. As by-products, we obtain a series of furry black hole. The Vaidya solution and its generalization in dRGT massive gravity describe the black holes with a variable horizon.

Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulationsof the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal component that is continuously being injected into the domain. This flux injection mimics the poloidal currents driven by the anode-cathode voltage drop in the experiment. The injected toroidal field stretches the poloidal fields to large distances, while forming a collimated jet along with several other key features. Detailed comparisons between 3D MHD simulations and experimental measurements provide a comprehensive description of the interplay among magnetic force, pressure and flow effects. In particular, we delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms. With suitably chosen parameters that are derived from experiments, the jet in the simulation agrees quantitatively with the experimental jet in terms of magnetic/kinetic/inertial energy, total poloidal current, voltage, jet radius, and jet propagation velocity. Specifically, the jet velocity in the simulation is proportional to the poloidal current divided by the square root of the jet density, in agreement with both the experiment and analytical theory. This work provides a new and quantitative method for relating experiments, numerical simulations and astrophysical observation, and demonstrates the possibility of using terrestrial laboratory experiments to study astrophysical jets.

We obtain the Schwarzschild solution from thermodynamic considerations using the assumptions of a quasi local mass form (the Misner-Sharp mass) and geometric surface gravity in a spherically symmetric spacetime. The deduction is extended to other cases such as the de Sitter, anti-de Sitter, Reissner-Nordstr\"om and higher dimensional spacetimes. This paper demonstrates the simple hypotheses to obtain these known spherically symmetric solutions via thermodynamics, where essentially the Misner-Sharp mass is the mass for an adiabatic system.

Depending on the choice of its potential, the scalar phantom field $\phi$ (the equation of state parameter $w<-1$) leads to various catastrophic fates of the universe including big rip, little rip and other future singularity. For example, big rip results from the evolution of the phantom field with an exponential potential and little rip stems from a quadratic potential in general relativity (GR). By choosing the same potential as in GR, we suggest a new mechanism to avoid these unexpected fates (big and little rip) in the inverse-\textit{R} gravity. As a pedagogical illustration, we give an exact solution where phantom field leads to a power-law evolution of the scale factor in an exponential type potential. We also find the sufficient condition for a universe in which the equation of state parameter crosses $w=-1$ divide. The phantom field with different potentials, including quadratic, cubic, quantic, exponential and logarithmic potentials are studied via numerical calculation in the inverse-\textit{R} gravity with $R^{2}$ correction. The singularity is avoidable under all these potentials. Hence, we conclude that the avoidance of big or little rip is hardly dependent on special potential.

We present theoretical results of a piecewise isothermal shock wind model devised for predicting the luminosity and surface brightness profile of diffuse X-ray emissions primarily from the inner shocked downstream wind zone of a planetary nebula (PN) surrounded by self-similar shocked dense shell and outer slow AGB wind envelope involving self-gravity and compare/fit our computational model results with available observations of a few grossly spherical X-ray emitting PNe. Matching shocked piecewise isothermal self-similar void (ISSV) solutions with self-gravity of Lou & Zhai (LZ) for the outer zone and a stationary isothermal fast tenuous wind with a reverse shock for the inner zone across an expanding contact discontinuity, we can consistently construct dynamic evolution models of PNe with diffuse X-ray emissions. On the basis of such a chosen dynamic wind interaction model, both X-ray luminosity and radial X-ray brightness profile are determined by three key parameters, namely the so-called X-ray parameter $X$, two radii $R_{rs}$ and $R_c$ of the reverse shock and the contact discontinuity. We find that morphologies of X-ray emissions would appear in the forms of either a central luminous sphere or a bright ring embedded within optically bright shells. In contrast to previous adiabatic models, the X-ray brightness peaks around the reverse shock, instead of the contact discontinuity surface just inside the outer shocked dense shell. Diffuse X-ray emissions of a few observed PNe appear to support this wind-wind dynamic interaction scenario.

We explore self-similar hydrodynamic evolution of central voids embedded in an isothermal gas of spherical symmetry under the self-gravity. More specifically, we study voids expanding at constant radial speeds in an isothermal gas and construct all types of possible void solutions without or with shocks in surrounding envelopes. We examine properties of void boundaries and outer envelopes. Voids without shocks are all bounded by overdense shells and either inflows or outflows in the outer envelope may occur. These solutions, referred to as type $\mathcal{X}$ void solutions, are further divided into subtypes $\mathcal{X}_{\rm I}$ and $\mathcal{X}_{\rm II}$ according to their characteristic behaviours across the sonic critical line (SCL). Void solutions with shocks in envelopes are referred to as type $\mathcal{Z}$ voids and can have both dense and quasi-smooth edges. Asymptotically, outflows, breezes, inflows, accretions and static outer envelopes may all surround such type $\mathcal{Z}$ voids. Both cases of constant and varying temperatures across isothermal shock fronts are analyzed; they are referred to as types $\mathcal{Z}_{\rm I}$ and $\mathcal{Z}_{\rm II}$ void shock solutions. We apply the `phase net matching procedure' to construct various self-similar void solutions. We also present analysis on void generation mechanisms and describe several astrophysical applications. By including self-gravity, gas pressure and shocks, our isothermal self-similar void (ISSV) model is adaptable to various astrophysical systems such as planetary nebulae, hot bubbles and superbubbles in the interstellar medium as well as supernova remnants.

We present new static spherically-symmetric solutions of Einstein equations with the quintessence-like matter surrounding a global monopole. These new solutions of the coupling scalar-Einstein equations are more complicated, which depend on the parameter of equation of state $-1 < w_{q} <-{1/3}$. A gravitating global monopole produces a gravitational field of de Sitter kind outside the core in addition to a solid angular deficit. In the $w_{q} = -{1/3}$ case, we have proved that the solution cannot exist since the density of quintessence-like tends to zero if $w_{q} \to -{1/3}$. As a new feature, these monopoles have the outer horizon depending on both Goldstone field and quintessence-like. Since current observations constrain $-1.14 < w_{q} < -0.93$, new global monopoles have interesting astrophysical applications.

Viscous GCG(generalized Chaplygin gas) cosmology is discussed, assuming that there is bulk viscosity in the barotropic fluid and GCG. The dynamical analysis indicates that the phase $w_g=-1+\sqrt{3}\gamma\kappa\tau_{g}/(\gamma-\sqrt{3}\kappa\tau_{\gamma})$ is a dynamical attractor and the equation of state of GCG approaches it from either $w_g>-1$ or $w_g<-1$ depending on the choice of its initial cosmic density parameter and the ratio of pressure to critical energy density. Obviously, the equation of state $w_g$ can cross the boundary $w_g=-1$. Also, from the point of view of dynamics, the parameters of viscous GCG should be in the range of $\gamma>\sqrt{3}\kappa\tau_{\gamma}/(1-\sqrt{3}\kappa\tau_{g})$ and $0<\alpha<1+\sqrt{3} \kappa\gamma\tau_{g}/(\gamma-\sqrt{3}\kappa\tau_{\gamma}-\sqrt{3}\kappa\gamma\t au_{g})$.

In this paper, we study a realistic model of quintessential inflation with radiation and matter. By the analysis of the dynamical system and numerical work about the evolution of the equation of state and cosmic density parameter, we show that this model is a good match for the current astronomical observation. The conclusion we obtain is in favor of the model where the modular part of complex field plays the role of the inflaton whereas the argument part is the quintessence field. The numerical calculation shows that a heteroclinic orbit (solution of dynamical system) which interpolates between early-time de Sitter phase (an unstable critical point) and a late-time de Sitter attractor.

We study the perturbations on both "radial" and "angular" components of the quintessence with an internal abelian and non-abelian symmetry. The properties of the perturbation on the "radial" component depend on the specific potential of the model and is similiar for both abelian and non-abelian case. We show that the consine-type potential is very interesting for the O(\textit{N}) quintessence model and also give a critical condition of instability for the potential. While the properties of perturbations on "angular" components depend on whether the internal symmetry is abelian or non-abelian, which we have discussed respectively. In the non-abelian case, the fluctuation of the "angular" component will increase rapidly with time while in the abelian case it will not.

We show two analytical solutions of the tachyon inflation for which the spectrum of curvature (density) perturbations can be calculated exactly to linear order, ignoring both gravity and the self-interactions of the tachyon field . The main feature of these solutions is that the spectral indices are independent with scale.