Showing 1–34 of 34 results for author: Matsumoto, J

Three-dimensional Magneto-hydrodynamic Simulations of Core-collapse Supernovae: I. Hydrodynamic evolution and protoneutron star properties

Authors: Ko Nakamura, Tomoya Takiwaki, Jin Matsumoto, Kei Kotake

Abstract: We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of sixteen progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 $M_{\odot}$. The examined progenitors cover a wide range of the compactness parameter including a peak around $23 M_{\odot}$. We find that neutri… ▽ More We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of sixteen progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 $M_{\odot}$. The examined progenitors cover a wide range of the compactness parameter including a peak around $23 M_{\odot}$. We find that neutrino-driven explosions occur for all models within 0.3 s after bounce. We also find that the properties of the explosions and the central remnants are well correlated with the compactness. Early shock evolution is sensitive to the mass accretion rate onto the central core, reflecting the density profile of the progenitor stars. The most powerful explosions with diagnostic explosion energy $E_{\rm exp} \sim 0.75 \times 10^{51}$ erg are obtained by 23 and 24 $M_{\odot}$ models, which have the highest compactness among the examined models. These two models exhibit spiral SASI motions during 150-230 ms after bounce preceding a runaway shock expansion and leave a rapidly rotating neutron star with spin periods $\sim 50$ ms. Our models predict the gravitational masses of the neutron star ranging between $1.22 M_{\odot}$ and $1.67 M_{\odot}$ and their spin periods 0.04-4 s. The number distribution of these values roughly matches observation. On the other hand, our models predict small hydrodynamic kick velocity (15-260 km/s), although they are still growing at the end of our simulations. Further systematic studies, including rotation and binary effects, as well as long-term simulations up to several seconds, will enable us to explore the origin of various core-collapse supernova explosions. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 15 pages, 13+5 figures, 1 appendix

Numerical simulation of photospheric emission in long gamma-ray bursts: prompt correlations, spectral shapes, and polarizations

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov, Daisuke Yonetoku

Abstract: We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the… ▽ More We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the observed Yonetoku, Amati, and Golenetskii correlations. It is also shown that the spectral peak energy and luminosity correlate with the bulk Lorentz factor, as indicated in the literature. On the other hand, synthetic spectral shapes tend to be narrower than those of the observations. The result indicates that an additional physical process that can provide non-thermal broadening is needed to reproduce the spectral features. Furthermore, the polarization analysis finds that, while the degree of polarization is low for the emission from the jet core ($Π< 4~\%$), it tends to increase with the viewing angle outside the core and can be as high as $Π\sim 20-40~\%$ in an extreme case. This suggests that the typical GRBs show systematically low polarization compared to softer, dimmer counterparts (X-ray-rich GRBs and X-ray flashes). Interestingly, our simulations indicate that photospheric emission exhibits large temporal variation in the polarization position angle ($Δψ\sim 90^{\circ}$), which may be compatible with those inferred in observations. A notable energy dependence of the polarization property is another characteristic feature found in the current study. Particularly, the difference in the position angle among different energy bands can be as large as $\sim 90^{\circ}$. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 20 pages, 12 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-23

Neutrino-driven massive stellar explosions in 3D fostered by magnetic fields via turbulent $α$-effect

Authors: Jin Matsumoto, Tomoya Takiwaki, Kei Kotake

Abstract: We investigate the influence of magnetic field amplification on the core-collapse supernovae in highly magnetized progenitors through three-dimensional simulations. By considering rotating models, we observe a strong correlation between the exponential growth of the magnetic field in the gain region and the initiation of shock revival, with a faster onset compared to the non-rotating model. We hig… ▽ More We investigate the influence of magnetic field amplification on the core-collapse supernovae in highly magnetized progenitors through three-dimensional simulations. By considering rotating models, we observe a strong correlation between the exponential growth of the magnetic field in the gain region and the initiation of shock revival, with a faster onset compared to the non-rotating model. We highlight that the mean magnetic field experiences exponential amplification as a result of $α$-effect in the dynamo process, which works efficiently with the increasing kinetic helicity of the turbulence within the gain region. Our findings indicate that the significant amplification of the mean magnetic fields leads to the development of locally intense turbulent magnetic fields, particularly in the vicinity of the poles, thereby promoting the revival of the shock by neutrino heating. △ Less

Submitted 14 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

Chiral plasma instability and inverse cascade from nonequilibrium left-handed neutrinos in core-collapse supernovae

Authors: Jin Matsumoto, Naoki Yamamoto, Di-Lun Yang

Abstract: We show that the backreaction of left-handed neutrinos out of equilibrium on the matter sector induces an electric current proportional to a magnetic field even without a chiral imbalance for electrons in core-collapse supernovae. We derive the transport coefficient of this effect based on the recently formulated chiral radiation transport theory for neutrinos. This chiral electric current generat… ▽ More We show that the backreaction of left-handed neutrinos out of equilibrium on the matter sector induces an electric current proportional to a magnetic field even without a chiral imbalance for electrons in core-collapse supernovae. We derive the transport coefficient of this effect based on the recently formulated chiral radiation transport theory for neutrinos. This chiral electric current generates a strong magnetic field via the so-called chiral plasma instability, which could provide a new mechanism for the strong and stable magnetic field of magnetars. We also numerically study the physical origin of the inverse cascade of the magnetic energy in the magnetohydrodynamics including this current. Our results indicate that incorporating the chiral effects of neutrinos would drastically modify the hydrodynamic evolutions of supernovae, which may also be relevant to the explosion dynamics. △ Less

Submitted 26 June, 2022; v1 submitted 18 February, 2022; originally announced February 2022.

Comments: 8 pages, 6 figures; v2: added figure 2 on the time evolution of magnetic helicity, minor revisions, published version

Journal ref: Phys. Rev. D 105, 123029 (2022)

Magnetic support for neutrino-driven explosion of 3D non-rotating core-collapse supernova models

Authors: Jin Matsumoto, Yuta Asahina, Tomoya Takiwaki, Kei Kotake, Hiroyuki R. Takahashi

Abstract: The impact of the magnetic field on postbounce supernova dynamics of non-rotating stellar cores is studied by performing three-dimensional magnetohydrodynamics simulations with spectral neutrino transport. The explodability of strongly and weakly magnetized models of $20$ and $27$ $M_{\odot}$ pre-supernova progenitors are compared. We find that although the efficiency for the conversion of the neu… ▽ More The impact of the magnetic field on postbounce supernova dynamics of non-rotating stellar cores is studied by performing three-dimensional magnetohydrodynamics simulations with spectral neutrino transport. The explodability of strongly and weakly magnetized models of $20$ and $27$ $M_{\odot}$ pre-supernova progenitors are compared. We find that although the efficiency for the conversion of the neutrino heating into turbulent energy including magnetic fields in the gain region is not significantly different between the strong and weak field models, the amplified magnetic field due to the neutrino-driven convection on large hot bubbles just behind stalled shock results in a faster and more energetic explosion in the strongly magnetized models. In addition, by comparing the difference between the 2nd- and 5th-order spatial accuracy of the simulation in the strong field model for $27$ $M_{\odot}$ progenitor, we also find that the higher order accuracy in space is beneficial to the explosion because it enhances the growth of neutrino-driven convection in the gain region. Based on our results of core-collapse supernova simulations for the non-rotating model, a new possibility for the origin of the magnetic field of the protoneutron star (PNS) is proposed. The magnetic field is accumulated and amplified to magnetar level, that is, $\mathcal{O}(10^{14})$ G, in the convectively stable shell near the PNS surface. △ Less

Submitted 18 August, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 16 pages, 21 figures, accepted for publication in MNRAS

Magnetic inhibition of the recollimation instability in relativistic jets

Authors: Jin Matsumoto, Serguei S. Komissarov, Konstantinos N. Gourgouliatos

Abstract: In this paper, we describe the results of three-dimensional relativistic magnetohydrodynamic simulations aimed at probing the role of regular magnetic field on the development of the instability that accompanies recollimation of relativistic jets. In particular, we studied the recollimation driven by the reconfinement of jets from active galactic nuclei (AGN) by the thermal pressure of galactic co… ▽ More In this paper, we describe the results of three-dimensional relativistic magnetohydrodynamic simulations aimed at probing the role of regular magnetic field on the development of the instability that accompanies recollimation of relativistic jets. In particular, we studied the recollimation driven by the reconfinement of jets from active galactic nuclei (AGN) by the thermal pressure of galactic coronas. We find that a relatively weak azimuthal magnetic field can completely suppress the recollimation instability in such jets, with the critical magnetisation parameter $σ_{cr} < 0.01$. We argue that the recollimation instability is a variant of the centrifugal instability (CFI) and show that our results are consistent with the predictions based on the study of magnetic CFI in rotating fluids. The results are discussed in the context of AGN jets in general and the nature of the Fanaroff-Riley morphological division of extragalactic radio sources in particular. △ Less

Submitted 18 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

Two-dimensional numerical study for magnetic field dependence of neutrino-driven core-collapse supernova models

Authors: Jin Matsumoto, Tomoya Takiwaki, Kei Kotake, Yuta Asahina, Hiroyuki R. Takahashi

Abstract: We study the effects of the magnetic field on the dynamics of non-rotating stellar cores by performing two-dimensional (2D), magnetohydrodynamics (MHD) simulations. To this end, we have updated our neutrino-radiation-hydrodynamics supernova code to include MHD employing a divergence cleaning method with both careful treatments of finite volume and area reconstructions. By changing the initial stre… ▽ More We study the effects of the magnetic field on the dynamics of non-rotating stellar cores by performing two-dimensional (2D), magnetohydrodynamics (MHD) simulations. To this end, we have updated our neutrino-radiation-hydrodynamics supernova code to include MHD employing a divergence cleaning method with both careful treatments of finite volume and area reconstructions. By changing the initial strength of the magnetic field, the evolution of $15.0$, $18.4$ and $27.0$ $M_\odot$ presupernova progenitors is investigated. An intriguing finding in our study is that the neutrino-driven explosion occurs regardless of the strength of the initial magnetic field. For the 2D models presented in this work, the neutrino heating is the main driver for the explosion, whereas the magnetic field secondary contributes to the pre-explosion dynamics. Our results show that the strong magnetic field weakens the growth of the neutrino-driven turbulence in the small scale compared to the weak magnetic field. This results in the slower increase of the turbulent kinetic energy in the postshock region, leading to the slightly delayed onset of the shock revival for models with the stronger initial magnetic field. △ Less

Submitted 6 October, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: 23 pages, 10 figures, accepted for publication in MNRAS

New measures to test modified gravity cosmologies

Authors: Jiro Matsumoto, Teppei Okumura, Misao Sasaki

Abstract: 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 $Λ$CDM model but the density perturbation evolves differently. In this paper, we investigate more general classes… ▽ More 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 $Λ$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 $Λ$CDM model. We introduce two phase diagrams, $α{\rm-}fσ_8$ and $H{\rm-}fσ_8$ diagrams; $H$ is the expansion rate, $fσ_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 $α$ 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 $Λ$CDM model appears in the $α{\rm-}fσ_8$ and $H{\rm-}fσ_8$ diagrams. The predicted deviations will be useful for future large-scale structure observations to exclude some of the modified gravity models. △ Less

Submitted 1 July, 2020; v1 submitted 19 May, 2020; originally announced May 2020.

Comments: 18 pages, 8 figures, accepted to JCAP

Report number: YITP-20-74

Journal ref: JCAP 07 (2020) 059

Propagation, cocoon formation, and resultant destabilization of relativistic jets

Authors: Jin Matsumoto, Youhei Masada

Abstract: 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 co… ▽ More 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. △ Less

Submitted 25 October, 2019; originally announced October 2019.

Comments: 11 pages, 6 figures. Accepted for publication in MNRAS

Magnetic Inhibition of Centrifugal Instability in Astrophysical Jets

Authors: Serguei S. Komissarov, Konstantinos N. Gourgouliatos, Jin Matsumoto

Abstract: 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 s… ▽ More 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. △ Less

Submitted 20 June, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: submitted to MNRAS

Oscillating solutions of the matter density contrast in Horndeski's theory

Authors: Jiro Matsumoto

Abstract: 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 mo… ▽ More 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. △ Less

Submitted 6 January, 2019; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: 21 pages, 6 figures. References added, sentences and expressions polished, version to appear in JCAP

Journal ref: JCAP01(2019)054

The photospheric origin of the Yonetoku relation in gamma-ray bursts

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov, Daisuke Yonetoku

Abstract: 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 discove… ▽ More 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. △ Less

Submitted 3 April, 2019; v1 submitted 2 June, 2018; originally announced June 2018.

Comments: 18 pages, 5 figures, The original submitted version of an article published in Nature Communications. The final authenticated version is available online at: https://www.nature.com/articles/s41467-019-09281-z

Report number: RIKEN-iTHEMS-Report-18

Journal ref: Nature Communications volume 10, Article number: 1504 (2019)

Phantom crossing dark energy in Horndeski's theory

Authors: Jiro Matsumoto

Abstract: The $Λ$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 int… ▽ More The $Λ$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. △ Less

Submitted 24 May, 2018; v1 submitted 28 December, 2017; originally announced December 2017.

Comments: 16 pages, 5 figures. References added, sentences and expressions polished, version to appear in PRD

Journal ref: Phys. Rev. D 97, 123538 (2018)

Linear Theory of the Rayleigh-Taylor Instability at a Discontinuous Surface of a Relativistic Flow

Authors: Jin Matsumoto, Miguel A. Aloy, Manel Perucho

Abstract: 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… ▽ More 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. △ Less

Submitted 2 August, 2017; v1 submitted 15 July, 2017; originally announced July 2017.

Comments: 11 pages, 6 figures. Accepted for publication in MNRAS

General dynamical properties of cosmological models with nonminimal kinetic coupling

Authors: Jiro Matsumoto, Sergey V. Sushkov

Abstract: We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $ηG^{μν}φ_{,μ}φ_{,ν}$, where $η$ is an arbitrary coupling parameter, and the scalar potential $V(φ)$ which assumed to be as general as possible. With an appropriate dimensionless parametrization we represent the field equations as an autonomous dynamical… ▽ More We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $ηG^{μν}φ_{,μ}φ_{,ν}$, where $η$ is an arbitrary coupling parameter, and the scalar potential $V(φ)$ 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 $χ(x)= 8 π\vert η\vert V(x/\sqrt{8 π})$ with $x=\sqrt{8 π}φ$. Then, assuming the rather general properties of $χ(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 $χ(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(φ)=M^4(φ/φ_0)^σ$, Higgs-like potential $V(φ)=\fracλ{4}(φ^2-φ_0^2)^2$, and exponential potential $V(φ)=M^4 e^{-φ/φ_0}$. △ Less

Submitted 7 January, 2018; v1 submitted 15 March, 2017; originally announced March 2017.

Comments: 31 pages, 5 figures, 3 tables. New section "Slow-roll regime" added, references added, sentences polished. Version accepted for publication in JCAP

Journal ref: JCAP01(2018)040

Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications

Authors: Yosuke Matsumoto, Yuta Asahina, Yuki Kudoh, Tomohisa Kawashima, Jin Matsumoto, Hiroyuki R. Takahashi, Takashi Minoshima, Seiji Zenitani, Takahiro Miyoshi, Ryoji Matsumoto

Abstract: 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 magne… ▽ More 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. △ Less

Submitted 18 August, 2019; v1 submitted 6 November, 2016; originally announced November 2016.

Comments: 45 pages, 18 figures, published in Publ. Astron. Soc. Japan

Journal ref: Publications of the Astronomical Society of Japan, Volume 71, Issue 4, 2019

Unified description of dark energy and dark matter in mimetic matter model

Authors: Jiro Matsumoto

Abstract: 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, t… ▽ More 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. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Comments: 5 pages, 3 figures

Scalar field dark energy with a minimal coupling in a spherically symmetric background

Authors: Jiro Matsumoto

Abstract: 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 t… ▽ More 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. △ Less

Submitted 9 August, 2017; v1 submitted 11 July, 2016; originally announced July 2016.

Comments: 15 pages, 3 figures; v2: equations corrected and a figure added, conclusions are changed a little, v3: references added, explanations expanded, version to appear in IJMPD

Journal ref: Int. J. Mod. Phys. D, 27, 1750173 (2018)

Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 27 November, 2015; v1 submitted 11 November, 2015; originally announced November 2015.

Comments: 7 pages, 4 figures, accepted for publication in ApJ Letters

Cosmology with nonminimal kinetic coupling and a Higgs-like potential

Authors: Jiro Matsumoto, Sergey V. Sushkov

Abstract: We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $κG^{μν}φ_{,μ}φ_{,ν}$, and the Higgs-like potential $V(φ)=\fracλ{4}(φ^2-φ_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… ▽ More We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as $κG^{μν}φ_{,μ}φ_{,ν}$, and the Higgs-like potential $V(φ)=\fracλ{4}(φ^2-φ_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 πGλφ_0^4)^{1/2}$ as $t\to\infty$, while the scalar field tends to zero, $φ(t)\to 0$, so that the Higgs potential reaches its local maximum $V(0)=\frac14 λφ_0^4$. (ii) {\em The Big Rip} when $H(t)\sim(t_*-t)^{-1}\to\infty$ and $φ(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 $φ(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. △ Less

Submitted 10 November, 2015; v1 submitted 12 October, 2015; originally announced October 2015.

Comments: 29 pages, 11 figures, discussions and references added, to be published on JCAP

Journal ref: JCAP 1511, no. 11, 047 (2015)

Matter Mixing in Core-collapse Supernova Ejecta: Large Density Perturbations in the Progenitor Star?

Authors: J. Mao, M. Ono, S. Nagataki, M. Hashimoto, H. Ito, J. Matsumoto, M. G. Dainotti, S. -H. Lee

Abstract: 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 mi… ▽ More 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. △ Less

Submitted 24 July, 2015; originally announced July 2015.

Comments: ApJ accepted

Cosmological perturbations in mimetic matter model

Authors: Jiro Matsumoto, Sergei D. Odintsov, Sergey V. Sushkov

Abstract: 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 dy… ▽ More 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 δ+ 2 H \dot δ- κ^2 ρδ/2 = 0$ is investigated. The possibility of peculiar evolution of $δ$ 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. △ Less

Submitted 7 March, 2015; v1 submitted 9 January, 2015; originally announced January 2015.

Comments: 20 pages, 7 figures, version to appear in PRD

Journal ref: Phys. Rev. D 91, 064062 (2015)

Spectral and Polarization Properties of Photospheric Emission From Stratified Jets

Authors: Hirotaka Ito, Shigehiro Nagataki, Jin Matsumoto, Shiu-Hang Lee, Alexey Tolstov, Jirong Mao, Maria Dainotti, Akira Mizuta

Abstract: 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 je… ▽ More 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 ($θ$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%. △ Less

Submitted 24 May, 2014; originally announced May 2014.

Comments: 21 pages, 12 figures, accepted for publication in ApJ

Cosmological perturbations in the models of dark energy and modified gravity

Authors: Jiro Matsumoto

Abstract: 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 expa… ▽ More 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. △ Less

Submitted 14 January, 2014; originally announced January 2014.

Comments: 7 pages, to appear in the proceedings of The 10th International Symposium on Cosmology and Particle Astrophysics (CosPA2013)

Journal ref: Universe 1, no. 1, 17 (2015)

Two-dimensional Numerical Study for Rayleigh-Taylor and Richtmyer-Meshkov Instabilities in Relativistic Jets

Authors: Jin Matsumoto, Youhei Masada

Abstract: 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… ▽ More 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 $η$ determines a threshold for the onset of instabilities. The condition $η< 1$ should be satisfied for the transverse structure of the jet being persisted. △ Less

Submitted 5 June, 2013; originally announced June 2013.

Comments: 6 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

Cosmological perturbations in $F(R)$ gravity

Authors: Jiro Matsumoto

Abstract: 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 mod… ▽ More 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. △ Less

Submitted 2 May, 2013; v1 submitted 27 March, 2013; originally announced March 2013.

Comments: 15 pages, note and reference added and typos corrected

Journal ref: Phys.Rev.D87:104002,2013

The Casimir effect as a candidate of dark energy

Authors: Jiro Matsumoto

Abstract: 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… ▽ More 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. △ Less

Submitted 25 December, 2013; v1 submitted 17 March, 2013; originally announced March 2013.

Comments: 14 pages, 4 figures. v2: the expression of the energy density revised, figures added, v3; figures and typos are corrected

Effect of Interacting Rarefaction Waves on Relativistically Hot Jets

Authors: Jin Matsumoto, Youhei Masada, Kazunari Shibata

Abstract: 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… ▽ More 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 $τ_{\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 $λ$ evolves according to a self-similar relation $λ\propto t^{α/2}$ where $α$ is the power-law index of the pressure distribution. △ Less

Submitted 25 April, 2012; originally announced April 2012.

Comments: 38 pages, 21 figures. Accepted for publication in ApJ

Generalized Galileon Model -- cosmological reconstruction and the Vainshtein mechanism --

Authors: Norihito Shirai, Kazuharu Bamba, Shota Kumekawa, Jiro Matsumoto, Shin'ichi Nojiri

Abstract: 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.… ▽ More 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. △ Less

Submitted 20 June, 2012; v1 submitted 22 March, 2012; originally announced March 2012.

Comments: 27 pages, no figure, version to appear in Physical Review D

Cosmological perturbations in $k$-essence model

Authors: Kazuharu Bamba, Jiro Matsumoto, Shin'ichi Nojiri

Abstract: 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 cla… ▽ More 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 $Λ$CDM model and $k$-essence model. △ Less

Submitted 16 March, 2012; v1 submitted 6 September, 2011; originally announced September 2011.

Comments: 21 pages, sentences and equations are corrected, conclusions are changed a little

Journal ref: Phys.Rev.D85:084026,2012

Scaling Law of Relativistic Sweet-Parker Type Magnetic Reconnection

Authors: Hiroyuki R. Takahashi, Takahiro Kudoh, Youhei Masada, Jin Matsumoto

Abstract: 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 inflow… ▽ More 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. △ Less

Submitted 19 August, 2011; originally announced August 2011.

Comments: accepted for publication in ApJL, 6 pages, 4 figures

On the necessity of the revisions for the cosmological matter perturbations from the general relativity

Authors: Jiro Matsumoto

Abstract: 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 $Λ$CDM model, we derive the general relativistic corrections exactly for the equations and show that the differential equations derived from the Newton gravity a… ▽ More 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 $Λ$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. △ Less

Submitted 1 June, 2011; v1 submitted 7 May, 2011; originally announced May 2011.

Comments: 8 pages, LaTeX; the arguments in section 3 are corrected and the results are changed

Journal ref: Phys.Rev.D83:124040,2011

Special Relativistic Magnetohydrodynamic Simulation of Two-Component Outflow Powered by Magnetic Explosion on Compact Stars

Authors: Jin Matsumoto, Youhei Masada, Eiji Asano, Kazunari Shibata

Abstract: 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 $ρ(r) \propto r^{- α}$ and is threaded by a dipole magnetic field. The injection of magnetic energy at… ▽ More 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 $ρ(r) \propto r^{- α}$ 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én 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 $Γ_{\rm sh} \propto r_{\rm sh}$, where $Γ_{\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. △ Less

Submitted 10 March, 2011; originally announced March 2011.

Comments: Accepted for publication in ApJ. 15 pages, 2 tables, 17 figures

Dark energy from modified gravity with Lagrange multipliers

Authors: Salvatore Capozziello, Jiro Matsumoto, Shin'ichi Nojiri, Sergei D. Odintsov

Abstract: 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 ($Λ$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 bro… ▽ More 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 ($Λ$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. △ Less

Submitted 11 August, 2010; v1 submitted 21 April, 2010; originally announced April 2010.

Comments: LaTeX 12 pages, discussion is improved

Journal ref: Phys.Lett.B693:198-208,2010