We study the post-merger mass ejection of low-mass binary neutron stars (NSs) with the system mass of $2.5\, M_\odot$, and subsequent nucleosynthesis by performing general-relativistic, neutrino-radiation viscous-hydrodynamics simulations in axial symmetry. We find that the merger remnants are long-lived massive NSs surviving more than several seconds, irrespective of the nuclear equations of state (EOSs) adopted. The ejecta masses of our fiducial models are $\sim 0.06$-$0.1\, M_\odot$ (depending on the EOS), being $\sim 30\%$ of the initial disk masses ($\sim 0.15$-$0.3\, M_\odot$). Post-processing nucleosynthesis calculations indicate that the ejecta is composed mainly of light $r$-process nuclei with small amounts of lanthanides (mass fraction $\sim 0.002$-$0.004$) and heavier species due to the modest average electron fraction ($\sim 0.32$-$0.34$) for a reasonable value of the viscous coefficient. Such abundance distributions are incompatible with the solar $r$-process-like abundance patterns found in all measured $r$-process-enhanced metal-poor stars. Therefore, low-mass binary NS mergers should be rare. If such low-mass NS mergers occur, their electromagnetic counterparts, kilonovae, will be characterized by an early bright blue emission because of the large ejecta mass as well as the small lanthanide fraction. We also show, however, that if the effective turbulent viscosity is very high, or there is an efficient mass ejection working in the early post-merger phase, the electron fraction of the ejecta could be low enough that the solar $r$-process-like abundance pattern is reproduced and the lanthanide fraction becomes so high that the kilonova would be characterized by early bright blue and late bright red emissions.

We reanalyze gravitational waves from binary-neutron-star mergers GW170817 and GW190425 using a numerical-relativity (NR) calibrated waveform model, the TF2+_Kyoto model, which includes nonlinear tidal terms. For GW170817, by imposing a uniform prior on the binary tidal deformability $\tilde{\Lambda}$, the symmetric $90\%$ credible interval of $\tilde{\Lambda}$ is estimated to be $481^{+436}_{-359}$ ($402^{+465}_{-279}$) for the case of $f_\mathrm{max}=1000~\mathrm{Hz}$ ($2048~\mathrm{Hz}$), where $f_\mathrm{max}$ is the maximum frequency in the analysis. We also reanalyze the event with other waveform models: two post-Newtonian waveform models (TF2_PNTidal and TF2+_PNTidal), the TF2+_NRTidal model that is another NR calibrated waveform model, and its upgrade, the TF2+_NRTidalv2 model. While estimates of parameters other than $\tilde{\Lambda}$ are broadly consistent among various waveform models, our results indicate that estimates of $\tilde{\Lambda}$ depend on waveform models. However, the difference is smaller than the statistical error. For GW190425, we can only obtain little information on the binary tidal deformability. The systematic error among the NR calibrated waveform models will become significant to measure $\tilde{\Lambda}$ as the number of detectors and events increase and sensitivities of detectors are improved.

We derive the upper limit to the ejecta mass of S190814bv, a black hole-neutron star merger candidate, through the radiative transfer simulations for kilonovae with the realistic ejecta density profile as well as the detailed opacity and heating rate models. The limits to the ejecta mass strongly depend on the viewing angle. For the face-on observations ($\le45^\circ$), the total ejecta mass should be smaller than $0.1\,M_\odot$ for the average distance of S190814bv ($D=267$ Mpc), while larger mass is allowed for the edge-on observations. We also derive the conservative upper limits of the dynamical ejecta mass to be $0.02\,M_\odot$, $0.03\,M_\odot$, and $0.05\,M_\odot$ for the viewing angle $\le 20^\circ$, $\le 50^\circ$, and for $\le 90^\circ$, respectively. We show that the {\it iz}-band observation deeper than $22$ mag within 2 d after the GW trigger is crucial to detect the kilonova with the total ejecta mass of $0.06\,M_\odot$ at the distance of $D=300$ Mpc. We also show that a strong constraint on the NS mass-radius relation can be obtained if the future observations put the upper limit of $0.03\,M_\odot$ to the dynamical ejecta mass for a BH-NS event with the chirp mass smaller than $\lesssim 3\,M_\odot$ and effective spin larger than $\gtrsim 0.5$.

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

We report results of numerical relativity simulations for {\it new} 26 non-spinning binary neutron star systems with 6 grid resolutions using an adaptive mesh refinement numerical re\ lativity code {\tt SACRA-MPI}. The finest grid spacing is $\approx 64$--$85$ m, depending on the systems. First, we derive long-term high-precision inspiral gravitational waveforms and show that the accumulated gravitational-wave phase error due to the finite grid resolution is less than $0.5$ rad during more than $200$ rad phase evolution irrespective of the systems. We also find that the gravitational-wave phase error for a binary system with a tabulated equation of state (EOS) is comparable to that for a piecewise polytropic EOS. Then we validate the SACRA inspiral gravitational waveform template, which will be used to extract tidal deformability from gravitational wave observation, and find that accuracy of \ our waveform modeling is $\lesssim 0.1$ rad in the gravitational-wave phase and $\lesssim 20 \%$ in the gravitational-wave amplitude up to the gravitational-wave frequency $1000$ Hz.\ Finally, we calibrate the proposed universal relations between a post-merger gravitational wave signal and tidal deformability/neutron star radius in the literature and show that th\ ey suffer from systematics and many relations proposed as universal are not very universal. Improved fitting formulae are also proposed.

We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed within its surrounding magnetosphere. A detailed study of the resulting magnetospheric configurations in our setting is presented, incorporating as well the effects due to neutron star spin and the misalignment of the magnetic and orbital axes. We find many features resembling those of pulsar magnetospheres for the orbiting neutron star --even in the absence of spin--, being of particular interest the development of a spiral current sheet that extends beyond the light cylinder. Then, we use recent advances in pulsar theory to estimate electromagnetic emissions produced at the reconnection regions of such current sheets.

We argue that the kilonova/macronova associated with the gravitational-wave event GW190425 could have been bright enough to be detected if it was caused by the merger of a low-mass black hole and a neutron star. Although tidal disruption occurs for such a low-mass black hole as is generally expected, the masses of the dynamical ejecta are limited to <~10^{-3}M_sun, which is consistent with previous work in the literature. The remnant disk could be as massive as 0.05--0.1M_sun, and the disk outflow of ~0.01--0.03M_sun is likely to be driven by viscous or magnetohydrodynamic effects. The disk outflow may not be neutron-rich enough to synthesize an abundance of lanthanide elements, even in the absence of strong neutrino emitter, if the ejection is driven on the viscous time scale of >~0.3s. If this is the case, the opacity of the disk outflow is kept moderate, and a kilonova/macronova at the distance of GW190425 reaches a detectable brightness of 20--21mag at 1day after merger for most viewing angles. If some disk activity ejects the mass within ~0.1s, instead, lanthanide-rich outflows would be launched and the detection of emission becomes challenging. Future possible detections of kilonovae/macronovae from GW190425-like systems will disfavor the prompt collapse of binary neutron stars and a non-disruptive low-mass black hole--neutron star binary associated with a small neutron-star radius, whose mass ejection is negligible. The host-galaxy distance will constrain the viewing angle and deliver further information about the mass ejection.

New viscous neutrino-radiation hydrodynamics simulations are performed for accretion disks surrounding a spinning black hole with low mass $3M_\odot$ and dimensionless spin 0.8 or 0.6 in full general relativity, aiming at modeling the evolution of a merger remnant of massive binary neutron stars or low-mass black hole-neutron star binaries. We reconfirm the following results found by previous studies of other groups: 15-30% of the disk mass is ejected from the system with the average velocity of $\sim $5-10% of the speed of light for the plausible profile of the disk as merger remnants. In addition, we find that for the not extremely high viscous coefficient case, the neutron richness of the ejecta does not become very high, because weak interaction processes enhance the electron fraction during the viscous expansion of the disk before the onset of the mass ejection, resulting in the suppression of the lanthanide synthesis. For high-mass disks, the viscous expansion timescale is increased by a longer-term neutrino emission, and hence, the electron fraction of the ejecta becomes even higher. We also confirm that the mass distribution of the electron fraction depends strongly on the magnitude of the given viscous coefficient. This demonstrates that a first-principle magnetohydrodynamics simulation is necessary for black hole-disk systems with sufficient grid resolution and with sufficiently long timescale (longer than seconds) to clarify the nucleosynthesis and electromagnetic signals from them.

We perform radiative transfer simulations for kilonova in various situations, including the cases of prompt collapse to a black hole from neutron-star mergers, high-velocity ejecta possibly accelerated by magnetars, and a black hole-neutron star merger. Our calculations are done employing ejecta profiles predicted by numerical-relativity simulations and a new line list for all the r-process elements. We found that (i) the optical emission for binary neutron stars promptly collapsing to a black hole would be fainter by $\gtrsim1$-$2\,{\rm mag}$ than that found in GW170817, while the infrared emission could be as bright as that in GW170817 if the post-merger ejecta is as massive as $\approx0.01\,M_\odot$; (ii) the kilonova would be brighter than that observed in GW170817 for the case that the ejecta is highly accelerated by the electromagnetic energy injection from the remnant, but it would decline rapidly and the magnitude would become fainter than in GW170817 within a few days; (iii) the optical emission from a black hole-neutron star merger ejecta could be as bright as that observed in GW170817 for the case that sufficiently large amount of matter is ejected ($\gtrsim0.02\,M_\odot$), while the infrared brightness would be brighter by $1$-$2\,{\rm mag}$ at the same time. We show that the difference in the ejecta properties would be imprinted in the differences in the peak brightness and time of peak. This indicates that we may be able to infer the type of the central engine for kilonovae by observation of the peak in the multiple band.

We find that the Hanford and Livingston detectors of Advanced LIGO derive a distinct posterior probability distribution of binary tidal deformability tilde{Lambda} of the first binary-neutron-star merger GW170817. By analyzing public data of GW170817 with a nested-sampling engine and the default TaylorF2 waveform provided by the LALInference package, the probability distribution of the binary tidal deformability derived by the LIGO-Virgo detector network turns out to be determined dominantly by the Hanford detector. Specifically, by imposing the flat prior on tidal deformability of individual stars, symmetric 90% credible intervals of tilde{Lambda} are estimated to be 527^{+619}_{-345} with the Hanford detector, 927^{+522}_{-619} with the Livingston detector, and 455^{+668}_{-281} with the LIGO-Virgo detector network. Furthermore, the distribution derived by the Livingston detector changes irregularly when we vary the maximum frequency of the data used in the analysis. This feature is not observed for the Hanford detector. While they are all consistent, the discrepancy and irregular behavior suggest that an in-depth study of noise properties might improve our understanding of GW170817 and future events.

Mergers of binary neutron stars and black hole-neutron star binaries are one of the most promising sources for the ground-based gravitational-wave (GW) detectors and also a high-energy astrophysical phenomenon as illustrated by the observations of gravitational waves and electromagnetic (EM) waves in the event of GW170817. Mergers of these neutron-star binaries are also the most promising site for r-process nucleosynthesis. Numerical simulation in full general relativity (numerical relativity) is a unique approach to the theoretical prediction of the merger process, GWs emitted, mass ejection process, and resulting EM emission. We summarize our current understanding for the processes of neutron star mergers and subsequent mass ejection based on the results of the latest numerical-relativity simulations. We emphasize that the predictions of the numerical-relativity simulations agrees broadly with the optical and infrared observations of GW170817.

We report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date.

We revisit the constraint on the maximum mass of cold spherical neutron stars coming from the observational results of GW170817. We develop a new framework for the analysis by employing both energy and angular momentum conservation laws as well as solid results of latest numerical-relativity simulations and of neutron stars in equilibrium. The new analysis shows that the maximum mass of cold spherical neutron stars can be only weakly constrained as $M_{\rm max} \alt 2.3M_\odot$. Our present result illustrates that the merger remnant neutron star at the onset of collapse to a black hole is not necessarily rapidly rotating and shows that we have to take into account the angular momentum conservation law to impose the constraint on the maximum mass of neutron stars.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

KAGRA is a new gravitational wave detector which aims to begin joint observation with Advanced LIGO and Advanced Virgo from late 2019. Here, we present KAGRA's possible upgrade plans to improve the sensitivity in the decade ahead. Unlike other state-of-the-art detectors, KAGRA requires different investigations for the upgrade since it is the only detector which employs cryogenic cooling of the test mass mirrors. In this paper, investigations on the upgrade plans which can be realized by changing the input laser power, increasing the mirror mass, and injecting frequency dependent squeezed vacuum are presented. We show how each upgrade affects to the detector frequency bands and also discuss impacts on gravitational-wave science. We then propose an effective progression of upgrades based on technical feasibility and scientific scenarios.

We revisit the lower bound on binary tidal deformability tilde{Lambda} imposed by a luminous kilonova/macronova, AT 2017gfo, by numerical-relativity simulations of models that are consistent with gravitational waves from the binary neutron star merger GW170817. Contrary to the claim made in the literature, we find that binaries with tilde{Lambda}<~400 can explain the luminosity of AT 2017gfo, as long as moderate mass ejection from the remnant is assumed as had been done in previous work. The reason is that the maximum mass of a neutron star is not strongly correlated with the tidal deformability of neutron stars with a typical mass of ~1.4M_sun. If the maximum mass is so large that the binary does not collapse into a black hole immediately after merger, the mass of the ejecta can be sufficiently large irrespective of the binary tidal deformability. We present models of binary mergers with tilde{Lambda} down to 242 that satisfy the requirement on the mass of the ejecta from the luminosity of AT 2017gfo. We further find that the luminosity of AT 2017gfo could be explained by models that do not experience bounce after merger. We conclude that the luminosity of AT 2017gfo is not very useful for constraining the binary tidal deformability. Accurate estimation of the mass ratio will be necessary to establish a lower bound using electromagnetic counterparts in the future. We also caution that merger simulations that employ a limited class of tabulated equations of state could be severely biased due to the lack of generality.

We investigate remnant neutron star masses (in particular, the minimum allowed mass) by performing advanced stellar evolution calculations and neutrino-radiation hydrodynamics simulations for core-collapse supernova explosions. We find that, based on standard astrophysical scenarios, low-mass carbon-oxygen cores can have sufficiently massive iron cores that eventually collapse, explode as supernovae, and give rise to remnant neutron stars that have a minimum mass of 1.17 M$_\odot$ --- compatible with the lowest mass of the neutron star precisely measured in a binary system of PSR J0453+1559.

Rapidly and differentially rotating compact stars are believed to be formed in binary neutron star merger events, according to both numerical simulations and the multi-messenger observation of GW170817. The lifetime and evolution of such a differentially rotating star, is tightly related to the observations in the post-merger phase. Various studies on the maximum mass of differentially rotating neutron stars have been done in the past, most of which assume the so-called $j$-const law as the rotation profile inside the star and consider only neutron star equations of state. In this paper, we extend the studies to strange star models, as well as to a new rotation profile model. Significant differences are found between differentially rotating strange stars and neutron stars, with both differential rotation laws. A moderate differential rotation rate for neutron stars is found to be too large for strange stars, resulting in a rapid drop in the maximum mass as the differential rotation degree is increased further from $\hat{A}\sim2.0$, where $\hat{A}$ is a parameter characterizing the differential rotation rate for $j$-const law. As a result the maximum mass of a differentially rotating self-bound star drops below the uniformly rotating mass shedding limit for a reasonable degree of differential rotation. The continuous transition to the toroidal sequence is also found to happen at a much smaller differential rotation rate and angular momentum than for neutron stars. In spite of those differences, $\hat{A}$-insensitive relation between the maximum mass for a given angular momentum is still found to hold, even for the new differential rotation law. Astrophysical consequences of these differences and how to distinguish between strange star and neutron star models with future observations are also discussed.

We compute gravitational waves emitted by the collapse of a rotating very massive star (VMS) core leading directly to a black hole in axisymmetric numerical-relativity simulations. The evolved rotating VMS is derived by a stellar evolution calculation and its initial mass and the final carbon-oxygen core mass are $320M_\odot$ and $\approx 150M_\odot$, respectively. We find that for the moderately rapidly rotating cases, the peak strain amplitude and the corresponding frequency of gravitational waves are $\sim 10^{-22}$ and $f \approx 300$--600\,Hz for an event at the distance of $D=50$~Mpc. Such gravitational waves will be detectable only for $D \lesssim 10$~Mpc by second generation detectors, advanced LIGO, advanced VIRGO, and KAGRA, even if the designed sensitivity for these detectors is achieved. However, third-generation detectors will be able to detect such gravitational waves for an event up to $D \sim 100$~Mpc. The detection of the gravitational-wave signal will provide a potential opportunity for verifying the presence of VMSs with mass $\gtrsim 300M_\odot$ and their pair-unstable collapse in the universe.

KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.

The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.

We explore the formation process of a black hole (BH) through the pair-instability collapse of a rotating Population III very massive star in axisymmetric numerical relativity. As the initial condition, we employ a progenitor star which is obtained by evolving a rapidly rotating zero-age main sequence (ZAMS) star with mass $320M_\odot$ until it reaches a pair instability region. We find that for such rapidly rotating model, a fraction of the mass, $\sim 10M_\odot$, forms a torus surrounding the remnant BH of mass $\sim 130M_\odot$ and an outflow is driven by a hydrodynamical effect. We also perform simulations, artificially reducing the initial angular velocity of the progenitor star, and find that only a small or no torus is formed and no outflow is driven. We discuss the possible evolution scenario of the remnant torus for the rapidly rotating model by considering the viscous and recombination effects and show that if the energy of $\sim 10^{52}$ erg is injected from the torus to the envelope, the luminosity and timescale of the explosion could be of the orders of $10^{43}$ erg/s and yrs, respectively. We also point out the possibility for observing gravitational waves associated with the BH formation for the rapidly rotating model by ground-based gravitational-wave detectors.

We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

Recent detection of gravitational waves from a binary-neutron star merger (GW170817) and the subsequent observations of electromagnetic counterparts provide a great opportunity to study the physics of compact binary mergers. The optical and near-infrared counterparts to GW170817 (SSS17a, also known as AT 2017gfo or DLT17ck) are found to be consistent with a kilonova/macronova scenario with red and blue components. However, in most of previous studies in which contribution from each ejecta component to the lightcurves is separately calculated and composited, the red component is too massive as dynamical ejecta and the blue component is too fast as post-merger ejecta. In this letter, we perform a 2-dimensional radiative transfer simulation for a kilonova/macronova consistently taking the interplay of multiple ejecta components into account. We show that the lightcurves and photospheric velocity of SSS17a can be reproduced naturally by a setup consistent with the prediction of the numerical-relativity simulations.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we employ the azimuthally averaged data of the MNS-torus system derived in a three-dimensional, numerical-relativity simulation for the binary neutron star merger. The viscous effect plays key roles for the remnant evolution and mass ejection from it in two phases of the evolution. In the first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to a rigidly rotating state, and as a result, a sound wave, which subsequently becomes a shock wave, is formed in the vicinity of the MNS due to the variation of the quasi-equilibrium state of the MNS. The shock wave induces significant mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim 0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the viscosity-driven ejecta could dominate that of the dynamical ejecta of mass $\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$ times lower than that of the dynamical ejecta. This indicates that the electromagnetic signal from the ejecta would be rapidly evolving, bright, and blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for which the effect of the dynamical ejecta is minor.

We perform a general-relativistic magnetohydrodynamics simulation for $\approx 30$ ms after merger of a binary neutron star to a remnant massive neutron star (RMNS) with a high spatial resolution of the finest grid resolution $12.5$ m. First, we estimate that the Kelvin-Helmholtz instability at merger could amplify the magnetic-field energy up to $\sim 1\%$ of the thermal energy. Second, we find that the magnetorotational instability in the RMNS envelope and torus with $\rho < 10^{13}~{\rm g~cm^{-3}}$ sustains magneto-turbulent state and the effective viscous parameter in these regions is likely to converge to $\approx 0.01$--$0.02$ with respect to the grid resolution. Third, the current grid resolution is not still fine enough to sustain magneto-turbulent state in the RMNS with $\rho \ge 10^{13}~{\rm g~cm^{-3}}$.

We find, using high resolution numerical relativistic simulations, that the tail of the dynamical ejecta of neutron star mergers extends to mildly relativistic velocities faster than $0.7c$. The kinetic energy of this fast tail is $\sim 10^{47}$--$10^{49}$ erg, depending on the neutron star equation of state and on the binary masses. The synchrotron flare arising from the interaction of this fast tail with the surrounding ISM can power the observed non-thermal emission that followed GW170817, provided that the ISM density is $\sim 10^{-2}\,{\rm cm^{-3}}$, the two neutron stars had roughly equal masses and the neutron star equation of state is soft (small neutron star radii). One of the generic predictions of this scenario is that the cooling frequency crosses the X-ray band on a time scale of a few months to a year, leading to a cooling break in the X-ray light curve. If this dynamical ejecta scenario is correct, we expect that the synchrotron radio flare from the ejecta that have produced the macronova/kilonova emission will be observable on time scales of $10^3$ to $10^5$ days. Further multi-frequency observations will confirm or rule out this dynamical ejecta scenario.

We develop a model for frequency-domain gravitational waveforms from inspiraling binary neutron stars. Our waveform model is calibrated by comparison with hybrid waveforms constructed from our latest high-precision numerical-relativity waveforms and the SEOBNRv2T waveforms in the frequency range of $10$--$1000\,{\rm Hz}$. We show that the phase difference between our waveform model and the hybrid waveforms is always smaller than $0.1\, {\rm rad}$ for the binary tidal deformability, ${\tilde \Lambda}$, in the range $300\lesssim{\tilde \Lambda}\lesssim1900$ and for the mass ratio between 0.73 and 1. We show that, for $10$--$1000\,{\rm Hz}$, the distinguishability for the signal-to-noise ratio $\lesssim50$ and the mismatch between our waveform model and the hybrid waveforms are always smaller than 0.25 and $1.1\times10^{-5}$, respectively. The systematic error of our waveform model in the measurement of ${\tilde \Lambda}$ is always smaller than $20$ with respect to the hybrid waveforms for $300\lesssim{\tilde \Lambda}\lesssim1900$. The statistical error in the measurement of binary parameters is computed employing our waveform model, and we obtain results consistent with the previous studies. We show that the systematic error of our waveform model is always smaller than $20\%$ (typically smaller than $10\%$) of the statistical error for events with the signal-to-noise ratio of $50$.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We study the merger of black hole-neutron star binaries by fully general-relativistic neutrino-radiation-hydrodynamics simulations throughout the coalescence, particularly focusing on the role of neutrino irradiation in dynamical mass ejection. Neutrino transport is incorporated by an approximate transfer scheme based on the truncated moment formalism. While we fix the mass ratio of the black hole to the neutron star to be 4 and the dimensionless spin parameter of the black hole to be 0.75, the equations of state for finite-temperature neutron-star matter are varied. The hot accretion disk formed after tidal disruption of the neutron star emits a copious amount of neutrinos with the peak total luminosity ~1--3x10^53 erg s^(-1) via thermal pair production and subsequent electron/positron captures on free nucleons. Nevertheless, the neutrino irradiation does not modify significantly the electron fraction of the dynamical ejecta from the neutrinoless beta-equilibrium value at zero temperature of initial neutron stars. The mass of the wind component driven by neutrinos from the remnant disk is negligible compared to the very neutron-rich dynamical component, throughout our simulations performed until a few tens milliseconds after the onset of merger, for the models considered in this study. These facts suggest that the ejecta from black hole-neutron star binaries are very neutron rich and are expected to accommodate strong r-process nucleosynthesis, unless magnetic or viscous processes contribute substantially to the mass ejection from the disk. We also find that the peak neutrino luminosity does not necessarily increase as the disk mass increases, because tidal disruption of a compact neutron star can result in a remnant disk with a small mass but high temperature.

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

Gravitational-wave observation together with a large number of electromagnetic observations shows that the source of the latest gravitational-wave event, GW170817, detected primarily by advanced LIGO, is the merger of a binary neutron star. We attempt to interpret this observational event based on our results of numerical-relativity simulations performed so far paying particular attention to the optical and infra-red observations. We finally reach a conclusion that this event is described consistently by the presence of a long-lived hypermassive or supramassive neutron star as the merger remnant, because (i) significant contamination by lanthanide elements along our line of sight to this source can be avoided by the strong neutrino irradiation from it and (ii) it could play a crucial role to produce an ejecta component of appreciable mass with fast motion in the post-merger phase. We also point out that (I) the neutron-star equation of state has to be sufficiently stiff (i.e., the maximum mass of cold spherical neutron stars, M_max, has to be appreciably higher than 2M_sun in order that a long-lived massive neutron star can be formed as the merger remnant for the binary systems of GW170817, for which the initial total mass is >~ 2.73M_sun and (II) no detection of relativistic optical counterpart suggests a not-extremely high value of M_max approximately as 2.15-2.25M_sun.

KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.

The location of the cutoff in the gravitational-wave spectrum of black hole-neutron star mergers is correlated strongly with the neutron-star radius for the case that the neutron star is disrupted by the black hole during the merger. However, the modulation which appears in the spectrum due to the mode mixing makes it difficult to measure the cutoff frequency if gravitational waves are observed from inclined direction or the binary is precessing. In this letter, we show that the cutoff frequency can be measured even in such situations with a method we have recently proposed to reconstruct the face-on waveforms only from the strain observed from a particular direction. We show that the systematic error in the measurement of the neutron-star radius can be reduced to $\lesssim 5\%$ for the case that tidal disruption of the neutron star occurs significantly.

Extending our previous studies, we perform high-resolution simulations of inspiraling binary neutron stars in numerical relativity. We thoroughly carry through a convergence study in our currently available computational resources with the smallest grid spacing of $\approx 63$--86~meter for the neutron-star radius 10.9--13.7\,km. The estimated total error in the gravitational-wave phase is of order 0.1~rad for the total phase of $\gtrsim 210$\,rad in the last $\sim 15$--16 inspiral orbits. We then compare the waveforms (without resolution extrapolation) with those calculated by the latest effective-one-body formalism (tidal SEOBv2 model referred to as TEOB model). We find that for any of our models of binary neutron stars, the waveforms calculated by the TEOB formalism agree with the numerical-relativity waveforms up to $\approx 3$\,ms before the peak of the gravitational-wave amplitude is reached: For this late inspiral stage, the total phase error is $\lesssim 0.1$\,rad. Although the gravitational waveforms have an inspiral-type feature for the last $\sim 3$\,ms, this stage cannot be well reproduced by the current TEOB formalism, in particular, for neutron stars with large tidal deformability (i.e., lager radius). The reason for this is described.

We performed general relativistic, long-term, axisymmetric neutrino radiation hydrodynamics simulations for the remnant formed after the binary neutron star merger, which consist of a massive neutron star and a torus surrounding it. As an initial condition, we employ the result derived in a three-dimensional, numerical relativity simulation for the binary neutron star merger. We investigate the properties of neutrino-driven ejecta. Due to the pair-annihilation heating, the dynamics of the neutrino-driven ejecta is significantly modified. The kinetic energy of the ejecta is about two times larger than that in the absence of the pair-annihilation heating. This suggests that the pair-annihilation heating plays an important role in the evolution of the merger remnants. The relativistic outflow, which is required for driving gamma-ray bursts, is not observed because the specific heating rate around the rotational axis is not sufficiently high due to the baryon loading caused by the neutrino-driven ejecta from the massive neutron star. We discuss the condition for launching the relativistic outflow and the nucleosynthesis in the ejecta.

We investigate the explosive nucleosynthesis during two dimensional neutrino-driven explosion of ultra-stripped Type Ic supernovae evolved from 1.45 and 1.5 M$_\odot$ CO stars. These supernovae explode with the explosion energy of $\sim 10^{50}$ erg and release $\sim 0.1$ M$_\odot$ ejecta. The light trans-iron elements Ga-Zr are produced in the neutrino-irradiated ejecta. The abundance distribution of these elements has a large uncertainty because of the uncertainty of the electron fraction of the neutrino-irradiated ejecta. The yield of the elements will be less than 0.01 M$_\odot$. Ultra-stripped supernova and core-collapse supernova evolved from a light CO core can be main sources of the light trans-iron elements. They could also produce neutron-rich nuclei $^{48}$Ca. The ultra-stripped supernovae eject $^{56}$Ni of $\sim$ 0.006 - 0.01 M$_\odot$. If most of neutrino-irradiated ejecta is proton-rich, $^{56}$Ni will be produced more abundantly. The light curves of these supernovae indicate sub-luminous fast decaying explosion with the peak magnitude of about $-15$ - $-16$. Future observations of ultra-stripped supernovae could give a constraint to the event rate of a class of neutron star mergers.

Employing a simplified version of the Israel-Stewart formalism of general-relativistic shear-viscous hydrodynamics, we explore the evolution of a remnant massive neutron star of binary neutron star merger and pay special attention to the resulting gravitational waveforms. We find that for the plausible values of the so-called viscous alpha parameter of the order $10^{-2}$, the degree of the differential rotation in the remnant massive neutron star is significantly reduced in the viscous timescale, $\alt 5$\,ms. Associated with this, the degree of non-axisymmetric deformation is also reduced quickly, and as a consequence, the amplitude of quasi-periodic gravitational waves emitted also decays in the viscous timescale. Our results indicate that for modeling the evolution of the merger remnants of binary neutron stars, we would have to take into account magnetohydrodynamics effects, which in nature could provide the viscous effects.

Supermassive stars (SMSs) of mass $\gtrsim 10^5 M_\odot$ are candidates for seeds of supermassive black holes found in the center of many massive galaxies. We simulate the gravitational collapse of a rigidly rotating SMS core including nuclear burning effects in axisymmetric numerical-relativity simulation. We find that for realistic initial conditions, the nuclear burning does not play an important role. After the collapse, a torus surrounding a rotating black hole is formed and a fraction of the torus material is ejected. We quantitatively study the relation between the properties of these objects and rotation. We find that if a SMS core is sufficiently rapidly rotating, the torus and outflow mass have approximately $6\%$ and $1\%$ of the initial mass, respectively. The typical average velocity and the total kinetic energy of the outflow are $0.2~c$ and $10^{54-56}$ erg where $c$ is the speed of light. Finally, we briefly discuss the possibility for observing the torus and outflow.

Employing a simplified version of the Israel-Stewart formalism for general-relativistic shear-viscous hydrodynamics, we perform axisymmetric general-relativistic simulations for a rotating neutron star surrounded by a massive torus, which can be formed from differentially rotating stars. We show that with our choice of a shear-viscous hydrodynamics formalism, the simulations can be stably performed for a long time scale. We also demonstrate that with a possibly high shear-viscous coefficient, not only viscous angular momentum transport works but also an outflow could be driven from a hot envelope around the neutron star for a time scale $\gtrsim 100\,$ms with the ejecta mass $\gtrsim 10^{-2}M_\odot$ which is comparable to the typical mass for dynamical ejecta of binary neutron star mergers. This suggests that massive neutron stars surrounded by a massive torus, which are typical outcomes formed after the merger of binary neutron stars, could be the dominant source for providing neutron-rich ejecta, if the effective shear viscosity is sufficiently high, i.e., if the viscous $\alpha$ parameter is $\gtrsim 10^{-2}$. The present numerical result indicates the importance of a future high-resolution magnetohydrodynamics simulation that is the unique approach to clarify the viscous effect in the merger remnants of binary neutron stars by the first-principle manner.

Black hole-neutron star mergers are among the promising gravitational-wave sources for ground-based detectors, and gravitational waves from black hole-neutron star mergers are expected to be detected in the next few years. Simultaneous detection of electromagnetic counterparts with gravitational-wave detection provides rich information about the merger events. Among the possible electromagnetic counterparts from the black hole-neutron star merger, the emission powered by the decay of radioactive r-process nuclei, so called kilonova/macronova, is one of the best targets for follow-up observation. We derive fitting formulas for the mass and the velocity of ejecta from a generic black hole-neutron star merger based on recently performed numerical relativity simulations. We combined these fitting formulas with a new semi-analytic model for a black hole-neutron star kilonova/macronova lightcurve which reproduces the results of radiation-transfer simulations. Specifically, the semi-analytic model reproduces the result of each band magnitude obtained by the radiation transfer simulations within ~1 mag. By using this semi-analytic model, we found that, at 400 Mpc, the kilonova/macronova is as bright as 22-24 mag for the cases with a small chirp mass and a high black hole spin, and >28 mag for a large chirp mass and a low black hole spin. We also apply our model to GRB130603B as an illustration, and show that a black hole-neutron star merger with a rapidly spinning black hole and a large neutron star radius is favored.

We derive the gravitational waveform from the collapse of a rapidly rotating supermassive star (SMS) core leading directly to a seed of a supermassive black hole (SMBH) in axisymmetric numerical-relativity simulations. We find that the peak strain amplitude of gravitational waves emitted during the black-hole formation is $\approx 5 \times 10^{-21}$ at the frequency $f \approx 5$\,mHz for an event at the cosmological redshift $z=3$, if the collapsing SMS core is in the hydrogen-burning phase. Such gravitational waves will be detectable by space laser interferometric detectors like eLISA with signal-to-noise ratio $\approx 10$, if the sensitivity is as high as LISA for $f=1$--10\,mHz. The detection of the gravitational-wave signal will provide a potential opportunity for testing the direct-collapse scenario for the formation of a seed of SMBHs.

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

We revisit secular stability against quasi-radial collapse for rigidly rotating supermassive stars (SMSs) in general relativity. We suppose that the SMSs are in a nuclear-burning phase and can be modeled by polytropic equations of state with the polytropic index $n_p$ slightly smaller than $3$. The stability is determined in terms of the turning-point method. We find a fitting formula of the stability condition for the plausible range of $n_p$ ($2.95 \alt n_p \alt 3$) for SMSs. This condition reconfirms that, while non-rotating SMSs with mass $\sim 10^5M_\odot$--$10^6M_\odot$ may undergo a general-relativistically induced quasi-radial collapse, rigidly rotating SMSs with a ratio of rotational to gravitational potential energy ($\beta$) of $\sim 10^{-2}$ are likely to be stable against collapse unless they are able to accrete $\sim 5$ times more mass during the (relatively brief) hydrogen-burning phase of their evolution. We discuss implications of our results.

We perform neutrino radiation-hydrodynamics simulations for the merger of asymmetric binary neutron stars in numerical relativity. Neutron stars are modeled by soft and moderately stiff finite-temperature equations of state (EOS). We find that the properties of the dynamical ejecta such as the total mass, neutron richness profile, and specific entropy profile depend on the mass ratio of the binary systems for a given EOS in a unique manner. For the soft EOS (SFHo), the total ejecta mass depends weakly on the mass ratio, but the average of electron number per baryon ($Y_e$) and specific entropy ($s$) of the ejecta decreases significantly with the increase of the degree of mass asymmetry. For the stiff EOS (DD2), with the increase of the mass asymmetry degree, the total ejecta mass significantly increases while the average of $Y_e$ and $s$ moderately decreases. We find again that only for the soft EOS (SFHo), the total ejecta mass exceeds $0.01M_\odot$ irrespective of the mass ratio chosen in this paper. The ejecta have a variety of electron number per baryon with its average approximately between $Y_e \sim 0.2$ and $\sim 0.3$ irrespective of the EOS employed, which is well-suited for the production of the r-process heavy elements (second and third peaks), although its averaged value decreases with the increase of the degree of mass asymmetry.

A new air-shower core-detector array (YAC : Yangbajing Air-shower Core-detector array) has been developed to measure the primary cosmic-ray composition at the "knee" energies in Tibet, China, focusing mainly on the light components. The prototype experiment (YAC-I) consisting of 16 detectors has been constructed and operated at Yangbajing (4300 m a.s.l.) in Tibet since May 2009. YAC-I is installed in the Tibet-III AS array and operates together. In this paper, we performed a Monte Carlo simulation to check the sensitivity of YAC-I+Tibet-III array to the cosmic-ray light component of cosmic rays around the knee energies, taking account of the observation conditions of actual YAC-I+Tibet-III array. The selection of light component from others was made by use of an artificial neural network (ANN). The simulation shows that the light-component spectrum estimated by our methods can well reproduce the input ones within 10\% error, and there will be about 30\% systematic errors mostly induced by the primary and interaction models used. It is found that the full-scale YAC and the Tibet-III array is powerful to study the cosmic-ray composition, in particular, to obtain the energy spectra of protons and helium nuclei around the knee energies.

We explore magnetic-field amplification due to the Kelvin-Helmholtz instability during binary neutron star mergers. By performing high-resolution general relativistic magnetohydrodynamics simulations with a resolution of $17.5$ m for $4$--$5$ ms after the onset of the merger on the Japanese supercomputer "K", we find that an initial magnetic field of moderate maximum strength $10^{13}$ G is amplified at least by a factor of $\approx 10^3$. We also explore the saturation of the magnetic-field energy and our result shows that it is likely to be $\gtrsim 4 \times 10^{50}$ erg, which is $\gtrsim 0.1\%$ of the bulk kinetic energy of the merging binary neutron stars.

Tidal disruption has a dramatic impact on the outcome of neutron star-black hole mergers. The phenomenology of these systems can be divided in three classes: nondisruptive, mildly disruptive or disruptive. The cutoff frequency of the gravitational radiation produced during the merger (which is potentially measurable by interferometric detectors) is very different in each regime, and when the merger is disuptive it carries information on the neutron star equation of state. Here we use semianalytical tools to derive a formula for the critical binary mass ratio $Q=M_{\rm BH}/M_{\rm NS}$ below which mergers are disruptive as a function of the stellar compactness $\mathcal{C}=M_{\rm NS}/R_{\rm NS}$ and the dimensionless black hole spin $\chi$. We then employ a new gravitational waveform amplitude model, calibrated to $134$ general relativistic numerical simulations of binaries with black hole spin (anti-)aligned with the orbital angular momentum, to obtain a fit to the gravitational-wave cutoff frequency in the disruptive regime as a function of $\mathcal{C}$, $Q$ and $\chi$. Our findings are important to build gravitational wave template banks, to determine whether neutron star-black hole mergers can emit electromagnetic radiation (thus helping multimessenger searches), and to improve event rate calculations for these systems.

We study explosion characteristics of ultra-stripped supernovae (SNe), which are candidates of SNe generating binary neutron stars (NSs). As a first step, we perform stellar evolutionary simulations of bare carbon-oxygen cores of mass from 1.45 to 2.0 $M_\odot$ until the iron cores become unstable and start collapsing. We then perform axisymmetric hydrodynamics simulations with spectral neutrino transport using these stellar evolution outcomes as initial conditions. All models exhibit successful explosions driven by neutrino heating. The diagnostic explosion energy, ejecta mass, Ni mass, and NS mass are typically $\sim 10^{50}$ erg, $\sim 0.1 M_\odot$, $\sim 0.01M_\odot$, and $\approx 1.3 M_\odot$, which are compatible with observations of rapidly-evolving and luminous transient such as SN 2005ek. We also find that the ultra-stripped SN is a candidate for producing the secondary low-mass NS in the observed compact binary NSs like PSR J0737-3039.

We report results of a high-resolution numerical-relativity simulation for the merger of black hole-magnetized neutron star binaries on Japanese supercomputer "K". We focus on a binary that is subject to tidal disruption and subsequent formation of a massive accretion torus. We find the launch of thermally driven torus wind, subsequent formation of a funnel wall above the torus and a magnetosphere with collimated poloidal magnetic field, and high Blandford-Znajek luminosity. We show for the first time this picture in a self-consistent simulation. The turbulence-like motion induced by the non-axisymmetric magnetorotational instability as well as the Kelvin-Helmholtz instability inside the accretion torus works as an agent to drive the mass accretion and converts the accretion energy to thermal energy, which results in the generation of a strong wind. By an in-depth resolution study, we reveal that high resolution is essential to draw such a picture. We also discuss the implication for the r-process nucleosynthesis, the radioactively-powered transient emission, and short gamma-ray bursts.

The gravitational radiation emitted during the merger of a black hole with a neutron star is rather similar to the radiation from the merger of two black holes when the neutron star is not tidally disrupted. When tidal disruption occurs, gravitational waveforms can be broadly classified in two groups, depending on the spatial extent of the disrupted material. Extending previous work by some of us, here we present a phenomenological model for the gravitational waveform amplitude in the frequency domain encompassing the three possible outcomes of the merger: no tidal disruption, "mild" and "strong" tidal disruption. The model is calibrated to 134 general-relativistic numerical simulations of binaries where the black hole spin is either aligned or antialigned with the orbital angular momentum. All simulations were produced using the SACRA code and piecewise polytropic neutron star equations of state. The present model can be used to determine when black-hole binary waveforms are sufficient for gravitational-wave detection, to extract information on the equation of state from future gravitational-wave observations, to obtain more accurate estimates of black hole-neutron star merger event rates, and to determine the conditions under which these systems are plausible candidates as central engines of gamma-ray bursts, macronovae and kilonovae.

We investigate properties of material ejected dynamically in the merger of black hole-neutron star binaries by numerical-relativity simulations. We systematically study the dependence of ejecta properties on the mass ratio of the binary, spin of the black hole, and equation of state of the neutron-star matter. Dynamical mass ejection is driven primarily by tidal torque, and the ejecta is much more anisotropic than that from binary neutron star mergers. In particular, the dynamical ejecta is concentrated around the orbital plane with a half opening angle of 10--20deg and often sweeps out only a half of the plane. The ejecta mass can be as large as ~0.1M_sun, and the velocity is subrelativistic with ~0.2--0.3c for typical cases. The ratio of the ejecta mass to the bound mass (disk and fallback components) is larger, and the ejecta velocity is larger, for larger values of the binary mass ratio, i.e., for larger values of the black-hole mass. The remnant black hole-disk system receives a kick velocity of O(100)km/s due to the ejecta linear momentum, and this easily dominates the kick velocity due to gravitational radiation. Structures of postmerger material, velocity distribution of the dynamical ejecta, fallback rates, and gravitational waves are also investigated. We also discuss the effect of ejecta anisotropy on electromagnetic counterparts, specifically a macronova/kilonova and synchrotron radio emission, developing analytic models.

We systematically performed numerical-relativity simulations for black hole (BH) - neutron star (NS) binary mergers with a variety of the BH spin orientation and equations of state (EOS) of the NS. The initial misalignment angles of the BH spin are chosen in the range of i_tilt,0 = 30--90[deg.]. We employed four models of NS EOS with which the compactness of the NS is in the range of C = M_NS/R_NS = 0.138--0.180, where M_NS and R_NS are the mass and the radius of the NS, respectively. The mass ratio of the BH to the NS, Q = M_BH/M_NS, and the dimensionless spin parameter of the BH, chi, are chosen to be Q = 5 and chi = 0.75, together with M_NS = 1.35 M_sun. We obtain the following results: (i) The inclination angle of i_tilt,0 < 70[deg.] and i_tilt,0 < 50[deg.] are required for the formation of a remnant disk with its mass larger than 0.1 M_sun for the case C = 0.140 and C = 0.160, respectively, while the disk mass is always smaller than 0.1M_sun for C = 0.175. The ejecta with its mass larger than 0.01 M_sun is obtained for i_tilt,0 < 85[deg.] with C = 0.140, for i_tilt,0 < 65[deg.] with C = 0.160, and for i_tilt,0 < 30[deg.] with C = 0.175. (ii) The rotational axis of the dense part of the remnant disk is approximately aligned with the remnant BH spin for i_tilt,0 = 30[deg.]. On the other hand, the disk axis is misaligned initially with ~ 30[deg.] for i_tilt,0 = 60[deg.], and the alignment with the remnant BH spin is achieved at ~ 50--60 ms after the onset of merger. The accretion time scale of the remnant disk is typically ~ 100 ms and depends only weakly on the misalignment angle and the EOS. (iii) The ejecta velocity is typically ~ 0.2--0.3c and depends only weakly on i_tilt,0 and the EOS of the NS, while the morphology of the ejecta depends on its mass. (iv) The gravitational-wave spectra contains the information of the NS compactness in the cutoff frequency for i_tilt,0 < 60[deg.].

We perform radiation-hydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese "K" supercomputer, taking into account neutrino cooling and heating by an updated leakage-plus-transfer scheme for the first time. Neutron stars are modeled by three modern finite-temperature equations of state (EOS) developed by Hempel and his collaborators. We find that the properties of the dynamical ejecta of the merger such as total mass, average electron fraction, and thermal energy depend strongly on the EOS. Only for a soft EOS (the so-called SFHo), the ejecta mass exceeds $0.01M_{\odot}$. In this case, the distribution of the electron fraction of the ejecta becomes broad due to the shock heating during the merger. These properties are well-suited for the production of the solar-like $r$-process abundance. For the other stiff EOS (DD2 and TM1), for which a long-lived massive neutron star is formed after the merger, the ejecta mass is smaller than $0.01M_{\odot}$, although broad electron-fraction distributions are achieved by the positron capture and the neutrino heating.

We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.

We calculate quasiequilibrium sequences of equal-mass, irrotational binary neutron stars (BNSs) in a scalar-tensor (ST) theory of gravity that admits dynamical scalarization. We model neutron stars with realistic equations of state (notably through piecewise polytropic equations of state). Using these quasiequilibrium sequences we compute the binary's scalar charge and binding energy versus orbital angular frequency. We find that the absolute value of the binding energy is smaller than in general relativity (GR), differing at most by ~14% at high frequencies for the cases considered. We use the newly computed binding energy and the balance equation to estimate the number of gravitational-wave (GW) cycles during the adiabatic, quasicircular inspiral stage up to the end of the sequence, which is the last stable orbit or the mass-shedding point, depending on which comes first. We find that, depending on the ST parameters, the number of GW cycles can be substantially smaller than in GR. In particular, we obtain that when dynamical scalarization sets in around a GW frequency of ~130 Hz, the sole inclusion of the ST binding energy causes a reduction of GW cycles from ~120 Hz up to the end of the sequence (~1200 Hz) of ~11% with respect to the GR case. We estimate that when the ST energy flux is also included the reduction in GW cycles becomes of ~24%. Quite interestingly, dynamical scalarization can produce a difference in the number of GW cycles with respect to the GR point-particle case that is much larger than the effect due to tidal interactions, which is on the order of only a few GW cycles. These results further clarify and confirm recent studies that have evolved BNSs either in full numerical relativity or in post-Newtonian theory, and point out the importance of developing accurate ST-theory waveforms for systems composed of strongly self-gravitating objects, such as BNSs.

We develop a method to compute low-eccentricity initial data of binary neutron stars required to perform realistic simulations in numerical relativity. The orbital eccentricity is controlled by adjusting the orbital angular velocity of a binary and incorporating an approaching relative velocity of the neutron stars. These modifications improve the solution primarily through the hydrostatic equilibrium equation for the binary initial data. The orbital angular velocity and approaching velocity of initial data are updated iteratively by performing time evolutions over ~3 orbits. We find that the eccentricity can be reduced by an order of magnitude compared to standard quasicircular initial data, specifically from ~0.01 to <~0.001, by three successive iterations for equal-mass binaries leaving ~10 orbits before the merger.

We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70\,m, which yields the highest-resolution results among those derived so far. By an in-depth resolution study, we clarify several amplification mechanisms of magnetic fields during the binary neutron star merger for the first time. First, the Kelvin-Helmholtz instability developed in the shear layer at the onset of the merger significantly amplifies the magnetic fields. A hypermassive neutron star (HMNS) formed after the merger is then subject to the nonaxisymmetric magnetorotational instability, which amplifies the magnetic field in the HMNS. These two amplification mechanisms cannot be found with insufficient-resolution runs. We also show that the HMNS eventually collapses to a black hole surrounded by an accretion torus which is strongly magnetized at birth.

Recent studies suggest that binary neutron star (NS-NS) mergers robustly produce the heavy r-process nuclei above the atomic mass number A ~ 130 because of their ejecta consisting of almost pure neutrons (electron fraction of Y_e < 0.1). However, little production of the lighter r-process nuclei (A = 90-120) conflicts with the spectroscopic results of r-process-enhanced Galactic halo stars. We present, for the first time, the result of nucleosynthesis calculations based on the fully general-relativistic simulation of a NS-NS merger with approximate neutrino transport. It is found that the bulk of the dynamical ejecta are appreciably shock-heated and neutrino-processed, resulting in a wide range of Y_e (= 0.09-0.45). The mass-averaged abundance distribution of calculated nucleosynthesis yields is in reasonable agreement with the full-mass range (A = 90-240) of the solar r-process curve. This implies, if our model is representative of such events, that the dynamical ejecta of NS-NS mergers can be the origin of the Galactic r-process nuclei. Our result also shows that the radioactive heating after ~ 1 day from the merging, giving rise to r-process-powered transient emission, is dominated by the beta-decays of several species close to stability with precisely measured half-lives. This implies that the total radioactive heating rate for such an event can be well constrained within about a factor of two if the ejected material has a solar-like r-process pattern.

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10{\deg}. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

We propose a possibility of ultrarelativistic electromagnetic counterparts to gravitational waves from binary neutron star mergers at nearly all the viewing angles. Our proposed mechanism relies on the merger-shock propagation accelerating a smaller mass in the outer parts of the neutron star crust to a larger Lorentz factor $\Gamma$ with smaller energy $\sim 10^{47} \Gamma^{-1}$ erg. This mechanism is difficult to resolve by current 3D numerical simulations. The outflows emit synchrotron flares for seconds to days by shocking the ambient medium. Ultrarelativistic flares shine at an early time and in high-energy bands, potentially detectable by current X-ray to radio instruments, such as Swift XRT and Pan-STARRS, and even in low ambient density $\sim 10^{-2}$ cm$^{-3}$ by EVLA. The flares probe the merger position and time, and the merger types as black hole--neutron star outflows would be non-/mildly relativistic.

Detection of electromagnetic counterparts of gravitational wave (GW) sources is important to unveil the nature of compact binary coalescences. We perform three-dimensional, time-dependent, multi-frequency radiative transfer simulations for radioactively powered emission from the ejecta of black hole (BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of the BH, and equations of state of dense matter, BH-NS mergers can eject more material than NS-NS mergers. In such cases, radioactively powered emission from the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We show that, in spite of the expected larger distances to BH-NS merger events, observed brightness of BH-NS mergers can be comparable to or even higher than that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger ejecta are confined to a small solid angle, the emission from BH-NS merger ejecta tends to be bluer than that from NS-NS merger ejecta for a given total luminosity. Thanks to this property, we might be able to distinguish BH-NS merger events from NS-NS merger events by multi-band observations of the radioactively powered emission. In addition to the GW observations, such electromagnetic observations can potentially provide independent information on the nature of compact binary coalescences.

Black hole-neutron star binary mergers display a much richer phenomenology than black hole-black hole mergers, even in the relatively simple case - considered in this paper - in which both the black hole and the neutron star are nonspinning. When the neutron star is tidally disrupted, the gravitational wave emission is radically different from the black hole-black hole case and it can be broadly classified in two groups, depending on the spatial extent of the disrupted material. We present a phenomenological model for the gravitational waveform amplitude in the frequency domain that encompasses the three possible outcomes of the merger: no tidal disruption, "mild", and "strong" tidal disruption. The model is calibrated to general relativistic numerical simulations using piecewise polytropic neutron star equations of state. It should prove useful to extract information on the nuclear equation of state from future gravitational-wave observations, and also to obtain more accurate estimates of black hole-neutron star merger event rates in second- and third-generation interferometric gravitational-wave detectors. We plan to extend and improve the model as longer and more accurate gravitational waveforms become available, and we will make it publicly available online as a Mathematica package. We also present in appendix analytical fits of the projected KAGRA noise spectral density, that should be useful in data analysis applications.

An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.

We carry out numerical-relativity simulations of coalescing binary neutron stars in a scalar-tensor theory that admits spontaneous scalarization. We model neutron stars with realistic equations of state. We choose the free parameters of the theory taking into account the constraints imposed by the latest observations of neutron-star-- white-dwarf binaries with pulsar timing. We show that even within those severe constraints, scalarization can still affect the evolution of the binary neutron stars not only during the late inspiral, but also during the merger stage. We also confirm that even when both neutron stars have quite small scalar charge at large separations, they can be strongly scalarized dynamically during the final stages of the inspiral. In particular, we identify the binary parameters for which scalarization occurs either during the late inspiral or only after the onset of the merger when a remnant, supramassive or hypermassive neutron star is formed. We also discuss how those results can impact the extraction of physical information on gravitational waves once they are detected.

The merger of black hole-neutron star binaries can eject substantial material with the mass ~0.01-0.1M_sun when the neutron star is disrupted prior to the merger. The ejecta shows significant anisotropy, and travels in a particular direction with the bulk velocity ~0.2c. This is drastically different from the binary neutron star merger, for which ejecta is nearly isotropic. Anisotropic ejecta brings electromagnetic-counterpart diversity which is unique to black hole-neutron star binaries, such as viewing-angle dependence, polarization, and proper motion. The kick velocity of the black hole, gravitational-wave memory emission, and cosmic-ray acceleration are also discussed.

Massive (hypermassive and supramassive) neutron stars are likely to be often formed after the merger of binary neutron stars. We explore the evolution process of the remnant massive neutron stars and gravitational waves emitted by them, based on numerical-relativity simulations for binary neutron star mergers employing a variety of equations of state and choosing a plausible range of the neutron-star mass of binaries. We show that the lifetime of remnant hypermassive neutron stars depends strongly on the total binary mass and also on the equations of state. Gravitational waves emitted by the remnant massive neutron stars universally have a quasiperiodic nature of an approximately constant frequency although the frequency varies with time. We also show that the frequency and time-variation feature of gravitational waves depend strongly on the equations of state. We derive a fitting formula for the quasiperiodic gravitational waveforms, which may be used for the data analysis of a gravitational-wave signal.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter {\Lambda} that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with {\Lambda}, and changes in the EOS that did not change {\Lambda} are not measurable. We estimate with two methods the minimal and expected measurability of {\Lambda} in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.

A new EAS hybrid experiment has been designed by constructing a YAC (Yangbajing Air shower Core) detector array inside the existing Tibet-III air shower array. The first step of YAC, called "YAC-I", consists of 16 plastic scintillator units (4 rows times 4 columns) each with an area of 40 cm * 50 cm which is used to check hadronic interaction models used in AS simulations. A Monte Carlo study shows that YAC-I can record high energy electromagnetic component in the core region of air showers induced by primary particles of several tens TeV energies where the primary composition is directly measured by space experiments. It may provide a direct check of the hadronic interaction models currently used in the air shower simulations in the corresponding energy region. In present paper, the method of the observation and the sensitivity of the characteristics of the observed events to the different interaction models are discussed.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We study gravitational waves emitted in the late inspiral stage of binary neutron stars by analyzing the waveform obtained in numerical-relativity simulations. For deriving the physical gravitational waveforms from the numerical results, the resolution extrapolation plays an essential role for our simulations. The extrapolated gravitational-wave phases are compared with those calculated in the post-Newtonian (PN) and effective-one-body (EOB) formalisms including corrections of tidal effects. We show that the extrapolated gravitational-wave phases in numerical relativity agree well with those by the PN and EOB calculations for most of the inspiral stage except for a tidally-dominated, final inspiral stage, in which the PN and EOB results underestimate the tidal effects. Nevertheless, the accumulated phase difference between our extrapolated results and the results by the PN/EOB calculations is at most 1--3 radian in the last 15 cycles.

Numerical-relativity simulations for the merger of binary neutron stars are performed for a variety of equations of state (EOSs) and for a plausible range of the neutron-star mass, focusing primarily on the properties of the material ejected from the system. We find that a fraction of the material is ejected as a mildly relativistic and mildly anisotropic outflow with the typical and maximum velocities $\sim 0.15$ -- $0.25c$ and $\sim 0.5$ -- $0.8c$ (where $c$ is the speed of light), respectively, and that the total ejected rest mass is in a wide range $10^{-4}$ -- $10^{-2}M_{\odot}$, which depends strongly on the EOS, the total mass, and the mass ratio. The total kinetic energy ejected is also in a wide range between $10^{49}$ and $10^{51} {\rm ergs}$. The numerical results suggest that for a binary of canonical total mass $2.7M_{\odot}$, the outflow could generate an electromagnetic signal observable by the planned telescopes through the production of heavy-element unstable nuclei via the $r$-process or through the formation of blast waves during the interaction with the interstellar matter, if the EOS and mass of the binary are favorable ones.

We construct a new three-dimensional general relativistic magnetohydrodynamics code, in which a fixed mesh refinement technique is implemented. To ensure the divergence-free condition as well as the magnetic flux conservation, we employ the method by Balsara (2001). Using this new code, we evolve differentially rotating magnetized neutron stars, and find that a magnetically driven outflow is launched from the star exhibiting a kink instability. The matter ejection rate and Poynting flux are still consistent with our previous finding (Shibata et al., 2011) obtained in axisymmetric simulations.

Using our new numerical-relativity code SACRA, long-term simulations for inspiral and merger of black hole (BH)-neutron star (NS) binaries are performed, focusing particularly on gravitational waveforms. As the initial conditions, BH-NS binaries in a quasiequilibrium state are prepared in a modified version of the moving-puncture approach. The BH is modeled by a nonspinning moving puncture and for the NS, a polytropic equation of state with $\Gamma=2$ and the irrotational velocity field are employed. The mass ratio of the BH to the NS, $Q=M_{\rm BH}/M_{\rm NS}$, is chosen in the range between 1.5 and 5. The compactness of the NS, defined by ${\cal C}=GM_{\rm NS}/c^2R_{\rm NS}$, is chosen to be between 0.145 and 0.178. For a large value of $Q$ for which the NS is not tidally disrupted and is simply swallowed by the BH, gravitational waves are characterized by inspiral, merger, and ringdown waveforms. In this case, the waveforms are qualitatively the same as that from BH-BH binaries. For a sufficiently small value of $Q \alt 2$, the NS may be tidally disrupted before it is swallowed by the BH. In this case, the amplitude of the merger and ringdown waveforms is very low, and thus, gravitational waves are characterized by the inspiral waveform and subsequent quick damping. The difference in the merger and ringdown waveforms is clearly reflected in the spectrum shape and in the "cut-off" frequency above which the spectrum amplitude steeply decreases. When an NS is not tidally disrupted (e.g., for Q=5), kick velocity, induced by asymmetric gravitational wave emission, agrees approximately with that derived for the merger of BH-BH binaries, whereas for the case that the tidal disruption occurs, the kick velocity is significantly suppressed.

We construct magnetized stars composed of a fluid stably stratified by entropy gradients in the framework of general relativity, assuming ideal magnetohydrodynamics and employing a barotropic equation of state. We first revisit basic equations for describing stably-stratified stationary axisymmetric stars containing both poloidal and toroidal magnetic fields. As sample models, the magnetized stars considered by Ioka and Sasaki (2004), inside which the magnetic fields are confined, are modified to the ones stably stratified. The magnetized stars newly constructed in this study are believed to be more stable than the existing relativistic models because they have both poloidal and toroidal magnetic fields with comparable strength, and magnetic buoyancy instabilities near the surface of the star, which can be stabilized by the stratification, are suppressed.

We describe the current status of our numerical simulations for the collapse of a massive stellar core to a BH and the BNS mergers, performed in the framework of full general relativity incorporating finite-temperature EOS and neutrino cooling. For the stellar core collapse simulation, we present the latest numerical results. We employed a purely nucleonic EOS (Shen-EOS). As an initial condition, we adopted a 100 $M_{\odot}$ presupernova model calculated by Umeda and Nomoto. Changing the degree of rotation for the initial condition, we clarify the strong dependence of the outcome of the collapse on this. When the rotation is rapid enough, the shock wave formed at the core bounce is deformed to be a torus-like shape. Then, the infalling matter is accumulated in the central region due to the oblique shock at the torus surface, hitting the PNS and dissipating the kinetic energy there. As a result, outflows can be launched. The PNS eventually collapses to a BH and an accretion torus is formed around it. We also found that the evolution of the BH and torus depends strongly on the rotation initially given. In the BNS merger simulations, we in addition employ an EOS incorporating a degree of freedom for hyperons. The numerical simulations show that for the purely nucleonic EOS, a HMNS with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $M \lesssim 3.0M_{\odot}$. By contrast, the formed HMNS collapses to a BH in a shorter time scale with the hyperonic EOS for $M \gtrsim 2.7M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$10\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 2--$6 \times 10^{-22}$ at $f\approx 2$--2.5 kHz for a source distance of 100 Mpc.

A radiation-magnetohydrodynamic simulation for the black hole-torus system is performed in the framework of full general relativity for the first time. A truncated moment formalism is employed for a general relativistic neutrino radiation transport. Several systems in which the black hole mass is $M_{\rm BH}=3$ or $6M_{\odot}$, the black hole spin is zero, and the torus mass is $\approx 0.14$--$0.38M_{\odot}$ are evolved as models of the remnant formed after the merger of binary neutron stars or black hole-neutron star binaries. The equation of state and microphysics for the high-density and high-temperature matter are phenomenologically taken into account in a semi-quantitative manner. It is found that the temperature in the inner region of the torus reaches $\agt 10$ MeV which enhances a high luminosity of neutrinos $\sim 10^{51}$ ergs/s for $M_{\rm BH}=6M_{\odot}$ and $\sim 10^{52}$ ergs/s for $M_{\rm BH}=3M_{\odot}$. It is shown that neutrinos are likely to be emitted primarily toward the outward direction in the vicinity of the rotational axis and their energy density may be high enough to launch a low-energy short gamma-ray burst via the neutrino-antineutrino pair-annihilation process with the total energy deposition $\sim 10^{47}$--$10^{49}$ ergs. It is also shown in our model that for $M_{\rm BH}=3M_{\odot}$, the neutrino luminosity is larger than the electromagnetic luminosity while for $M_{\rm BH}=6M_{\odot}$, the neutrino luminosity is comparable to or slightly smaller than the electromagnetic luminosity.

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling. It is found that for the nucleonic and hyperonic EOS, a hyper massive neutron star (HMNS) with a long lifetime $(t_{\rm life}\gtrsim 10 {\rm ms})$ is the outcome for the total mass $\approx 2.7 M_\odot$. For the total mass $\approx 3 M_\odot$, a long-lived (short-lived with $t_{\rm life}\approx 3 {\rm ms}$) HMNS is the outcome for the nucleonic (hyperonic) EOS. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$ -- $6 \times 10^{53} {\rm erg /s}$ and the effective amplitude of gravitational waves from the HMNS is 1 -- $4\times 10^{-22}$ at $f\approx 2$ -- $3.2 {\rm kHz}$ for a source of distance of 100 Mpc. During the HMNS phase, characteristic frequencies of gravitational waves shift to a higher frequency for the hyperonic EOS in contrast to the nucleonic EOS in which they remain constant approximately. Our finding suggests that the effects of hyperons are well imprinted in gravitational wave and its detection will give us a potential opportunity to explore the composition of the neutron star matter. We present the neutrino luminosity curve when a black hole is formed as well.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling for the first time. It is found that even for the hyperonic EOS, a hypermassive neutron star is first formed after the merger for the typical total mass $\approx$ 2.7M\bigodot, and subsequently collapses to a black hole (BH). It is shown that hyperons play a substantial role in the post-merger dynamics, torus formation around the BH, and emission of gravitational waves (GWs). In particular, the existence of hyperons is imprinted in GWs. Therefore, GW observations will provide a potential opportunity to explore the composition of the neutron star matter.

Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.

We study the merger of black hole (BH)-neutron star (NS) binaries with a variety of BH spins aligned or anti-aligned with the orbital angular momentum, and with the mass ratio in the range MBH/MNS = 2--5, where MBH and MNS are the mass of the BH and NS, respectively. We model NS matter by systematically parametrized piecewise polytropic equations of state. The initial condition is computed in the puncture framework adopting an isolated horizon framework to estimate the BH spin and assuming an irrotational velocity field for the fluid inside the NS. Dynamical simulations are performed in full general relativity by an adaptive mesh refinement code, SACRA. The treatment of hydrodynamic equations and estimation of the disk mass are improved. We find that the NS is tidally disrupted irrespective of the mass ratio when the BH has a moderately large prograde spin, whereas only binaries with low mass ratios, MBH/MNS <~ 3 or small compactnesses of the NSs, bring the tidal disruption when the BH spin is zero or retrograde. The mass of the remnant disk is accordingly large as >~ 0.1 Msun, which is required by central engines of short gamma-ray bursts, if the BH spin is prograde. Information of the tidal disruption is reflected in a clear relation between the compactness of the NS and an appropriately defined "cutoff frequency" in the gravitational-wave spectrum, above which the spectrum damps exponentially. We find that the tidal disruption of the NS and excitation of the quasinormal mode of the remnant BH occur in a compatible manner in high mass-ratio binaries with the prograde BH spin. The correlation between the compactness and the cutoff frequency still holds for such cases. It is also suggested by extrapolation that the merger of an extremely spinning BH and an irrotational NS binary does not lead to formation of an overspinning BH.

The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters. To within the accuracy of the simulations, we find that, apart from the neutron-star mass, a single physical parameter Lambda, describing its deformability, can be extracted from the late inspiral, merger, and ringdown waveform. This parameter is related to the radius, mass, and l=2 Love number, k_2, of the neutron star by Lambda = 2k_2 R^5/3M_{NS}^5, and it is the same parameter that determines the departure from point-particle dynamics during the early inspiral. Observations of gravitational waves from BHNS inspiral thus restrict the EOS to a surface of constant Lambda in the parameter space, thickened by the measurement error. Using various configurations of a single Advanced LIGO detector, we find that neutron stars are distinguishable from black holes of the same mass and that Lambda^{1/5} or equivalently R can be extracted to 10-40% accuracy from single events for mass ratios of Q=2 and 3 at a distance of 100 Mpc, while with the proposed Einstein Telescope, EOS parameters can be extracted to accuracy an order of magnitude better.

In February 2011, a burst event of the H$_{2}$O maser in Orion KL (Kleinmann-Low object) has started after 13-year silence. This is the third time to detect such phenomena in Orion KL, followed by those in 1979-1985 and 1998. We have carried out astrometric observations of the bursting H$_{2}$O maser features in Orion KL with VERA (VLBI Exploration of Radio Astrometry), a Japanese VLBI network dedicated for astrometry. The total flux of the bursting feature at the LSR velocity of 7.58 km s$^{-1}$ reaches 4.4$\times10^{4}$ Jy in March 2011. The intensity of the bursting feature is three orders of magnitudes larger than that of the same velocity feature in the quiescent phase in 2006. Two months later, another new feature appears at the LSR velocity of 6.95 km s$^{-1}$ in May 2011, separated by 12 mas north of the 7.58 km s$^{-1}$ feature. Thus, the current burst occurs at two spatially different features. The bursting masers are elongated along the northwest-southeast direction as reported in the previous burst in 1998. We determine the absolute positions of the bursting features for the first time ever with a submilli-arcsecond (mas) accuracy. Their positions are coincident with the shocked molecular gas called the Orion Compact Ridge. We tentatively detect the absolute proper motions of the bursting features toward southwest direction. It is most likely that the outflow from the radio source I or another young stellar object interacting with the Compact Ridge is a possible origin of the H$_{2}$O maser burst.

We report results of a numerical-relativity simulation for the merger of a black hole-neutron star binary with a variety of equations of state (EOSs) modeled by piecewise polytropes. We focus in particular on the dependence of the gravitational waveform at the merger stage on the EOSs. The initial conditions are computed in the moving-puncture framework, assuming that the black hole is nonspinning and the neutron star has an irrotational velocity field. For a small mass ratio of the binaries (e.g., MBH/MNS = 2 where MBH and MNS are the masses of the black hole and neutron star, respectively), the neutron star is tidally disrupted before it is swallowed by the black hole irrespective of the EOS. Especially for less-compact neutron stars, the tidal disruption occurs at a more distant orbit. The tidal disruption is reflected in a cutoff frequency of the gravitational-wave spectrum, above which the spectrum amplitude exponentially decreases. A clear relation is found between the cutoff frequency of the gravitational-wave spectrum and the compactness of the neutron star. This relation also depends weakly on the stiffness of the EOS in the core region of the neutron star, suggesting that not only the compactness but also the EOS at high density is reflected in gravitational waveforms. The mass of the disk formed after the merger shows a similar correlation with the EOS, whereas the spin of the remnant black hole depends primarily on the mass ratio of the binary, and only weakly on the EOS. Properties of the remnant disks are also analyzed.

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $\alt 3.0M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$8\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 4--$6 \times 10^{-22}$ at $f=2.1$--2.5 kHz for a source distance of 100 Mpc. We also present the neutrino luminosity curve when a black hole is formed for the first time.

Black hole (BH)--torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observable by forthcoming ground-based and spacecraft detectors.

We present the first numerical result of fully general relativistic axisymmetric simulations for the collapse of a rotating high-entropy stellar core to a black hole and an accretion disk. The simulations are performed taking into account the relevant microphysics. We adopt as initial condition a spherical core with constant electron fraction ($Y_e = 0.5$) and entropy per baryon $s$ = 8 $k_B$, and angular velocity is superimposed. In the early phase, the core collapses in a homologous manner. Then, it experiences a weak bounce due to the gas pressure of free nucleons. Because the bounce is weak, the core collapses eventually to a black hole. Subsequent evolution depends on initial angular velocity. When the rotation is not fast, a geometrically thin (but optically thick) accretion disk is formed, and shock waves are formed in the inner part of the disk. For the moderately rotating case, the thin accretion disk expands eventually to be a geometrically thick torus after sufficient accumulation of the thermal energy generated at the shocks. Furthermore, convection occurs inside the torus. Neutrino luminosities vary violently with time because of the convective motion. For the rapidly rotating case, by contrast, a geometrically thick torus is formed soon after the black hole formation, and convective activity is weak due to the presence of epicyclic mode.

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one \sigma level. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.

The merger of two neutron stars results often in a rapidly and differentially rotating hypermassive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromagnetic radiation with a large luminosity ~ 0.1 B^2 R^3 Omega is emitted with baryon (B, R, and Omega are poloidal magnetic-field strength at stellar surface, stellar radius, and angular velocity of a HMNS). The predicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole direction, is ~ 10^{47} (B/10^{13} G)^2(R/10 km)^3(Omega/10^4 rad/s) ergs/s, that is comparable to the luminosity of quasars.

The aim of this paper is to clarify the stabilities of neutron stars with strong toroidal magnetic fields against non-axisymmetric perturbation. The motivation comes from the fact that super magnetized neutron stars of $\sim 10^{15}$G, magnetars, and magnetized proto-neutron stars born after the magnetically-driven supernovae are likely to have such strong toroidal magnetic fields. Long-term, three-dimensional general relativistic magneto-hydrodynamic simulations are performed, preparing isentropic neutron stars with toroidal magnetic fields in equilibrium as initial conditions. To explore the effects of rotations on the stability, simulations are done for both non-rotating and rigidly rotating models. We find the emergence of the Parker and/or Tayler instabilities in both the non-rotating and rotating models. For both non-rotating and rotating models, the Parker instability is the primary instability as predicted by the local linear perturbation analysis. The interchange instability also appears in the rotating models. It is found that rapid rotation is not enough to suppress the Parker instability, and this finding does not agree with the perturbation analysis. The reason for this is that rigidly and rapidly rotating stars are marginally stable, and hence, in the presence of stellar pulsations by which the rotational profile is deformed, unstable regions with negative gradient of angular momentum profile is developed. After the onset of the instabilities, a turbulence is excited. Contrary to the axisymmetric case, the magnetic fields never reach an equilibrium state after the development of the turbulence. This conclusion suggests that three-dimensional simulation is indispensable for exploring the formation of magnetars or prominence activities of magnetars such as giant flares.

A truncated moment formalism for general relativistic radiation hydrodynamics, based on the Thorne's moment formalism, is derived. The fluid rest frame is chosen to be the fiducial frame for defining the radiation moments. Then, zeroth-, first-, and second-rank radiation moments are defined from the distribution function with a physically reasonable assumption for it in the optically thin and thick limits. The source terms are written, focusing specifically on the neutrino transfer and neglecting higher harmonic angular dependence of the reaction angle. Finally, basic equations for a truncated moment formalism for general relativistic radiation hydrodynamics in a closed covariant form are derived assuming a closure relation among the radiation stress tensor, energy density, and energy flux, and a variable Eddington factor, which works well.

We present the first quantitative comparison of two independent general-relativistic hydrodynamics codes, the Whisky code and the SACRA code. We compare the output of simulations starting from the same initial data and carried out with the configuration (numerical methods, grid setup, resolution, gauges) which for each code has been found to give consistent and sufficiently accurate results, in particular in terms of cleanness of gravitational waveforms. We focus on the quantities that should be conserved during the evolution (rest mass, total mass energy, and total angular momentum) and on the gravitational-wave amplitude and frequency. We find that the results produced by the two codes agree at a reasonable level, with variations in the different quantities but always at better than about 10%.

The properties of neutron star matter above nuclear density are not precisely known. Gravitational waves emitted from binary neutron stars during their late stages of inspiral and merger contain imprints of the neutron-star equation of state. Measuring departures from the point-particle limit of the late inspiral waveform allows one to measure properties of the equation of state via gravitational wave observations. This and a companion talk by J. S. Read reports a comparison of numerical waveforms from simulations of inspiraling neutron-star binaries, computed for equations of state with varying stiffness. We calculate the signal strength of the difference between waveforms for various commissioned and proposed interferometric gravitational wave detectors and show that observations at frequencies around 1 kHz will be able to measure a compactness parameter and constrain the possible neutron-star equations of state.

Quasi-equilibrium sequences of binary neutron stars are constructed for a variety of equations of state in general relativity. Einstein's constraint equations in the Isenberg-Wilson-Mathews approximation are solved together with the relativistic equations of hydrostationary equilibrium under the assumption of irrotational flow. We focus on unequal-mass sequences as well as equal-mass sequences, and compare those results. We investigate the behavior of the binding energy and total angular momentum along a quasi-equilibrium sequence, the endpoint of sequences, and the orbital angular velocity as a function of time, changing the mass ratio, the total mass of the binary system, and the equation of state of a neutron star. It is found that the orbital angular velocity at the mass-shedding limit can be determined by an empirical formula derived from an analytic estimation. We also provide tables for 160 sequences which will be useful as a guideline of numerical simulations for the inspiral and merger performed in the near future.

Results from the first fully general relativistic numerical simulations in axisymmetry of a system formed by a black hole surrounded by a self-gravitating torus in equilibrium are presented, aiming to assess the influence of the torus self-gravity on the onset of the runaway instability. We consider several models with varying torus-to-black hole mass ratio and angular momentum distribution orbiting in equilibrium around a non-rotating black hole. The tori are perturbed to induce the mass transfer towards the black hole. Our numerical simulations show that all models exhibit a persistent phase of axisymmetric oscillations around their equilibria for several dynamical timescales without the appearance of the runaway instability, indicating that the self-gravity of the torus does not play a critical role favoring the onset of the instability, at least during the first few dynamical timescales.

We elucidate the feature of gravitational waves (GWs) from binary neutron star merger collapsing to a black hole by general relativistic simulation. We show that GW spectrum imprints the coalescence dynamics, formation process of disk, equation of state for neutron stars, total masses, and mass ratio. A formation mechanism of the central engine of short $\gamma$-ray bursts, which are likely to be composed of a black hole and surrounding disk, therefore could be constrained by GW observation.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

General relativistic simulations for the merger of binary neutron stars are performed as an extension of a previous work\cite{Shibata:2006nm}. We prepare binary neutron stars with a large initial orbital separation and employ the moving-puncture formulation, which enables to follow merger and ringdown phases for a long time, even after black hole formation. For modeling inspiraling neutron stars, which should be composed of cold neutron stars, the Akmal-Pandhalipande-Ravenhall (APR) equation of state (EOS) is adopted. After the onset of the merger, the hybrid-type EOS is used; i.e., the cold and thermal parts are given by the APR and $\Gamma$-law EOSs, respectively. Three equal-mass binaries each with mass $1.4M_\odot,1.45M_\odot,1.5M_\odot$ and two unequal-mass binaries with mass 1.3--$1.6M_\odot$, 1.35--$1.65M_\odot$ are prepared. We focus primarily on the black hole formation case, and explore mass and spin of the black hole, mass of disks which surround the black hole, and gravitational waves emitted during the black hole formation. We find that (i) for the systems of $m_0=2.9$--$3.0M_\odot$ and of mass ratio $\approx 0.8$, the mass of disks which surround the formed black hole is 0.006--$0.02M_{\odot}$; (ii) the spin of the formed black hole is $0.78 \pm 0.02$ when a black hole is formed after the merger in the dynamical time scale. This value depends weakly on the total mass and mass ratio, and is about 0.1 larger than that of a black hole formed from nonspinning binary black holes; (iii) for the black-hole formation case, Fourier spectrum shape of gravitational waves emitted in the merger and ringdown phases has a universal qualitative feature irrespective of the total mass and mass ratio, but quantitatively, the spectrum reflects the parameters of the binary neutron stars.

General relativistic quasiequilibrium states of black hole-neutron star binaries are computed in the moving-puncture framework. We propose three conditions for determining the quasiequilibrium states and compare the numerical results with those obtained in the excision framework. We find that the results obtained in the moving-puncture framework agree with those in the excision framework and with those in the third post-Newtonian approximation for the cases that (i) the mass ratio of the binary is close to unity irrespective of the orbital separation, and (ii) the orbital separation is large enough ($m_0\Omega \alt 0.02$ where $m_0$ and $\Omega$ are the total mass and the orbital angular velocity, respectively) irrespective of the mass ratio. For $m_0 \Omega \agt 0.03$, both of the results in the moving-puncture and excision frameworks deviate, more or less, from those in the third post-Newtonian approximation. Thus the numerical results do not provide a quasicircular state, rather they seem to have a nonnegligible eccentricity of order 0.01--0.1. We show by numerical simulation that a method in the moving-puncture framework can provide approximately quasicircular states in which the eccentricity is by a factor of $\sim 2$ smaller than those in quasiequilibrium given by other approaches.

We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.

We have observed cosmic-ray electrons from 10 GeV to 800 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an average altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates with 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of Iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. From the long duration balloon observations, we derived the energy spectrum of cosmic-ray electrons in the energy range from 100 GeV to 800 GeV. In addition, for the first time we derived the electron arrival directions above 100 GeV, which is consistent with the isotropic distribution.

We report our new code (named SACRA) for numerical relativity simulations in which an adaptive mesh refinement algorithm is implemented. In this code, the Einstein equations are solved in the BSSN formalism with a fourth-order finite differencing, and the hydrodynamic equations are solved by a third-order high-resolution central scheme. The fourth-order Runge-Kutta scheme is adopted for integration in time. To test the code, simulations for coalescence of black hole-black hole (BH-BH), neutron star-neutron star (NS-NS), and black hole-neutron star (BH-NS) binaries are performed, and also, properties of BHs formed after the merger and gravitational waveforms are compared among those three cases. For the simulations of BH-BH binaries, we adopt the same initial conditions as those by Buonanno et al. and compare numerical results. We find reasonable agreement except for a slight disagreement possibly associated with the difference in choice of gauge conditions and numerical schemes. For an NS-NS binary, we performed simulations employing both SACRA and Shibata's previous code, and find reasonable agreement. For a BH-NS binary, we compare numerical results with our previous ones, and find that gravitational waveforms and properties of the BH formed after the merger agree well with those of our previous ones, although the disk mass formed after the merger is less than 0.1% of the total rest mass, which disagrees with the previous result. We also report numerical results of a longterm simulation (with $\sim 4$ orbits) for a BH-NS binary for the first time. All these numerical results show behavior of convergence, and extrapolated numerical results for time spent in the inspiral phase agree with post-Newtonian predictions in a reasonable accuracy.

We study the stability of neutron stars with toroidal magnetic fields by magnetohydrodynamic simulation in full general relativity under assumption of axial symmetry. Nonrotating and rigidly rotating neutron stars are prepared for a variety of magnetic field configuration. For modeling the neutron stars, the polytropic equation of state with the adiabatic index $\Gamma=2$ is used for simplicity. It is found that nonrotating neutron stars are dynamically unstable for the case that toroidal magnetic field strength varies $\propto \varpi^{2k-1}$ with $k\geq 2$ (here $\varpi$ is the cylindrical radius), whereas for $k=1$ the neutron stars are stable. After the onset of the instability, unstable modes grow approximately in the Alfv\'en time scale and, as a result, a convective motion is excited to change the magnetic field profile until a new state, which is stable against axisymmetric perturbation, is reached. We also find that rotation plays a role in stabilization, although the instability still sets in in the Alfv\'en time scale when the ratio of magnetic energy to rotational kinetic energy is larger than a critical value $\sim 0.2$. Implication for the evolution of magnetized protoneutron stars is discussed.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

We systematically perform the merger simulation of black hole-neutron star (BH-NS) binaries in full general relativity, focusing on the case that the NS is tidally disrupted. We prepare BH-NS binaries in a quasicircular orbit as the initial condition in which the BH is modeled by a nonspinning moving puncture. For modeling the NS, we adopt the $\Gamma$-law equation of state with $\Gamma=2$ and the irrotational velocity field. We change the BH mass in the range $M_{\rm BH} \approx 3.3$--$4.6M_{\odot}$, while the rest mass of the NS is fixed to be $M_{*}=1.4 M_{\odot}$ (i.e., the NS mass $M_{\rm NS} \approx 1.3M_{\odot}$). The radius of the corresponding spherical NS is set in the range $R_{\rm NS} \approx 12$--15 km (i.e., the compactness $GM_{\rm NS}/R_{\rm NS}c^2 \approx 0.13$--0.16). We find for all the chosen initial conditions that the NS is tidally disrupted near the innermost stable circular orbit. For the model of $R_{\rm NS}=12$ km, more than 97 % of the rest mass is quickly swallowed into the BH and the resultant torus mass surrounding the BH is less than $0.04M_{\odot}$. For the model of $R_{\rm NS} \approx 14.7$ km, by contrast, the torus mass is about $0.16M_{\odot}$ for the BH mass $\approx 4M_{\odot}$. The thermal energy of the material in the torus increases by the shock heating occurred in the collision between the spiral arms, resulting in the temperature $10^{10}$--$10^{11}$ K. (.. omission ..) We also present gravitational waveforms during the inspiral, tidal disruption of the NS, and subsequent evolution of the disrupted material. (.. omission ..)

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

We present our first numerical results of axisymmetric magnetohydrodynamic simulations for neutrino-cooled accretion tori around rotating black holes in general relativity. We consider tori of mass $\sim 0.1$--0.4$M_{\odot}$ around a black hole of mass $M=4M_{\odot}$ and spin $a=0$--$0.9M$; such systems are candidates for the central engines of gamma-ray bursts (GRBs) formed after the collapse of massive rotating stellar cores and the merger of a black hole and a neutron star. In this paper, we consider the short-term evolution of a torus for a duration of $\approx 60$ ms, focusing on short-hard GRBs. Simulations were performed with a plausible microphysical equation of state that takes into account neutronization, the nuclear statistical equilibrium of a gas of free nucleons and $\alpha$-particles, black body radiation, and a relativistic Fermi gas (neutrinos, electrons, and positrons). Neutrino-emission processes, such as $e^{\pm}$ capture onto free nucleons, $e^{\pm}$ pair annihilation, plasmon decay, and nucleon-nucleon bremsstrahlung are taken into account as cooling processes. Magnetic braking and the magnetorotational instability in the accretion tori play a role in angular momentum redistribution, which causes turbulent motion, resultant shock heating, and mass accretion onto the black hole. The mass accretion rate is found to be $\dot M_* \sim 1$--$10 M_{\odot}$/s, and the shock heating increases the temperature to $\sim 10^{11}$ K. This results in a maximum neutrino emission rate of $L_{\nu}=$ several $\times 10^{53}$ ergs/s and a conversion efficiency $L_{\nu}/\dot M_* c^2$ on the order of a few percent for tori with mass $M_{\rm t} \approx 0.1$--0.4$M_{\odot}$ and for moderately high black hole spins.

The latest studies of massive star evolution indicate that an initially rapidly rotating star with sufficiently low metallicity can produce a rapidly rotating, massive stellar core that could be a progenitor of long-soft gamma-ray bursts (LGRBs). Motivated by these studies, we follow the collapse of a rapidly rotating massive stellar core to a 'collapsar' black hole (BH) surrounded by a massive, hot accretion disk performing fully general relativistic simulations. We focus on the general relativistic dynamics of the collapse, and the relevant microphysics is treated in a qualitative manner. The simulations are performed until the system consisting of the BH and the disk has relaxed to a quasi-stationary state. A novel mechanism found in this study is that strong shock waves are formed at the inner part of the disk after the formation of the BH. These shock waves propagate mainly along the rotation axis, heating the disk and sweeping materials around the rotational axis, and thereby forming a low density region. The temperature of the disk is high enough for copious neutrino emission. All these features indicate that the direct formation of a rapidly rotating BH is a promising source of LGRBs even in the absence of strong magnetic fields.

The merger of two neutron stars usually produces a remnant with a mass significantly above the single (nonrotating) neutron star maximum mass. In some cases, the remnant will be stabilized against collapse by rapid, differential rotation. MHD-driven angular momentum transport eventually leads to the collapse of the remnant's core, resulting in a black hole surrounded by a massive accretion torus. Here we present simulations of this process. The plausibility of generating short duration gamma ray bursts through this scenario is discussed.

We present our latest results for simulation for merger of black hole (BH)-neutron star (NS) binaries in full general relativity which is performed preparing a quasicircular state as initial condition. The BH is modeled by a moving puncture with no spin and the NS by the $\Gamma$-law equation of state with $\Gamma=2$ and corotating velocity field as a first step. The mass of the BH is chosen to be $\approx 3.2 M_{\odot}$ or $4.0M_{\odot}$, and the rest-mass of the NS $\approx 1.4 M_{\odot}$ with relatively large radius of the NS $\approx 13$--14 km. The NS is tidally disrupted near the innermost stable orbit but $\sim 80$--90% of the material is swallowed into the BH and resulting disk mass is not very large as $\sim 0.3M_{\odot}$ even for small BH mass $\sim 3.2M_{\odot}$. The result indicates that the system of a BH and a massive disk of $\sim M_{\odot}$ is not formed from nonspinning BH-NS binaries irrespective of BH mass, although a disk of mass $\sim 0.1M_{\odot}$ is a possible outcome for this relatively small BH mass range as $\sim 3$--4$M_{\odot}$. Our results indicate that the merger of low-mass BH and NS may form a central engine of short-gamma-ray bursts.

We study magnetohydrodynamic (MHD) effects arising in the collapse of magnetized, rotating, massive stellar cores to proto-neutron stars (PNSs). We perform axisymmetric numerical simulations in full general relativity with a hybrid equation of state. The formation and early evolution of a PNS are followed with a grid of 2500 x 2500 zones, which provides better resolution than in previous (Newtonian) studies. We confirm that significant differential rotation results even when the rotation of the progenitor is initially uniform. Consequently, the magnetic field is amplified both by magnetic winding and the magnetorotational instability (MRI). Even if the magnetic energy E_EM is much smaller than the rotational kinetic energy T_rot at the time of PNS formation, the ratio E_EM/T_rot increases to 0.1-0.2 by the magnetic winding. Following PNS formation, MHD outflows lead to losses of rest mass, energy, and angular momentum from the system. The earliest outflow is produced primarily by the increasing magnetic stress caused by magnetic winding. The MRI amplifies the poloidal field and increases the magnetic stress, causing further angular momentum transport and helping to drive the outflow. After the magnetic field saturates, a nearly stationary, collimated magnetic field forms near the rotation axis and a Blandford-Payne type outflow develops along the field lines. These outflows remove angular momentum from the PNS at a rate given by \dot{J} \sim \eta E_EM C_B, where \eta is a constant of order 0.1 and C_B is a typical ratio of poloidal to toroidal field strength. As a result, the rotation period quickly increases for a strongly magnetized PNS until the degree of differential rotation decreases. Our simulations suggest that rapidly rotating, magnetized PNSs may not give rise to rapidly rotating neutron stars.

Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences that model inspiraling binary neutron stars during the final several orbits are successfully computed. The binding energy of the system near its final orbit deviates from earlier results of third post-Newtonian and of spatially conformally flat calculations. The new solutions may serve as initial data for merger simulations and as members of quasiequilibrium sequences to generate gravitational wave templates, and may improve estimates of the gravitational-wave cutoff frequency set by the last inspiral orbit.

The capacity to model magnetohydrodynamical (MHD) flows in dynamical, strongly curved spacetimes significantly extends the reach of numerical relativity in addressing many problems at the forefront of theoretical astrophysics. We have developed and tested an evolution code for the coupled Einstein-Maxwell-MHD equations which combines a BSSN solver with a high resolution shock capturing scheme. As one application, we evolve magnetized, differentially rotating neutron stars under the influence of a small seed magnetic field. Of particular significance is the behavior found for hypermassive neutron stars (HMNSs), which have rest masses greater the mass limit allowed by uniform rotation for a given equation of state. The remnant of a binary neutron star merger is likely to be a HMNS. We find that magnetic braking and the magnetorotational instability lead to the collapse of HMNSs and the formation of rotating black holes surrounded by massive, hot accretion tori and collimated magnetic field lines. Such tori radiate strongly in neutrinos, and the resulting neutrino-antineutrino annihilation (possibly in concert with energy extraction by MHD effects) could provide enough energy to power short-hard gamma-ray bursts. To explore the range of outcomes, we also evolve differentially rotating neutron stars with lower masses and angular momenta than the HMNS models. Instead of collapsing, the non-hypermassive models form nearly uniformly rotating central objects which, in cases with significant angular momentum, are surrounded by massive tori.

We study the effects of magnetic fields on the evolution of differentially rotating neutron stars, which can form in stellar core collapse or binary neutron star coalescence. Magnetic braking and the magnetorotational instability (MRI) both redistribute angular momentum; the outcome of the evolution depends on the star's mass and spin. Simulations are carried out in axisymmetry using our recently developed codes which integrate the coupled Einstein-Maxwell-MHD equations. For initial data, we consider three categories of differentially rotating, equilibrium configurations, which we label normal, hypermassive and ultraspinning. Hypermassive stars have rest masses exceeding the mass limit for uniform rotation. Ultraspinning stars are not hypermassive, but have angular momentum exceeding the maximum for uniform rotation at the same rest mass. We show that a normal star will evolve to a uniformly rotating equilibrium configuration. An ultraspinning star evolves to an equilibrium state consisting of a nearly uniformly rotating central core, surrounded by a differentially rotating torus with constant angular velocity along magnetic field lines, so that differential rotation ceases to wind the magnetic field. In addition, the final state is stable against the MRI, although it has differential rotation. For a hypermassive neutron star, the MHD-driven angular momentum transport leads to catastrophic collapse of the core. The resulting rotating black hole is surrounded by a hot, massive, magnetized torus undergoing quasistationary accretion, and a magnetic field collimated along the spin axis--a promising candidate for the central engine of a short gamma-ray burst. (Abridged)

Three-dimensional simulations for the merger of binary neutron stars (BNSs) are performed in the framework of full general relativity. We pay particular attention to the black hole (BH) formation case and to the resulting mass of the surrounding disk for exploring possibility for formation of the central engine of short-duration gamma-ray bursts. Hybrid equations of state (EOSs) are adopted mimicking realistic, stiff nuclear EOSs, for which the maximum allowed gravitational mass of cold and spherical neutron stars (NSs), M_sph, is larger than 2M_sun. For the simulations, we focus on BNSs of the ADM mass M>2.6M_sun. For M>M_thr, the merger results in prompt formation of a BH irrespective of the mass ratio Q_M with 0.65<Q_M<1. The value of M_thr is approximately written as 1.3-1.35M_sph for the chosen EOSs. For the BH formation case, we evolve the spacetime using a BH excision technique and determine the mass of a quasistationary disk surrounding the BH. The disk mass steeply increases with decreasing the value of Q_M for given ADM mass and EOS. For M<M_thr, the outcome is a hypermassive neutron star (HMNS) of a large ellipticity. If the HMNS collapses to a BH after the longterm angular momentum transport, the disk mass may be >0.01M_sun. Gravitational waves (GWs) are computed in terms of a gauge-invariant wave extraction technique. In the formation of the HMNS, quasiperiodic GWs of frequency (3-3.5kHz) are emitted. The effective amplitude of GWs can be >5x10^{-21} at a distance of 50 Mpc. For the BH formation case, the BH excision technique enables a longterm computation and extraction of ring-down GWs associated with a BH quasinormal mode. It is found that the frequency and amplitude are 6.5-7kHz and 10^{-22} at a distance of 50Mpc for M=2.7-2.9M_sun.

Hypermassive neutron stars (HMNSs) -- equilibrium configurations supported against collapse by rapid differential rotation -- are possible transient remnants of binary neutron star mergers. Using newly developed codes for magnetohydrodynamic simulations in dynamical spacetimes, we are able to track the evolution of a magnetized HMNS in full general relativity for the first time. We find that secular angular momentum transport due to magnetic braking and the magnetorotational instability results in the collapse of an HMNS to a rotating black hole, accompanied by a gravitational wave burst. The nascent black hole is surrounded by a hot, massive torus undergoing quasistationary accretion and a collimated magnetic field. This scenario suggests that HMNS collapse is a possible candidate for the central engine of short gamma-ray bursts.

A hypermassive neutron star (HMNS) is a possible transient formed after the merger of a neutron star binary. In the latest magnetohydrodynamic simulations in full general relativity, we find that a magnetized HMNS undergoes `delayed' collapse to a rotating black hole (BH) as a result of angular momentum transport via magnetic braking and the magnetorotational instability. The outcome is a BH surrounded by a massive, hot torus with a collimated magnetic field. The torus accretes onto the BH at a quasi-steady accretion rate ~10 solar mass/s; the lifetime of the torus is ~10 ms. The torus has a temperature \sim 10^{12} K, leading to copious neutrino-antineutrino thermal radiation. Therefore, the collapse of an HMNS is a promising scenario for generating short-duration gamma-ray bursts and an accompanying burst of gravitational waves and neutrinos.

A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving {\it dynamical} spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accurate as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC scheme has increased computational efficiency as it avoids the costly solution of Riemann problems and has practical advantages in the modeling of neutron star spacetimes. Namely, it allows simulations with stiff equations of state by successfully dealing with very low-density unphysical atmospheres. These facts not only suggest that such a HRC scheme may be a desirable tool for hydrodynamical simulations in general relativity, but also open the possibility to perform accurate magnetohydrodynamical simulations in curved dynamic spacetimes.

A new implementation for magnetohydrodynamics (MHD) simulations in full general relativity (involving dynamical spacetimes) is presented. In our implementation, Einstein's evolution equations are evolved by a BSSN formalism, MHD equations by a high-resolution central scheme, and induction equation by a constraint transport method. We perform numerical simulations for standard test problems in relativistic MHD, including special relativistic magnetized shocks, general relativistic magnetized Bondi flow in stationary spacetime, and a longterm evolution for self-gravitating system composed of a neutron star and a magnetized disk in full general relativity. In the final test, we illustrate that our implementation can follow winding-up of the magnetic field lines of magnetized and differentially rotating accretion disks around a compact object until saturation, after which magnetically driven wind and angular momentum transport inside the disk turn on.

We perform axisymmetric simulations for gravitational collapse of a massive iron core to a black hole in full general relativity. The iron cores are modeled by $\Gamma = 4/3$ equilibrium polytrope for simplicity and a parametric equation of state is adopted during the collapse. Performing simulations for a wide variety of initial conditions changing the mass, the angular momentum, the rotational velocity profile of the core, and the parameters of the equations of state (which are chosen so that the maximum mass of the cold spherical polytrope is $\approx 1.6M_{\odot}$), a criterion for the prompt black hole formation is clarified in terms of them. It is found that (i) with the increase of the thermal energy generated by shocks, the threshold mass for the prompt black hole formation is increased by 20--40%, (ii) the rotational centrifugal force increases the threshold mass by $\alt 25%$, (iii) with the increase of the degree of differential rotation, the threshold mass is also increased, and (iv) the amplification factors shown in the results (i)--(iii) depend sensitively on the equation of state. We also find that the collapse dynamics and the structure of the shock formed at the bounce depend strongly on the stiffness of the adopted equation of state. In particular, as a new feature, a strong bipolar explosion is observed for the collapse of rapidly rotating iron cores with an equation of state which is stiff in subnuclear density and soft in supranuclear density. Gravitational waves are computed in terms of a quadrupole formula. It is also found that the waveform depends sensitively on the equations of state.

Latest general relativistic simulations for merger of binary neutron stars with realistic equations of state (EOSs) show that a hypermassive neutron star of an ellipsoidal figure is formed after the merger if the total mass is smaller than a threshold value which depends on the EOSs. The effective amplitude of quasiperiodic gravitational waves from such hypermassive neutron stars is $\sim 6$--$7 \times 10^{-21}$ at a distance of 50 Mpc, which may be large enough for detection by advanced laser interferometric gravitational wave detectors although the frequency is high $\sim 3$ kHz. We point out that the detection of such signal may lead to constraining the EOSs for neutron stars.

We present numerical results of three-dimensional simulations for the merger of binary neutron stars (BNSs) in full general relativity. Hybrid equations of state (EOSs) are adopted to mimic realistic nuclear EOSs. In this approach, we divide the EOSs into two parts, i.e., the thermal part and the cold part. For the cold part, we assign a fitting formula for realistic EOSs of cold nuclear matter slightly modifying the formula developed by Haensel and Potekhin. We adopt the SLy and FPS EOSs for which the maximum allowed ADM mass of cold and spherical neutron stars (NSs) is ~ 2.04Mo and 1.80Mo, respectively. Simulations are performed for BNSs of the total ADM mass in the range between 2.4Mo and 2.8Mo with the rest-mass ratio Q_M to be in the range 0.9 < Q_M < 1. It is found that if the total ADM mass of the system is larger than a threshold M_{thr}, a black hole (BH) is promptly formed in the merger irrespective of the mass ratios. In the other case, the outcome is a hypermassive NS of a large ellipticity, which results from the large adiabatic index of the realistic EOSs adopted. The value of M_{thr} depends on the EOS: ~ 2.7Mo and ~ 2.5Mo for the SLy and FPS EOSs, respectively. Gravitational waves are computed in terms of a gauge-invariant wave extraction technique. In the formation of the hypermassive NS, quasiperiodic gravitational waves of a large amplitude and of frequency between 3 and 4 kHz are emitted. The estimated emission time scale is < 100 ms, after which the hypermassive NS collapses to a BH. Because of the long emission time, the effective amplitude may be large enough to be detected by advanced laser interferometric gravitational wave detectors if the distance to the source is smaller than ~ 100 Mpc.

We derive a tidal potential for a self-gravitating fluid star orbiting Kerr black hole along a timelike geodesic extending previous works by Fishbone and Marck. In this paper, the tidal potential is calculated up to the third and fourth-order terms in $R/r$, where $R$ is the stellar radius and $r$ the orbital separation, in the Fermi-normal coordinate system following the framework developed by Manasse and Misner. The new formulation is applied for determining the tidal disruption limit (Roche limit) of corotating Newtonian stars in circular orbits moving on the equatorial plane of Kerr black holes. It is demonstrated that the third and fourth-order terms quantitatively play an important role in the Roche limit for close orbits with $R/r \agt 0.1$. It is also indicated that the Roche limit of neutron stars orbiting a stellar-mass black hole near the innermost stable circular orbit may depend sensitively on the equation of state of the neutron star.

We perform fully general relativistic simulations of rotating stellar core collapse in three spatial dimension. A parametric equation of state is adopted following Dimmelmeier et al. The early stage of the collapse is followed by an axisymmetric code. When the stellar core becomes compact enough, we start a 3-dimensional simulation adding a bar-mode nonaxisymmetric density perturbation. In the axisymmetric simulations, it is clarified that the maximum value of $\beta \equiv T/W$ achieved during the stellar collapse and depends sensitively on the velocity profile and total mass of the initial core, and equations of state. It is also found that for all the models with high degree of differential rotation, a funnel structure is formed around the rotational axis after the formation of neutron stars. For selected models in which the maximum value of $\beta$ is larger than $\sim 0.27$, 3-dimensional simulations are performed. It is found that the bar-mode dynamical instability sets in for the case that the following conditions are satisfied: (i) the progenitor should be rapidly rotating with the initial value of $0.01 \alt \beta \alt 0.02$, (ii) the degree of differential rotation of the initial condition should be sufficiently high, and (iii) a depletion factor of pressure in an early stage of collapse should be large enough to induce a significant contraction for which an efficient spin-up can be achieved. As a result of the onset of the bar-mode dynamical instabilities, the amplitude of gravitational waves can be by a factor of $\sim 10$ larger than that in the axisymmetric collapse. It is found that a dynamical instability with the $m=1$ mode is also induced, but the perturbation does not grow significantly and, hence, it does not contribute to an outstanding amplification of gravitational waves.

The evolution of a nonaxisymmetric bar-mode perturbation of rapidly rotating stars due to a secular instability induced by gravitational wave emission is studied in post-Newtonian simulations taking into account gravitational radiation reaction. A polytropic equation of state with the polytropic index $n=1$ is adopted. The ratio of the rotational kinetic energy to the gravitational potential energy $T/|W|$ is chosen in the range between 0.2 and 0.26. Numerical simulations were performed until the perturbation grows to the nonlinear regime, and illustrate that the outcome after the secular instability sets in is an ellipsoidal star of a moderately large ellipticity $\agt 0.7$. A rapidly rotating protoneutron star may form such an ellipsoid, which is a candidate for strong emitter of gravitational waves for ground-based laser interferometric detectors. A possibility that effects of magnetic fields neglected in this work may modify the growth of the secular instability is also mentioned.

We study gravitational collapse of rapidly rotating relativistic polytropes of the adiabatic index $\Gamma = 1.5$ and 2, in which the spin parameter $q \equiv J/M^{2} > 1$ where $J$ and $M$ are total angular momentum and gravitational mass, in full general relativity. First, analyzing initial distributions of the mass and the spin parameter inside stars, we predict the final outcome after the collapse. Then, we perform fully general relativistic simulations on assumption of axial and equatorial symmetries and confirm our predictions. As a result of simulations, we find that in contrast with the previous belief, even for stars with $q > 1$, the collapse proceeds to form a seed black hole at central region, and the seed black hole subsequently grows as the ambient fluids accrete onto it. We also find that growth of angular momentum and mass of the seed black hole can be approximately determined from the initial profiles of the density and the specific angular momentum. We define an effective spin parameter at the central region of the stars, $q_{c}$, and propose a new criterion for black hole formation as $q_{c} \alt 1$. Plausible reasons for the discrepancy between our and previous results are clarified.

We study secular stability against a quasi-radial oscillation for rigidly rotating stars with soft equations of state in general relativity. The polytropic equations of state with polytropic index $n$ between 3 and 3.05 are adopted for modeling the rotating stars. The stability is determined in terms of the turning-point method. It is found that (i) for $n \agt 3.04$, all the rigidly rotating stars are unstable against the quasi-radial oscillation and (ii) for $n \agt 3.01$, the nondimensional angular momentum parameter $q \equiv cJ/GM^2$ (where $J$, $M$, $G$, and $c$ denote the angular momentum, the gravitational mass, the gravitational constant, and the speed of light, respectively) for all marginally stable rotating stars is larger than unity. A semi-analytic calculation is also performed, and good agreement with the numerical results is confirmed. The final outcome after axisymmetric gravitational collapse of rigidly rotating and marginally stable massive stars with $q > 1$ is predicted, assuming that the rest-mass distribution as a function of the specific angular momentum is preserved and that the pressure never halt the collapse. It is found that even for $1 < q \alt 2.5$, a black hole may be formed as a result of the collapse, but for $q \agt 2.5$, the significant angular momentum will prevent the direct formation of a black hole.

Gravitational waves from oscillating neutron stars in axial symmetry are studied performing numerical simulations in full general relativity. Neutron stars are modeled by a polytropic equation of state for simplicity. A gauge-invariant wave extraction method as well as a quadrupole formula are adopted for computation of gravitational waves. It is found that the gauge-invariant variables systematically contain numerical errors generated near the outer boundaries in the present axisymmetric computation. We clarify their origin, and illustrate it possible to eliminate the dominant part of the systematic errors. The best corrected waveforms for oscillating and rotating stars currently contain errors of magnitude $\sim 10^{-3}$ in the local wave zone. Comparing the waveforms obtained by the gauge-invariant technique with those by the quadrupole formula, it is shown that the quadrupole formula yields approximate gravitational waveforms besides a systematic underestimation of the amplitude of $O(M/R)$ where $M$ and $R$ denote the mass and the radius of neutron stars. However, the wave phase and modulation of the amplitude can be computed accurately. This indicates that the quadrupole formula is a useful tool for studying gravitational waves from rotating stellar core collapse to a neutron star in fully general relativistic simulations. Properties of the gravitational waveforms from the oscillating and rigidly rotating neutron stars are also addressed paying attention to the oscillation associated with fundamental modes.

Axisymmetric numerical simulations of rotating stellar core collapse to a neutron star are performed in the framework of full general relativity. The so-called Cartoon method, in which the Einstein field equations are solved in the Cartesian coordinates and the axisymmetric condition is imposed around the $y=0$ plane, is adopted. The hydrodynamic equations are solved in the cylindrical coordinates (on the $y=0$ plane in the Cartesian coordinates) using a high-resolution shock-capturing scheme with the maximum grid size $(2500,2500)$. A parametric equation of state is adopted to model collapsing stellar cores and neutron stars following Dimmelmeier et al. It is found that the evolution of central density during the collapse, bounce, and formation of protoneutron stars agree well with those in the work of Dimmelmeier et al. in which an approximate general relativistic formulation is adopted. This indicates that such approximation is appropriate for following axisymmetric stellar core collapses and subsequent formation of protoneutron stars. Gravitational waves are computed using a quadrupole formula. It is found that the waveforms are qualitatively in good agreement with those by Dimmelmeier et al. However, quantitatively, two waveforms do not agree well. Possible reasons for the disagreement are discussed.

We present results of three dimensional numerical simulations of the merger of unequal-mass binary neutron stars in full general relativity. A $\Gamma$-law equation of state $P=(\Gamma-1)\rho\epsilon$ is adopted, where $P$, $\rho$, $\varep$, and $\Gamma$ are the pressure, rest mass density, specific internal energy, and the adiabatic constant, respectively. We take $\Gamma=2$ and the baryon rest-mass ratio $Q_M$ to be in the range 0.85--1. The typical grid size is $(633,633,317)$ for $(x,y,z)$ . We improve several implementations since the latest work. In the present code, the radiation reaction of gravitational waves is taken into account with a good accuracy. This fact enables us to follow the coalescence all the way from the late inspiral phase through the merger phase for which the transition is triggered by the radiation reaction. It is found that if the total rest-mass of the system is more than $\sim 1.7$ times of the maximum allowed rest-mass of spherical neutron stars, a black hole is formed after the merger irrespective of the mass ratios. The gravitational waveforms and outcomes in the merger of unequal-mass binaries are compared with those in equal-mass binaries. It is found that the disk mass around the so formed black holes increases with decreasing rest-mass ratios and decreases with increasing compactness of neutron stars. The merger process and the gravitational waveforms also depend strongly on the rest-mass ratios even for the range $Q_M= 0.85$--1.

We study the final state of the gravitational collapse of uniformly rotating supramassive neutron stars by axisymmetric simulations in full general relativity. The rotating stars provided as the initial condition are marginally stable against quasiradial gravitational collapse and its equatorial radius rotates with the Kepler velocity (i.e., the star is at the mass-shedding limit). To model the neutron stars, we adopt the polytropic equations of state for a wide range of the polytropic index as $n=2/3$, 4/5, 1, 3/2 and 2. We follow the formation and evolution of the black holes, and show that irrespective of the value of $n (2/3\leq n \leq 2)$, the final state is a Kerr black hole and the disk mass is very small ($< 10^{-3}$ of the initial stellar mass).

Nonspherical stellar collapse to a black hole is one of the most promising gravitational wave sources for gravitational wave detectors. We numerically study gravitational waves from a slightly nonspherical stellar collapse to a black hole in linearized Einstein theory. We adopt a spherically collapsing star as the zeroth-order solution and gravitational waves are computed using perturbation theory on the spherical background. In this paper we focus on the perturbation of odd-parity modes. Using the polytropic equations of state with polytropic indices $n_p=1$ and 3, we qualitatively study gravitational waves emitted during the collapse of neutron stars and supermassive stars to black holes from a marginally stable equilibrium configuration. Since the matter perturbation profiles can be chosen arbitrarily, we provide a few types for them. For $n_p=1$, the gravitational waveforms are mainly characterized by a black hole quasinormal mode ringing, irrespective of perturbation profiles given initially. However, for $n_p=3$, the waveforms depend strongly on the initial perturbation profiles. In other words, the gravitational waveforms strongly depend on the stellar configuration and, in turn, on the ad hoc choice of the functional form of the perturbation in the case of supermassive stars.

As an extension of our previous work, we investigate the dynamical instability against nonaxisymmetric bar-mode deformations of differentially rotating stars in Newtonian gravity varying the equations of state and velocity profiles. We performed the numerical simulation and the followup linear stability analysis adopting polytropic equations of state with the polytropic indices n=1, 3/2, and 5/2 and with two types of angular velocity profiles (the so-called j-constant-like and Kepler-like laws). It is confirmed that rotating stars of a high degree of differential rotation are dynamically unstable against the bar-mode deformation, even for the ratio of the kinetic energy to the gravitational potential energy $\beta$ of order 0.01. The criterion for onset of the bar-mode dynamical instability depends weakly on the polytropic index n and the angular velocity profile as long as the degree of differential rotation is high. Gravitational waves from the final nonaxisymmetric quasi-stationary states are calculated in the quadrupole formula. For proto-neutron stars of mass $1.4M_{\odot}$, radius $\sim 30$ km and $\beta \alt 0.1$, such gravitational waves have the frequency of $\sim$ 600--1,400 Hz, and the effective amplitude is larger than $10^{-22}$ at a distance of about 100 Mpc irrespective of n and the angular velocity profile.

We report a new implementation for axisymmetric simulation in full general relativity. In this implementation, the Einstein equations are solved using the Nakamura-Shibata formulation with the so-called cartoon method to impose an axisymmetric boundary condition, and the general relativistic hydrodynamic equations are solved using a high-resolution shock-capturing scheme based on an approximate Riemann solver. As tests, we performed the following simulations: (i) long-term evolution of non-rotating and rapidly rotating neutron stars, (ii) long-term evolution of neutron stars of a high-amplitude damping oscillation accompanied with shock formation, (iii) collapse of unstable neutron stars to black holes, and (iv) stellar collapses to neutron stars. The tests (i)--(iii) were carried out with the $\Gamma$-law equation of state, and the test (iv) with a more realistic parametric equation of state for high-density matter. We found that this new implementation works very well: It is possible to perform the simulations for stable neutron stars for more than 10 dynamical time scales, to capture strong shocks formed at stellar core collapses, and to accurately compute the mass of black holes formed after the collapse and subsequent accretion. In conclusion, this implementation is robust enough to apply to astrophysical problems such as stellar core collapse of massive stars to a neutron star and black hole, phase transition of a neutron star to a high-density star, and accretion-induced collapse of a neutron star to a black hole. The result for the first simulation of stellar core collapse to a neutron star started from a realistic initial condition is also presented.

The collapse of a uniformaly rotating, supermassive star (SMS) to a supermassive black hole (SMBH) has been followed recently by means of hydrodynamic simulations in full general relativity. The initial SMS of arbitrary mass M in these simulations rotates uniformly at the mass--shedding limit and is marginally unstable to radial collapse. The final black hole has mass M_h/M = 0.9 and and spin J_h/M_h^2 = 0.75, approximately. The remaining mass goes into a disk of mass M_disk/M = 0.1, also approximately. Here we show that these black hole and disk parameters can be calculated analytically from the initial stellar density and angular momentum distribution. The analytic calculation thereby corroborates and provides a simple physical explanation for the computational discovery that SMS collapse inevitably terminates in the simultaneous formation of a SMBH and a rather substantial ambient disk. This disk arises even though the total spin of the progenitor star, J/M^2 = 0.97, is safely below the Kerr limit. The calculation performed here applies to any marginally unstable n = 3 polytrope uniformly rotating at the break--up speed, independent of stellar mass or the source of internal pressure. It illustrates how the black hole and disk parameters can be determined for the collapse of other types of stars with different initial density and rotation profiles.

We study the dynamical instability against bar-mode deformation of differentially rotating stars. We performed numerical simulation and linear perturbation analysis adopting polytropic equations of state with the polytropic index $n=1$. It is found that rotating stars of a high degree of differential rotation are dynamically unstable even for the ratio of the kinetic energy to the gravitational potential energy of $O(0.01)$. Gravitational waves from the final nonaxisymmetric quasistationary states are calculated in the quadrupole formula. For rotating stars of mass $1.4M_{\odot}$ and radius several 10 km, gravitational waves have frequency several 100 Hz and effective amplitude $\sim 5 \times 10^{-22}$ at a distance of $\sim 100$ Mpc.

We follow the collapse in axisymmetry of a uniformly rotating, supermassive star (SMS) to a supermassive black hole (SMBH) in full general relativity. The initial SMS of arbitrary mass $M$ is marginally unstable to radial collapse and rotates at the mass-shedding limit. The collapse proceeds homologously early on and results in the appearance of an apparent horizon at the center. Although our integration terminates before final equilibrium is achieved, we determine that the final black hole will contain about 90% of the total mass of the system and have a spin parameter $J/M^2 \sim 0.75$. The remaining gas forms a rotating disk about the nascent hole.

We study the gravitational collapse of a rotating supermassive star (SMS) by means of a (3+1) hydrodynamical simulation in a post-Newtonian (PN) approxi- mation of general relativity. This problem is particularly challenging because of the vast dynamical range in space which must be covered in the course of col- lapse. We evolve a uniformly rotating SMS from the onset of radial instability at R_{p}/M=411, where R_{p} is the proper polar radius of the star and M is the total mass-energy, to the point at which the PN approximation breaks down. We introduce a scale factor and a "comoving" coordinate to handle the large varia- tion in radius during the collapse and focus on the central core. Since T/W, the ratio of the rotational kinetic energy to the gravitational binding energy, is nearly proportional to 1/R_{p} throughout the collapse, the imploding star may ultimately exceed the critical value of T/W for dynamical instability to bar-mode formation. However, for stars rotating uniformly at the onset of col- lapse, we do not find any unstable growth of bars prior to the termination of our simulation. We do find that the collapse is likely to form a supermassive black hole (BH) coherently, with almost all of the matter falling into the hole, leaving very little ejected matter to form a disk. In the absence of nonaxi- symmetric bar formation, the collapse of a uniformly rotating SMS does not lead to appreciable quasi-periodic gravitational wave (GW) emission by the time our integrations terminate. However, the coherent nature of the implosion suggests that rotating SMS collapse will be a promising source of GW bursts. We also expect that, following BH formation, long wavelength quasi-periodic waves will result from quasi-normal ringing. These waves may be detectable by LISA.

We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^\alpha, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes. Binary systems with M\Omega small have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation \delta M = \Omega \delta J, from which one can obtain a turning point criterion that governs the stability of orbits.

We model the adiabatic inspiral of relativistic binary neutron stars in a quasi-equilibrium (QE) approximation, and compute the gravitational wavetrain from the late phase of the inspiral. We compare corotational and irrotational sequences and find a significant difference in the inspiral rate, which is almost entirely caused by differences in the binding energy. We also compare our results with those of a point-mass post-Newtonian calculation. We illustrate how the late inspiral wavetrain computed with our QE numerical scheme can be matched to the subsequent plunge and merger waveform calculated with a fully relativistic hydrodynamics code.

We perform 3D numerical simulations for merger of equal mass binary neutron stars in full general relativity preparing irrotational binary neutron stars in a quasiequilibrium state as initial conditions. Simulations have been carried out for a wide range of stiffness of equations of state and compactness of neutron stars, paying particular attention to the final products and gravitational waves. We take a fixed uniform grid in Cartesian coordinates with typical grid size $(293,293,147)$ in $(x,y,z)$ assuming a plane symmetry with respect to the equatorial plane. A result of one new large-scale simulation performed with grid size $(505,505,253)$ is also presented. We find that the final product depends sensitively on the initial compactness of the neutron stars: In a merger between sufficiently compact neutron stars, a black hole is formed in a dynamical timescale. As the compactness is decreased, the formation timescale becomes longer and longer. For less compact cases, a differentially rotating massive neutron star is formed instead of a black hole. In the case of black hole formation, the disk mass around the black hole appears to be smaller than 1% of the total rest mass. It is also indicated that waveforms of high-frequency gravitational waves after merger depend strongly on the compactness of neutron stars.

We investigate the equilibrium and stability of supermassive stars of mass $M \agt 10^5M_{\odot}$ in binary systems. We find that corotating binaries are secularly unstable for close, circular orbits with $r \alt 4R(M/10^6M_{\odot})^{1/6}$ where $r$ is the orbital separation and $R$ the stellar radius. We also show that corotation cannot be achieved for distant orbits with $r \agt 12 R (M/10^6M_{\odot})^{-11/24}$, since the timescale for viscous angular momentum transfer associated with tidal torques is longer than the evolution timescale due to emission of thermal radiation. These facts suggest that the allowed mass range and orbital separation for corotating supermassive binary stars is severely restricted. In particular, for supermassive binary stars of large mass $M \agt 6\times 10^6M_{\odot}$, corotation cannot be achieved, as viscosity is not adequate to mediate the transfer between orbital and spin angular momentum. One possible outcome for binary supermassive stars is the onset of quasi-radial, relativistic instability which drives each star to collapse prior to merger: We discuss alternative outcomes of collapse and possible spin states of the resulting black holes. We estimate the frequency and amplitude of gravitational waves emitted during several inspiral and collapse scenarios.

We study by computational means the dynamical stability against bar-mode deformation of rapidly and differentially rotating stars in a post-Newtonian approximation of general relativity. We vary the compaction of the star $M/R$ (where $M$ is the gravitational mass and $R$ the equatorial circumferential radius) between 0.01 and 0.05 to isolate the influence of relativistic gravitation on the instability. For compactions in this moderate range, the critical value of $\beta = T/W$ for the onset of the dynamical instability (where $T$ is the rotational kinetic energy and W the gravitational binding energy) slightly decreases from $\sim$ 0.26 to $\sim$ 0.25 with increasing compaction for our choice of the differential rotational law. Combined with our earlier findings based on simulations in full general relativity for stars with higher compaction, we conclude that relativistic gravitation enhances the dynamical bar-mode instability, i.e. the onset of instability sets in for smaller values of $\beta$ in relativistic gravity than in Newtonian gravity. We also find that once a triaxial structure forms after the bar-mode perturbation saturates in dynamically unstable stars, the triaxial shape is maintained, at least for several rotational periods. To check the reliability of our numerical integrations, we verify that the general relativistic Kelvin-Helmholtz circulation is well-conserved, in addition to rest-mass and total mass-energy, linear and angular momentum. Conservation of circulation indicates that our code is not seriously affected by numerical viscosity. We determine the amplitude and frequency of the quasi-periodic gravitational waves emitted during the bar formation process using the quadrupole formula.

We study the dynamical stability against bar-mode deformation of rapidly and differentially rotating stars in the first post-Newtonian approximation of general relativity. We vary the compaction of the star $M/R$ (where $M$ is the gravitational mass and $R$ the equatorial circumferential radius) between 0.01 and 0.05 to isolate the influence of relativistic gravitation on the instability. For compactions in this moderate range, the critical value of $\beta \equiv T/W$ for the onset of the dynamical instability (where $T$ is the rotational kinetic energy and $W$ the gravitational binding energy) slightly decreases from $\sim 0.26$ to $\sim 0.25$ with increasing compaction for our choice of the differential rotational law. Combined with our earlier findings based on simulations in full general relativity for stars with higher compaction, we conclude that relativistic gravitation {\em enhances} the dynamical bar-mode instability, i.e. the onset of instability sets in for smaller values of $\beta$ in relativistic gravity than in Newtonian gravity. We also find that once a triaxial structure forms after the bar-mode perturbation saturates in dynamically unstable stars, the triaxial shape is maintained, at least for several rotational periods. To check the reliability of our numerical integrations, we verify that the general relativistic Kelvin-Helmholtz circulation is well-conserved, in addition to rest-mass energy, total mass-energy, linear and angular momentum. Conservation of circulation indicates that our code is not seriously affected by numerical viscosity. We determine the amplitude and frequency of the quasi-periodic gravitational waves emitted during the bar formation process using the quadrupole formula.

We numerically construct time-symmetric initial data sets of a black hole in the Randall-Sundrum brane world model, assuming that the black hole is spherical on the brane. We find that the apparent horizon is cigar-shaped in the 5D spacetime.

We investigate close binary neutron stars in quasiequilibrium states in a general relativistic framework. We assume conformal flatness for the spatial metric and irrotational velocity field for the neutron stars. We adopt the polytropic equation of state. The computation is performed for the polytropic index n(=0.5, 0.66667, 0.8, 1, 1.25), and compactness of neutron stars M/R(=0.03 - 0.3). Results of this paper are as follows. (i) The sequences of the irrotational binary are always terminated at an innermost orbit where a cusp (inner Lagrange point)appears at the inner edges of the stellar surface. The binaries with cusps are found to be dynamically unstable for n=0.5 and stable for n > 0.8 irrespective of M/R < 0.2. For n=0.66667, the stability changes depending on M/R. (ii) The gravitational wave frequency at the innermost orbit turns out to be between 800 and 1500 Hz for realistic compactness 0.14 < M/R < 0.2. (iii) The ISCO for n=0.5 appears to be determined by a hydrodynamic instability for M/R < 0.2. We derive fitting formulae for the relation between the orbital angular velocity at the ISCO and the compactness to clarify it. (iv) The maximum density of neutron stars in binary systems slightly decreases with decreasing the orbital separation and hence they are stable against individual radial collapse during the inspiral. (v) q = J_tot/M_ADM^2 at the innermost orbits is always less than unity for M/R > 0.13 irrespective of n, which indicates that the realistic binary neutron stars satisfy a necessary condition (q<1) for formation of a black hole before the merger. (vi) The specific angular momentum of any mass element in irrotational binary neutron stars at the innermost orbit appears to be too small to form a disk around black holes formed after the merger.

We study the dynamical stability against bar-mode deformation of rapidly spinning neutron stars with differential rotation. We perform fully relativistic 3D simulations of compact stars with $M/R \geq 0.1$, where $M$ is the total gravitational mass and $R$ the equatorial circumferential radius. We adopt an adiabatic equation of state with adiabatic index $\Gamma=2$. As in Newtonian theory, we find that stars above a critical value of $\beta \equiv T/W$ (where $T$ is the rotational kinetic energy and $W$ the gravitational binding energy) are dynamically unstable to bar formation. For our adopted choices of stellar compaction and rotation profile, the critical value of $\beta = \beta_{dGR}$ is $\sim 0.24-0.25$, only slightly smaller than the well-known Newtonian value $\sim 0.27$ for incompressible Maclaurin spheroids. The critical value depends only very weakly on the degree of differential rotation for the moderate range we surveyed. All unstable stars form bars on a dynamical timescale. Models with sufficiently large $\beta$ subsequently form spiral arms and eject mass, driving the remnant to a dynamically stable state. Models with moderately large $\beta \gtrsim \beta_{dGR}$ do not develop spiral arms or eject mass but adjust to form dynamically stable ellipsoidal-like configurations. If the bar-mode instability is triggered in supernovae collapse or binary neutron star mergers, it could be a strong and observable source of gravitational waves. We determine characteristic wave amplitudes and frequencies.

We have performed 3D numerical simulations for merger of equal mass binary neutron stars in full general relativity. We adopt a $\Gamma$-law equation of state in the form $P=(\Gamma-1)\rho\epsilon$ where P, $\rho$, $\varep$ and $\Gamma$ are the pressure, rest mass density, specific internal energy, and the adiabatic constant with $\Gamma=2$. As initial conditions, we adopt models of corotational and irrotational binary neutron stars in a quasi-equilibrium state which are obtained using the conformal flatness approximation for the three geometry as well as an assumption that a helicoidal Killing vector exists. In this paper, we pay particular attention to the final product of the coalescence. We find that the final product depends sensitively on the initial compactness parameter of the neutron stars : In a merger between sufficiently compact neutron stars, a black hole is formed in a dynamical timescale. As the compactness is decreased, the formation timescale becomes longer and longer. It is also found that a differentially rotating massive neutron star is formed instead of a black hole for less compact binary cases, in which the rest mass of each star is less than 70-80% of the maximum allowed mass of a spherical star. In the case of black hole formation, we roughly evaluate the mass of the disk around the black hole. For the merger of corotational binaries, a disk of mass $\sim 0.05-0.1M_*$ may be formed, where M_* is the total rest mass of the system. On the other hand, for the merger of irrotational binaries, the disk mass appears to be very small : < 0.01M_*.

We perform 3D numerical simulations in full general relativity to study the stability of rapidly rotating, supramassive neutron stars at the mass-shedding limit to dynamical collapse. We adopt an adiabatic equation of state with $\Gamma = 2$ and focus on uniformly rotating stars. We find that the onset of dynamical instability along mass-shedding sequences nearly coincides with the onset of secular instability. Unstable stars collapse to rotating black holes within about one rotation period. We also study the collapse of stable stars which have been destabilized by pressure depletion (e.g. via a phase transition) or mass accretion. In no case do we find evidence for the formation of massive disks or any ejecta around the newly formed Kerr black holes, even though the progenitors are rapidly rotating.

We construct relativistic equilibrium models of differentially rotating neutron stars and show that they can support significantly more mass than their nonrotating or uniformly rotating counterparts. We dynamically evolve such ``hypermassive'' models in full general relativity and show that there do exist configurations which are dynamically stable against radial collapse and bar formation. Our results suggest that the remnant of binary neutron star coalescence may be temporarily stabilized by differential rotation, leading to delayed collapse and a delayed gravitational wave burst.

We present our first successful numerical results of 3D general relativistic simulations in which the Einstein equation as well as the hydrodynamic equations are fully solved. This paper is especially devoted to simulations of test problems such as spherical dust collapse, stability test of perturbed spherical stars, and preservation of (approximate) equilibrium states of rapidly rotating neutron star and/or corotating binary neutron stars. These test simulations confirm that simulations of coalescing binary neutron stars are feasible in a numerical relativity code. It is illustrated that using our numerical code, simulations of these problems, in particular those of corotating binary neutron stars, can be performed stably and fairly accurately for a couple of dynamical timescales. These numerical results indicate that our formulation for solving the Einstein field equation and hydrodynamic equations are robust and make it possible to perform a realistic simulation of coalescing binary neutron stars for a long time from the innermost circular orbit up to formation of a black hole or neutron star.

We report on our numerical implementation of fully relativistic hydrodynamics coupled to Einstein's field equations in three spatial dimensions. We briefly review several steps in our code development, including our recasting of Einstein's equations and several tests which demonstrate its advantages for numerical integrations. We outline our implementation of relativistic hydrodynamics, and present numerical results for the evolution of both stable and unstable Oppenheimer-Volkov equilibrium stars, which represent a very promising first test of our code.

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.

We investigate the innermost stable circular orbit (ISCO) of a test particle moving on the equatorial plane around rotating relativistic stars such as neutron stars. First, we derive approximate analytic formulas for the angular velocity and circumferential radius at the ISCO making use of an approximate relativistic solution which is characterized by arbitrary mass, spin, mass quadrupole, current octapole and mass $2^4$-pole moments. Then, we show that the analytic formulas are accurate enough by comparing them with numerical results, which are obtained by analyzing the vacuum exterior around numerically computed geometries for rotating stars of polytropic equation of state. We demonstrate that contribution of mass quadrupole moment for determining the angular velocity and, in particular, the circumferential radius at the ISCO around a rapidly rotating star is as important as that of spin.

We perform simulations of relativistic binary stars in post-Newtonian gravity to investigate their dynamical stability prior to merger against gravitational collapse in a tidal field. In general, our equations are only strictly accurate to first post-Newtonian order, but they recover full general relativity for spherical, static stars. We study both corotational and irrotational binary configurations of identical stars in circular orbits. We adopt a soft, adiabatic equation of state with $\Gamma = 1.4$, for which the onset of instability occurs at a sufficiently small value of the compaction $M/R$ that a post-Newtonian approximation is quite accurate. For such a soft equation of state there is no innermost stable circular orbit, so that we can study arbitrarily close binaries. This choice still allows us to study all the qualitative features exhibited by any adiabatic equation of state regarding stability against gravitational collapse. We demonstrate that, independent of the internal stellar velocity profile, the tidal field from a binary companion stabilizes a star against gravitational collapse.

We present a formalism to obtain equilibrium configurations of uniformly rotating fluid in the second post-Newtonian approximation of general relativity. In our formalism, we need to solve 29 Poisson equations, but their source terms decrease rapidly enough at the external region of the matter(i.e., at worst $O(r^{-4})$). Hence these Poisson equations can be solved accurately as the boundary value problem using standard numerical methods.This formalism will be useful to obtain nonaxisymmetric uniformly rotating equilibrium configurations such as synchronized binary neutron stars just before merging and the Jacobi ellipsoid.

Using the (3+1) formalism in general relativity, we perform the post-Newtonian(PN) approximation to clarify what sort of gauge condition is suitable for numerical analysis of coalescing compact binary neutron stars and gravitational waves from them. We adopt a kind of transverse gauge condition to determine the shift vector. On the other hand, for determination of the time slice, we adopt three slice conditions(conformal slice, maximal slice and harmonic slice) and discuss their properties. Using these conditions, the PN hydrodynamic equations are obtained up through the 2.5PN order including the quadrupole gravitational radiation reaction. In particular, we describe methods to solve the 2PN tensor potential which arises from the spatial 3-metric. It is found that the conformal slice seems appropriate for analysis of gravitational waves in the wave zone and the maximal slice will be useful for describing the equilibrium configurations. The PN approximation in the (3+1) formalism will be also useful to perform numerical simulations using various slice conditions and, as a result, to provide an initial data for the final merging phase of coalescing binary neutron stars which can be treated only by fully general relativistic simulations.

We investigate the distance-redshift relation in the simple void model. As discussed by Moffat and Tatarski, if the observer stays at the center of the void, the observed Hubble parameter is not so different from the background Hubble parameter. However, if the position of observer is off center of the void, we must consider the peculiar velocity correction which is measured by the observed dipole anisotropy of cosmic microwave background. This peculiar velocity correction for the redshift is crucial to determine the Hubble parameter and we shall discuss this effect. Further the results of Turner et al by the N-body simulation will be also considered.

Using the Teukolsky and Sasaki-Nakamura equations for the gravitational perturbation of the Kerr spacetime, we calculate the post-Newtonian expansion of the energy and angular momentum luminosities of gravitational waves from a test particle orbiting around a rotating black hole up through ${\rm P^{5/2}N}$ order beyond the quadrupole formula. We apply a method recently developed by Sasaki to the case of a rotating black hole. We take into account a small inclination of the orbital plane to the lowest order of the Carter constant. The result to $ P^{3/2}N}$ order is in agreement with a similar calculation by Poisson as well as with the standard post-Newtonian calculation by Kidder et al. Using our result, we calculate the integrated phase of gravitational waves from a neutron star-neutron star binary and a black hole-neutron star binary during their inspiral stage. We find that, in both cases, spin-dependent terms in the P$^2$N and P$^{5/2}$N corrections are important to construct effective template waveforms which will be used for future laser-interferometric gravitational wave detectors.