We have derived new physical quantities for several High-Mass X-ray Binaries (HMXBs) with supergiant (SG) companions through their cyclotron lines. The parameters are: the terminal velocity of the wind, the mass loss rate of the donor, the effective temperature and the magnetic fields. These parameters influence significantly the improvement of the model of accretion. In spite of the variety of their observational properties, the corresponding magnetic field is around B ~ 10^12 G. This result can be constrained by the effects on stellar evolution. In addition, we have performed a segmentation in the parameter space of donors intended for several SG-HMXB listed in our sample set. The parameter space can be categorized into five regimes depending on the possibility of disk formation associated with accretion from the stellar wind. This can give a quantitative clarification of the observed variability and the properties of these objects. We show that, when these systems come into the direct accretion region, systems with corresponding parameters can emit X-rays.

Recent X-ray observations have revealed the complexity and diversity of high-mass X-ray binaries (HMXBs). This diversity challenges a classical understanding of the accretion process onto the compact objects. In this study, we reinforce the conventional concept of the nature of wind-fed accretion onto a neutron star considering the geometrical effect of radiatively accelerated wind, and re-evaluate the transported angular momentum by using a simple wind model. Our results suggest that even in an OB-type HMXB fed by stellar wind, a large amount of angular momentum could be transported to form an accretion disk due to the wind-inhomogeneity, if the binary separation is tight enough and/or stellar wind is slow. We apply our model into actual systems such as LMC X-4 and OAO 1657-415, and discuss the possibility of disk formations in these systems.

We have considered the existence of neutron star magnetic field given by the cyclotron lines. We collected the data of 9 sources of high-mass X-ray binaries with supergiant companions as a case of testing our model, to demonstrate their distribution and evolution. The wind velocity, spin period and magnetic field strength are studied under different mass loss rate. In our model, correlations between mass-loss rate and wind velocity are found and can be tested in further observations. We examined the parameter space where wind accretion is allowed, avoiding barrier of rotating magnetic fields, with robust data of magnetic field of neutron stars. Our model shows that most of sources (6 of 9 systems) can be fed by the wind with relatively slow velocity, and this result is consistent with previous predictions. In a few sources, our model cannot fit under the standard wind accretion scenario. In these peculiar cases, other scenarios (disk formation, partial Roche lobe overflow) should be considered. This would provide information about the evolutionary tracks of various types of binaries, and thus show a clear dichotomy behavior in wind-fed X-ray binary systems.

In this study, we examine the properties of donor stars of the three recently discovered ultraluminous X-ray sources (ULXs) powered by rotating neutron stars. For this purpose, we constructed a theoretical relation between the X-ray luminosity ($L_{\rm{X}}$) and the orbital period ($P_{\rm{orb}}$) suitable for ULXs with neutron stars. By using this new $L_{\rm{X}} - P_{\rm{orb}}$ relation, we attempt to determine the currently unknown nature of donor stars in ULXs associated with neutron stars. Especially, comparing the observed properties with the stellar evolution tracks, we suggest that the donor star in the NGC5907 ULX-1 system is a moderately massive star with $6 - 12 \rm{M}_{\odot}$, just departing from the main sequence phase. The results of our models for the other two ULX systems (M82 X-2 and NGC7793 P-13) are consistent with those in previous studies. Although there are only a few samples, observed ULX systems with neutron stars seems to involve relatively massive donors.

Recently the first finding of a spin period of an accreting neutron star in M31 is reported. The observed spin period is 1.2 s and it shows 1.27 d modulations due to orbital motion. From the orbital information, the mass donor could not be a giant massive star. On the other hand, the observed properties are quite odd for typical low mass X-ray binaries. In this study, we compare observed binary parameters with theoretical models given by a stellar evolution track and make a restriction on the possible mass range of the donor. According to the standard stellar evolution model, the donor star should be larger than 1.5 solar mass, and this suggests that this system is a new member of a rare category, intermediate mass X-ray binary. The magnetic field strength of the neutron star suggested by spin-up/down tendency in this system supports the possibility of intermediate mass donor.

Periodic pulsations have been found in emission from the ultra-luminous X-ray source (ULX) M82 X-2, strongly suggesting that the emitter is a rotating neutron star rather than a black hole. However, the radiation mechanisms and accretion mode involved have not yet been clearly established. In this paper, we examine the applicability to this object of standard accretion modes for high mass X-ray binaries (HMXBs). We find that spherical wind accretion, which drives OB-type HMXBs, cannot apply here but that there is a natural explanation in terms of an extension of the picture for standard Be-type HMXBs. We show that a neutron star with a moderately strong magnetic field, accreting from a disc-shaped wind emitted by a Be-companion, could be compatible with the observed relation between spin and orbital period. A Roche lobe overflow picture is also possible under certain conditions.

We study an influence of X-ray photo-ionization from an accreting neutron star in a high mass X-ray binary. Our aim is to unveil a new principle governing X-ray luminosities of X-ray binaries, with a simple analysis of fluid equations simulating line-driven wind flow under influence of X-ray irradiation. In this study, we solve equation of motion of the accretion flow taking into account the line-driven acceleration and X-ray photo-ionization. Under the influence of X-ray irradiation, we find the flow equations take two types of solutions. The first solution is characterized by a slow wind velocity which causes a large accretion rate. The second solution is a fast wind flow which results in a small accretion rate. We find that only the solution with a fast wind and faint X-ray luminosity is a steady solution. On the other hand, slow wind solution with a large X-ray luminosity is not a realizable solution. In bright X-ray binary systems, X-ray luminosity would increase until strong X-ray reduces the line-driven acceleration and causes a stagnation of the wind. This result implies an important consequence that the X-ray luminosity of the wind-fed, X-ray emitting binary will be settled by the wind stagnation limit. At the same time, the fast wind solution with small X-ray luminosity also can be a steady state. Bright X-ray sources, such as Vela X-1, would have limiting luminosities of wind stagnation, while faint systems such a quiescent Super giant Fast X-ray Transients could be on the faint solutions.

Supergiant fast X-ray transient (SFXT) is a new class of the high mass X-ray binary that shows short X-ray flares. The physical mechanism of SFXT short flares is still open for discussion. The accretion process of dense clumps in stellar wind onto neutron star (NS) has been proposed as the origin of such short flares. In order to examine the applicability of the clumpy wind scenario, we focus on the accretion mode that depends on orbital parameters. Our goal is to impose restrictions on the orbital parameters of SFXT. Assuming a simple analytic model of clumpy wind, we investigate the condition where the size of accretion cylinder overcomes the clump size. The allowed parameter region for SFXT is restricted in a relatively narrow window in $P_{\rm{orb}} - e$ diagram. Binary systems with large eccentricities ($e \ga 0.4$) and moderate orbital periods ($P_{\rm{orb}} \sim 10$ d) are prone to show periodic X-ray outbursts which are characteristic for SFXT. We confirm that systems with a long orbital period of more than 100 days cannot produce bright X-ray flares in the simple clumpy wind scenario.

In this paper, we initiate a new study of steady funnel-flow accretion onto strongly magnetized neutron stars, including a full treatment of shock generation. As a first step, we adopt a simplified model considering the flow within Newtonian theory and neglecting radiative pressure and cooling. The flow is taken to start from an accretion disc and then to follow magnetic field lines, forming a transonic funnel flow onto the magnetic poles. A standing shock occurs at a certain point in the flow and beyond this material accretes subsonically onto the star with high pressure and density. We calculate the location of the standing shock and all other features of the flow within the assumptions of our model. Applications to observed X-ray pulsars are discussed.

We examine the published light curves (LCs) of 117 X-ray binary pulsars, focusing on the dependence of their light curves on the observed energy bands. It is found that the energy dependence of the LCs appears only when the X-ray luminosity is larger than ~ 5 x 10^36 erg/s. Assuming that the behavior of light curve is related to the radiative accretion column on the neutron star surface, this energy threshold can be considered as the observational proof of the accretion column formation proposed by Basko and Sunyaev. Once we can grasp the existence of radiative column, we can also obtain several useful informations on the neutron star properties. As an instance, we perform the statistical analysis of the orientation angle of the magnetic axis, and we find that the inclination angle of magnetic axis should be small in order to explain the observed statistics.

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

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.

It is usually believed that the Cowling approximation can give satisfactory solutions if the stars %are highly compressible have soft equations of state and/or if the strongly general relativistic stars are treated in the case of rigid rotation. Since, however, there have been no systematic studies about the accuracy of the Cowling approximation for differentially rotating compressible stars, we investigate eigenfrequencies and eigenfunctions of the oscillation modes in {\it rapidly} and {\it differentially} rotating compressible stars by employing the exact method including full-perturbations and the Cowling approximation. We have found that the Cowling approximation for f-mode oscillations is not a good approximation in rapidly and differentially rotating stars, although rapid rotation makes this approximation better for rigidly rotating stars. This result suggests that we must be careful when we apply the Cowling approximation to differentially rotating stars even in the framework of general relativity. On the other hand, the approximation will work well for r-modes even if the star is rotating differentially. Therefore, the Cowling approximation can be used as a strong tool for the investigation of r-mode oscillations in the general relativistic framework that it is difficult to compute including the perturbations of gravity.

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.

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.

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.

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.