We have analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) cycle 2 data of band 6 toward the G345.4938+01.4677 massive young protostellar object (G345.5+1.47 MYSO) in the IRAS 16562--3959 high-mass star-forming region with an angular resolution of $\sim 0.3"$, corresponding to $\sim 760$ au. We spatially resolve the central region which consists of three prominent molecular emission cores. A hypercompact (HC) H$_{\rm {II}}$ region (Core A) and two molecule-rich cores (Core B and Core C) are identified using the moment zero images of the H30$\alpha$ line and a CH$_{3}$OH line, respectively. Various oxygen-bearing complex organic molecules (COMs), such as (CH$_{3}$)$_{2}$CO and CH$_{3}$OCHO, have been detected toward the positions of Core B and Core C, while nitrogen-bearing species, CH$_{3}$CN, HC$_{3}$N and its $^{13}$C isotopologues, have been detected toward all of the cores. We discuss the formation mechanisms of H$_{2}$CO by comparing the spatial distribution of C$^{18}$O with that of H$_{2}$CO. The $^{33}$SO emission, on the other hand, shows a ring-like structure surrounding Core A, and it peaks on the outer edge of the H30$\alpha$ emission region. These results imply that SO is enhanced in a shock produced by the expanding motion of the ionized region.

We have carried out mapping observations of molecular emission lines of HC$_{3}$N and CH$_{3}$OH toward two massive cluster-forming clumps, NGC2264-C and NGC2264-D, using the Nobeyama 45-m radio telescope. We derive an $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) integrated intensity ratio map, showing a higher value at clumps including 2MASS point sources at the northern part of NGC2264-D. Possible interpretations of the $I$(HC$_{3}$N)/$I$(CH$_{3}$OH) ratio are discussed. We have also observed molecular emission lines from CCS and N$_{2}$H$^{+}$ toward five positions in each clump. We investigate the $N$(N$_{2}$H$^{+}$)/$N$(CCS) and $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) column density ratios among the ten positions in order to test whether they can be used as chemical evolutionary indicators in these clumps. The $N$(N$_{2}$H$^{+}$)/$N$(CCS) ratio shows a very high value toward a bright embedded IR source (IRS1), whereas the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) ratio at IRS1 is comparable with those at the other positions. These results suggest that UV radiation affects the chemistry around IRS1. We find that there are positive correlations between these column density ratios and the excitation temperatures of N$_{2}$H$^{+}$, which implies the chemical evolution of clumps. These chemical evolutionary indicators likely reflect the combination of evolution along the filamentary structure and evolution of each clump.

We report a first clear detection of the Zeeman splitting of a CCS emission line at 45 GHz toward a nearby prestellar dense filament, Taurus Molecular Cloud-1. We observed HC$_3$N non-Zeeman line simultaneously as the CCS line, and did not detect any significant splitting of HC$_3$N line. Thus, we conclude that our detection of the CCS Zeeman splitting is robust. The derived \textcolor{black}{line-of-sight} magnetic field strength is about 117 $\pm$ 21 $\mu$G, which corresponds to the normalized mass-to-magnetic flux ratio of 2.2 if we adopt the inclination angle of 45$^\circ$. Thus, we conclude that the TMC-1 filament is magnetically supercritical. Recent radiative transfer calculations of CCS and HC$_3$N lines along the line of sight suggest that the filament is collapsing with a speed of $\sim$ 0.6 km s$^{-1}$, which is comparable to three times the isothermal sound speed. This infall velocity appears to be consistent with the evolution of a gravitationally-infalling core.

We report the statistical physical properties of the C$^{18}$O($J=1-0$) clumps present in a prominent cluster-forming region, Cygnus X, using the dataset obtained by the Nobeyama 45-m radio telescope. This survey covers 9 deg$^2$ of the north and south regions of Cygnus X, and totally 174 C$^{18}$O clumps are identified using the dendrogram method. Assuming a distance of 1.4 kpc, these clumps have radii of 0.2-1 pc, velocity dispersions of $<2.2~\mathrm{km~s^{-1}}$, gas masses of 30-3000 $M_\odot$, and H$_2$ densities of (0.2-5.5)$\times10^4~\mathrm{cm^{-3}}$. We confirm that the C$^{18}$O clumps in the north region have a higher H$_2$ density than those in the south region, supporting the existence of a difference in the evolution stages, consistent with the star formation activity of these regions. The difference in the clump properties of the star-forming and starless clumps is also confirmed by the radius, velocity dispersion, gas mass, and H$_2$ density. The average virial ratio of 0.3 supports that these clumps are gravitationally bound. The C$^{18}$O clump mass function shows two spectral index components, $\alpha=-1.4$ in 55-140 $M_\odot$ and $\alpha=-2.1$ in $>140~M_\odot$, which are consistent with the low- and intermediate-mass parts of the Kroupa's initial mass function. The spectral index in the star-forming clumps in $>140~M_\odot$ is consistent with that of the starless clumps in 55-140 $M_\odot$, suggesting that the latter will evolve into star-forming clumps while retaining the gas accretion. Assuming a typical star formation efficiency of molecular clumps (10%), about ten C$^{18}$O clumps having a gas mass of $>10^3~M_\odot$ will evolve into open clusters containing one or more OB stars.

We have carried out observations of CCH and its two $^{13}$C isotopologues, $^{13}$CCH and C$^{13}$CH, in the 84 - 88 GHz band toward two starless cores, L1521B and L134N (L183), using the Nobeyama 45 m radio telescope. We have detected C$^{13}$CH with a signal-to-noise (S/N) ratio of 4, whereas no line of $^{13}$CCH was detected in either the dark clouds. The column densities of the normal species were derived to be ($1.66 \pm 0.18$)$\times 10^{14}$ cm$^{-2}$ and ($7.3 \pm 0.9$)$\times 10^{13}$ cm$^{-2}$ ($1 \sigma$) in L1521B and L134N, respectively. The column density ratios of $N$(C$^{13}$CH)/$N$($^{13}$CCH) were calculated to be $>1.1$ and $>1.4$ in L1521B and L134N, respectively. The characteristic that $^{13}$CCH is less abundant than C$^{13}$CH is likely common for dark clouds. Moreover, we find that the $^{12}$C/$^{13}$C ratios of CCH are much higher than those of HC$_{3}$N in L1521B by more than a factor of 2, as well as in Taurus Molecular Cloud-1 (TMC-1). In L134N, the differences in the $^{12}$C/$^{13}$C ratios between CCH and HC$_{3}$N seem to be smaller than those in L1521B and TMC-1. We discuss the origins of the $^{13}$C isotopic fractionation of CCH and investigate possible routes that cause the significantly high $^{12}$C/$^{13}$C ratio of CCH especially in young dark clouds, with the help of chemical simulations. The high $^{12}$C/$^{13}$C ratios of CCH seem to be caused by reactions between hydrocarbons (e.g., CCH, C$_{2}$H$_{2}$, $l,c$-C$_{3}$H) and C$^{+}$.

Recent radio astronomical observations have revealed that HC$_{5}$N, the second shortest cyanopolyyne (HC$_{2n+1}$N), is abundant around some massive young stellar objects (MYSOs), which is not predicted by classical carbon-chain chemistry. For example, the observed HC$_{5}$N abundance toward the G28.28$-$0.36 MYSO is higher than that in L1527, which is one of the warm carbon chain chemistry (WCCC) sources, by more than one order of magnitude (Taniguchi et al., 2017). In this paper, we present chemical simulations of hot-core models with a warm-up period using the astrochemical code Nautilus. We find that the cyanopolyynes are formed initially in the gas phase and accreted onto the bulk and surface of granular ice mantles during the lukewarm phase, which occurs at $25 < T < 100$ K. In slow warm-up period models, the peak abundances occur as the cyanopolyynes desorb from dust grains after the temperature rises above 100 K. The lower limits of the abundances of HC$_{5}$N, CH$_{3}$CCH, and CH$_{3}$OH observed in the G28.28$-$0.36 MYSO can be reproduced in our hot-core models, after their desorption from dust grains. Moreover, previous observations suggested chemical diversity in envelopes around different MYSOs. We discuss possible interpretations of relationships between stages of the star-formation process and such chemical diversity, such as the different warm-up timescales. This timescale depends not only on the mass of central stars but also on the relationship between the size of warm regions and their infall velocity.

A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions associated with mixed poloidal and toroidal components of magnetic fields are successfully obtained in generic (non-circular) spacetimes. We introduce the formulation of the problem and the numerical method in detail, then present examples of extremely magnetized compact star solutions and their convergence tests. It is found that, in extremely magnetized stars, the stellar matter can be expelled from the region of strongest toroidal fields. Hence we conjecture that a toroidal electro-vacuum region may appear inside of the extremely magnetized compact stars, which may seem like the neutron star becoming the strongest toroidal solenoid coil in the universe.

We present the results of mapping observations toward a nearby starless filamentary cloud, the Taurus Molecular Cloud 1 (TMC-1), in the CCS(JN=43-32, 45.379033 GHz) emission line, using the Nobeyama 45-m telescope. The map shows that the TMC-1 filament has a diameter of ~0.1 pc and a length of ~0.5 pc at a distance of 140 pc. The position-velocity diagrams of CCS clearly indicate the existence of velocity-coherent substructures in the filament. We identify 21 substructures that are coherent in the position-position-velocity space by eye. Most of the substructures are elongated along the major axis of the TMC-1 filament. The line densities of the subfilaments are close to the critical line density for the equilibrium (~17 Mo/pc for the excitation temperature of 10 K), suggesting that self-gravity should play an important role in the dynamics of the subfilaments.

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.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

We have carried out survey observations of molecular emission lines from HC$_{3}$N, N$_{2}$H$^{+}$, CCS, and cyclic-C$_{3}$H$_{2}$ in the 81$-$94 GHz band toward 17 high-mass starless cores (HMSCs) and 28 high-mass protostellar objects (HMPOs) with the Nobeyama 45-m radio telescope. We have detected N$_{2}$H$^{+}$ in all of the target sources except one and HC$_{3}$N in 14 HMSCs and in 26 HMPOs. We investigate the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) column density ratio as a chemical evolutionary indicator of massive cores. Using the Kolmogorov-Smirnov (K-S) test and Welch's t test, we confirm that the $N$(N$_{2}$H$^{+}$)/$N$(HC$_{3}$N) ratio decreases from HMSCs to HMPOs. This tendency in high-mass star-forming regions is opposite to that in low-mass star-forming regions. Furthermore, we found that the detection rates of carbon-chain species (HC$_{3}$N, HC$_{5}$N, and CCS) in HMPOs are different from those in low-mass protostars. The detection rates of cyanopolyynes (HC$_{3}$N and HC$_{5}$N) are higher and that of CCS is lower in high-mass protostars, compared to low-mass protostars. We discuss a possible interpretation for these differences.

We have carried out observations in the 42$-$46 and 82$-$103 GHz bands with the Nobeyama 45-m radio telescope, and in the 338.2$-$339.2 and 348.45$-$349.45 GHz bands with the ASTE 10-m telescope toward three high-mass star-forming regions containing massive young stellar objects (MYSOs), G12.89+0.49, G16.86$-$2.16, and G28.28$-$0.36. We have detected HC$_{3}$N including its $^{13}$C and D isotopologues, CH$_{3}$OH, CH$_{3}$CCH, and several complex organic molecules (COMs). Combining our previous results of HC$_{5}$N in these sources, we compare the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratios in the three observed sources. The ratio in G28.28$-$0.36 is derived to be $0.091^{+0.109}_{-0.039}$, which is higher than that in G12.89+0.49 by one order of magnitude, and that in G16.86$-$2.16 by a factor of $\sim 5$. We investigate the relationship between the $N$(HC$_{5}$N)/$N$(CH$_{3}$OH) ratio and the $N$(CH$_{3}$CCH)/$N$(CH$_{3}$OH) ratio. The relationships of the two column density ratios in G28.28$-$0.36 and G16.86$-$2.16 are similar to each other, while HC$_{5}$N is less abundant when compared to CH$_{3}$CCH in G12.89+0.49. These results imply a chemical diversity in the lukewarm ($T \sim 20-30$ K) envelope around MYSOs. Besides, several spectral lines from complex organic molecules, including very-high-excitation energy lines, have been detected toward G12.89+0.49, while the line density is significantly low in G28.28$-$0.36. These results suggest that organic-poor MYSOs are surrounded by a carbon-chain-rich lukewarm envelope (G28.28$-$0.36), while organic-rich MYSOs, namely hot cores, are surrounded by a CH$_{3}$OH-rich lukewarm envelope (G12.89+0.49 and G16.86$-$2.16).

We have carried out interferometric observations of cyanopolyynes, HC$_{3}$N, HC$_{5}$N, and HC$_{7}$N, in the 36 GHz band toward the G28.28$-$0.36 high-mass star-forming region using the Karl G. Jansky Very Large Array (VLA) Ka-band receiver. The spatial distributions of HC$_{3}$N and HC$_{5}$N are obtained. HC$_{5}$N emission is coincident with a 450 $\mu$m dust continuum emission and this clump with a diameter of $\sim 0.2$ pc is located at the east position from the 6.7 GHz methanol maser by $\sim 0.15$ pc. HC$_{7}$N is tentatively detected toward the clump. The HC$_{3}$N : HC$_{5}$N : HC$_{7}$N column density ratios are estimated at 1.0 : $\sim 0.3$ : $\sim 0.2$ at an HC$_{7}$N peak position. We discuss possible natures of the 450 $\mu$m continuum clump associated with the cyanopolyynes. The 450 $\mu$m continuum clump seems to contain deeply embedded low- or intermediate-mass protostellar cores, and the most possible formation mechanism of the cyanopolyynes is the warm carbon chain chemistry (WCCC) mechanism. In addition, HC$_{3}$N and compact HC$_{5}$N emission is detected at the edge of the 4.5 $\mu$m emission, which possibly implies that such emission is the shock origin.

We clarify the line-of-sight structure of the Taurus Molecular Cloud 1 (TMC-1) on the basis of the CCS($J_N=4_3-3_2$) and HC$_3$N($J=5-4$) spectral data observed at a very high velocity resolution and sensitivity of $\Delta V \simeq 0.0004$ km s$^{-1}$ ($=61$ Hz) and $\Delta T_{\rm mb} \simeq 40$ mK. The data were obtained toward the cyanopolyyne peak with $\sim$30 hours integration using the Z45 receiver and the PolariS spectrometer installed in the Nobeyama 45m telescope. Analyses of the optically thin $F=4-4$ and $5-5$ hyperfine lines of the HC$_3$N emission show that the spectra consist of four distinct velocity components with a small line width ($\lesssim 0.1$ km s$^{-1}$) at $V_{\rm LSR}=$5.727, 5.901, 6.064, and 6.160 km s$^{-1}$, which we call A, B, C, and D, respectively, in the order of increasing LSR velocities. Utilizing the velocity information of the four velocity components, we further analyzed the optically thicker CCS spectrum and the other hyperfine lines of the HC$_3$N emission by solving the radiative transfer to investigate how the four velocity components overlap along the line of sight. Results indicate that they are located in the order of A, B, C, and D from far side to near side to the observer, indicating that TMC-1 is shrinking, moving inward as a whole.

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 carried out survey observations of HC$_{3}$N and HC$_{5}$N in the 42$-$45 GHz band toward 17 high-mass starless cores (HMSCs) and 35 high-mass protostellar objects (HMPOs) with the Nobeyama 45 m radio telescope. We have detected HC$_{3}$N from 15 HMSCs and 28 HMPOs, and HC$_{5}$N from 5 HMSCs and 14 HMPOs, respectively. The average values of the column density of HC$_{3}$N are found to be ($5.7 \pm 0.7$)$\times 10^{12}$ and ($1.03 \pm 0.12$)$\times 10^{13}$ cm$^{-2}$ in HMSCs and HMPOs, respectively. The average values of the fractional abundance of HC$_{3}$N are derived to be ($6.6 \pm 0.8$)$\times 10^{-11}$ and ($3.6 \pm 0.5$)$\times 10^{-11}$ in HMSCs and HMPOs, respectively. We find that the fractional abundance of HC$_{3}$N decreases from HMSCs to HMPOs using the Kolmogorov-Smirnov test. On the other hand, its average value of the column density slightly increases from HMSCs to HMPOs. This may imply that HC$_{3}$N is newly formed in dense gas in HMPO regions. We also investigate the relationship between the column density of HC$_{3}$N in HMPOs and the luminosity-to-mass ratio ($L/M$), a physical evolutional indicator. The column density of HC$_{3}$N tends to decrease with increasing the $L/M$ ratio, which suggests that HC$_{3}$N is destroyed by the stellar activities.

We present an unbiased large-scale (9 deg$^2$) CN ($N$=1-0) and C$^{18}$O ($J$=1-0) survey of Cygnus-X conducted with the Nobeyama 45m Cygnus-X CO survey. CN and C$^{18}$O are detected in various objects towards the Cygnus-X North and South (e.g., DR17, DR18, DR21, DR22, DR23, and W75N). We find that CN/C$^{18}$O integrated intensity ratios are systematically different from region to region, and are especially enhanced in DR17 and DR18 which are irradiated by the nearby OB stars. This result suggests that CN/C$^{18}$O ratios are enhanced via photodissociation reactions. We investigate the relation between the CN/C$^{18}$O ratio and strength of the UV radiation field. As a result, we find that CN/C$^{18}$O ratios correlate with the far-UV intensities, $G_0$. We also find that CN/C$^{18}$O ratios decrease inside molecular clouds, where the interstellar UV radiation is reduced due to the interstellar dust extinction. We conclude that the CN/C$^{18}$O ratio is controlled by the UV radiation, and is a good probe of photon-dominated regions.

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.

Outcomes of numerical relativity simulations of massive core collapses or binary neutron star mergers with moderate masses suggest formations of rapidly and differentially rotating neutron stars. Subsequent fall back accretion may also amplify the degree of differential rotations. We propose new formulations for modeling differential rotations of those compact stars, and present selected solutions of differentially rotating, stationary, and axisymmetric compact stars in equilibriums. For the cases when rotating stars reach break-up velocities, the maximum masses of such rotating models are obtained.

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 observed the $J=5-4$ rotational lines of the normal species and three $^{13}$C isotopologues of HC$_{3}$N at the 45 GHz band toward two low-mass starless cores, L1521B and L134N (L183), using the Nobeyama 45 m radio telescope in order to study the main formation pathways of HC$_{3}$N in each core. The abundance ratios of the three $^{13}$C isotopologues in L1521B are derived to be [H$^{13}$CCCN]:[HC$^{13}$CCN]:[HCC$^{13}$CN] = 0.98 (+- 0.14) : 1.00 : 1.52 (+- 0.16) (1 sigma). The fractionation pattern is consistent with that at the cyanopolyyne peak in Taurus Molecular Cloud-1. This fractionation pattern suggests that the main formation pathway of HC$_{3}$N is the neutral-neutral reaction between C$_{2}$H$_{2}$ and CN. On the other hand, their abundance ratios in L134N are found to be [H$^{13}$CCCN]:[HC$^{13}$CCN]:[HCC$^{13}$CN]= 1.5 (+- 0.2) : 1.0 : 2.1 (+- 0.4) (1 sigma), which are different from those in L1521B. From this fractionation pattern, we propose that the reaction between HNC and CCH is a possible main formation pathway of HC$_{3}$N in L134N. We find out that the main formation pathways of the same molecule are not common even in the similar physical conditions. We discuss the possible factors to make a difference in fractionation pattern between L134N and L1521B/TMC-1.

We report the first detection of HC$_{5}$$^{15}$N with the $J=9-8$ rotational line from the cyanopolyyne peak in Taurus Molecular Cloud-1 (TMC-1 CP) using the 45-m radio telescope of the Nobeyama Radio Observatory. The column density of HC$_{5}$$^{15}$N is derived to be (1.9 +- 0.5)*$10^{11}$ cm$^{-2}$ (1 sigma). We apply the double isotope method to derive the $^{14}$N/$^{15}$N ratios of HC$_{5}$N and HC$_{3}$N in TMC-1 CP. The $^{14}$N/$^{15}$N ratios are calculated to be 344 +- 53 and 257 +- 54 for HC$_{5}$N and HC$_{3}$N, respectively. The $^{14}$N/$^{15}$N ratio of HC$_{5}$N is lower than the elemental ratio in the local interstellar medium (~440) and slightly higher than that of HC$_{3}$N in TMC-1 CP. Since HC$_{3}$N is formed via the neutral-neutral reaction between C$_{2}$H$_{2}$ and CN, the slightly higher $^{14}$N/$^{15}$N ratio of HC$_{5}$N may support our previous suggestions that the main formation mechanism of HC$_{5}$N is the ion-molecule reactions between hydrocarbon ions (C$_{5}$H$_{n}^{+}$) and nitrogen atoms.

We carried out line survey observations at the 26-30 GHz band toward the four high-mass star-forming regions containing hot cores, G10.30-0.15, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Robert C. Byrd Green Bank Telescope. We have detected HC5N from all of the sources, and HC7N from the three sources, except for G10.30-0.15. We further conducted observations of HC5N at the 42-46 GHz and 82-103 GHz bands toward the three sources, G12.89+0.49, G16.86-2.16, and G28.28-0.36, with the Nobeyama 45 m radio telescope. The rotational lines of HC5N with the high excitation energies (Eu/k=63-100 K), which are hardly excited in the cold dark clouds, have been detected from the three sources. The rotational temperatures of HC5N are found to be 13-20 K in the three sources. The detection of the lines with the high excitation energies and the derived rotational temperatures indicate that HC5N exists in the warm gas within 0.07-0.1 pc radii around massive young stellar objects. The column densities of HC5N in the three sources are derived to be (2.0+-2.8)*10^(13) cm-2 . We compare the ratios between N(HC5N) the column density of HC5N and W(CH3OH) the integrated intensity of the thermal CH3OH emission line among the three high-mass star-forming regions. We found a possibility of the chemical differentiation in the three high-mass star-forming regions; G28.28-0.36 shows the largest N (HC5N)/W (CH3OH) ratio of > 8.0*10^(14) in units of (K km s-1 )-1 cm-2 , while G12.89+0.49 and G16.86-2.16 show the smaller values (~2*10^(13) ).

Recently a new method to discover obscured active galactic nuclei (AGNs) by utilizing X-ray and Infrared data has been developed. We carried out a survey of H2O maser emission toward ten obscured AGNs with the Nobeyama 45-m telescope. We newly detected the maser emission with the signal-noise-ratio (SNR) of above 4 from two AGNs, NGC 1402 and NGC 7738. We also found a tentative detection with its SNR > 3 in NGC 5037. The detection rate of 20% is higher than those of previous surveys (usually several percents).

We aim to reveal the chemistry of carbon-chain molecules in high-mass star-forming regions from molecular level to chemical evolution, applying the methods established in low-mass star-forming regions. In this proceeding, we summarize the topic about the main formation pathway of HC3N. We carried out observations of three 13C isotopologues of HC3N toward a hot core G28.28-0.36 with the 45-m radio telescope of the Nobeyama Radio Observatory. From the observational results, we propose that the main formation pathway of HC3N is the neutral-neutral reaction between C2H2 and CN. We also compare the results among the different star-forming regions, from a low-mass starless core to a high-mass star-forming core.

We study quasiequilibrium solutions of triaxially deformed rotating compact stars -- a generalization of Jacobi ellipsoids under relativistic gravity and compressible equations of state (EOS). For relatively stiff (piecewise) polytropic EOSs, we find supramassive triaxial solutions whose masses exceed the maximum mass of the spherical solution, but are always lower than those of axisymmetric equilibriums. The difference in the maximum masses of triaxial and axisymmetric solutions depends sensitively on the EOS. If the difference turns out to be only about 10%, it will be strong evidence that the EOS of high density matter becomes substantially softer in the core of neutron stars. This finding opens a novel way to probe phase transitions of high density nuclear matter using detections of gravitational waves from new born neutron stars or magnetars under fallback accretion.

We observed the J=9-8 and 10-9 rotational lines of three13C isotopologues of HC3N in L1527 and G28.28-0.36 with the 45-m radio telescope of the Nobeyama Radio Observatory in order to constrain the main formation mechanisms of HC3N in each source. The abundance ratios of the three 13C isotopologues of HC3N are found to be 0.9 (0.2) : 1.00 : 1.29 (0.19) (1sigma) and 1.0 (0.2) : 1.00 : 1.47 (0.17) (1sigma) for [H13CCCN]: [HC13CCN]: [HCC13CN] in L1527 and G28.28-0.36, respectively. We recognize a similar 13C isotopic fractionation pattern that the abundances of H13CCCN and HC13CCN are comparable, and HCC13CN is more abundant than the others. Based on the results, we discuss the main formation pathway of HC3N. The 13C isotopic fractionation pattern derived from our observations can be explained by the neutral-neutral reaction between C2H2 and CN in both the low-mass (L1527) and high-mass (28.28-0.36) star forming regions.

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.

We observed the J = 9-8 and 16-15 rotational transitions of the normal species and five 13C isotopologues of HC5N to study its formation mechanisms toward the cyanopolyyne peak in Taurus Molecular Cloud-1, with the 45-m radio telescope of Nobeyama Radio Observatory. We detected the five 13C isotopologues with high signal-to-noise ratios between 12 and 20, as well as the normal species. The abundance ratios of the five 13C isotopologues of HC5N are found to be 1.00:0.97:1.03:1.05:1.16 (0.19) for [H13CCCCCN]:[HC13CCCCN]:[HCC13CCCN]:[HCCC13CCN]:[HCCCC13CN]. We do not find any significant differences among the five 13C isotopologues. The averaged [HC5N]/[13C isotopologues] abundance ratio is determined to be 94 (6), which is slightly higher than the local interstellar elemental 12C/13C ratio of 60-70. Possible formation pathways are discussed on the basis of these results.

We developed a dual-linear-polarization HEMT (High Electron Mobility Transistor) amplifier receiver system of the 45-GHz band (hereafter Z45), and installed it in the Nobeyama 45-m radio telescope. The receiver system is designed to conduct polarization observations by taking the cross correlation of two linearly-polarized components, from which we process full-Stokes spectroscopy. We aim to measure the magnetic field strength through the Zeeman effect of the emission line of CCS ($J_N=4_3-3_2$) toward pre-protostellar cores. A linear-polarization receiver system has a smaller contribution of instrumental polarization components to the Stokes $V$ spectra than that of the circular polarization system, so that it is easier to obtain the Stokes $V$ spectra. The receiver has an RF frequency of 42 $-$ 46 GHz and an intermediate frequency (IF) band of 4$-$8 GHz. The typical noise temperature is about 50 K, and the system noise temperature ranges from 100 K to 150K over the frequency of 42 $-$ 46 GHz. The receiver system is connected to two spectrometers, SAM45 and PolariS. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. PolariS is a newly-developed digital spectrometer with a finest frequency resolution of 60 Hz, having a capability to process the full-Stokes spectroscopy. The Half Power Beam Width (HPBW) of the beam was measured to be 37$"$ at 43 GHz. The main beam efficiency of the Gaussian main beam was derived to be 0.72 at 43 GHz. The SiO maser observations show that the beam pattern is reasonably round at about 10 \% of the peak intensity and the side-lobe level was less than 3 \% of the peak intensity. Finally, we present some examples of astronomical observations using Z45.

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.

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 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 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.

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.

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.

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.

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 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.

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.

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.

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.

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 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.

Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss our new fully general relativistic calculations of merging BHNS binaries, which use high-accuracy, low-eccentricity, conformal thin-sandwich configurations as initial data. Our evolutions are performed using the moving puncture method and include a fully relativistic, high-resolution shock-capturing hydrodynamics treatment. Focusing on systems in which the neutron star is irrotational and the black hole is nonspinning with a 3:1 mass ratio, we investigate the inspiral, merger, and disk formation in the system. We find that the vast majority of material is promptly accreted and no more than 3% of the neutron star's rest mass is ejected into a tenuous, gravitationally bound disk. We find similar results for mass ratios of 2:1 and 1:1, even when we reduce the NS compaction in the 2:1 mass ratio case. These ambient disks reach temperatures suitable for triggering SGRBs, but their masses may be too small to produce the required total energy output. We measure gravitational waveforms and compute the effective strain in frequency space, finding measurable differences between our waveforms and those produced by binary black hole mergers within the advanced LIGO band. These differences appear at frequencies corresponding to the emission that occurs when the NS is tidally disrupted and accreted by the black hole. The resulting information about the radius of the neutron star may be used to constrain the neutron star equation of state.

Binary neutron stars (NSNS) are expected to be among the leading sources of gravitational waves observable by ground-based laser interferometers and may be the progenitors of short-hard gamma ray bursts. We present a series of general relativistic NSNS coalescence simulations both for unmagnetized and magnetized stars. We adopt quasiequilibrium initial data for circular, irrotational binaries constructed in the conformal thin-sandwich (CTS) framework. We adopt the BSSN formulation for evolving the metric and a high-resolution shock-capturing scheme to handle the magnetohydrodynamics. Our simulations of unmagnetized binaries confirm the results of Shibata, Taniguchi and Uryu (2003). In cases in which the mergers result in a prompt collapse to a black hole, we are able to use puncture gauge conditions to extend the evolution and determine the mass of the material that forms a disk. We find that the disk mass is less than 2% of the total mass in all cases studied. We then add a small poloidal magnetic field to the initial configurations and study the subsequent evolution. For cases in which the remnant is a hypermassive neutron star, we see measurable differences in both the amplitude and phase of the gravitational waveforms following the merger. For cases in which the remnant is a black hole surrounded by a disk, the disk mass and the gravitational waveforms are about the same as the unmagnetized cases. Magnetic fields substantially affect the long-term, secular evolution of a hypermassive neutron star (driving `delayed collapse') and an accretion disk around a nascent black hole.

We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together with the relativistic equations of hydrostationary equilibrium. We adopt maximal slicing, assume spatial conformal flatness, and impose equilibrium boundary conditions on an excision surface (i.e., the apparent horizon) to model the black hole. In our previous treatment we adopted a "leading-order" approximation for a parameter related to the black-hole spin in these boundary conditions to construct approximately nonspinning black holes. Here we improve on the models by computing the black hole's quasilocal spin angular momentum and setting it to zero. As before, we adopt a polytropic equation of state with adiabatic index Gamma=2 and assume the neutron star to be irrotational. In addition to recomputing several sequences for comparison with our earlier results, we study a wider range of neutron star masses and binary mass ratios. To locate the innermost stable circular orbit we search for turning points along both the binding energy and total angular momentum curves for these sequences. Unlike for our previous approximate boundary condition, these two minima now coincide. We also identify the formation of cusps on the neutron star surface, indicating the onset of tidal disruption. Comparing these two critical binary separations for different mass ratios and neutron star compactions we distinguish those regions that will lead to a tidal disruption of the neutron star from those that will result in the plunge into the black hole of a neutron star more or less intact, albeit distorted by tidal forces.

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 ..)

We present new sequences of general relativistic, quasiequilibrium black hole-neutron star binaries. We solve for the gravitational field in the conformal thin-sandwich decomposition of Einstein's field equations, coupled to the equations of relativistic hydrostatic equilibrium for a perfect fluid. We account for the black hole by solving these equations in the background metric of a Schwarzschild black hole whose mass is much larger than that of the neutron star. The background metric is treated in Kerr-Schild as well as isotropic coordinates. For the neutron star, we assume a polytropic equation of state with adiabatic index Gamma=2, and solve for both irrotational and corotational configurations. By comparing the results of irrotational and synchronized configurations with the same background metric, we conclude that the effect of the rotation on the location of tidal break-up is only on the order of a few percent. The different choices in the background also lead to differences of order a few percent, which may be an indication of the level to which these configurations approximate quasiequilibrium.

We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spatial background geometry and do not assume extreme mass ratios. We adopt a Gamma=2 polytropic equation of state and focus on irrotational neutron star configurations as well as approximately nonspinning black holes. We present numerical results for ratios of the black hole's irreducible mass to the neutron star's ADM mass in isolation of M_{irr}^{BH}/M_{ADM,0}^{NS} = 1, 2, 3, 5, and 10. We consider neutron stars of baryon rest mass M_B^{NS}/M_B^{max} = 83% and 56%, where M_B^{max} is the maximum allowed rest mass of a spherical star in isolation for our equation of state. For these sequences, we locate the onset of tidal disruption and, in cases with sufficiently large mass ratios and neutron star compactions, the innermost stable circular orbit. We compare with previous results for black hole-neutron star binaries and find excellent agreement with third-order post-Newtonian results, especially for large binary separations. We also use our results to estimate the energy spectrum of the outgoing gravitational radiation emitted during the inspiral phase for these binaries.

Many of the recent numerical simulations of binary black holes in vacuum adopt the moving puncture approach. This successful approach avoids the need to impose numerical excision of the black hole interior and is easy to implement. Here we wish to explore how well the same approach can be applied to moving black hole punctures in the presence of relativistic hydrodynamic matter. First, we evolve single black hole punctures in vacuum to calibrate our BSSN (Baumgarte-Shapiro-Shibata-Nakamura) implementation and to confirm that the numerical solution for the exterior spacetime is invariant to any ``junk'' (i.e., constraint-violating) initial data employed in the black hole interior. Then we focus on relativistic Bondi accretion onto a moving puncture Schwarzschild black hole as a numerical testbed for our high-resolution shock-capturing relativistic hydrodynamics scheme. We find that the hydrodynamical equations can be evolved successfully in the interior without imposing numerical excision. These results help motivate the adoption of the moving puncture approach to treat the binary black hole-neutron star problem using conformal thin-sandwich initial data.

We construct quasiequilibrium sequences of black hole-neutron star binaries for arbitrary mass ratios by solving the constraint equations of general relativity in the conformal thin-sandwich decomposition. We model the neutron star as a stationary polytrope satisfying the relativistic equations of hydrodynamics, and account for the black hole by imposing equilibrium boundary conditions on the surface of an excised sphere (the apparent horizon). In this paper we focus on irrotational configurations, meaning that both the neutron star and the black hole are approximately nonspinning in an inertial frame. We present results for a binary with polytropic index n=1, mass ratio M_{irr}^{BH}/M_{B}^{NS}=5 and neutron star compaction M_{ADM,0}^{NS}/R_0=0.0879, where M_{irr}^{BH} is the irreducible mass of the black hole, M_{B}^{NS} the neutron star baryon rest-mass, and M_{ADM,0}^{NS} and R_0 the neutron star Arnowitt-Deser-Misner mass and areal radius in isolation, respectively. Our models represent valid solutions to Einstein's constraint equations and may therefore be employed as initial data for dynamical simulations of black hole-neutron star binaries.

We have calculated the first dynamical evolutions of merging black hole-neutron star binaries that treat the combined spacetime in a nonperturbative general relativistic framework. Using the conformal flatness approximation, we have studied how the location of the tidal disruption radius with respect to the the black hole horizon and innermost stable circular orbit (ISCO) affects the qualitative evolution of the system. Based on simple arguments, we show that for a binary mass ratio q>~0.24, tidal disruption occurs outside the ISCO, while the opposite is true for q<~0.24. When tidal disruption occurs sufficiently far outside the ISCO, mass is transferred unstably from the neutron star to the black hole, resulting in the complete disruption of the neutron star. When tidal disruption occurs slightly within the ISCO, we find that some of the mass forms an extremely hot disk around the black hole. The resulting configurations in this case are excellent candidates for the progenitors of short-hard gamma ray bursts.

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.

The recent localization of some short-hard gamma ray bursts (GRBs) in galaxies with low star formation rates has lent support to the suggestion that these events result from compact object binary mergers. We discuss how new simulations in general relativity are helping to identify the central engine of short-hard GRBs. Motivated by our latest relativistic black hole-neutron star merger calculations, we discuss a scenario in which these events may trigger short-hard GRBs, and compare this model to competing relativistic models involving binary neutron star mergers and the delayed collapse of hypermassive neutron stars. Distinguishing features of these models may help guide future GRB and gravitational wave observations to identify the nature of the sources.

We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness approximation, and assume that the black hole mass is sufficiently large compared to that of the neutron star so that the black hole remains fixed in space. Using a spheroidal spectral methods solver, we solve the resulting field equations for a neutron star orbiting a Schwarzschild black hole. The matter is evolved using a relativistic, Lagrangian, smoothed particle hydrodynamics (SPH) treatment. We take as our initial data recent quasiequilibrium models for synchronized neutron star polytropes generated as solutions of the conformal thin-sandwich (CTS) decomposition of the Einstein field equations. We are able to construct from these models relaxed SPH configurations whose profiles show good agreement with CTS solutions. Our adiabatic evolution calculations for neutron stars with low compactness show that mass transfer, when it begins while the neutron star orbit is still outside the innermost stable circular orbit, is more unstable than is typically predicted by analytical formalisms. This dynamical mass loss is found to be the driving force in determining the subsequent evolution of the binary orbit and the neutron star, which typically disrupts completely within a few orbital periods. The majority of the mass transferred onto the black hole is accreted promptly; a significant fraction (~30%) of the mass is shed outward as well, some of which will become gravitationally unbound and ejected completely from the system. The remaining portion forms an accretion disk around the black hole, and could provide the energy source for short-duration gamma ray bursts.

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 present calculations of quasiequilibrium sequences of irrotational binary neutron stars based on realistic equations of state (EOS) for the whole neutron star interior. Three realistic nuclear EOSs of various softness and based on different microscopic models have been joined with a recent realistic EOS of the crust, giving in this way three different EOSs of neutron-star interior. Computations of quasiequilibrium sequences are performed within the Isenberg-Wilson-Mathews approximation to general relativity. For all evolutionary sequences, the innermost stable circular orbit (ISCO) is found to be given by mass-shedding limit (Roche lobe overflow). The EOS dependence on the last orbits is found to be quite important: for two 1.35 M_sol neutron stars, the gravitational wave frequency at the ISCO (which marks the end of the inspiral phase) ranges from 800 Hz to 1200 Hz, depending upon the EOS. Detailed comparisons with 3rd order post-Newtonian results for point-mass binaries reveals a very good agreement until hydrodynamical effects (dominated by high-order functions of frequency) become important, which occurs at a frequency ranging from 500 Hz to 1050 Hz, depending upon the EOS.

Quasiequilibrium sequences of binary neutron stars are numerically calculated in the framework of the Isenberg-Wilson-Mathews (IWM) approximation of general relativity. The results are presented for both rotation states of synchronized spins and irrotational motion, the latter being considered as the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation of state and compute several evolutionary sequences of binary systems composed of different-mass stars as well as identical-mass stars with adiabatic indices gamma=2.5, 2.25, 2, and 1.8. From our results, we propose as a conjecture that if the turning point of binding energy (and total angular momentum) locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of the adiabatic index gamma_0, that of the ADM mass (and total angular momentum) should exist in the IWM approximation of general relativity for the same value of the adiabatic index.

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.

Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the determination of M/R will allow very strong constraints to be placed on the equation of state of nuclear matter at high densities.

We present the first computations of quasiequilibrium binary neutron stars with different mass components in general relativity, within the Isenberg-Wilson-Mathews approximation. We consider both cases of synchronized rotation and irrotational motion. A polytropic equation of state is used with the adiabatic index gamma=2. The computations have been performed for the following combinations of stars: (M/R)_{star 1} vs. (M/R)_{star 2} = 0.12 vs. (0.12, 0.13, 0.14), 0.14 vs. (0.14, 0.15, 0.16), 0.16 vs. (0.16, 0.17, 0.18), and 0.18 vs. 0.18, where (M/R) denotes the compactness parameter of infinitely separated stars of the same baryon number. It is found that for identical mass binary systems there is no turning point of the binding energy (ADM mass) before the end point of the sequence (mass shedding point) in the irrotational case, while there is one before the end point of the sequence (contact point) in the synchronized case. On the other hand, in the different mass case, the sequence ends by the tidal disruption of the less massive star (mass shedding point). It is then more difficult to find a turning point in the ADM mass. Furthermore, we find that the deformation of each star depends mainly on the orbital separation and the mass ratio and very weakly on its compactness. On the other side, the decrease of the central energy density depends on the compactness of the star and not on that of the companion.

We study equilibrium sequences of close binary systems on circular orbits and composed of different mass stars with polytropic equation of state in Newtonian gravity. The solving method is a multidomain spectral method which we have recently developed. The computations are performed for both cases of synchronized and irrotational binary systems with adiabatic indices $\gamma = 3, 2.5, 2.25, 2$ and 1.8, and for three mass ratios: $M_1/M_2 =0.5, 0.2$ and 0.1. It is found that the equilibrium sequences always terminate at a mass shedding limit (appearance of a cusp on the surface of the less massive star). For synchronized binaries, this contrasts with the equal mass case, where the sequences terminate instead by a contact configuration. Regarding the turning point of the total angular momentum (or total energy) along a sequence, we find that it is difficult to get it for small mass ratios. Indeed, we do not find any turning points for $M_1/M_2 \le 0.5$ in the irrotational case. However, in the synchronized case, it becomes easier again to find one for mass ratios much smaller than $M_1/M_2 \sim 0.2$.

We study equilibrium sequences of close binary systems composed of identical polytropic stars in Newtonian gravity. The solving method is a multi-domain spectral method which we have recently developed. An improvement is introduced here for accurate computations of binary systems with stiff equation of state ($\gamma > 2$). The computations are performed for both cases of synchronized and irrotational binary systems with adiabatic indices $\gamma=3,~2.5,~2.25,~2$ and 1.8. It is found that the turning points of total energy along a constant-mass sequence appear only for $\gamma \ge 1.8$ for synchronized binary systems and $\gamma \ge 2.3$ for irrotational ones. In the synchronized case, the equilibrium sequences terminate by the contact between the two stars. On the other hand, for irrotational binaries, it is found that the sequences terminate at a mass shedding limit which corresponds to a detached configuration.

We present a numerical method to compute quasiequilibrium configurations of close binary neutron stars in the pre-coalescing stage. A hydrodynamical treatment is performed under the assumption that the flow is either rigidly rotating or irrotational. The latter state is technically more complicated to treat than the former one (synchronized binary), but is expected to represent fairly well the late evolutionary stages of a binary neutron star system. As regards the gravitational field, an approximation of general relativity is used, which amounts to solving five of the ten Einstein equations (conformally flat spatial metric). The obtained system of partial differential equations is solved by means of a multi-domain spectral method. Two spherical coordinate systems are introduced, one centered on each star; this results in a precise description of the stellar interiors. Thanks to the multi-domain approach, this high precision is extended to the strong field regions. The computational domain covers the whole space so that exact boundary conditions are set to infinity. Extensive tests of the numerical code are performed, including comparisons with recent analytical solutions. Finally a constant baryon number sequence (evolutionary sequence) is presented in details for a polytropic equation of state with gamma=2.

Solutions to equilibrium sequences of irrotational binary polytropic stars in Newtonian gravity are expanded in a power of $\epsilon=a_0/R$, where R and $a_0$ are the orbital separation of the binary system and the radius of each star for $R=\infty$. For each order of $\epsilon$, we should solve ordinary differential equations for arbitrary polytropic indices n. We show solutions for polytropic indices n= 0.5, 1, 1.5 and 2 up to $\epsilon^6$ orders. Our semi-analytic solutions can be used to check the validity of numerical solutions.

We investigate the effects of the magnetic dipole-dipole coupling and the electromagnetic radiation on the frequency evolution of gravitational waves from inspiralling binary neutron stars with magnetic dipole moments. This study is motivated by the discovery of the superstrongly magnetized neutron stars, i.e., magnetar. We derive the contributions of the magnetic fields to the accumulated cycles in gravitational waves as $N_{mag} \sim 6 \times 10^{-3} (H/10^{16}{\rm G})^{2}$, where $H$ denotes the strength of the polar magnetic fields of each neutron star in the binary system. It is found that the effects of the magnetic fields will be negligible for the detection and the parameter estimation of gravitational waves, if the upper limit for magnetic fields of neutron stars are less than $\sim 10^{16}$G, which is the maximum magnetic field observed in the soft gamma repeaters and the anomalous X-ray pulsars up to date. We also discuss the implications of electromagnetic radiation from the inspiralling binary neutron stars for the precursory X-ray emission prior to the gamma ray burst observed by the Ginga satellite.

A solution to an equilibrium of irrotational binary polytropic stars in Newtonian gravity is expanded in a power of \epsilon=a_0/R, where R and a_0 are the separation of the binary system and the radius of each star for R=\infty. For the polytropic index n=1, the solutions are given almost analytically up to order \epsilon^6. We have found that in general an equilibrium solution should have the velocity component along the orbital axis and that the central density should decrease when R decreases. Our almost analytic solutions can be used to check the validity of numerical solutions.

The first post-Newtonian (PN) hydrostatic equations for an irrotational fluid are solved for an incompressible binary system. The equilibrium configuration of the binary system is given by a small deformation from the irrotational Darwin-Riemann ellipsoid which is the solution at Newtonian order. It is found that the orbital separation at the innermost stable circular orbit (ISCO) decreases when one increases the compactness parameter $M_{\ast}/c^2 a_{\ast}$, in which $M_{\ast}$ and $a_{\ast}$ denote the mass and the radius of a star, respectively. If we compare the 1PN angular velocity of the binary system at the ISCO in units of $\sqrt{M_{\ast}/a_{\ast}^3}$ with that of Newtonian order, the angular velocity at the ISCO is almost the same value as that at Newtonian order when one increases the compactness parameter. Also, we do not find the instability point driven by the deformation at 1PN order, where a new sequence bifurcates throughout the equilibrium sequence of the binary system until the ISCO. We also investigate the validity of an ellipsoidal approximation, in which a 1PN solution is obtained assuming an ellipsoidal figure and neglecting the deformation. It is found that the ellipsoidal approximation gives a fairly accurate result for the total energy, total angular momentum and angular velocity. However, if we neglect the velocity potential of 1PN order, we tend to overestimate the angular velocity at the ISCO regardless of the shape of the star (ellipsoidal figure or deformed figure).

The innermost stable circular orbits (ISCO) of coalescing neutron star-black hole binary are studied taking into account both the tidal and relativistic effects. We adopt the generalized pseudo-Newtonian potential to mimic the general relativistic effects of gravitation. It is found that the separation of the neutron star and the black hole at the ISCO is greater than that obtained by using either the Newtonian potential or the second post-Newtonian equation of motion of point mass. In equal mass cases, it is found that for $\bar{a}/m_1 > 3.5$ with $\bar{a}$ and $m_1$ being the mean radius and the mass of the neutron star, the tidal effect dominates the stability of the binary system while for $\bar{a}/m_1 < 3.5$, the relativistic effect, i.e. the fact that the interaction potential has unstable orbit, does.