the be / x - ray binaries represent the largest subclass of high - mass x - ray binaries .
about two - thirds of the identified systems fall into this category .
these systems consist of , generally , a neutron star and a be star with a cool ( @xmath2 ) equatorial disc , which is geometrically thin and nearly keplerian .
the orbit is wide @xmath3 and usually eccentric ( @xmath4 ) ( e.g. , @xcite ) .
most of the be / x - ray binaries show only transient activity in the x - ray emission and are termed be / x - ray transients .
be / x - ray transients show periodical ( type i ) outbursts , which are separated by the orbital period and have the lumiocity of @xmath5 , and giant ( type ii ) outbursts of @xmath6 with no orbital modulation .
the transient nature in the x - ray activity of be / x - ray binaries is considered to result from the interaction between the accreted material from the circumsteller disc of the be star and the rotating magnetized neutron star @xcite . if the accreted material is dense enough to make the magnetospheric radius smaller than the corotation radius , the accretion on to the neutron star causes a bright x - ray emission ( the direct accretion regime ) . otherwise , the magnetospheric radius is larger than the corotation radius , and the accretion is prevented by the centrifugal inhibition .
the system is then in quiescence ( the propeller regime ) .    during the outbursts , the neutron stars in
most be / x - ray binaries show spin - up ( e.g. , 4u0115 + 63 , gs1843 - 02 , exo2030 + 375 and xtej1946 + 26 ) ( @xcite ; @xcite ; @xcite ; @xcite ) , and some even show quasi - periodic oscillations ( qpos ) ( e.g. , a0535 + 26 and exo2030 + 375 ) ( @xcite ; @xcite ) .
these observational features strongly suggest the presence of an accretion disc during the outbursts .
the accretion disc is also indicated to exist between the outbursts in several systems ( e.g. , a0535 + 26 ) ( @xcite ; see also @xcite ) . on the other hand
, @xcite searched the optical / infrared contribution from the accretion disc in a0535 + 26 during the x - ray quiescent phase , and found no signature of the accretion disc in propeller regime .
obviously , more work is needed to determine whether the accretion disc also exists in the quiescent state .
recently , @xcite studied the interaction between the be - star disc and the neutron star , using a three - dimensional ( 3d ) smoothed particle hydrodynamics ( sph ) code ( @xcite ; @xcite ) .
they also qualitatively discussed the relation between the x - ray luminosity and the accretion rate using the simulated mass - capture rate by the neutron star . in their simulations , however , the neutron star was modelled by a sink particle @xcite with the size of the roche lobe .
hence , no information was available how the captured material accretes on to the neutron star .
moreover , no direct comparison with the observed x - ray data was possible .    in this paper , we study the accretion on to the neutron star in be / x - ray binaries , using a 3d sph code
. the simulation should ideally take account of all the processes at work in a be / x - ray binary , such as the disc evolution around the be star , the mass transfer from the be - star disc to the neutron star , and the accretion on to the neutron star .
such a simulation , however , would require very intense computations .
therefore , we confine ourselves to simulate only the accretion flow around the neutron star , with an outer boundary condition constructed by the data of the particles captured by the neutron star from a high - resolution simulation by @xcite . this situation is illustrated in fig .
1 , where the right panel taken from their simulation shows the mass transfer from the be - star disc to the neutron star at periastron , while the left panel shows the accretion disc formed around the neutron star in one of our simulations .
simulations presented here were performed with a 3d sph code .
the sph code is basically the same as that used by @xcite , which is based on a version originally developed by benz ( @xcite ; @xcite ) .
the smoothing length is variable in time and space .
the sph equation with the standard cubic - spline kernel are integrated using a second - order runge - kutta - fehlberg integrator with individual time steps for each particle @xcite , which results in an enormous computational saving when a large range of dynamical time - scales are involved .
in our code , the accretion flow on to the neutron star is modelled by an ensemble of gas particles , each of which has an arbitrary and negligible mass chosen to be @xmath7 with a variable smoothing length .
note that the accretion flow is nonself - gravitating in our simulations . on the other hand , the be star and the neutron star
are modelled by two sink particles @xcite with corresponding masses . the gas particles which fall within a specified accretion radius
are accreted by the sink particles .
we assume that the neutron star has the fixed accretion radius of @xmath8 , @xmath9 being the semi - major axis of the binary .
for the be star , we adopt a variable accretion radius of @xmath10 , where @xmath11 is the roche - lobe radius for a circular binary .
this is because our interest is only in the flow around the neutron star .
an approximate formula for @xmath11 is given by @xmath12 ( see , for example , @xcite ; @xcite ) with the mass ratio @xmath13 , where @xmath14 and @xmath15 are the masses of the neutron star and the be star , respectively , and @xmath16 is the distance between the stars .    in disc simulations using a three dimensional sph code with the standard cubic spline kernel
, one has an approximate relation between the shakura - sunyaev viscosity parameter @xmath17 and the linear sph artificial viscosity parameter @xmath18 as @xmath19 if the non - linear sph artifical viscosity parameter @xmath20 is set equal to zero @xcite . here
@xmath21 is the sph smoothing length and @xmath22 is the scale - height of the disc .    in a simulation
we report in this paper , we adopt @xmath23 and @xmath24 in which case @xmath17 is variable in time and space . in the other simulations , however , we adopt a constant value of @xmath17 in order to roughly model the @xmath25 disc . in these simulations , @xmath26 was variable in time and space and @xmath27 .    in our simulations ,
we adopt a polytropic equation of state specified by the exponent @xmath28 , in which the effects of an energy dissipation resulting from viscosity and a radiative cooling are approximately boxed .
our binary configuration is the same as that of @xcite .
we set the binary orbit on the @xmath29-@xmath30 plane with the semi - major axis along the @xmath29-axis .
initially , the neutron star is at the periastron .
it orbits about the be star with the orbital period @xmath31 and the eccentricity @xmath32 .
the unit of time adopted here and throughout the paper is @xmath33 .
the be disc is coplanar with the binary orbital plane . as the be star , we take a b0v star with @xmath34 , @xmath35 and @xmath36 , where @xmath37 and @xmath38 are the stellar radius and effective temprature , respectively . for the neutron star , we take @xmath39 and @xmath40 .
these parameters are adopted to target 4u0115 + 63 , which has shown only type  ii x - ray outbursts followed by a short series of type  i x - ray outbursts @xcite .
as mentioned earlier , we concentrate only on the simulations of the accretion flow on to the neutron star .
this is done by constructing an outer boundary condition from the data of particles captured by the neutron star in a high - resolution simulation by @xcite , which was run for @xmath41 and finally had @xmath42 particles with @xmath43 .
we first fold the data on the orbital period over @xmath44 to reduce the fluctuation noise .
the phase dependence of the mass - capture rate by the neutron star obtained by this procedure is shown in fig .
[ fig : mdot_cap ] .
note that there is a strong phase - dependence in the mass - transfer rate .
the peak comes slightly after the periastron passage and decreases by two orders of magnitude by apastron passage .
we then fit the spatial and velocity distribution of the captured particles in each phase bin with separate gaussian functions .
this is for running simulations with an arbitrary mass - transfer rate from the be disc .
we finally model the mass - transfer process from the be disc to the neutron star by injecting gas particles at a given phase - dependent rate shown in fig .
2 with the spatial and velocity distributions specified by these gaussian functions . the injected particles are assumed to have the initial temperature of @xmath45 .    as described earlier
, the mass of each gas particle is arbitrary and negligible in our simulations .
all the mass - related quantities discussed later , such as the disc mass and the mass - accretion rate , are normalised by using the mass - capture rate shown in fig .  2 .    .
the mass - capture rate is measured in units of @xmath46 , where @xmath47 is highest local density in the be - star disc normalized by @xmath48 , a typical value for be stars .
the right axis shows the x - ray luminosity corresponding to the mass - capture rate with the x - ray emission efficiency @xmath49 , where @xmath50 is defined by @xmath51.,width=302,height=226 ]    .
, width=264 ]    . ,
width=264 ]    [ cols="<,^,^,^,^,^,^,^,^ " , ]      before describing the results of our simulations , we discuss below the possibility of the formation of an accretion disc around the neutron star .
we first consider the radius @xmath52 at which the material transferred from the be disc is circularized via the interaction with the previously transferred gas .
we assume the circularization radius @xmath52 to be approximately given by @xmath53 , where @xmath54 is the initial specific angular momentum of the transferred material .    in the upper panel of fig .
[ fig : rcirc ] , we show the dependence of @xmath52 on the orbital phase , based on the data taken from a simulation by @xcite .
it is noted that there is a strong phase - dependence of @xmath52 .
it has a minimum of @xmath55 slightly before periastron and then rapidly increases with phase . for comparison
, we also show the corotation radius , @xmath56 , for 4u0115 + 63 by the dotted line . here
, @xmath56 is given by @xmath57 where @xmath58 is the spin period of the neutron star in 4u0115 + 63 .
note that @xmath59 for most of the orbital phase .
therefore , we expect that an accretion disc is formed around the neutron star in 4u0115 + 63 .    next , we discuss whether the accretion disc formed around the neutron star is persistent or transient .
this is done by comparing the viscous accretion time - scale with the orbital period . for simplicity
, we assume the accretion disc to be geometrically thin and isothermal with the shakura - sunyaev viscosity parameter @xmath17 and evaluate the accretion time - scale at @xmath60 .
the ratio of these two time - scales for 4u0115 + 63 is then given by @xmath61^{1/2 } \left (    \frac{t_{\rm{d}}}{10^{4}\rm{k } } \right)^{-1 } ,    \label{eq : tratio}\end{aligned}\ ] ] where @xmath62 is the sound speed and @xmath63 is the disc temprature .
the orbital - phase dependence of @xmath64 is shown in the lower panel of fig .
[ fig : rcirc ] .
it is immediately observed that the viscous accretion time - scale is much longer than the orbital period for most of the orbital phase .
therefore , once formed , the accretion disc around the neutron star in 4u0115 + 63 will survive over the whole orbital phase .
note that this result depends little on the orbital period , because @xmath65 increases almost linearly with increasing @xmath9 .
based on the procedure described in the previous section , we have performed several simulations with different parameters .
details of the models are given in table  [ tbl : models ] .
the polytropic exponents considered are @xmath66 ( model  3 ) , 1.2 ( models  1 , 2 , 5 and 6 ) and 5/3 ( model  4 ) .
we have taken the inner radius of the simulation region , @xmath67 , to be @xmath68 in all simulations but model  1 with a smaller inner radius of @xmath69 , which was run to study the effect of the inner radius on the simulation results . in models  1 - 4 and 6 ,
we have assumed @xmath43 , where @xmath18 varies in time and space and @xmath27 . for comparison purpose
, we have also run model  5 with @xmath23 and @xmath24 , where @xmath17 varies in time and space .
model  6 was run solely for clearly illustrating the initial build - up phase of an accretion disc . for this purpose
, we have taken the injection rate of the sph particles in model  6 four times as high as in other models and run only for one orbital period .
analysing multi - wavelength long term monitoring observations of 4u0115 + 63 , @xcite found that the be star undergoes quasi - cyclic @xmath70 activity , losing and reforming its circumsteller disc .
then , the mass transfer to the neutron star should vary , depending on the state of the be disc . by the time the be disc is fully developed and becomes capable of supplying mass to the neutron star , even a large accretion disc would be exhausted .
therefore , to investigate the evolution of an accretion flow around the neutron star after the be disc is fully developed , we performed 3d sph simulations based on the numerical models described above , where there is intially no gas particle around the neutron star .
since all models have shown similar structure and evolution of the accretion flow , we will mainly show the results from model  1 with @xmath71 and @xmath72 .
+     +    we consider the initial phase of accretion flow around the neutron star . as discussed in section  [ sec : possibility ] , it is expected that the material transferred from the be disc near the periastron does not directly fall on to the neutron star and forms a disc - like structure with a size corresponding to its specific angular momentum .    for illustration purpose ,
we give in fig .
[ fig : model6_ss ] the snapshots of the initial accretion flow around the neutron star from model  6 , which covers @xmath73 . at @xmath74 ,
the be star is at periastron .
each panel shows the surface density in the range of four orders of magnitude in the logarithmic scale .
we note that the accreting matter first forms an elliptical ring around the neutron star , which is then relaxed to an elliptical disc via the viscous diffusion .
because the viscous time - scale is longer than the orbital period , the orbits of the accreting particles are not circularized and the flow is highly eccentric during the initial phase of accretion .
comparing the number of particles , @xmath75 , with the number of particles inside the effective roche radius of the neutron star , @xmath76 , we also note that most of the material transferred from the be disc stays inside the effective roche radius of the neutron star .
it should be noted that these dynamical features are common to all models studied in this paper .
[ fig : model6_rad ] shows the time - dependence of the radial structure of the accretion flow around the neutron star at the same times as in fig .
[ fig : model6_ss ] . in each panel ,
the solid line , the dashed line and the dash - dotted line denote the surface density @xmath77 , the azimuthal velocity @xmath78 normalized by the keplerian velocity at the inner boundary and the radial mach number @xmath79 , respectively . in the figure ,
the density is integrated vertically and averaged azimuthally , while the velocity components are averaged vertically and azimuthally . from fig .
[ fig : model6_rad ] we note that while the azimuthal velocity distribution gradually approaches the keplerian one , the behaviour of the radial velocity is very complex , reflecting the dynamical process of the disc formation shown in fig .
[ fig : model6_ss ] .
the surface density distribution also shows the initial evolution of the accretion flow described above . in the panel at @xmath80
, there appears a hump at @xmath81 caused by the infall of injected particles .
outside the hump the density decreases rapidly , making a ring - like structure with a sharp outer edge .
this rapid decrease in the density is due to the ram pressure of the infalling particles .
as the hump moves inward , the density gradient outside the hump becomes less steep . by the next periastron passage ( @xmath82 ) ,
a disc - like structure is formed , the inner radius of which is at @xmath83 .    in our simulations ,
the presence of the inner simulation boundary causes an artificial decrease in the density near the boundary , because the gas particles which pass through the boundary are removed from the simulation . as a result
, the density distribution has a break near the inner boundary , as seen in fig .
[ fig : model6_rad ] . we have found that our simulations are reliable for @xmath84 .
+     +      many be / x - ray binaries with known orbital parameters have orbital eccentricities in the range of 0.3 to 0.9 .
in such systems , the mass - transfer rate from the be disc is expected to have the strong phase dependence @xcite .
accordingly , the structure of an accretion disc formed around the neutron star in these systems is also likely to be strongly phase dependent . in this section ,
we first discuss the phase - dependent structure of a developed accretion disc in model  1 and then compare that with those in other models .      fig .
[ fig : model1_ss ] gives a sequence of snapshots of the accretion flow around the neutron star covering one orbital period for @xmath85 .
the format of the figure is the same as that of fig .
[ fig : model6_ss ] , except that the surface density is shown in the range of three orders of magnitude .
we note from the figure that the disc shrinks at periastron passage by a negative torque exerted by the be star and restores its radius afterwards by the viscous diffusion .
we also note that the disc is significantly eccentric even after it is fully developed , although the eccentricity is much smaller than that in the initial phase of the disc formation .
[ fig : model1_rad ] shows the radial disc structure at the orbital phases corresponding to fig .
[ fig : model1_ss ] .
the format of the figure is the same as that of fig .
is noted from fig .
[ fig : model1_rad ] that the disc is nearly keplerian at any phase and the accretion flow is highly subsonic except in the outermost part for a short period of time when the material transferred from the be disc falls towards the neutron star .
we call this state of the accretion disc developed .
we have found that these features of a developed disc hold for @xmath86 in model  1 .
as in the initial phase of the disc formation , the ram pressure by the supersonic infall of matter transferred from the be disc around periastron enhances the density in the outermost part of the accretion disc , making the disc outer edge sharp ( see panels for @xmath87 in figs .  [
fig : model1_ss ] and [ fig : model1_rad ] ) .
afterwards , the disc gradually expands until the next periastron passage of the be star .
[ fig : chardisc ] shows how hot the disc is at @xmath88 in models  1 and 3 - 5 . in each panel , the solid line , the dash - dotted line and the dashed line denote the disc temprature normalized by the initial temprature @xmath89 , the ratio of the smoothing length to the harf - thickness of the disc @xmath90 and the relative disc thickness @xmath91 , respectively . from the upper - left panel of fig .
[ fig : chardisc ] , it is noted that the disc in model  1 is vertically resolved except for the innermost region of the disc ( @xmath92 ) .
the disc temprature varies from @xmath93 at @xmath94 to @xmath95 at the inner boundary .
the disc is geometrically thin with the relative thickness of about 0.1 - 0.2 .      in this section
, we describe the effects of the simulation parameters in table  [ tbl : models ] on the structure of the accretion flow around the neutron star .
we ran model  2 to study the effect of the inner boundary on the disc strcture .
[ fig : model2 ] shows the snapshots and the radial structures of the accretion flow around the neutron star at @xmath96 ( periastron ) and @xmath97 ( apastron ) for model  2 , which has the same simulation parameters as model  1 , except for a larger inner boundary of @xmath98 . comparing fig .
[ fig : model2 ] with the panels at the corresponding times in figs .
[ fig : model1_ss ] and [ fig : model1_rad ] for model  1 , we note that the presence of the inner boundary artificially reduces the density in a narrow region ( @xmath99 ) near the boundary , but outside of it , the disc structure is similar for both models .    in order to study the dependence of the disc structure on the polytropic exponent , we also ran models  3 ( @xmath100 ; isothermal case ) and 4 ( @xmath101 ; adiabatic case ) , for which the polytropic exponent of 1.2 for models  1 , 2 and 6 is in between .
[ fig : model3 ] shows the snapshots and the radial structures of the accretion disc for model  3 ( isothermal case ) at the same times as in fig .
[ fig : model2 ] . from the figure
, we note that the distributions of the surface density and both of the azimuthal and radial velocity components are similar to those in the initial phase of the disc formation for @xmath102 ( see figs .
[ fig : model6_ss ] and [ fig : model6_rad ] ) . consequently , in model  3 the accretion disc is still developing at @xmath103 .
this is because the viscous time - scale of the isothermal disc , particularly in an inner region , is longer than those in other models ( see also fig .
[ fig : rcirc ] ) . from the upper - right panel of fig .
[ fig : chardisc ] , it is noted that the disc in model  3 is not vertically resolved in most of the disc region .
[ fig : model4 ] shows the disc structure for model  4 ( adiabatic case ) .
we note that the disc is nearly keplerian at both orbital phases and that the surface density is lower and the disc is bigger than those for model  2 with @xmath102 .
this is due to a higher viscous stress in model  4 , resulting from a higher pressure .
as expected , in model  4 , the accretion rate is much higher than the models with lower polytropic exponents , which we will see in a later section . from the lower - left panel of fig .
[ fig : chardisc ] , it is noted that the disc in model  4 is much hotter than in other models , reaching the temprature of @xmath104 near the inner boundary .
the disc is vertically resolved except for the innermost region .
we also studied a model with the standard , artificial viscosity parameters .
model  5 has the same simulation parameters as those in model  2 , except that it has constant articficial viscosity parameters @xmath105 and @xmath106 , for which @xmath17 varies in time and space .
[ fig : model5 ] shows the disc structure for model  5 . from the figure and the lower - right panel of fig .
[ fig : chardisc ] , it is immediately noted that the inner boundary affets the disc structure much more widely in this model .
this is mainly because the inner disc region is not resolved vertically for this number of sph particles .
since the ratio of the smoothing length to the harf - thickness of the disc @xmath90 in the inner region is much higher than unity , the shear viscosity parameter @xmath17 in this model is much higher in the inner region than that in other models .    except for the above - mentioned dependence of the disc structure on the simulation parameters , the accretion discs in our models share the common features :
they are nearly keplerian after developed and their structures modulate with the orbital phase .
they are also significantly eccentric .
we will analyse the non - axisymmentry of the accretion discs in be / x - ray binaries in a subsequent paper .
+     +     +     +     +      @xcite investigated the tidal / resonant truncation of circumsteller and circumbinary discs in eccentric binaries and found that a gap is always opened between the disc and the binary orbit .
following their formulation , @xcite and @xcite have shown that the be disc is truncated via the resonant interaction with the neutron star , as long as @xmath107 , which has recently been confirmed by numerical simulations by @xcite .
it is interesting to study whether the tidal / resonant torque by the be star also truncates the accretion disc around the neutron star . in order to have a measure of the accretion disc radius
at which the disc density has a major break , we have applied a non - linear least squares fitting method to the radial distribution of the azimuthally - averaged surface density @xmath77 , adopting the following simple fitting function , @xmath108 where @xmath109 and @xmath110 are constants and @xmath111 is the radius of accretion disc .
this method is the same as that adopted by @xcite .
[ fig : radius ] shows the orbital - phase dependence of the accretion disc radius in model  1 .
to reduce the fluctuation noise , we folded the data on the orbital period over @xmath112 .
it is noted that the disc radius varies between @xmath113 and @xmath114 .
the disc shirinks at periastron due to a negative torque by the be star and gradually restores its radius towards apastron due to the viscous diffusion .
the rapid increase of the disc radius before periastron is due to the infall of matter from the be star .
next , we examine whether the disc radius shown in fig .
13 is determined by the resonant truncation as is the case in the be disc in be / x - ray binaries .
fig .  8 of @xcite gives the truncation radius of a circum - secondary disc as a function of the orbital eccentricity @xmath115 and the disc reynolds number @xmath116 for systems with the mass ratio of 0.1 , which is similar to that of be / x - ray binaries .
since the accretion disc in our model is geometrically thin and nearly keplerian , the reynolds number @xmath116 is given by @xmath117 @xmath118 , where @xmath119 is the keplerian velocity . in model  1 ,
@xmath120 for @xmath121 . from fig .  8 of @xcite , with @xmath122 and @xmath123
, we have the truncation radius of @xmath124 , which corresponds to the 6:1 resonance radius . note
that this radius is much larger than the disc radius shown in fig .
[ fig : radius ] .
it is likely that the interaction with the be star is significantly non - linear owing to its large mass .
then , the accretion disc would be truncated well inside the 6:1 radius given above .
even such a strong interaction , however , is unlikely to truncate the disc at a radius as small as that shown in fig .
[ fig : radius ] .
it is suggestive that the disc radius is about the same as the circularization radius @xmath65 of the matter transferred from the be disc after the periastron passage ( see fig .
[ fig : rcirc ] ) .
since the viscous time - scale at the circularization radius is much longer than the run time of our simulations , the matter once circularized dose not have time enough to change its radius significantly .
therefore , we conclude that the disc radius shown in fig .  [
fig : radius ] is determined by the specific angular momentum of the matter transferred from the be disc .
it is probable , however , that the disc expands well beyond the circularization radius and is truncated by a tidal torque in a time - scale longer than the viscous one .
as mentioned in section  1 , most of the be / x - ray binaries show transient x - ray activities , of which nature is considered to result from the interactions between the accreted material and the rotating magnetized neutron star @xcite . in this section ,
we first describe the growth of the accretion disc in mass and radius .
we then examine the evolution of the mass - accretion rate obtained by our simulations , from which we finally discuss the x - ray activity expected in be / x - ray binaries .
[ fig : masrad](a ) shows the evolution of the disk mass in model  1 over the period of @xmath125 . in the long term , the disc mass increases with time , because the the gas transferred from the be star accumulates in the accretion disc owing to the viscous time - scale longer than the orbital period . on the other hand , in the short term
, it modulates with the orbital phase .
the disk mass starts increasing , when the mass - supply from the be disc begins slightly before every periastron passage .
it shows a rapid increase for a short while , followed by a gradual decrease by accretion with no mass - supply from the be disc , which lasts until the next periastron passage . at first glance
, the disc mass in model  1 seems to increase secularly and approach no steady average state . on closer inspection
, however , we note that the increasing rate of the disc mass gradually decrases with time after the disc is developed .
this suggests an interesting possibility that the disc reachs a steady average state in a time - scale much longer than the term shown in fig .
14(a ) .    in fig .
[ fig : masrad](b ) , we show the evolution of the disc radius over @xmath125 .
one may expect that the disc radius gradually increases with the disc mass . from the figure
, however , we note that this is the case only for @xmath86 .
in fact , for @xmath126 , the disc radius decreases as the disc grows .
this counter - intuitive decrease in radius is related to the change in the disc eccentricity . in the early stage of disc formation , the disc / ring is highly eccentric ( see fig .  [
fig : model6_ss ] ) , for which the radius @xmath111 calculated by equation  ( [ eq : fitting ] ) is larger than that for a circular disc with the same disc mass . since the disc eccentricity decreases with time ,
so does the disk radius . in model  1 ,
the effect of the disc circularization is balanced with that of the disc growth at @xmath127 .
after this epoch , the disc radius increases with the disc mass .
thus , the evolution of the accretion flow around the neutron star in be / x - ray binaries involves a two - stage process , which consists of the initial developing stage and the later developed stage .
[ fig : mdot_acc](a ) shows the evolution of the mass - accretion rate and the corresponding x - ray luminosity in model  1 for @xmath128 .
here we calculated the x - ray luminosity by @xmath129 with the x - ray emission effciency @xmath49 for the neutron star .
the thick and thin lines denote the mass - accretion rate in models  1 and 2 , respectively .
note that the mass - accretion rate in our simulations generally has double peaks per orbit .
the first peak , which occurs at the periastron , is due to the contraction of the disc by the tidal torque of the be star , except that the peaks at periastron during the initial phase of disc formation are mainly due to the direct accretion of gas with low specific angular momentum .
the second peak , which lags behind the peak in the mass - transfer rate from the be disc shown in fig .
[ fig : mdot_cap ] , results from the viscous accretion of matter transferred from the be disc .
the first peak could be artificial since it is related to the presence of the inner simulation boundary .
the orbital modulation in the mass - accretion rate is significantly different between the developing disc and the developed disc .
the first peak becomes weak as the disc developes .
after the disc is fully developed , the second peak is higher than the first peak .
[ fig : mdot_acc](a ) also shows that , once developed , the accretion disc in 4u0115 + 63 shows the peak x - ray luminosity of @xmath130 .
this level of x - ray emission enters the transition regime ( see @xcite ) and is consistent with the observed x - ray luminosity in the quiescent state .    for comparison purpose , we present in fig .  [ fig : mdot_234 ] the evolution of the mass - accretion rate and the corresponding x - ray luminosity in models  2 - 4 for @xmath125 . the thick and thin solid lines and the thin dotted line are for models  2 , 3 and 4 , respectively .
we see from the figure that the accretion rate increases with increasing polytropic exponent . in models  2 ( @xmath131 ) and 3 ( @xmath100 )
, the accretion rate profile has the double peaks as described above .
the long - term change in the peak profiles is also similar to that in model  1 .
we note that the disc in model  2 entered the stage of the developed disc at @xmath132 . on the other hand ,
the disc in model  3 is at the developing stage even at @xmath133 , which is consistent with the disc structure shown in fig .
[ fig : model3 ] .
in contrast to the models with low @xmath28 , the accretin rate in model  4 has a single and periodic peak per orbit , even after the disc is developed . with @xmath134 ,
the pressure in the disc in model  4 is much higher than that in other models .
this means that much larger viscous torque exerts on the disc in model  4 than in other models , because the viscous stress in the @xmath25-viscosity prescription is proportional to the pressure . as a result
, the viscous torque dominates the tidal torque in model  4 even at periastron .
this is why there is no peak at periastron and the accretion rate shows the regular orbital modulation with the amplitude much higher than in other models .
it is important to note that the disc evolution and its orbital modulation is sensitive to the polytropic exponent @xmath28 .
this strongly implies that the cooling effciency in the accretion disc plays an important role in the x - ray activity in be / x - ray binaries .
changes in the mass flux on the magnetosphere of the neutron star determine whether the system is in the direct accretion regime or in the propeller regime . while the matter falls directly on to the neutron star in the direct accretion regime , no material accretes on to the neutron star in the propeller regime by the centrifugal inhibitation .
if the magnetospheric radius @xmath135 is smaller than the corotation radius @xmath136 , the system is in the direct accretion regime , otherwise , in the propeller regime .    assuming the neutron star magnetic field as the dipole - like field , the magnetospheric radius @xmath135 , obtained by equating the magnetic field pressure of the neutron star to the ram pressure of the accreting matter from the accretion disc , is written by @xmath137
( @xcite ; @xcite ; see also @xcite ) , where @xmath138 is the magnetic moment . here
, @xmath139 is a constant model parameter of order unity , depending on the physics and geometry of accretion . for spherical accretion , @xmath140 .
the original model calculations by @xcite for disc accretion on to the aligned rotating magnetized neutron star obtained @xmath141 , whereas calculations by @xcite , which included an induced wind not considered by @xcite , resulted in the value of @xmath142 for disc accretion .    fig .
[ fig : mdot_acc](b ) shows the evolution of the magnetospheric radius @xmath135 normalized by the corotation radius @xmath136 , where @xmath143 was adopted . in the figure , the thick and thin lines are for models  1 and 2 , respectively .
recall that these models have the same parameters except for the radius of the inner simulation boundary @xmath67 ; @xmath144 in model  1 , while @xmath145 in model  2 .
it should be noted that the magnetospheric radii in both models evolve in a very similar way , despite the different inner radii .
our result on the magnetospheric radius is robust in this sense .
we also note from fig .
[ fig : mdot_acc](b ) that the ratio of the magnetospheric radius to the corotation radius , @xmath146 , varies between @xmath147 and @xmath148 .
consequently , for a plausible range of @xmath139 ( @xmath149 ) , the magnetospheric radius in models with @xmath150 is almost always larger than the corotation radius so that the accretion on to the neutron star is inhibited , which is consistent with the oberved x - ray luminosity of 4u0115 + 63 in quiescence ( @xcite ) . on the other hand
, we can see from equation  ( [ eq : magnetosphere ] ) and fig .
[ fig : mdot_234 ] that model  4 with @xmath101 is consistent with the observed x - ray behaviour of 4u0115 + 63 only if @xmath151 . for @xmath152 ,
the neutron star in model  4 enters the direct accretion regime after every peristron passage of the be star .
then , the system would show regular , periodic x - ray outbursts at a level of @xmath153 , which has not been observed for 4u0115 + 63 .
for the first time , we have performed numerical simulations of accretion on to the neutron star in be / x - ray binaries , taking into account the phase - dependent mass transfer from the be disc .
we have adopted the mass - transfer rate from a high - resolution simulation by @xcite for a coplanar system with a short period ( @xmath31 ) and moderate eccentricity @xmath1 , which targeted 4u0115 + 63 , one of the best studied be / x - ray binaries .
we have found that a nearly keplerian , non - steady accretion disc is formed around the neutron star regardless of the parameters adopted .
the formation of the accretion disc involves a two - stage process , which consists of the initial developing stage and the later developed stage .
the disc structure shows a strong dependence on the orbital phase , because both the tidal potential and the mass - transfer rate change periodically .
the disc is also significantly eccentric , because the gas particles from the be disc originally have elliptical orbits .    in 4u0115 + 63 , type  i x - ray outbursts have occurred only as a short series after occasional type  ii x - ray outbursts @xcite . except for such occasions ,
the system is in the quiescent state ( @xmath154 ) @xcite .
the peak x - ray luminosity estimated from our simulations corresponds to the transition regime between the direct accretion regime and the propeller regime and is consistent with the observed x - ray behaviour .
note , however , that the discussion based on the x - ray luminosity estimated from the simulations only gives the necessary condition for the quiescent state .    in order to see whether our model satisfies a more stringent condition ,
we also have estimated the magnetospheric radius by using a general formula with the mass - accretion rate from our simulations and compared with the corotation radius for 4u0115 + 63 .
the result was that our models with @xmath150 are consistent with the observed x - ray behaviour .
other models with @xmath101 and @xmath155 were ruled out , because in these models the system would show an x - ray outburst after every periastron passage , which has never been observed .
we have to admit , however , that the radius of our inner simulation boundary is about 60 - 100 times as large as the corotation radius of the neutron star in 4u0115 + 63 . to make the constraint on the system parameters more reliable
, we need numerical simulations with a much smaller inner boundary .
it is important to note that our model is not only consistent with the quiescent state of 4u0115 + 63 .
it can also give a natural explanation for a short series of type  i x - ray outbursts after occasional type  ii x - ray outbursts .
the mass - transfer rate from the be disc in the x - ray outbursts is likely much higher than in the quiescent state .
suppose that the mass - transfer rate from the be disc in 4u0115 + 63 is temporarily increased by an order of magnitude .
then , the estimated x - ray luminosity in our model increases by an order of magnitude to a level of several @xmath156 , a typical x - ray luminosity of type  i x - ray outbursts . at the same time , the magnetosphric radius @xmath135 decreases by a factor of about two , for which @xmath135 is almost always smaller than the corotation radius , @xmath136 , if @xmath152 ( see fig .  [
fig : mdot_acc](b ) ) .
this means that the system is in the direct accretion regime .
thus , our model explains the type  i x - ray outbursts in 4u0115 + 63 as a temporal phenomena with an enhanced mass - transfer from the be disc .
be / x - ray binaries are distributed over a wide range of orbital periods @xmath157 and eccentricities @xmath158 .
it is interesting to consider here the effects of orbital parameters . in this paper
, we have seen that the main ingredients which determines the x - ray behaviour of a system are the viscous time - scale relative to the orbital period and the mass - transfer rate from the be disc .
since the former is insensitive to the orbital period itself ( see equation  ( [ eq : tratio ] ) ) , the latter is most important if the viscosity parameter and the equation of state are similar in accretion discs of be / x - ray binaries .    in systems with @xmath159 , the viscous time - scale in an outer part of the accretion disc is much longer than the orbital period unless the cooling is ineffective and the gas obeys an adiabatic equation of state . in these systems , persistent accretion discs are expected to form .
then , the mass - transfer rate from the be disc determines whether the system shows type  i x - ray outbursts .
for example , for a system with a moderate orbital eccentricity of @xmath160 and a typical density of be disc , the peak mass - transfer rate is @xmath161 , for which the system stays in quiescence .
it is obvious , however , that the situation should be different in systems with much higher rate of mass - transfer .
a high mass - transfer rate is expected to result from a high orbital eccentricity and/or an intrinsically dense be disc .
therefore , it is not strange that many systems which have regularly shown type  i x - ray outbursts have orbital eccentricities higher than that of 4u0115 + 63 studied in this paper . in a subsequent paper , we will perform numerical simulations of accretion in systems with high orbital eccentricities .
above we have assumed that the viscous time - scale is much longer than the orbital period in all be / x - ray binaries .
however , it is possible that an accretion disc has a relatively short viscous time - scale by some reason .
one possibility is that the disc has a low cooling efficiency and becomes hot , particularly in an inner region .
our model  4 with @xmath101 gives an extreme example of such discs ( see also fig .  [
fig : chardisc ] ) .
as shown in fig .
[ fig : mdot_234 ] , in a system with a short viscous time - scale , the mass - accretion rate has a high and sharp peak shortly after every periastron passage .
the resemblance between the accretion rate profile in model  4 and the x - ray light curves of regular type  i x - ray outbursts from systems such as exo2030 + 375 @xcite implies that these are systems with relatively short viscous time - scales .
be / x - ray binaries are an ideal group of objects for studying the physics of accretion under the circumstances which have been studied little . unlike close binaries ,
many be / x - ray binaries have eccentric orbits .
they can also be highly inclined systems .
such systems provide us a valuable oppotunity to study the effects of the periodically - changing tidal potential and mass - transfer rate and the inclination angle on the structure and evolution of the accretion flow .
another important characteristic of be / x - ray binaries is that they are systems with double circumstellar discs ( the be disc and the accretion disc ) .
the interaction in such systems is threefold : the interactions between the neutron star and the be disc , the be star and the accretion disc and the be disc and the accretion disc ( mainly via the mass transfer ) .
much more work is desirable both theoretically and observationally in order to understand this interesting and important group of objects .
we acknowlege the anonymous refree for constructive comments .
we would like to thank matthew r.  bate for allowing us to use his sph code .
kh acknowledes james r. murray for useful comments during the short stay at the swinburne university of technology , australia .
kh also thanks noboru kaneko for his continuous encouragement .
the simulations reported here were performed using the facility at the hokkaido university information initiative center , japan .
this work was supported in part by grant - in - aid for scientific reserch ( 13640244 ) of japan society for the promotion of science .
artymowicz p & lubow s.h . , 1994 ,
apj , 421 , 651 bate m.r . , bonnell i.a.,price n.m . , 1995 , mnras , 285 , 33 benz w. , 1990 , in buchler j. r.,ed.,the numerical modelling of nonlinear stellar pulsations .
kluwer , dordrecht , p.269 benz w. , bowers r.l . ,
cameron a.g.w . ,
press w.h . , 1990 ,
apj , 348 , 647 bildsten , l. , chakrabarty , d. , chiu , j. , et al .
1997 , apjs , 113 , 367 campana s. , gastaldello f. , stella l. , israel g.l . , colpi m. , pizzolato f. , orlandini m. , dal fiume d. , 2001 , apj , 561 , 924 clark j.s . , lyuty v.m .
, zaitseva g.v . , larionov v.m . ,
larionova l.v . , finger m. , tarasov a.e . , roche p. , coe m.j
. , 1999 , mnras , 302 , 167 finger m.h . , wilson r.b .
, harmon b.a . , 1996 , apj , 459 , 288 finger m.h . , bildsten l. , chakrabarty d. , prince t.a . , scott d.m . ,
wilson c.a . , wilson r.b & zhang s.n . , 1999 , apj , 517 , 449 frank j. , king a.r . ,
raine d.j . , 2002 ,
accretion power in astrophysics , 3rd edn .
cambridge univ . press ,
cambridge ghosh p & lamb f.k .
, 1979 , apj , 234 , 296 lamb f.k .
, pethick c.j & pines d. , 1973 , apj , 184 , 271 negueruela i. , okazaki a.t .
, fabregat j. , coe m.j .
, munari u & tomov . , 2001 , a&a , 369 , 108 negueruela i & okazaki a.t . , 2001 , a&a , 369 , 108 okazaki a.t & negueruela i. , 2001 , a&a , 377 , 161 okazaki a.t . , bate m.r .
, ogilvie g.i & pringle j.e . , 2002 , mnras , 337 , 967 ostriker e.c . , shu f.h . , 1995 , 447 , 813 pringle j.e & rees m.j . , 1972 ,
a&a , 21 , 1 stella l. , white n.e . , rosner r. , 1986 , apj , 208 , 669 warner b. , 1995 , cataclysmic variable stars . cambridge univ . press , cambridge wilson c.a & finger m.h .
, coe m.j . , laycock s. , fabregat j. , 2002 , apj , 570 , 287 wilson c.a & finger m.h . , coe m.j .
, negueruela i. , 2003 , apj , 584 , 996 ziolkowski j. , 2002 , mem .
ital , 73 , 1038