The observed accelerated expansion of the Universe may be explained by dark energy or the breakdown of general relativity (GR) on cosmological scales. When the latter case, a modified gravity scenario, is considered, it is often assumed that the background evolution is the same as the $\Lambda$CDM model but the density perturbation evolves differently. In this paper, we investigate more general classes of modified gravity, where both the background and perturbation evolutions are deviated from those in the $\Lambda$CDM model. We introduce two phase diagrams, $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams; $H$ is the expansion rate, $f\sigma_8$ is a combination of the growth rate of the Universe and the normalization of the density fluctuation which is directly constrained by redshift-space distortions, and $\alpha$ is a parameter which characterizes the deviation of gravity from GR and can be probed by gravitational lensing. We consider several specific examples of Horndeski's theory, which is a general scalar-tensor theory, and demonstrate how deviations from the $\Lambda$CDM model appears in the $\alpha{\rm-}f\sigma _8$ and $H{\rm-}f\sigma _8$ diagrams. The predicted deviations will be useful for future large-scale structure observations to exclude some of the modified gravity models.

A cocoon is a by-product of a propagating jet that results from shock heating at the jet head. Herein, considering simultaneous cocoon formation, we study the stability of relativistic jets propagating through the uniform ambient medium. Using a simple analytic argument, we demonstrate that independent from the jet launching condition, the effective inertia of the jet is larger than that of the cocoon when the fully relativistic jet oscillates radially owing to the pressure mismatch between jet and cocoon. In such situations, it is expected that the onset condition for the oscillation-induced Rayleigh-Taylor instability is satisfied at the jet interface, resulting in the destabilization of the relativistic jet during its propagation. We have quantitatively verified and confirmed our prior expectation by performing relativistic hydrodynamic simulations in three dimensions. The possible occurrences of the Richtmyer-Meshkov instability, oscillation-induced centrifugal instability, and Kelvin-Helmholtz instability are also discussed.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

Recently it was shown that the centrifugal instability may be important in the dynamics of astrophysical jets undergoing reconfinement by external pressure. However, these studies were limited to the case of unmagnetised flows. Here we explore the role of the magnetic field within both the Newtonian and relativistic frameworks. Since the jet problem is rather complicated, we focus instead on the simpler problem of cylindrical rotation and axial magnetic field, which shares significant similarity with the jet problem, and consider only axisymmetric perturbations. The studied equilibrium configurations involve a cylindrical interface and they are stable to non-magnetic centrifugal and magneto-rotational instabilities everywhere except this interface. We use a heuristic approach to derive the local stability criterion for the interface in the magnetic case and numerical simulations to verify the role of the magnetic field. The theory and simulations agree quite well for Newtonian models but indicate a potential discrepancy for the relativistic models in the limit of high Lorentz factor of the rotational motion at the interface. In general, the magnetic field sets a critical wavelength below which the centrifugal modes are stabilised. We discuss the implication of our findings for the astrophysical jets, which suggest that the centrifugal instability develops only in jets with relatively low magnetisation. Namely, the magnetic pressure has to be below the thermal one and for the relativistic case the jets have to be kinetic-energy dominated.

Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discovered in observations. This Yonetoku relation is the tightest correlation found in the properties of the prompt phase of GRB emission, providing the best diagnostic for the radiation mechanism. Here we present 3D hydrodynamical simulations, and post-process radiation transfer calculations, of photospheric emission from a relativistic jet. Our simulations reproduce the Yonetoku relation as a natural consequence of viewing angle. Although jet dynamics depend sensitively on luminosity, the Ep-Lp correlation holds regardless. This result strongly suggests that photospheric emission is the dominant component in the prompt phase of GRBs.

Evolution of the large scale structure of the Universe is described by galaxy bias and the growth rate function of the matter density contrast. The behavior of the growth rate function, in general, depends on the types of dark energy or modified gravity theories, because they modify Poisson's equation for gravity by changing the Newtonian constant in the quasi-static solutions. Whereas, another mode of the matter density contrast "oscillating mode" has not been studied well, although there is a possibility that it drastically changes the gravitational law. In this paper, we will construct the way to investigate oscillating solutions and show the concrete forms of the solutions by using the WKB approximation.

The $\Lambda$CDM model is a remarkably successful model which is consistent with the observations of cosmic microwave background radiation (CMB), baryon acoustic oscillation (BAO), and large scale structure of the Universe. However, the discrepancy in the value of $H_0$ between the local observations and PLANCK observation of CMB was recently pointed out. One of the way to ease the discrepancy is to introduce phantom dark energy instead of the cosmological constant. While, phantom dark energy often suffer from the instabilities. We will investigate the general solution to overcome the difficulty of phantom dark energy and construct some particular models which have a phantom crossing and can be consistent with the observations.

We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $\eta G^{\mu\nu}\phi_{,\mu}\phi_{,\nu}$, where $\eta$ is an arbitrary coupling parameter, and the scalar potential $V(\phi)$ which assumed to be as general as possible. With an appropriate dimensionless parametrization we represent the field equations as an autonomous dynamical system which contains ultimately only one arbitrary function $\chi (x)= 8 \pi \vert \eta \vert V(x/\sqrt{8 \pi})$ with $x=\sqrt{8 \pi}\phi$. Then, assuming the rather general properties of $\chi(x)$, we analyze stationary points and their stability, as well as all possible asymptotical regimes of the dynamical system. It has been shown that for a broad class of $\chi(x)$ there exist attractors representing three accelerated regimes of the Universe evolution, including de Sitter expansion (or late-time inflation), the Little Rip scenario, and the Big Rip scenario. As the specific examples, we consider a power-law potential $V(\phi)=M^4(\phi/\phi_0)^\sigma$, Higgs-like potential $V(\phi)=\frac{\lambda}{4}(\phi^2-\phi_0^2)^2$, and exponential potential $V(\phi)=M^4 e^{-\phi/\phi_0}$.

Dark energy models and modified gravity theories have been actively studied and the behaviors in the solar system have been also carefully investigated in a part of the models. However, the isotropic solutions of the field equations in the simple models of dark energy, e.g. quintessence model without matter coupling, have not been well investigated. One of the reason would be the nonlinearity of the field equations. In this paper, a method to evaluate the solution of the field equations is constructed, and it is shown that there is a model that can easily pass the solar system tests, whereas, there is also a model that is constrained from the solar system tests.

We address the linear stability of a discontinuous surface of a relativistic flow in the context of a jet that oscillates radially as it propagates. The restoring force of the oscillation is expected to drive a Rayleigh-Taylor instability (RTI) at the interface between the jet and its cocoon. We perform a linear analysis and numerical simulations of the growth of the RTI in the transverse plane to the jet flow with a uniform acceleration. In this system, an inertia force due to the uniform acceleration acts as the restoring force for the oscillation. We find that not only the difference in the inertia between the two fluids separated by the interface but also the pressure at the interface helps to drive the RTI because of a difference in the Lorenz factor across the discontinuous surface of the jet. The dispersion relation indicates that the linear growth rate of each mode becomes maximum when the Lorentz factor of the jet is much larger than that of the cocoon and the pressure at the jet interface is relativistic. By comparing the linear growth rates of the RTI in the analytical model and the numerical simulations, the validity of our analytically derived dispersion relation for the relativistic RTI is confirmed.

The existence of dark matter and dark energy in cosmology is implied by various observations, however, they are still unclear because they have not been directly detected. In this Letter, an unified model of dark energy and dark matter that can explain the evolution history of the Universe later than inflationary era, the time evolution of the growth rate function of the matter density contrast, the flat rotation curves of the spiral galaxies, and the gravitational experiments in the solar system is proposed in mimetic matter model.

We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra, even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.

We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $\kappa G^{\mu\nu}\phi_{,\mu}\phi_{,\nu}$, and the Higgs-like potential $V(\phi)=\frac{\lambda}{4}(\phi^2-\phi_0^2)^2$. Using the dynamical system method, we analyze stationary points, their stability, and all possible asymptotical regimes of the model under consideration. We show that the Higgs field with the kinetic coupling provides an existence of accelerated regimes of the Universe evolution. There are three possible cosmological scenarios with acceleration: (i) {\em The late-time inflation} when the Hubble parameter tends to the constant value, $H(t)\to H_\infty=(\frac23 \pi G\lambda\phi_0^4)^{1/2}$ as $t\to\infty$, while the scalar field tends to zero, $\phi(t)\to 0$, so that the Higgs potential reaches its local maximum $V(0)=\frac14 \lambda\phi_0^4$. (ii) {\em The Big Rip} when $H(t)\sim(t_*-t)^{-1}\to\infty$ and $\phi(t)\sim(t_*-t)^{-2}\to\infty$ as $t\to t_*$. (iii) {\em The Little Rip} when $H(t)\sim t^{1/2}\to\infty$ and $\phi(t)\sim t^{1/4}\to\infty$ as $t\to\infty$. Also, we derive modified slow-roll conditions for the Higgs field and demonstrate that they lead to the Little Rip scenario.

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

We investigate the cosmological evolution of mimetic matter model with arbitrary scalar potential. The cosmological reconstruction, which is the way to construct a model for arbitrary evolutions of the scale factor, is explicitly done for different choices of potential. The cases that mimetic matter model shows the evolution as Cold Dark Matter(CDM), wCDM model, dark matter and dark energy with dynamical $Om(z)$, where $Om(z) \equiv [ (H(z)/H_0)^2-1 ]/[(1+z)^3 -1],$ or phantom dark energy with phantom-non-phantom crossing are presented in detail. The cosmological perturbations for such evolution are studied in mimetic matter model. For instance, the evolution behavior of the matter density contrast which is different from usual one, i.e. $\ddot \delta + 2 H \dot \delta - \kappa ^2 \rho \delta /2 = 0$ is investigated. The possibility of peculiar evolution of $\delta$ in the model under consideration is shown. Special attention is paid to the behavior of matter density contrast near to future singularity where decay of perturbations may occur much earlier the singularity.

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($\theta$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

The quasi-static solutions of the matter density perturbation in various dark energy models and modified gravity models have been investigated in numerous papers. However, the oscillating solutions in those models have not been investigated enough so far. In this paper, the oscillating solutions, which have a possibility to unveil the difference between the models of the late-time accelerated expansion of the Universe, are also mentioned by using appropriate approximations.

It is known that the simply evaluated value of the zero point energy of quantum fields is extremely deviated from the observed value of dark energy density. In this paper, we consider whether the Casimir energy, which is the zero point energy brought from boundary conditions, can cause the accelerating expansion of the Universe by using proper renormalization method and introducing the fermions of finite temperature living in $3+n+1$ space-time. We show that the zero temperature Casimir energy and the finite temperature corrections can explain the current accelerated expansion of the Universe.

We study the stability of a non-rotating single-component jet using two-dimensional special relativistic hydrodynamic simulations. By assuming translational invariance along the jet axis, we exclude the destabilization effect by Kelvin-Helmholtz mode. The nonlinear evolution of the transverse structure of the jet with a normal jet velocity is highlighted. An intriguing finding in our study is that Rayleigh-Taylor and Richtmeier-Meshkov type instabilities can destroy cylindrical jet configuration as a result of spontaneously induced radial oscillating motion. This is powered by in-situ energy conversion between the thermal and bulk kinetic energies. The effective inertia ratio of the jet to the surrounding medium $\eta$ determines a threshold for the onset of instabilities. The condition $\eta < 1$ should be satisfied for the transverse structure of the jet being persisted.

The quasi-static solutions of the matter density perturbation in $F(R)$ gravity models have been investigated in numerous papers. However, the oscillating solutions in $F(R)$ gravity models have not been investigated enough so far. In this paper, the oscillating solutions are also examined by using appropriate approximations. And the behaviors of the matter density perturbation in F(R) gravity models with singular evolutions of the physical parameters are shortly investigated as applications of the approximated calculations.

Explicit formulae of the equations in the generalized Galileon models are given. We also develop the formulation of the reconstruction. By using the formulation, we can explicitly construct an action which reproduces an arbitrary development of the expansion of the universe. The conditions how the reconstructed solution becomes stable and therefore it becomes an attractor solution are also given. Working in the static and spherically symmetric space-time, we investigate how the Vainshtein mechanism works in the generalized Galileon model and the correction to the Newton law becomes small. It is also shown that any spherically symmetric and static geometry can be realized by properly choosing the form of the action, which may tell that the solution could have fourth hair corresponding to the scalar field.

The effect of rarefaction acceleration on the propagation dynamics and structure of relativistically hot jets is studied through relativistic hydrodynamic simulations. We emphasize the nonlinear interaction of rarefaction waves excited at the interface between a cylindrical jet and the surrounding medium. From simplified one-dimensional models with radial jet structure, we find that a decrease in the relativistic pressure due to the interacting rarefaction waves in the central zone of the jet transiently yields a more powerful boost of the bulk jet than that expected from single rarefaction acceleration. This leads to a cyclic in-situ energy conversion between thermal and bulk kinetic energies which induces radial oscillating motion of the jet. The oscillation timescale is characterized by the initial pressure ratio of the jet to the ambient medium, and follows a simple scaling relation $\tau_{\rm oscillation} \propto (P_{\rm jet,0}/P_{\rm amb,0})^{1/2}$. It is confirmed from extended two-dimensional simulations that this radial oscillating motion in the one-dimensional system manifests as modulation of the structure of the jet in a more realistic situation where a relativistically hot jet propagates through an ambient medium. It is found that when the ambient medium has a power law pressure distribution, the size of the reconfinement region along the propagation direction of the jet in the modulation structure $\lambda$ evolves according to a self-similar relation $\lambda \propto t^{\alpha/2}$ where $\alpha$ is the power-law index of the pressure distribution.

Subhorizon approximation is often used in cosmological perturbation theory. In this paper, however, it is shown that the subhorizon approximation is not always a good approximation at least in case of $k$-essence model. We also show that the sound speed given by $k$-essence model exerts a huge influence on the time evolution of the matter density perturbation, and the future observations could clarify the differences between the $\Lambda$CDM model and $k$-essence model.

Relativistic Sweet-Parker type magnetic reconnection is investigated by relativistic resistive magnetohydrodynamic (RRMHD) simulations. As an initial setting, we assume anti-parallel magnetic fields and a spatially uniform resistivity. A perturbation imposed on the magnetic fields triggers magnetic reconnection around a current sheet, and the plasma inflows into the reconnection region. The inflows are then heated due to ohmic dissipation in the diffusion region, and finally become relativistically hot outflows. The outflows are not accelerated to ultra-relativistic speeds (i.e., Lorentz factor ~ 1), even when the magnetic energy dominates the thermal and rest mass energies in the inflow region. Most of the magnetic energy in the inflow region is converted into the thermal energy of the outflow during the reconnection process. The energy conversion from magnetic to thermal energy in the diffusion region results in an increase in the plasma inertia. This prevents the outflows from being accelerated to ultra-relativistic speeds. We find that the reconnection rate R obeys the scaling relation R S^{-0.5}, where S is the Lundquist number. This feature is the same as that of non-relativistic reconnection. Our results are consistent with the theoretical predictions of Lyubarsky (2005) for Sweet-Parker type magnetic reconnection.

The differential equations, which are used for matter perturbations, are usually derived from the Newton gravity and the Euler equation in the expanded universe. In this paper, by the explicit calculations of metric perturbation theory in $\Lambda$CDM model, we derive the general relativistic corrections exactly for the equations and show that the differential equations derived from the Newton gravity and the Euler equation are valid in the small scale region from non-relativistic matter dominant era onwards. However the corrections are not negligible in radiation dominant era.

The nonlinear dynamics of outflows driven by magnetic explosion on the surface of a compact star is investigated through special relativistic magnetohydrodynamic simulations. We adopt, as the initial equilibrium state, a spherical stellar object embedded in hydrostatic plasma which has a density $\rho(r) \propto r^{- \alpha}$ and is threaded by a dipole magnetic field. The injection of magnetic energy at the surface of compact star breaks the equilibrium and triggers a two-component outflow. At the early evolutionary stage, the magnetic pressure increases rapidly around the stellar surface, initiating a magnetically driven outflow. A strong forward shock driven outflow is then excited. The expansion velocity of the magnetically driven outflow is characterized by the Alfv\'en velocity on the stellar surface, and follows a simple scaling relation $v_{\rm mag} \propto {v_{\rm A}}^{1/2}$. When the initial density profile declines steeply with radius, the strong shock is accelerated self-similarly to relativistic velocity ahead of the magnetically driven component. We find that it evolves according to a self-similar relation $\Gamma_{\rm sh} \propto r_{\rm sh}$, where $\Gamma_{\rm sh}$ is the Lorentz factor of the plasma measured at the shock surface $r_{\rm sh}$. Purely hydrodynamic process would be responsible for the acceleration mechanism of the shock driven outflow. Our two-component outflow model, which is the natural outcome of the magnetic explosion, can provide a better understanding of the magnetic active phenomena on various magnetized compact stars.

We study scalar-tensor theory, k-essence and modified gravity with Lagrange multiplier constraint which role is to reduce the number of degrees of freedom. Dark Energy cosmology of different types ($\Lambda$CDM, unified inflation with DE, smooth non-phantom/phantom transition epoch) is reconstructed in such models. It is shown that mathematical equivalence between scalar theory and $F(R)$ gravity is broken due to presence of constraint. The cosmological dynamics of $F(R)$ gravity is modified by the second $F_2(R)$ function dictated by the constraint. Dark Energy cosmology is defined by this function while standard $F_1(R)$ function is relevant for local tests (modification of newton regime). A general discussion on the role of Lagrange multipliers to make higher-derivative gravity canonical is developed.