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 have analyzed magnetized equilibrium states and showed a condition for appearance of the prolate and the toroidal magnetic field dominated stars by analytic approaches. Both observations and numerical stability analysis support that the magnetized star would have the prolate and the large internal toroidal magnetic fields. In this context, many investigations concerning magnetized equilibrium states have tried to obtain the prolate and the toroidal dominant solutions, but many of them have failed to obtain such configurations. Since the Lorentz force is a cross product of current density and magnetic field, the prolate shaped configurations and the large toroidal magnetic fields in stars require a special relation between current density and the Lorentz force. We have analyzed simple analytical solutions and found that the prolate and the toroidal dominant configuration require non force-free toroidal current density that flows in the opposite direction with respect to the bulk current within the star. Such current density results in the Lorentz force which makes the stellar shape prolate. Satisfying this special relation between the current density and the Lorentz force is a key for appearance of the prolate and the toroidal magnetic field dominated magnetized star.

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

We have shown analytically that shapes of incompressible stars could be prolate if appropriate meridional flows exist. Although this result is strictly valid only if either the meridional flow or the rotation is absent and the vorticity is associated uniformly with meridional flow, this implies that perpendicular forces against centrifugal and/or magnetic forces might play important roles within stars. A consequence of the presence of meridional flows might be to decrease stellar oblateness due to centrifugal and/or magnetic fields.

We have succeeded in obtaining magnetized star models that have extremely strong magnetic fields in the interior of the star. In our formulation, arbitrary functions of the magnetic flux function appear in the expression for the current density. By appropriately choosing the functional form for one of the arbitrary functions that corresponds to the distribution of the toroidal current density, we have obtained configurations with magnetic field distributions that are highly localized within the central part and near the magnetic axis region. The absolute values of the central magnetic fields are stronger than those of the surface region by two orders of magnitude. By applying our results to magnetars, the internal magnetic poloidal fields could be 10^17 G, although the surface magnetic fields are about 10^15 G in the case of magnetars. For white dwarfs, the internal magnetic poloidal fields could be 10^12 G, when the surface magnetic fields are 10^9-10^10 G .

We obtain the general forms for the current density and the vorticity from the integrability conditions of the basic equations which govern the stationary states of axisymmetric magnetized self-gravitating barotropic objects with meridional flows under the ideal magnetohydrodynamics (MHD) approximation. As seen from the stationary condition equations for such bodies, the presence of the meridional flows and that of the poloidal magnetic fields act oppositely on the internal structures. The different actions of these two physical quantities, the meridional flows and the poloidal magnetic fields, could be clearly seen through stationary structures of the toroidal gaseous configurations around central point masses in the framework of Newtonian gravity because the effects of the two physical quantities can be seen in an amplified way for toroidal systems compared to those for spheroidal stars. The meridional flows make the structures more compact, i.e. the widths of toroids thinner, while the poloidal magnetic fields are apt to elongate the density contours in a certain direction depending on the situation. Therefore, the simultaneous presence of the internal flows and the magnetic fields would work as if there were no such different actions within and around the stationary gaseous objects such as axisymmetric magnetized toroids with internal motions around central compact objects under the ideal MHD approximation, although these two quantities might exist in real systems.

We will show the importance of coexistence of oppositely flowing $\varphi$-currents for magnetized stars to sustain strong toroidal magnetic fields within the stars by an- alyzing stationary states of magnetized stars with surface currents which flow in the opposite direction with respect to the bulk currents within the stars. We have imposed boundary conditions for currents and toroidal magnetic fields to vanish outside the stars. It is important to note that these boundary conditions set an upper limit for the total current within the stars. This upper limit for the total current results in the presence of an upper limit for the magnitude of the energy for the toroidal magnetic fields of the stars. If the stars could have the toroidal surface currents which flow in the opposite directions to the internal toroidal currents, the positively flowing inter- nal toroidal currents can become stronger than the upper limit value of the current for configurations without surface toroidal currents. Thus the energies for the toroidal magnetic fields can become much larger than those for the magnetized stars with- out surface toroidal currents. We have also analyzed the same phenomena appearing in spherical incompressible stars for dipole-like magnetic fields with or without sur- face toroidal currents by employing the zero-flux-boundary method. We have applied those configurations with surface toroidal currents to magnetars and discussed their flares through which magnetic helicities could arise outside the stellar surfaces.

A fully geometrical treatment of general relativistic magnetohydrodynamics (GRMHD) is developed under the hypotheses of perfect conductivity, stationarity and axisymmetry. The spacetime is not assumed to be circular, which allows for greater generality than the Kerr-type spacetimes usually considered in GRMHD. Expressing the electromagnetic field tensor solely in terms of three scalar fields related to the spacetime symmetries, we generalize previously obtained results in various directions. In particular, we present the first relativistic version of the Soloviev transfield equation, subcases of which lead to fully covariant versions of the Grad-Shafranov equation and of the Stokes equation in the hydrodynamical limit. We have also derived, as another subcase of the relativistic Soloviev equation, the equation governing magnetohydrodynamical equilibria with purely toroidal magnetic fields in stationary and axisymmetric spacetimes.

We analyze the impact of the choice rotation law on equilibrium sequences of relativistic differentially-rotating neutron stars in axisymmetry. The maximum allowed mass for each model is strongly affected by the distribution of angular velocity along the radial direction and by the consequent degree of differential rotation. In order to study the wide parameter space implied by the choice of rotation law, we introduce a functional form that generalizes the so called "j-const. law" adopted in all previous work. Using this new rotation law we reproduce the angular velocity profile of differentially-rotating remnants from the coalescence of binary neutron stars in various 3-dimensional dynamical simulations. We compute equilibrium sequences of differentially rotating stars with a polytropic equation of state starting from the spherically symmetric static case. By analyzing the sequences at constant ratio, T/|W|, of rotational kinetic energy to gravitational binding energy, we find that the parameters that best describe the binary neutron star remnants cannot produce equilibrium configurations with values of T/|W| that exceed 0.14, the criterion for the onset of the secular instability.

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

We develop a numerical scheme for obtaining the $r$-mode oscillations of rapidly rotating relativistic stars. In the present scheme, we neglect all metric perturbations and only take account of the dynamics of the fluid in the background spacetime of the unperturbed star (the relativistic Cowling approximation). We also assume the star is barotropic, i.e. neutrally stable against convection under the assumption of adiabatic oscillations. Our numerical scheme is based on the Yoshida-Eriguchi formulation for the analysis of the general relativistic f-mode oscillations in the Cowling approximation and a general relativistic generalization of the Karino-Yoshida-Yoshida-Eriguchi's numerical scheme for obtaining oscillations of rapidly rotating Newtonian stars. By this new numerical scheme, frequencies of the $r$-mode oscillations are obtained as functions of the ratio of the rotational energy to the absolute value of the gravitational energy $T/|W|$ along sequences of polytropic equilibrium stars whose ratio of the pressure to the total energy density at the center of the star and the polytropic index are kept constant. It is found that the dimensionless oscillation frequency $\sigma/\Omega$ is a linearly decreasing function of $T/|W|$, where $\sigma$ and $\Omega$ are the oscillation frequency and the angular velocity of the star measured in an inertial frame at spatial infinity. We also find that oscillation frequencies of the $r$-modes are highly dependent on the relativistic factor $M/R$ of the star as already found in previous studies in which the slow rotation approximation has been used. Here $M$ and $R$ denote the mass and radius of the star, respectively.

We develop a formulation for constructing and examining rapidly rotating Newtonian neutron star models that contain two superfluids, taking account of the effect of the rotation velocity difference between two superfluids. We assume neutron stars to be composed of the superfluid neutrons and the mixture of the superfluid protons and the normal fluid electrons. To describe Newtonian dynamics of the two superfluids, the Newtonian version of the so-called two-fluid formalism is employed. The effect of the rotation velocity difference on the structure of equilibrium state is treated as a small perturbation to rapidly rotating superfluid stars whose angular velocities of two superfluids are assumed to be exactly the same. We derive basic equations for the perturbed structures of rapidly rotating superfluid stars due to the rotation velocity difference between two superfluids. Assuming the superfluids to obey a simple analytical equation of state proposed by Prix, Comer, and Andersson, we obtain numerical solutions for the perturbations and find that the density distributions of the superfluids are strongly dependent on the parameter $\sigma$ which appears in the analytical equation of state and characterizes the so-called symmetry energy. It is also found that if Prix et al.'s analytical equation of state is assumed, the perturbations can be represented in terms of the universal functions that are independent of the parameters of the equation of state.

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

We have studied the $f$-mode oscillations of differentially rotating polytropes by making use of the linear stability analysis. We found that the critical values of $T/|W|$ where the dynamical instability against the $m = 2$ $f$-mode oscillations sets in decrease down to $T/|W| \sim 0.20$ as the degree of differential rotation becomes higher. Here $m$ is an azimuthal mode number and $T$ and $W$ are the rotational energy and the gravitational potential energy, respectively. This tendency is almost independent of the compressibility of the polytropes. These are the {\it first exact results} of the linear stability analysis for the occurrence of the dynamical instability against the $m = 2$ $f$-modes.

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

We have computed the eigenfrequencies of f-modes for a constant-rest-mass sequences of rapidly rotating relativistic inviscid stars in differential rotation. The frequencies have been calculated neglecting the metric perturbations (the relativistic Cowling approximation) and expressed as a function of the ratio between the rotational kinetic energy and the absolute value of the gravitational energy of the stellar model beta=T/|W|. The zeros and the end-points of these sequences mark respectively the onset of the secular instability driven by gravitational radiation-reaction and the maximum value of beta at which an equilibrium model exists. In differentially rotating stars the secular stability limits appear at a beta larger than those found for uniformly rotating stars. Differential rotation, on the other hand, also allows for the existence of equilibrium models at values of beta larger than those for uniformly rotating stars, moving the end-point of the sequences to larger beta. As a result, for some degrees of differential rotation, the onset of the secular instability for f-modes is generally favoured by the presence of differential rotation.

We have computed quasiequilibrium sequences of synchronously rotating compact binary star systems with constant rest masses. This computation has been carried out by using the numerical scheme which is different from the scheme based on the conformally flat assumption about the space. Stars are assumed to be polytropes with polytropic indices of N=0.5, N=1.0, and N=1.5. Since we have computed binary star sequences with a constant rest mass, they provide approximate evolutionary tracks of binary star systems. For relatively stiff equations of state (N < 1.0), there appear turning points along the quasiequilibrium sequences plotted in the angular momentum -- angular velocity plane. Consequently secular instability against exciting internal motion sets in at those points. Qualitatively, these results agree with those of Baumgarte et al. who employed the conformally flat condition. We further discuss the effect of different equations of state and different strength of gravity by introducing two kinds of dimensionless quantities which represent the angular momentum and the angular velocity. Strength of gravity is renormalized in these quantities so that the quantities are transformed to values around unity. Therefore we can clearly see relations among quasiequilibrium sequences for a wide variety of strength of gravity and for different compressibility.

For the analysis of the r-mode oscillation of hot young neutron stars, it is necessary to consider the effect of it differential rotation, because viscosity is not strong enough for differentially rotating young neutron stars to be lead to uniformly rotating configurations on a very short time scale after their birth. In this paper, we have developed a numerical scheme to solve r-mode oscillations of differentially rotating polytropic inviscid stars. This is the extended version of the method which was applied to compute r-mode oscillations of uniformly rotating Newtonian polytropic stars. By using this new method, we have succeeded in obtaining eigenvalues and eigenfunctions of r-mode oscillations of differentially rotating polytropic stars. Our numerical results show that as the degree of differential rotation is increased, it becomes more difficult to solve r-mode oscillations for slightly deformed configurations from sphere compared to solving r-mode oscillations of considerably deformed stars. One reason for it seems that for slightly deformed stars corotation points appear near the surface region if the degree of differential rotation is strong enough. This is similar to the situation that the perturbational approach of r-mode oscillations for it slowly rotating stars in general relativity results in a singular eigenvalue problem.

By using the Cowling approximation, quasi-radial modes of rotating general relativistic stars are computed along equilibrium sequences from non-rotating to maximally rotating models. The eigenfrequencies of these modes are decreasing functions of the rotational frequency. The eigenfrequency curve of each mode as a function of the rotational frequency has discontinuities, which arise from the avoided crossing with other curves of axisymmetric modes.

We have developed a new numerical scheme to solve r-mode oscillations of {\it rapidly rotating polytropic stars} in Newtonian gravity. In this scheme, Euler perturbations of the density, three components of the velocity are treated as four unknown quantities together with the oscillation frequency. For the basic equations of oscillations, the compatibility equations are used instead of the linearized equations of motion. By using this scheme, we have solved the classical r-mode oscillations of rotational equilibrium sequences of polytropes with the polytropic indices $N = 0.5, 1.0$ and 1.5 for $m = 2, 3$ and 4 modes. Here $m$ is the rank of the spherical harmonics $Y_l^m$. These results have been applied to investigate evolution of uniformly rotating hot young neutron stars by considering the effect of gravitational radiation and viscosity. We have found that the maximum angular velocities of neutron stars are around 10-20% of the Keplerian angular velocity irrespective of the softness of matter. This confirms the results obtained from the analysis of r-modes with the slow rotation approximation employed by many authors.

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

We have investigated an influence of the r-mode instability on hypercritically accreting ($\dot{M}\sim 1M_\odot {y}^{-1}$) neutron stars in close binary systems during their common envelope phases based on the scenario proposed by Bethe et al. \shortcite{bethe-brown-lee}. On the one hand neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency. The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning-up of the star, whose final rotational frequency is well below the mass-shedding limit, typically as small as 10% of that of the mass-shedding state. On a very short time scale the rotational frequency tends to approach a certain constant value and saturates there as far as the amount of the accreted mass does not exceed a certain limit to collapse to a black hole. This implies that the similar mechanism of gravitational radiation as the so-called Wagoner star may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by the gravitational wave detectors such as LIGO II.

The first results of numerical analysis of classical r-modes of {\it rapidly} rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity are numerically solved without the slow rotation approximation. A critical curve of gravitational wave emission induced instability which restricts the rotational frequencies of hot young neutron stars is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the `evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning-down proceed. Rotational frequencies of $1.4M_{\odot}$ stars suffering from this instability decrease to around 100Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors who adopted the slow rotation approximation.

We propose a new numerical method to compute quasi-equilibrium sequences of general relativistic irrotational binary neutron star systems. It is a good approximation to assume that (1) the binary star system is irrotational, i.e. the vorticity of the flow field inside component stars vanishes everywhere (irrotational flow), and (2) the binary star system is in quasi-equilibrium, for an inspiraling binary neutron star system just before the coalescence as a result of gravitational wave emission. We can introduce the velocity potential for such an irrotational flow field, which satisfies an elliptic partial differential equation (PDE) with a Neumann type boundary condition at the stellar surface. For a treatment of general relativistic gravity, we use the Wilson--Mathews formulation, which assumes conformal flatness for spatial components of metric. In this formulation, the basic equations are expressed by a system of elliptic PDEs. We have developed a method to solve these PDEs with appropriate boundary conditions. The method is based on the established prescription for computing equilibrium states of rapidly rotating axisymmetric neutron stars or Newtonian binary systems. We have checked the reliability of our new code by comparing our results with those of other computations available. We have also performed several convergence tests. By using this code, we have obtained quasi-equilibrium sequences of irrotational binary star systems with strong gravity as models for final states of real evolution of binary neutron star systems just before coalescence. Analysis of our quasi-equilibrium sequences of binary star systems shows that the systems may not suffer from dynamical instability of the orbital motion and that the maximum density does not increase as the binary separation decreases.

A numerical method of mode analysis of rapidly rotating relativistic stellar models by the Cowling approximation is applied to rotating neutron stars with realistic equations of state. For selected equations of state, eigenvalues and eigenfunctions of f-modes are numerically solved for stellar models from non-rotating to maximally rotating states. Neutral points of the lower order f-modes are determined as a function of the stellar rotational frequency. As in the polytropic case, we find that the bar mode can have neutral points for models with relatively strong gravity. The rotational frequency at the neutral point increases as the gravitational mass of the model becomes larger. As astrophysical applications of our analysis, we discuss the time scales of gravitational radiation induced instability and the possibility of the resonant excitation of f-modes during inspiraling motion of compact binary systems.

We propose a new numerical method to calculate irrotational binary systems composed of compressible gaseous stars in Newtonian gravity. Assuming irrotationality, i.e. vanishing of the vorticity vector everywhere in the star in the inertial frame, we can introduce the velocity potential for the flow field. Using this velocity potential we can derive a set of basic equations for stationary states which consist of (i) the generalized Bernoulli equation, (ii) the Poisson equation for the Newtonian gravitational potential and (iii) the equation for the velocity potential with the Neumann type boundary condition. We succeeded in developing a new code to compute numerically exact solutions to these equations for the first time. Such irrotational configurations of binary systems are appropriate models for realistic neutron star binaries composed of inviscid gases, just prior to coalescence of two stars caused by emission of gravitational waves. Accuracies of our numerical solutions are so high that we can compute reliable models for fully deformed final stationary configurations and hence determine the inner most stable circular orbit of binary neutron star systems under the approximations of weak gravity and inviscid limit.

Circularly orbiting black hole-gaseous star close binary systems are examined by using numerically exact stationary configurations in the framework of Newtonian gravity. We have chosen a polytropic star for the fluid component of the binary system and considered two ideal situations: 1) a synchronously rotating star and 2) an irrotationally rotating star. They correspond to a rotating star under the influence of viscosity and to that in the inviscid limit, respectively. By analyzing the stationary sequences of binary systems with small separations, we can discuss the final stages of black hole-gaseous star close binary systems. Our computational results show that the binary systems reach the Roche(-Riemann) limit states or the Roche lobe filling states without suffering from hydrodynamical instability due to tidal force for a large parameter range of the mass ratio and the polytropic index. It is very likely that such stable Roche(-Riemann) limits or Roche lobe filling states survive even under the general relativistic effect. Therefore, at the final stage of the evolution which is caused by the emission of gravitational waves, the Roche overflow will occur instead of merging of a black hole and a star.

In this paper we have examined numerically exact configurations of close binary systems composed of incompressible fluids with internal flows. Component stars of binary systems are assumed to be circularly orbiting each other but rotating nonsynchronously with the orbital motion, i.e. stars in binary systems have steady motions seen from a rotational frame of reference. We have computed several equilibrium sequences by taking fully into account the tidal effect of Newtonian gravity without approximation. We consider two binary systems consisting of either 1) a point mass and a fluid star or 2) a fluid star and a fluid star. Each of them corresponds to generalization of the Roche--Riemann or the Darwin--Riemann problem, respectively. Our code can be applied to various types of incompressible binary systems with various mass ratios and spin as long as the vorticity is constant. We compare these equilibrium sequences of binaries with approximate solutions which were studied extensively by Lai, Rasio and Shapiro (LRS) as models for neutron star--neutron star (NS--NS) binary systems or black hole--neutron star (BH--NS) binary systems. Our results coincide qualitatively with those of LRS but are different from theirs for configurations with small separations. For these binary systems, our sequences show that dynamical or secular instability sets in as the separation decreases. The quantitative errors of the ellipsoidal approximation amount to 2\sim 10% for configurations near the instability point. Compared to the results of LRS, the separation of the stars at the point where the instability sets in is larger and correspondingly the orbital frequency at the critical state is smaller for our models.

We conduct a direct comparison of three different representative numerical codes for constructing models of rapidly rotating neutron stars in general relativity. Our aim is to evaluate the accuracy of the codes and to investigate how the accuracy is affected by the choice of interpolation, domain of integration and equation of state. In all three codes, the same physical parameters, equations of state and interpolation method are used. We construct 25 selected models for polytropic equations of state and 22 models with realistic neutron star matter equations of state. The three codes agree well with each other (typical agreement is better than 0.1% to 0.01%) for most models, except for the extreme assumption of uniform density stars. We conclude that the codes can be used for the construction of highly accurate initial data configurations for polytropes of index N > 0.5 (which typically correspond to realistic neutron stars), when the domain of integration includes all space and for realistic equations with no phase transitions. With the exception of the uniform density case, the obtained values of physical parameters for the models considered in this paper can be regarded as ``standard'' and we display them in detail for all models.

We have computed line profiles from self-gravitating toroids around black holes. The specific angular momentum of the toroids is assumed to be constant in space. The images of the toroids show peculiar feature in the rear side of the black holes. Concerning the line profiles, the red wing extends to the very small frequency region because the location of the inner edge is rather near the event horizon of the black hole and consequently the velocity of the inner edge of toroids can be faster than that of the Kepler disks.

We have succeeded in obtaining exact configurations of irrotational binary systems for compressible (polytropic) equations of state. Our models correspond to binaries of equal mass neutron star systems in the viscosity free limit. By using the obtained sequences of stationary states, the evolution of binary systems of irrotational neutron stars due to gravitational wave emission has been examined. For inviscid binary systems the spin angular velocity of each component in a detached phase is smaller than the orbital angular velocity at a contact phase. Thus the irrotational approximation during evolution of binary neutron stars due to gravitational wave emission can be justified. Our computational results show that the binary will never reach a dynamically unstable state before a contact phase even for rather stiff polytropes with the index N > 0.7, as the separation of two components decreases due to gravitational wave emission. This conclusion is quantitatively different from that of Lai, Rasio & Shapiro who employed approximate solutions for polytropic binary systems.

Relativistic Cowling approximation in which all metric perturbations are omitted is applied to non-axisymmetric infinitesimal oscillations of uniformly rotating general relativistic polytropes. Frequencies of lower order f-modes which are important in analysis of secular instability driven by gravitational radiation are investigated and neutral points of the mode along equilibrium sequences of rotating polytropes are determined. Since this approximation becomes more accurate as stars are more relativistic and/or as they rotates more rapidly, we will be able to analyze how a rotation period of a neutron star may be limited by this instability. Possible errors in determining neutral points caused by omitting metric perturbations are also estimated.