The effect of pasta phases on the quark-hadron phase transition is investigated for a set of relativistic mean-field equations of state for both hadron and quark matter. The results of the full numerical solution with pasta phases are compared with those of an interpolating construction used in previous works, for which we demonstrate an adequate description of the numerical results. A one-to-one mapping of the free parameter of the construction to the physical surface tension of the quark-hadron interface is obtained for which a fit formula is given. For each pair of quark and hadron matter models the critical value of the surface tension is determined, above which the phase transition becomes close to the Maxwell construction. This result agrees well with earlier theoretical estimates. The study is extended to neutron star matter in beta equilibrium with electrons and muons and is applied to investigate the effect of pasta phases on the structure of hybrid compact stars and the robustness of a possible third family solution.

We study the interface effects of quark-hadron mixed phase in compact stars. The properties of nuclear matter are obtained based on the relativistic-mean-field model. For the quark phase, we adopt perturbation model with running quark masses and coupling constant. At certain choices of parameter sets, it is found that varying the quark-hadron interface tension will have sizable effects on the radii ($\Delta R \approx 600$ m) and tidal deformabilities ($\Delta \Lambda/\Lambda \approx 50\%$) of hybrid stars. These provide possibilities for us to constrain the quark-hadron interface tension with future gravitational wave observations as well as the ongoing NICER mission.

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method.

The existence of neutron stars with $2M_\odot$ requires the strong stiffness of the equation of state (EoS) of neutron-star matter. We introduce a multi-pomeron exchange potential (MPP) working universally among 3- and 4-baryons to stiffen the EoS. Its strength is restricted by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoSs are derived using the Brueckner-Hartree-Fock (BHF) and the cluster variational method (CVM) with the nuclear interactions ESC and AV18. The mass-radius relations are derived by solving the Tolmann-Oppenheimer-Volkoff (TOV) equation, where the maximum masses over $2M_\odot$ are obtained on the basis of the terrestrial data. Neutron-star radii $R$ at a typical mass $1.5M_\odot$ are predicted to be $12.3\!\sim\!13.0$ km. The uncertainty of calculated radii is mainly from the ratio of 3- and 4-pomeron coupling constants, which cannot be fixed by any terrestrial experiment. Though values of $R(1.5M_\odot)$ are not influenced by hyperon-mixing effects, finely-observed values for them indicate degrees of EoS softening by hyperon mixing in the region of $M\!\sim\!2M_\odot$. If $R(1.5M_\odot)$ is less than about 12.4 km, the softening of EoS by hyperon mixing has to be weak. Useful information can be expected by the space mission NICER offering precise measurements for neutron-star radii within $\pm 5\%$.

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme for the barotropic case, and those by Fujisawa's self-consistent field scheme for the baroclinic case. Included in these rotational configurations are those with shellular-type rotations, which are commonly assumed in the evolution calculation of rotating stars. Although radiation processes, convections and meridional flows have not been taken into account in this study, we have in mind the application of this method to the two-dimensional evolution calculations of rotating stars, for which the Lagrangian formulation is best suited.

We show a scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star with high effective temperature, and it is consistent with the cooling without exotic phases. The Cas A observation also gives the mass range of $M \geq 1.5 M_\odot$. It may conflict with the current cooling scenarios of compact stars that heavy stars show rapid cooling. We include the effect of the color superconducting (CSC) quark matter phase on the thermal evolution of compact stars. We assume the gap energy of CSC quark phase is large ($\Delta \gtrsim \mathrm{10 MeV}$), and we simulate the cooling of compact stars. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency.

Neutron stars with large masses $\sim 2M_{\odot}$ require the hard stiffness of equation of state (EoS) of neutron-star matter. On the other hand, hyperon mixing brings about remarkable softening of EoS. In order to solve this problem, a multi-pomeron exchange potential (MPP) is introduced as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) baryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interactions in $\Lambda\!N$, $\Sigma\!N$ and $\Xi\!N$ channels are shown to be consistent with experimental indications. The EoS in neutron-star matter with hyperon mixing is obtained from ESC in addition of MPP, and mass-radius relations of neutron stars are derived. The maximum mass is shown to reach $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

We study the quark-hadron phase transition with the finite-size effects in neutron stars. The finite-size effects should be, generally, taken into account in the phase transition of multi-component system. The behavior of the phase transition, however, strongly depends on the models for quark and hadron matter, surface tension, neutrino fraction, and temperature. We find that, if the surface tension is strong, the EOS becomes similar to the case of a Maxwell construction for any hadron and/or quark model, though we adopt the Gibbs conditions. We also find that the mass-radius relations for that EOS are consistent with the observations, and our model is then applicable to realistic astrophysical phenomena such as the thermal evolution of compact stars.

We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle, and treats not only barotropic but also baroclinic equations of state, for which angular momentum distributions are not necessarily cylindrical. We adopt a Monte Carlo technique, which is analogous to those employed in other fields, e.g. nuclear physics, in minimizing the energy functional, which is evaluated on a triangulated mesh. This letter is a proof of principle and detailed comparisons with existing results will be reported in the sequel, but some test calculations are presented, in which we have achieved an error of $O(10^{-4})$ in the Virial relation. We have in mind the application of this method to two-dimensional calculations of the evolutions of rotating stars, for which the Lagrangian formulation is best suited.

We study the role of finite-size effects at the hadron-quark phase transition in a new hybrid equation of state constructed from an ab-initio Br\"uckner-Hartree-Fock equation of state with the realistic Bonn-B potential for the hadronic phase and a covariant non-local Nambu--Jona-Lasinio model for the quark phase. We construct static hybrid star sequences and find that our model can support stable hybrid stars with an onset of quark matter below $2 M_\odot$ and a maximum mass above $2.17 M_\odot$ in agreement with recent observations. If the finite-size effects are taken into account the core is composed of pure quark matter. Provided that the quark vector channel interaction is small, and the finite size effects are taken into account, quark matter appears at densities 2-3 times the nuclear saturation density. In that case the proton fraction in the hadronic phase remains below the value required by the onset of the direct URCA process, so that the early onset of quark matter shall affect on the rapid cooling of the star.

A multi-pomeron exchange potential (MPP) is proposed as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) bryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interaction in $\Lambda N$ channels is shown to reproduce well the experimental $\Lambda$ binding energies. The equation of state (EoS) in neutron matter with hyperon mixing is obtained including the MPP contribution, and mass-radius relations of neutron stars are derived. It is shown that the maximum mass can be larger than the observed one $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

We calculate the spin deceleration of proto-neutron stars by asymmetric neutrino absorption in the context of a fully relativistic mean field theory. We calculate for the first time the spin deceleration of neutron stars due to asymmetric neutrino absorption in a toroidal magnetic field configurations. We find a surprising effect that the deceleration can much larger for asymmetric neutrino absorption in a toroidal magnetic field than the usually presumed braking due to magnetic dipole radiation. This may explain the observation that magnetars appear to have had a more rapid deceleration in the past.

Bearing in mind the application to high-magnetic-field (high-B) radio pulsars, we investigate two-dimensional (2D) thermal evolutions of neutron stars (NSs). We pay particular attention to the influence of different equilibrium configurations on the surface temperature distributions. The equilibrium configurations are constructed in a systematic manner, in which both toroidal and poloidal magnetic fields are determined self-consistently with the inclusion of general relativistic effects. To solve the 2D heat transfer inside the NS interior out to the crust, we have developed an implicit code based on a finite-difference scheme that deals with anisotropic thermal conductivity and relevant cooling processes in the context of a standard cooling scenario. In agreement with previous studies, the surface temperatures near the pole become higher than those in the vicinity of the equator as a result of anisotropic heat transfer. Our results show that the ratio of the highest to the lowest surface temperatures changes maximally by one order of magnitude, depending on the equilibrium configurations. Despite such difference, we find that the area of such hot and cold spots is so small that the simulated X-ray spectrum could be well reproduced by a single temperature blackbody fitting.

We study the hadron-quark phase transition, taking into account the finite-size effects for neutron star matter. For the hadron phase, we adopt a realistic equation of state within the framework of the Brueckner-Hartree-Fock theory. For the quark phase, we apply the Dyson-Schwinger method. The properties of the mixed phase are clarified by considering the finite-size effects. We find that, if the surface tension is strong enough, the equation of state becomes to be close the one with the Maxwell condition, though we properly adopt the Gibbs conditions. This result is qualitatively the same with the one by the use of the simple bag model. We also find that the mass-radius relation by the EoS is consistent with the observations of massive neutron stars.

A multi-pomeron exchange potential (MPP) is proposed as a model for the three-body repulsion indicated in neutron-star matter, which works universally among three- and four-baryons. Its strength is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoS in neutron matter is obtained including the MPP contribution. The neutron-star mass is calculated by solving the TOV equation. The maximum mass is obtained to be larger than the observed one $1.97 M_{solar}$ on the basis of the experimental data.

We determine the quark-hadron transition line in the whole region of temperature (T) and baryon-number chemical potential (\mu_B) from lattice QCD results and neutron-star mass measurements, making the quark-hadron hybrid model that is consistent with the two solid constraints. The quark part of the hybrid model is the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with entanglement vertex that reproduces lattice QCD results at mu_B/T=0, while the hadron part is the hadron resonance gas model with volume-exclusion effect that reproduces neutron-star mass measurements and the neutron-matter equation of state calculated from two- and three-nucleon forces based on the chiral effective field theory. The lower bound of the critical mu_B of the quark-hadron transition at zero T is mu_B = 1.6 GeV. The interplay between the heavy-ion collision physics around mu_B/T =6 and the neutron-star physics where mu_B/T is infinity is discussed.

We study the hadron-quark phase transition in the interior of hot protoneutron stars, combining the Brueckner-Hartree-Fock approach for hadronic matter with the MIT bag model or the Dyson-Schwinger model for quark matter. We examine the structure of the mixed phase constructed according to different prescriptions for the phase transition, and the resulting consequences for stellar properties. We find important effects for the internal composition, but only very small influence on the global stellar properties.

We calculate asymmetric neutrino absorption and scattering cross sections on hot and dense magnetized neutron-star matter including hyperons in fully relativistic mean field theory. The absorption/scattering cross sections are suppressed/enhanced incoherently in the direction of the magnetic fielld, B. The asymmetry is 2-4% at the matter density \rho_0 < \rho_B < 3\rho_0 and temperature T < 40 MeV for B=2 X 10^{17} G. This asymmetry is comparable to the effects owing to parity violation or asymmetric magnetic field topology proposed for the origin of pulsar kicks.

We present a new scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star that has high effective temperature, and it is consistent with the cooling without exotic phases. The observation also gives the mass range of $M \geqslant 1.5 M_\odot$, which may conflict with the current plausible cooling scenario of compact stars. There are some cooled compact stars such as Vela or 3C58, which can be barely explained by the minimal cooling scenario, which includes the neutrino emission by nucleon superfluidity (PBF). Therefore, we invoke the exotic cooling processes, where a heavier star cools faster than lighter one. However, the scenario seems to be inconsistent with the observation of Cas A. Therefore, we present a new cooling scenario to explain the observation of Cas A by constructing models that include a quark color superconducting (CSC) phase with a large energy gap; this phase appears at ultrahigh density region and reduces neutrino emissivity. In our model, a compact star has CSC quark core with a low neutrino emissivity surrounded by high emissivity region made by normal quarks. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency without considering nucleon superfluidity. Furthermore, we show that our scenario is consistent with the recent observations of the effective temperature of Cas A during the last 10 years, including nucleon superfluidity.

We study the hadron-quark mixed phase in protoneutron stars, where neutrinos are trapped and lepton number becomes a conserved quantity besides the baryon number and electric charge. Considering protoneutron-star matter as a ternary system, the Gibbs conditions are applied together with the Coulomb interaction. We find that there no crystalline ("pasta") structure appears in the regime of high lepton-number fraction; the size of pasta becomes very large and the geometrical structure becomes mechanically unstable due to the charge screening effect. Consequently the whole system is separated into two bulk regions like an amorphous state, where the surface effect is safely neglected. There, the local charge neutrality is approximately attained, so that the equation of state is effectively reduced to the one for a binary system. Hence, we conclude that there is no possibility for the density discontinuity to appear in protoneutron-star matter, which is a specific feature in a pure system. These features are important when considering astrophysical phenomena such as supernova explosions or radiation of the gravitational wave from protoneutron stars.

One of the most promising possibilities may be the appearance of quark matter in astrophysical phenomena in the light of recent progress in observations. The mechanism of deconfinement is not well understood, but the thermodynamical aspects of the hadron-quark (HQ) phase transition have been extensively studied in recent years. Then the mixed phase of hadron and quark matter becomes important; the proper treatment is needed to describe the HQ phase transition and derive the equation of state (EOS) for the HQ matter, based on the Gibbs conditions for phase equilibrium. We here adopt a EOS based on the baryon-baryon interactions including hyperons for the hadron phase, while we use rather simple EOS within the MIT bag model in the quark phase. For quark matter we further try to improve the previous EOS by considering other effective models of QCD. One of the interesting consequences may be the appearance of the inhomogeneous structures called "pasta", which are brought about by the surface and the Coulomb interaction effects. We present here a comprehensive review of our recent works about the HQ phase transition in various astrophysical situations: cold catalyzed matter, hot matter and neutrino-trapped matter. We show how the pasta structure becomes unstable by the charge screening of the Coulomb interaction, thermal effect or the neutrino trapping effect. Such inhomogeneous structure may affect astrophysical phenomena through its elasticity or thermal properties. Here we also discuss some implications on supernova explosion, gravitational wave and cooling of compact stars.

We make a perturbative calculation of neutrino scattering and absorption in hot and dense hyperonic neutron-star matter in the presence of a strong magnetic field. We find that the absorption cross-sections show a remarkable angular dependence in that the neutrino absorption strength is reduced in a direction parallel to the magnetic field and enhanced in the opposite direction. This asymmetry in the neutrino absorbtion can be as much as 2.2 % of the entire neutrino momentum for an interior magnetic field of \sim 2 x 10^{17} G. We estimate the pulsar kick velocities associated with this asymmetry in a fully relativistic mean-field theory formulation. We show that the kick velocities calculated here are comparable to observed pulsar velocities.

Quest for a new form of matter inside compact stars compels us to examine the thermodynamical properties of the phase transitions. We closely consider the first-order phase transitions and the phase equilibrium on the basis of the Gibbs conditions, taking the liquid-gas phase transition in asymmetric nuclear matter as an example. Characteristic features of the mixed phase are figured out by solving the coupled equations for mean-fields and densities of constituent particles self-consistently within the Thomas-Fermi approximation. The mixed phase is inhomogeneous matter composed of two phases in equilibrium; it takes a crystalline structure with a unit of various geometrical shapes, inside of which one phase with a characteristic shape, called "pasta", is embedded in another phase by some volume fraction. This framework enables us to properly take into account the Coulomb interaction and the interface energy, and thereby sometimes we see the mechanical instability of the geometric structures of the mixed phase. Thermal effect on the liquid-gas phase transition is also elucidated. Similarly hadron-quark deconfinement transition is studied in hyperonic matter, where the neutrino-trapping effect as well as the thermal effect is discussed in relation to the properties of the mixed phase. Specific features of the mixed phase are elucidated and the equation of state is presented.

Bearing in mind the application to the outer crust of the neutron stars (NSs), we investigate the magnetic field decay by means of the fully relativistic Particle-In-Cell simulations. Numerical computations are carried out in 2-dimensions, in which the initial magnetic fields are set to be composed both of the uniform magnetic fields that model the global fields penetrating the NS and of the turbulent magnetic fields that would be originated from the Hall cascade of the large-scale turbulence. Our results show that the whistler cascade of the turbulence transports the magnetic energy preferentially in the direction perpendicular to the uniform magnetic fields. It is also found that the distribution function of electrons becomes anisotropic because electrons with lower energies are predominantly heated in the direction parallel to the uniform magnetic fields due to the Landau resonance, while electrons with higher energies are heated mainly by the cyclotron resonance that makes the distribution function isotropic for the high energy tails. Furthermore we point out that the degree of anisotropy takes maximum as a function of the initial turbulent magnetic energy. As an alternative to the conventional ohmic dissipation, we propose that the magnetic fields in the outer crust of NSs, cascading down to the electron inertial scale via the whistler turbulence, would decay predominantly by the dissipation processes through the Landau damping and the cyclotron resonance.

We calculate stellar oscillations including the hadron-quark mixed phase considering the finite size effects. We find that it is possible to distinguish whether the density discontinuity exists or not in the stars, even if one will observe the gravitational waves of the fundamental mode. Additionally, the normalized eigenfrequencies of pressure modes depend strongly on the stellar mass and on the adopted equation of state. Especially, in spite of the fact that the radius of the neutron star with $1.4M_\odot$, which is standard mass, is almost independent from the equation of state with quark matter, the frequencies of pressure modes depend on the adopted equation of state. Thus, via observing the many kinds of gravitational waves, it will be possible to make a restriction on the equation of state.

We compare the Gibbs and Maxwell constructions for the hadron-quark phase transition in neutron and protoneutron stars, including interacting hyperons in the confined phase. We find that the hyperon populations are suppressed, and that neutrino trapping shifts the onset of the phase transition. The effects on the (proto)neutron star maximum mass are explored.

We study the quark-hadron mixed phase in proto-neutron stars with the finite-size effects. In the calculations of pasta structures appeared in the mixed phase, the Gibbs conditions require the pressure balance and chemical equilibrium between two phases besides the thermal equilibrium. We find that the region of the mixed phase is limited due to thermal instability. Moreover, we study the effects of neutrinos to the pasta structures. As a result, we find that the existence of neutrinos make the pasta structures unstable, too. These characteristic features of the hadron-quark mixed phase should be important for the middle stage of the evolutions of proto-neutron stars.

We investigate the properties of nuclear matter at the first-order phase transitions such as liquid-gas phase transition and hadron-quark phase transition. As a general feature of the first-order phase transitions of matter consisting of many species of charged particles, there appears a mixed phases with geometrical structures called ``pasta'' due to the balance of the Coulomb repulsion and the surface tension between two phases. The equation of state (EOS) of mixed phase is different from the one obtained by a bulk application of the Gibbs conditions or by the Maxwell construction due to the effects of the non-uniform structure. We show that the charge screening and strong surface tension make the EOS close to that of the Maxwell construction. The thermal effects are elucidated as well as the above finite-size effects.

We study the hadron-quark phase transition with the finite size effects at finite temperature. For the hadron phase, we adopt a realistic equation of state in the framework of the Brueckner-Hartree-Fock theory including hyperons. The properties of the mixed phase are clarified by considering the finite size effects under the Gibbs conditions. We find that the equation of state becomes softer than that at zero-temperature for some density region. We also find that the equation of state gets closer to that given by the Maxwell construction. Moreover, the number of hyperons is suppressed by the presence of quarks. These are characteristic features of the hadron-quark mixed phase, and should be important for many astrophysical phenomena such as mergers of binary neutron stars.

We report the effects of quark-hadron phase transition on the structures of general relativistic stars with purely toroidal magnetic field. For the mixed phase, we take into account of the finite-size effects, which lead to non-uniform "Pasta" structures. Our study is based on axisymetric and stationary formalism including purely toroidal magnetic field. For hybrid stars, we find the characteristic distribution of magnetic field, which has a discontinuity originated in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical phenomena, such as cooling processes.

We investigate structures of hybrid stars, which feature quark core surrounded by a hadronic matter mantle, with super-strong toroidal magnetic fields in full general relativity. Modeling the equation of state (EOS) with a first order transition by bridging the MIT bag model for the description of quark matter and the nuclear EOS by Shen et al., we numerically construct thousands of the equilibrium configurations to study the effects of the phase transition. It is found that the appearance of the quark phase can affect distributions of the magnetic fields inside the hybrid stars, making the maximum field strength about up to 30 % larger than for the normal neutron stars. Using the equilibrium configurations, we explore the possible evolutionary paths to the formation of hybrid stars due to the spin-down of magnetized rotating neutron stars. We find that the energy release by the phase transition to the hybrid stars is quite large ($\la 10^{52} \rm erg$) even for super strongly magnetized compact stars. Our results suggest that the strong gravitational-wave emission and the sudden spin-up signature could be observable signals of the QCD phase transition, possibly for a source out to Megaparsec distances.

We study structures of general relativistic compact stars with exotic matter. Our study is based on axisymetric and stationary formalism including purely toroidal magnetic field. We also study the finite size effects of quark-hadron mixed phase on structures of magnetars. For hybrid stars, we find a characteristic distribution of magnetic field, which has a discontinuity originated in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical phenomena, such as cooling processes.

We study the structures of hybrid stars with leptons at finite temperature under beta equilibrium. For the quark phase, we use the three flavor Nambu-Jona-Lasinio (NJL) model. For the hadron phase, we adopt nuclear equation of state (EOS) by Shen et al.. This EOS is in the framework of the relativistic mean field theory including the tree body effects. For the hadron-quark phase transition, we impose the bulk Gibbs construction or the Maxwell construction to take into account uncertainties by {\it finite size effects}. We find that the pure quark phase does not appear in stable star cores in all cases. With the phase transition, the maximum masses increase $\sim 10 %$ for high lepton fraction. On the contrary, without the transition, they decrease $\sim 10 %$. We also find that, in the NJL model, the lepton fraction is more important for structures of unstable stars than the temperature. This result is important for many astrophysical phenomena such as the core collapse of massive stars.

We perform two-dimensional, magnetohydrodynamical core-collapse simulations of massive stars accompanying the QCD phase transition. We study how the phase transition affects the gravitational waveforms near the epoch of core-bounce. As for initial models, we change the strength of rotation and magnetic fields. Particularly, the degree of differential rotation in the iron core (Fe-core) is changed parametrically. As for the microphysics, we adopt a phenomenological equation of state above the nuclear density, including two parameters to change the hardness before the transition. We assume the first order phase transition, where the conversion of bulk nuclear matter to a chirally symmetric quark-gluon phase is described by the MIT bag model. Based on these computations, we find that the phase transition can make the maximum amplitudes larger up to $\sim$ 10 percents than the ones without the phase transition. On the other hand, the maximum amplitudes become smaller up to $\sim$ 10 percents owing to the phase transition, when the degree of the differential rotation becomes larger. We find that even extremely strong magnetic fields $\sim 10^{17}$ G in the protoneutron star do not affect these results.

We present a scenario of formation of strange stars due to spin-down of {\it rapidly rotating} neutron stars left after supernova explosions . By assuming a process where the total baryon mass is conserved but the angular momentum is lost due to emission of gravitational waves and/or the magnetic braking, we find that the transition from rapidly rotating neutron stars to slowly rotating strange stars is possible; a large amount of energy $\sim 10^{53} ergs $ could be released. The liberated energy might become a new energy source for a delayed explosion of supernova. Furthermore, our scenario suggests that the supernova associated with gamma-ray bursts could become candidates for targets in the future observation of gravitational waves.