the accretion processes are recognized to power many astrophysical sources , which are observable due to the release of energy when matter is spiralling down into the gravity well .
the study of viscous accretion disks as underlying mechaisms for the extraction of potential gravitational energy began with the works of  @xcite , @xcite .
time - dependent problems for the accretion disks arise when describing a wide variety of outburst phenomena , observed in the accreting sources , i.e. binary systems with mass transfer and galactic nuclei , and systems with planet formation .
it is usually possible to consider separately the vertical and radial structure of an accretion disk due to significantly different characteristic time - scales .
in such a case the time - dependent radial structure of a disk can be described by a second - order partial differential equation , which is a consequence of the angular momentum and mass conservation .
the energy conservation in a geometrically thin optically thick disk is provided by the local balance of the viscous heating and radiative cooling .
the basic equation of the viscous evolution , eq .   or   below , describes the viscous angular momentum flux along the disk radius .
this is a diffusion type equation with a variable ` diffusion ' coefficient .
choice of an analytic approach to solve the equation depends on the nature of the viscosity involved . the preferential way of describing viscosity in astrophysical disks , although not exclusive , is to use  @xcite , i.e. , the proportionality of the viscous stress tensor to the total pressure in the disk .
one can generally assume that the kinematic viscosity @xmath4 is a product of two power functions : one function of radius and another of some local physical parameter , the latter usually represented by the surface density .
the boundary conditions imposed on the accretion disk are important .
while the outer boundary can be a freely expanding surface due to the spreading of the matter corresponding to the outward viscous angular momentum flux , in some situations certain conditions should be posed at some distance from the accreting object .
disks in the binary systems are the main focus for such models .    in a binary system , the tidal torques ,
acting inside the roche lobe of the accreting object , truncate the disk and provide a sink for the angular momentum .
@xcite show that the effective tidal radius is close to the roche lobe : @xmath5 of the mean radius of the roche lobe , well coinciding with the largest non - intersecting periodic orbits in the restricted three - body problem @xcite .
numerical models confirm that most of the tidal torque is applied in a narrow region at the edge of the disk , where perturbations become nonlinear and strong spiral shocks appear  ( see * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) . in the present study ,
we leave aside the outer disk radius variations and other manifestations of the tidal interactions in the binary except for the disk truncation .
we confine ourselves to searching for a solution to a problem of the viscous evolution of an accretion disk with the constant outer radius .
we investigate a geometrically thin , axisymmetric keplerian finite accretion disk with a viscosity in the form @xmath6 , which is time - independent .
we obtain green s functions ( [ eq : green_zerobi ] ) , ( [ eq : mdot_green_zerobi ] ) , ( [ eq : green_zeroai ] ) , and ( [ eq : mdot_green_zeroai ] ) , which are the kernels in the integrals ( [ eq : f_from_green ] ) , ( [ eq : mdot_from_green ] ) , and ( [ eq : sol_nonzeromdot_full ] ) to compute the viscous angular momentum flux of  @xcite . ] @xmath7 and the accretion rate @xmath8 .
the method of green s functions allow one to compute @xmath9 and @xmath10 for arbitrary initial surface density distributions .
different boundary conditions are considered : no stress or no accretion at the inner boundary , which is located at the zero radial coorinate , zero or time - dependent mass inflow at the outer boundary .
it is beleived that the mechanism of bursts in x - ray novae involves the thermal instability connected with opacity variations in the disks , as in dwarf novae . when a burst is ignited , and a sufficient portion of the disk is in the hot state , the viscous diffusion starts to govern the evolution , driving the characteristic rise and drop of the accretion rate onto the central object .
it is a question , which probably allow no universal answer , why x - ray novae light curves have different profiles and pecularities .
if the disk is not very large , as it happens in short - period x - ray transients , its whole body is likely to be engaged in the viscous evolution .
it has been shown that an episodic mass input at the outer radius causes the fast - rise - exponential - decay ( fred ) profile due to the viscous evolution  @xcite . using the green s function
, one can explore the situation with an arbitrary initial distribution of the matter . actually , as we show in the present study after obtaining the green function , for realistic initial density distributions and for the case with the zero inner viscous stress , the evolution is generally a fred .
it is interesting to compare the results of the model with what happens in the @xmath3-disks , the prevalent model for the disks turbulence in the stellar binaries .
the second case , when the viscous stress is non - zero at the center , applies to the disks with the central source of the angular momentum , e.g. , rotating stars with magnetosphers .
if the magnetosphere s radius is large enough , the matter has a super - keplerian angular momentum at its boundary .
different regimes for such disks have been proposed : propellers and ` dead ' disks .
using the green founction , the disks of constant mass that evolve to dead disks can be described .    in binary systems
the mass transfer rate is a key factor . at some stages
it can be considered negligible comparing to the accretion rate in the disk .
in other systems , mass transfer variations may cause the activity in the disk .
the case of time - dependent mass transfer , which includes the constant transfer rate as a particular case , can be also studied analytically if @xmath11 . using the green s functions ,
it is possible to express @xmath9 and @xmath10 for any outer boundary conditions ( for any @xmath12 ) as well as for any initial conditions .    in section 2
, we present the viscous evolution equation and review its known solutions for accretion disks . in section  [ sec : ev_fin_disk ] , the green s functions are found for different inner and outer boundary conditions .
an analytic way to calculate the disk evolution with the variable mass transfer at the outer radius is presented ; the features of a radiating dead disk are considered . in section  [ s.decay ]
, a stage of decaying accretion is dealt with and a comparison to the @xmath3-viscosity models is made .
application to the fred light curves of short - period x - ray transients is presented in section  [ s.xnovae ] .
the last sections are dedicated to the discussion and summary .
the equation of the viscous evolution in the accretion disk is the diffusion - type equation ( e.g. , * ? ? ?
* ) : @xmath13\ ,   , \label{eq : time - dep - acc1}\ ] ] where , among the standard designations of the time , radius , angular speed , and the surface density we have the component of the viscous stress tensor @xmath14 , integrated over the full thickness of the disk @xmath15 where @xmath16 is the half - thickness of the disk at radius @xmath17 , and the stress tensor is related to the differential rotation of gas masses by the kinematic viscosity @xmath4 as follows : @xmath18 for the keplerian disks @xmath19 , and we rewrite @xmath20 with the formula for the surface density @xmath21 assuming that @xmath4 is independent on @xmath22 , we have @xmath23    it is convenient to introduce a new independent variable  the specific angular momentum . then , for the keplerian disk , equation  ( [ eq : time - dep - acc1 ] ) can be rewritten as : @xmath24    the acknowledged approach is to solve the disk evolution equation for the variable proportional to the full viscous torque acting between the adjacent rings in the disk , @xmath25 .
this way , the boundary conditions are easy to express . from the above expressions , the following relation between @xmath26 and the accretion rate can be obtained : @xmath27 where the radial velocity @xmath28 is negative .
one chooses a procedure to solve ( [ eq : basic_tda_kepler ] ) provided the form of @xmath29 is known .
this form is determined by the physical conditions in the disk .
for example , the vertical structure can be solved for the @xmath3-disk , and the required relation between @xmath4 and @xmath30 can be found .
let us give a perspective of some studies ( mainly analytic ) of the viscously evolving accretion disks .
they differ in respect of the viscosity prescription and boundary conditions .
the analytic solutions are the benchmarks for more sophisticated and involved numerical models .      for the kinematic viscosity in the form @xmath6 ,
the equation of the viscous evolution   becomes a linear partial differential equation . in @xcite @xcite
obtained particular solutions to the equation proposed by his teacher von  weizscker  ( @xcite ) and described the principles to determine the general solution for the infinite and finite problems .
for the accretion disks that can extend _ infinitely _ , @xcite ( hereafter lp74 ) find , for the two types of the inner boundary conditions , the green s functions that provide means to describe the whole development of the disk .
the inner radius of the disk is zero in their exact solution .
the long - term self - similar evolution of the accretion rate goes as a power - law @xmath31 with a parameter @xmath32 ( see also appendix a ) .
@xcite considers the properties of an infinite accretion disk with a central source of the angular momentum using the green - function technique .
a related problem is solved by @xcite who finds the green s function for an infinite disk with the boundary conditions posed at the finite inner radius .
@xcite consider a _
finite _ disk with the viscosity constant over the radius and in time and derive the exponential decay of the accretion rate in the disk .
@xcite study the mass flow rate through a disk resulting from a varying mass - supply rate , derive the green s function for the accretion - rate and its fourier transform . by means of numerical simulations
, they investigate a _ finite _ disk , residing between non - zero @xmath33 and @xmath34 , with the non - zero accretion rate at the outer radius .
@xcite find general green s function for an explicit dependence of the inner boundary @xmath33 on time and a non - zero mass across @xmath33 .    in the present work ,
we study analytically a finite accretion disk with the zero @xmath33 and find an analytic solution describing its whole evolution , from the rise to decay . for the case of the zero stress tensor in the disk center and the mass deposition at the outer disk edge
, the green s function of the problem is found in @xcite .
similarity solutions of the first kind of the non - linear differential equation  @xcite have been found at the `` decaying - accretion stage , '' when the total angular momentum of an _ infinite _ disk is constant  @xcite .
accretion rate is found to decrease as @xmath36 for the thomson opacity , and as @xmath37 for the kramers opacity  ( see also * ? ? ?
* ; * ? ? ?
@xcite obtain self - similar solutions that describe separately three stages of the disk evolution .
the first two are the self - similar solutions of the second kind : the initial movement of the inner edge of a ring of matter toward the center and the rise of the accretion rate in the center ( @xmath38 for the kramers opacity and @xmath39 for the thomson opacity ) .
the third stage is the decay of the accretion and corresponds to the solution of @xcite .
the first two stages of @xcite may be applicable to the _ finite _ disks when the conditions at the distant outer boundary do not affect the behaviour of the disk near the center .
@xcite derive an analytic self - similar solution for a disk subject to the gravitational instability , which transfers mass and angular momentum in the disk , for which @xmath40 , and a numerical solution for the whole evolution .
@xcite found a self - similar solution for @xmath41 , if the viscosity is due to the action of convectively driven turbulent viscous stresses ( @xmath42 ) .
accretion disks with a central source of angular momentum and no mass transfer through the inner boundary were considered by @xcite who derived a self - similar solution ( see also @xcite ) .
@xcite in the extensive study of the circumbinary disks around supermassive black holes obtained self - similar solutions for the disks with a possible mass transfer across the orbit of the secondary black hole .    for a _ finite _ disk , the radial and temporal solution
is obtained by @xcite yielding @xmath43 for the kramers opacity and @xmath44 for the thomson opacity .    in the context of the self - similarity solutions ,
we would like to note that for an advection - dominated infinitely expanding accretion flow , @xcite obtains by similarity methods a time - dependent solution with the conserved total angular momentum .
if the kinematic viscosity @xmath4 is not a function of @xmath30 and is a function of the radius alone , we arrive at the linear differential equation  .
we write the kinematic viscosity as @xmath45 then eq .  (
[ eq : basic_tda_kepler ] ) for the dependent variable @xmath7 is as follows @xmath46 where @xmath47 we rewrite in the form similar to that in lp74 : @xmath48 where the constants are defined as follows : @xmath49 there is a degeneracy of the index @xmath50 at @xmath51 .
but it is only apparent as @xmath50 enters equation ( [ eq : diff_lyn - pringle ] ) in the denominators .
value @xmath51 corresponds to the change of the class of functions , which satisfy the differential equation , and is left out of the following consideration . in the case of @xmath51 the viscous time is constant over the disk radii @xmath52 .
we limit ourselves to the finite positive values of @xmath50 and the keplerian disks .    in the astrophysical hot accretion disks ,
the turbulent viscosity is at work . according to the prandtl concept  @xcite
, the kinematic turbulent viscosity is an averaged product of the two random values : @xmath53 , where @xmath54 and @xmath55 are the path and speed of the turbulent motion . for the component of the viscous stress tensor
we write @xmath56 .
following the prandtl consideration , , which leads to .
@xcite proposed the @xmath3-disks : @xmath57 . from the hydrostatic balance one
derives the sound speed @xmath58 and obtains for a keplerian disk @xmath59 .
if the relative half - thickness of the @xmath3-disk is invariable then @xmath60 and @xmath61 .
this value of @xmath62 takes place for the ` @xmath63-viscosity disks ' with @xmath64 , suggested by @xcite to act in the protoplanetary and galactic disks .
the standard model of the @xmath3-disk implies that @xmath65 varies with @xmath17 and @xmath30 .
let us find the green s function of the linear equation  ( [ eq : diff_lyn - pringle ] ) with the specific boundary conditions at the boundaries .
a green s function or a diffusion kernel is a reaction of the system to the delta impulse , in other words , the green s function is the solution to the differential equation with the initial condition @xmath66 , where @xmath67 is the dirac delta function .
we search a solution by separating the variables , , where @xmath68 is a constant and @xmath69 is the outer disk radius . substituting this into ( [ eq : diff_lyn - pringle ] ) , we get @xmath70 where @xmath71 .
this is a lommel s transformation of the bessel equation , which can be derived by introducing a new independent variable @xmath72 .
a particular solution is represented as @xmath73\,\ ] ] @xcite , where @xmath74 and @xmath75 are the bessel functions of a non - integer order .
if @xmath50 is an integer number , the bessel functions of the second kind should be used instead of @xmath75 . for @xmath51 ,
the solution is no more a bessel function  ( see appendix i of * ? ? ?
the value @xmath51 marks a situation when the characteristic viscosity time @xmath76 is constant over the radius . we do not consider this case in the present study .
a general solution for an infinite problem is a superposition of particular solutions and is expressed by lp74 as an integral over all positive values of @xmath77 : @xmath78 for a problem inside a finite interval , a general solution is a superposition of particular solutions and is expressed as an infinite sum  @xcite : @xmath79\ , , \label{eq : gen_sol_finite}\ ] ] where constants @xmath80 , @xmath81 .
and @xmath82 are to be defined from two boundary conditions and one initial condition .
the dimensionless coordinate @xmath83 lies in the interval @xmath84 $ ] , while the physical coordinate @xmath85 ( specific angular momentum ) lies inside @xmath86 $ ] .
the outer boundary condition , which sets the accretion rate at the outer radius of the disk , can be written as follows @xmath87 in the simplest case , the homogeneous outer - boundary condition corresponds to the zero mass - inflow rate at the outer boundary . as this happens , the viscous angular momentum flux @xmath7 takes some non - zero value , which implies withdrawal of the angular momentum at the outer boundary .
this corresponds to the situation of the disk in the binary system , when the angular momentum is removed from the outer edge of the disk by the tidal action of the secondary .
at , one trivial inner boundary condition approximates the structure for the so - called ` standard ' disk , for which the absence of the viscous stresses at the disk inner edge is assumed .
another simple case , @xmath88 at @xmath89 , describes the situation of no accretion through the inner disk edge .
this case is applicable for the stars with strong magnetosphere , when the alfvn radius is greater than the corotation radius , impeding the accretion of the matter .
we consider these two cases separately .
first , we consider the case of the zero viscous stress at the inner boundary .
approximately , this corresponds to the solution @xmath90 at @xmath89 and only the coefficients at @xmath74-terms in ( [ eq : gen_sol_finite ] ) remain non - zero .
thus , sum ( [ eq : gen_sol_finite ] ) at @xmath91 becomes @xmath92    we assume the zero accretion rate at @xmath34 .
such boundary condition approximates the situation of the low mass - transfer from the secondary in the binary system , for example , when the accretion rate in the disk is much higher than the mass transfer rate during an outburst . dropping the time part of the solution
, we arrive at condition ( [ eq : zero_dotm_out ] ) required for every term in the sum ( [ eq : gen_sol_finite_zerobi_t0 ] ) , which can be rewritten as follows : @xmath93    thus , the general solution for the case of zero viscous stress tensor at the center and zero accretion rate at the outer boundary is the sum ( [ eq : gen_sol_finite ] ) , where @xmath94 and @xmath80 are the positive roots of the transcendental equation ( [ eq : out_boundary ] ) . to define @xmath81 from the initial condition ,
the finite hankel transforms are used .
the series @xmath95 with condition ( [ eq : out_boundary ] ) are known as the dini s series ( see * ? ? ?
function @xmath96 can be expressed as the dini s expansion if the function satisfies the dirichlet conditions in the interval and the coefficients are defined as @xmath97 @xcite with the finite hankel transform @xmath98 where @xmath99 .
we look for a solution for the @xmath100-function as the initial condition : @xmath101 . using the properties of the dirac function , we obtain @xmath102    thus , for the specific boundary conditions
, we derive the green s function , which is the solution to ( [ eq : diff_lyn - pringle ] ) with @xmath100-function as the initial distribution @xmath103 where @xmath80 are the positive roots of eq .
( [ eq : out_boundary ] ) and @xmath104 . the function is plotted in fig .  [ fig.green_finite ] for consecutive times .
the curve for time @xmath105 corresponds to the maximum accretion rate in the center .
it is plotted at @xmath106 given by below .    , @xmath107 , @xmath108 , @xmath109 , @xmath110 .
initial spike was at position @xmath111 .
constants @xmath112 , @xmath61 .
[ fig.green_finite ] , scaledwidth=45.0% ]         given the dimensional initial distribution @xmath113 in the finite accretion disk , the distribution of @xmath7 at any @xmath114 is @xmath115 the accretion rate can be found from @xmath116 where the green s function for the accretion rate @xmath117 functions @xmath118 and @xmath119 are found by  @xcite for the case @xmath120 along with an analytical form obtained from the asymptotic forms for small and large @xmath121 .
one can express the initial @xmath7 , proportional to the disk mass , from the initial distribution of the surface density , using ( [ eq : stress_kepl_av ] ) and ( [ eq : consts_lb ] ) : @xmath122 where @xmath123 and @xmath124 .
let the initial mass be located at some radius @xmath125 where the keplerian specific angular momentum is @xmath126 , and the mass of the infinitely thin ring is @xmath127 , @xmath69 is the maximum possible specific angular momentum of the matter .
the surface density of the @xmath100-ring is @xmath128 we take into account that @xmath129 and @xmath130 .
applying ( [ eq : f_sigma_l ] ) , we obtain @xmath131 and the accretion rate evolution is as follows @xmath132 where we designate @xmath133    the time @xmath134 of the maximum and its accretion rate in the disk center @xmath135 are of order of those for the infinite disks . using the solution of lp74
, we find the time of the accretion rate peak for the infinite disk @xmath136 and value @xmath137  .
the numerical factor in is presented in table  [ tab : zeroes ] for various , physically motivated values of @xmath50 .        in fig .
[ fig.cmax1 ] the ratios of maximum times and accretion rates for the finite and infinite disk are shown for several cases of @xmath50 .
the ratios approach @xmath138 for @xmath139  the disk ` does not feel ' the boundary at the time when the accretion rate peaks .     in the disk , whose inner accretion rate evolution is shown in fig .
[ fig.mdotin ] .
top - hats are located at @xmath140 ( solid line ) and @xmath141 of @xmath34 ( dashed line ) and the wide distribution , corresponding to @xmath142 , goes in @xmath143 ( dotted line ) .
the disk mass is the same in all cases and equals 1 . the outer specific angular momentum @xmath144 .
the top - hat distributions are constructed using @xmath145 .
parameter @xmath61 .
[ fig.g_init ] , scaledwidth=45.0% ]    if the initial surface density is distributed as a power - law with radius @xmath146 , we use the following expression in eq .
( [ eq : f_sigma_l ] ) : @xmath147 assuming the mass of the disk is enclosed between the radii @xmath148 and @xmath149 .    in fig
.  [ fig.mdotin ] the variation of the accretion rate at the disk center is presented for three different initial distributions of @xmath7  ( see fig .
[ fig.g_init ] ) : two narrow top - hats ( @xmath150 ) with different @xmath148 and @xmath149 and the one corresponding to @xmath151 .
all three initial distributions are constructed for the same initial mass of the disk .
the evolution of a thin ring located far from the outer disk edge has a characteristic power - law interval , when the disk behaves as if it was infinite .
the interval lasts about @xmath152 $ ] .
if @xmath151 , the mass of the disk at @xmath91 is enclosed mainly in the outer parts , because it grows as the cube of the radius . as expected , the accretion rate s evolution after the peak depends weakly on the peculiarities of the initial distribution of mass in the disk ( compare the dotted and solid line in fig .
[ fig.mdotin ] ) , as long as its bulk was similarly positioned .
the viscous heating in the disk ( its half ) is found from @xmath7 : @xmath153 which can be equated to @xmath154 assuming the local energy balance .
in the center of the disk , the quasi - stationary distribution develops soon and the accretion rate is constant over radii : @xmath155 .
the zone of constant @xmath8 expands over time .
the viscous heating with the quasi - stationary distribution of @xmath7 diverges at the zero coordinate . in the real astrophysical disks ,
the inner radius of the disk is a finite value @xmath156 .
note that integrating over the disk surface from @xmath157 to @xmath158 also yields the incorrect result for the total power released in the disk .
( as a remark , the integration of gives the correct result for the case of the zero accretion rate at the center , see
[ s.deaddisk ] ) .
in fact , setting @xmath159 at the finite radius @xmath157 leads to a time - independent correction of @xmath7 near @xmath157 .
this correction is found for the infinite disks by @xcite and represents the classical factor @xmath160 .
naturally , the correction holds for a finite disk .
it follows that the bolometric luminosity can be calculated via the classical formula @xmath161 at times @xmath162 with the first term in the correction of order of @xmath163 for earlier times .
for the disks in the stellar binary systems , usually , ratio @xmath164 , and we are safe to calculate @xmath165 by the classical formula with @xmath166 found from .
0.5 cm      condition @xmath88 at necessitates that all @xmath167 in , which becomes for @xmath91 : @xmath168 where @xmath80 are the positive roots . in the previous section ,
the zero root contributed nothing to ( [ eq : gen_sol_finite_zerobi_t0 ] ) . here
the zero root provides a term like @xmath169 , which must contribute something , because the bessel function @xmath75 behaves as @xmath170 in the vicinity of zero .
the method to find @xmath82 for @xmath171 is analogous to that for @xmath81 , since the dini s series are convergent for the order of the bessel function @xmath172  @xcite  @xcite . ] .
we consider only non - integer @xmath173 and obtain @xmath174    it is known about the dini s series of the negative order that an additional first term is to be added to the sum , if the equation defining roots @xmath80 has particular coefficients=0 , the the additional term appears if  @xcite . ] .
just this case occurs for the properly rewritten outer boundary condition ( [ eq : zero_dotm_out ] ) : @xmath175 the additional term in the dini s series is expressed as @xmath176 where @xmath22 is a free variable  @xcite .
thus , adopting @xmath177 , from we obtain the green s function for the equation of the non - stationary accretion with the zero inner accretion rate and the finite outer radius : @xmath178 where @xmath80 are the positive roots of eq .
( [ eq : out_boundary_neg ] ) and @xmath179 . formulas ( [ eq : f_from_green ] ) and ( [ eq : mdot_from_green ] ) can be applied without changes .
this green function is illustrated in fig .
[ fig.green_finiteb1 ] . at the late stages ,
the accretion disks forget all information about the initial distribution of the viscous stresses .
for the disks without accretion on to the center , the uniform distribution of @xmath7 develops , with its magnitude being proportional to the mass stored in the disk .
we take the derivative of ( [ eq : green_zeroai ] ) and obtain the green function for the accretion rate : @xmath180 it approaches zero when @xmath181 for any time @xmath121 .
this corresponds to our boundary condition at the center : @xmath182 .
this disk , without sink of mass and additional mass supply , emits radiation even when it reaches the ` end ' of its evolution , as the viscous heating does not stop when the flow of matter ceases . if fig .
[ fig.lumin_zeroai ] its bolometric luminosity is plotted for the two initial positions of the ring ( the details are given in the next subsection ) .      the stable configuration of the confined ` dead ' disk ( we adopt the name from the work by @xcite ) has the following viscous angular momentum flux at each @xmath17 : @xmath183 where is used and @xmath184 as before .
the total luminosity of the dead disk is a constant value @xmath185 .
it can be obtained by the integration of over the disk surface between @xmath157 and @xmath158 .
unlikely to the case of the zero torque at @xmath186 , the late - time value of @xmath7 does not depend on the location at which condition is set and solely depends on the evolution at the large radii .
this was shown for the infinite disks by @xcite .
value of @xmath157 affects the limunosity only by a change in the size of the emitting surface .
there is still a time - dependent correction obtained by @xcite having an effect at early times @xmath187 .
the dead - disk state has to end sometime , as the disk can not receive infinitely the angular momentum from the central object .
if the central star has a magnetosphere , its radius @xmath188 can be estimated from the equality of the viscous and magnetic torque at the magnetosphere boundary  ( lp74 ) : @xmath189 where @xmath190 is the magnetic dipole of the star , @xmath191 is a dimensionless factor of order of unity  @xcite . within some factor ,
this condition is equivalent to the equality of the gas and magnetic pressure at the inner disk boundary  @xcite .
the accretion is inhibited if the inner disk radius , @xmath188 , is greater than the corotation radius , because the drag exerted by the magnetic field is super - keplerian : @xmath192 where @xmath193 , @xmath194 and @xmath195 are the mass , angular speed and period of the revolution of the central star . in the last expression
the mass is normalized by @xmath196  m@xmath197 .
thus we obtain for a final steady configuration @xmath198 the ratio can be rather close to unity .
the normalizations are @xmath199  g for the disk mass , @xmath200  [ g@xmath201 for the magnetic dipole , 100  days for the viscous time , r@xmath197 for the outer disk radius .
situation of @xmath202 is studied by  @xcite , who argue that no considerable expulsion of matter from the disk is expected in such regime , as contrasted to the propeller scenario . using the magnetosphere radius as the inner disk radius , we find the steady disk luminosity : @xmath203 for @xmath204 from the both sides , for typical @xmath205 ( see table  [ tab : zeroes ] ) .
the upper limit on the luminosity of the disk with conserved mass is estimated when the inner radius equals the corotation radius : @xmath206 the regime can be sustained for sufficiently strong magnetic dipole and fast rotation of the neutron star : @xmath207    .
the disk starts from the outer radius ( the solid line ) or from @xmath208 ( the dashed line ) , corresponding to the green s functions in fig .
[ fig.green_finiteb1 ] .
[ fig.lumin_zeroai ] , scaledwidth=43.0% ]    the flux of the angular momentum @xmath7 at small radii and the inner radius of the disk always manage to adjust themselves to the outer variations that proceed on longer viscous times .
the steady value @xmath209 could be low due to a degradation of @xmath210 in the outer disk with low luminosity .
then the normalized viscous time in the formulae above may be greater accordingly to @xmath211 .
luminosity during the spreading of the ring ( fig .
[ fig.green_finiteb1 ] ) is calculated by integrating expression over both surfaces of the disk confined between the outer and inner radius ; it is plotted in fig .
[ fig.lumin_zeroai ] .
the inner radius of the disk is varying according to eq .   and
depends on the star s magnetic dipole .
its typical value @xmath212r@xmath197 .    the torque transfers the angular momentum of the central star to the disk and , via the disk , to the orbital motion of the binary .
the corotation radius increases gradually and , eventually , accretion on the center begins  @xcite .
the characteristic brake time is the angular momentum of the star @xmath213 divided by the torque @xmath214 : @xmath215 where @xmath216  g @xmath217 and a typical value of @xmath50 is substituted .
if the mass of the disk grows due to the matter income , the braking time of the star shortens , while the magnetosphere radius decreases at the same time .
we return now to the standard accretion disks with zero stress tensor at the center . for large @xmath121 ,
the first term of the sum dominates , and a single exponential dependence on time prevails : @xmath222 this corresponds to the well - known exponential dependence on time of the accretion rate in the disk with the time - independent viscosity @xmath4 .
the characteristic time of the exponential decay is @xmath223 where we use @xmath224 . in table
[ tab : zeroes ] we give locations @xmath225 of the first zeroes of equation ( [ eq : out_boundary ] ) for several values of @xmath50 .
the table also lists the numerical factors in the expressions for the rise and decay time , and .    at late times , the profile of @xmath226 is self - similar , because @xmath227 .
the distribution @xmath228 is also self - similar .
actually , this profile develops already at the very beginning of the exponential decay .
this is illustrated in fig .
[ fig.mdisk_mdot ] by the fact that the ratio of the inner accretion rate to the disk mass , which is the integral of the surface density , approaches a constant value that depends only on the form of the radial distribution . in table
[ tab : zeroes ] , the values of the parameter @xmath229 can be found .
parameter @xmath230 describes the self - similar profile of @xmath231 and was calculated by  @xcite for the non - stationary finite @xmath3-disks . for the stationary disks with non - zero accretion rate , @xmath232 .
what does the presence of this profile mean for the released energy in the disk ? using the energy balance equation for the optically thick disks
, we obtain that at the outer radius @xmath233 , i.e. the outer ring luminosity is by factor @xmath234 less for the viscously evolving disk ( @xmath235 ) than in the case of the stationary disk . the factor @xmath236 for @xmath237 .      .
[ fig.mdisk_mdot ] , scaledwidth=45.0% ]    let us now compare the finite - disk solution for the time - independent @xmath4 and that obtained for the turbulent viscosity in the general form @xmath238 .
for the latter form of viscosity , eq .   acquires the following view : @xmath239 embracing a dimension constant @xmath240 and dimensionless parameters @xmath241 and @xmath242 : @xmath243 for a finite disk with such type of viscosity the inner accretion rate decays as a power law  @xcite : @xmath244 where @xmath245 is the accretion rate at the time @xmath91 , which can be attributed to any moment of the decay stage .
.   describes the viscous evolution of the @xmath3-disks .
the characteristic time of the solution @xmath246 , where @xmath247 is a numerical constant that can be calculated for specific @xmath248 and @xmath62 . using the expression for @xmath249 and relation @xmath250
, we obtain @xmath251 for the @xmath3-disk in the kramers regime of the opacity , parameters @xmath252 , @xmath253 , and @xmath254  @xcite should be used in eqs .  . in this case , the dimensionless factor in   is close to unity .
value of @xmath255 depends on the choice of the zero time .
ccccccl 0 & @xmath256 & 1.0585 & @xmath257 & @xmath258 & 1.267 & constant @xmath4 + @xmath259 & @xmath260 & 1.2430 & @xmath261 & @xmath262 & 1.363 & @xmath263disks with constant @xmath65 ; @xmath264disks + @xmath265 & @xmath266 & 1.2927 & @xmath267 & @xmath268 & 1.392 & @xmath263disks , thomson opacity ( @xmath269 ) + @xmath270 & @xmath271 & 1.3793 & @xmath272 & @xmath273 & 1.444 & @xmath263disks , kramers opacity ( @xmath274 ) + 1 & @xmath259 & 1.5708 & @xmath275 & @xmath276 & 1.571 & @xmath277 + 2 & @xmath278 &  &  &  &  & @xmath279 does not depend on @xmath17 +   + @xmath271 & @xmath280 & 3.4045 & @xmath281 &  & 0 & dead @xmath263disks , thomson opacity + @xmath265 & @xmath266 & 3.3425 & @xmath282 &  & 0 & dead @xmath263disks , kramers opacity    in the shakura  sunyaev disks ,
the viscous stress is proportional to the total pressure in the disk , and this proportionality is parametrized by the @xmath3-parameter  @xcite : @xmath283 then the kinematic viscosity can be expressed using ( [ eq : stress_kepl_av ] ) as @xmath284 or @xmath285    to approximate the viscous evolution of the @xmath3-disk by the disk with steady viscosity , one has to estimate the appropriate parameter @xmath62 for eq .   and .
this can be effectively done using the relation @xmath286 .
the solution for the stationary pressure - dominated stationary regions  @xcite gives the relative half - thickness @xmath287 , consequently , @xmath288 or @xmath289 in , for the free - free opacity regime .
table  [ tab : zeroes ] lists the parameters approximating in a similar way the standard @xmath3-disk with the thomson opacity and the ` dead ' disks models  @xcite .
( @xmath290 , the solid line ) and @xmath291 ( kramers opacity , @xmath252 , @xmath253 , the dashed line ) .
disk parameters @xmath34 and @xmath292 are the same in the both models .
[ fig.teordep ] , scaledwidth=45.0% ]    in fig .
[ fig.teordep ] the two laws are plotted : exponential @xmath293 for @xmath235 and power - law for @xmath294 .
viscosity at the outer radius and the outer radius itself are the same .
if applied to the real bursts during the time spans of order of @xmath295 , the models can not be easily discriminated  ( this was also pointed at by * ? ? ?
this fact provides a possibility to estimate @xmath3 using the @xmath1-folding time of the light curve . in the next section ,
we show how it could be done using the observations of the disks in x - ray novae , whose light curves in many cases show exponential decays .
in the x - ray novae , binary systems with a compact star and , typically , a less - than - solar - mass optical star , transient outbursts are thought to originate due to a dwarf nova - type instability  ( e.g. , * ? ? ?
the temperature rises in a portion of the ` cold ' neutral disk , which becomes ionized , and the viscosity @xmath4 rises .
moreover , the turbulent parameter @xmath3 changes from @xmath296 to higher @xmath297 .
the heating front moves outward from some ignition radius with the speed @xmath298 .
it stops at a radius , where the surface density of the disk is too low to sustain the stable ( on the thermal time - scale ) hot state  @xcite .    in the dwarf novae ,
the hot state accretion is quenched swiftly ( comparing to the viscous time ) by the matter returning to the cold state .
it was suggested that in the x - ray novae the illumination from the center has a stabilizing effect on the outer disk , keeping the temperatures above the hydrogen ionization temperature @xcite . @xcite
propose that in order for the disk to undergo the viscous evolution an intensive irradiation must be included .
a canonical light curve of an x - ray nova is described as a fred  @xcite . theoretically ,
if the disk is truncated from outside , a steeply decaying light curve is produced  @xcite .
the main property is that the disk does not expand infinitely , which is a consequence of its being in a binary system , and its radius does not change on a viscous time scale .
scenario by @xcite produces exponential decays due to the steadiness of the parameter @xmath4 in time .
more realistic temporal dependence of @xmath4 leads to a steep power - law decay  @xcite , which is close to an exponential one during few viscous times ( see the previous section ) .
the possibility to model the x - ray novae bursts , which have the fred profile , by the viscously evolving disk with the steady viscosity and after a descrete mass - transfer event , was shown by  @xcite ( see also the modeling in  @xcite ) .
if a burst is due to an instability affecting a slab of matter extended along the radii , the green s function obtained in the present work should be used . effectively , after a heating front has passed
, the dwarf - nova type instability leads to some distribution of the hot - ionized matter with high viscosity . in the most trivial case , when the bulk mass is concentrated near @xmath158 , the evolution of the disk approaches the one found by @xcite .
let us compare the typical time of the heating front propagation @xmath299 and the time @xmath134 it takes the inner accretion rate to peak , which is of order of @xmath300  ( see and fig .
[ fig.cmax1 ] ) .
relating @xmath3 and @xmath4 by , we obtain : @xmath301 where the gas temperature is normalized by @xmath302  k , central star mass by a solar mass , and radius by @xmath303  cm , @xmath304 is the molecular weight of the hot matter .
radius @xmath125 is the characteristic radius where the bulk mass resides at the beginning of an outburst .
this ratio indicates that we are safe to model an x - ray nova outburst at time @xmath305 using the green s function obtained .    if one assumes the commonly suggested surface density distribution in the quiescent state in the x - ray novae , @xmath306  @xcite , then , as it is shown in fig .
[ fig.mdotin ] , after the peak of the burst , the only important parameters are the disk mass and the outer radius . for such density distribution , the bulk mass of the disk resides close to the outer radius . as we show in sect .
[ s.decay ] , already at the maximum of an outburst the disk acquires the self - similar radial distribution , which is not dependent on the initial distribution of matter .
the necessary condition for this is that the bulk mass of the viscously evolving disk is located near the outer radius .
this is consistent with the absence of a power - law interval before the exponential decline on the the light curves of the x - ray novae that show the fred behaviour .
times the square relative half - thickness of the disk at @xmath158 for the six bursts of x - ray novae .
lower panel : corresponding estimates of the viscosity parameter @xmath4 at @xmath307r@xmath197 .
[ fig.alpha_p ] , scaledwidth=45.0% ]    eqs .   and
yield the ratio @xmath308 taking into account that @xmath134 is very close to @xmath309 for @xmath310 .
the value of @xmath311 is depicted in fig .
[ fig.cmax1 ] , where @xmath134 is the peak time for the finite disk . for @xmath235
( see table  [ tab : zeroes ] ) we obtain @xmath312 .
this provides a rough test whether the particular light curve is produced by the viscously evolving disk made of hot ionized matter and with approximately constant outer radius .
the viscous - disk relation @xmath313 for the hot - disks with the kramers opacity in the outer parts leads to a conclusion that to obtain a secondary maximum on the x - ray light curve of an x - ray nova , the disk has to provide an extra mass input at the radius that is determined by the time of the input .
the secondary maximuma are usually observed around @xmath314 after the peak .
if the extra mass input happens at the time of the burst maximum then it should happen at least at a radius about @xmath315 ( from @xmath316 and @xmath317 ) , without mass input from the intermediate radii .
another possibility is that the secondary maximum is triggered later after the peak and then it can happen at a radius @xmath318 ( as in the model of @xcite ) .
if an x - ray nova has the viscously evolving accretion disk with constant @xmath158 during the interval of the ` exponential decay ' ( usually , first @xmath319  days after the peak ) , we can infer the value of the turbulent parameter @xmath3 , using the closeness of the model with steady visosity and the @xmath3-disk . as the irradiation , which provides the steadiness of the hot zone outer radius , does not change the vertical structure of the disk  @xcite
, we can adopt @xmath235 applicable for the hot ionized standard @xmath3-disks .
combining eqs .   and , we get @xmath320 which relates the @xmath3-parameter to the observed @xmath321 , the time of the bolometric luminosity or the inner accretion rate .
the last formula resembles quite a few others found in the literature .
similarly , the relation of @xmath3 to the half - thickness of the disk can not be eluded . however , an advatage is that contains the e - folding time which can be inferred quite accurately from observations .
we have sampled six fred light curves in different x - ray novae ( fig .
[ fig.srcs_lc ] ) and fitted them with the steady-@xmath4 viscously - evolving disk model .
for particular values of @xmath322 , the model curves agree well with the peak - normalized count - rates before the second maximum , which is observed in most of the cases . for a0620 - 00
, we use the value of @xmath321 , estimated by @xcite from the @xmath323  kev light curve during the decline .
variation of the initial inner radius of the hot zone leads to slight changes of the model curve for the particular initial distribution of @xmath228 , as can be seen in the panel for .    in fig .
[ fig.alpha_p ] we plot the estimates for @xmath3-parameter times the normalized square of the relative half - thikness at the outer radius of the disk .
theses estimates are obtained using with the values @xmath322 indicated in fig .
[ fig.srcs_lc ] . in the lower panel ,
the viscosity @xmath4 is plotted , found accordingly to .
as @xmath4 is not constant over radius , we choose to present its value at the same radius r@xmath197 for all cases .
the considered x - ray novae have known orbital periods .
the masses of the primary and secondary components , mass ratios are from @xcite . for black hole candidate xtej1753.5 - 0127 with @xmath324
hr  @xcite we set rather arbitrarily
@xmath325  m@xmath197 and @xmath326  m@xmath197 .
the roche lobes @xmath327 are calculated using the formula by  @xcite , and @xmath328 .
the vertical barrs in fig .
[ fig.alpha_p ] reflect the errors related to the uncertainties of the binary parameters and do not take into account uncertainty of @xmath322 .
@xcite give the half - thickness of the outer parts for the stationary @xmath3-disks when the hydrogen is completely ionized : @xmath329 with the parameters normalizations equal to @xmath330gs@xmath331 and @xmath332  cm .
it can be rewritten for the factor in fig .
[ fig.alpha_p ] as follows @xmath333 where @xmath3 is normalized to 0.1 and @xmath334 is the conventional rate of accretion producing the critical luminosity : @xmath335  ergs@xmath331 . for the black hole masses in the range @xmath336  m@xmath197 , accretion rates in @xmath337 , @xmath3 in @xmath338 , and @xmath158 in @xmath339r@xmath197 , the above factor changes between extremes @xmath340 and @xmath341 .
for 4u1543 - 47 with parameters @xmath342  @xcite , @xmath343  m@xmath197  @xcite , @xmath344  r@xmath197  @xmath345 and @xmath346 one obtains the factor @xmath347 .
it follows from fig .
[ fig.alpha_p ] that @xmath3 can be less than 1 in this case , but this is still not a self - consistent value . apparently , the disk in this longer - period x - ray nova was not entierly involved in the viscous evolution that produced the fred light curve in 2002 .
levels of @xmath3 of the other five bursts ( fig .
[ fig.alpha_p ] , the upper panel ) agree with the estimates for the fully ionized , time - dependent accretion disks obtained by various approaches in x - ray and dwarf novae  @xcite , in the decretion disks of be stars  @xcite .
the values of @xmath322 , indicated in fig .
[ fig.srcs_lc ] , are to be considered with caution .
note that the observed light curves are limited to an energy band .
spectral modeling might be required in many cases to obtain accurate @xmath322 .
if an x - ray nova during the burst is in the _ high / soft _ state , when the disk emission dominates the spectrum , the spectral modeling can provide @xmath348 in a straightforward way .
we expect @xmath322 not to alter by a factor larger than 2 for the x - ray flux dominated by the disk ( preliminary resluts for 2 - 10  kev ) .
on the other hand , some x - ray novae bursts proceed entirely in the _ low / hard state _
, such as the burst of xtej1753.5.-0127  @xcite .
still , the form of the corresponding x - ray light curves is informatve . even if the optically thick disk does not reach the innermost orbits , it is very probable that the radiated energy is proportional to the accretion rate , because the released potential gravitational energy is the ultimate source to power the source .
second , the innermost disk peculiarities do not influence the behavior of @xmath349 as the latter is defined by the process of the viscous matter redistribution in the outer disk .
equation is written without the tidal - stress term , following the proposition that it should be negligible everywhere in the disk , except in a thin ring near the tidal truncation radius  @xcite . according to @xcite
, effects near the tidal - torque radius can be approximated with an effective boundary condition @xmath182 .
numerical models confirm that most of the tidal torque is applied in a narrow region at the edge of the disk , where perturbations become non - linear and strong spiral shocks appear  ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) . as @xcite notes , in the case of dwarf novae , during outbursts , the observed values of the outer radius appear to be approximately consistent with the theoretical predictions
@xcite argue that the action of the tidal torques are important also well inside the tidal radius .
they notice at the same time that observations do not allow discriminating between the rapidly growing and smoother tidal stresses , as the model light curves are not strongly affected .
the tidal radius for all mass ratios can be approximated as @xmath350 of the average roche lobe radius of a given component in circular binaries  @xcite .
the disk can extend beyond this radius due to the high inclination of the binary orbit  @xcite .
effects of the high eccentricity of the binary orbit on the disk truncation radius have been studied by various methods ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) .    the gravitational influence of the secondary on the disk in a binary system leads to a number of interactions , non - resonant ( or tidal ) and resonant , any of those transferring angular momentum from the disk to the binary  @xcite .
the tidal distortions generally lead to the truncation of the disk , whereas the resonances of different strength impact or do not the disk structure , depending on their growth rate and other conditions  @xcite .
an example of a resonance action is the debated explanation of the superhamps in su uma stars  ( see , e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) .    in our mathematical set - up ,
we have presumed a fixed outer radius of the disk without going into details of its actual value .
thus the model is appplicable to the coplanar disks in low - eccentricity binaries .
the outer radius of the disk can also change due to intrinsic processes . reproducing the whole evolutionary cycles of the dwarf novae and x - ray novae , @xcite ,
@xcite should have taken into account the variations of the outer disk radius , but these changes are significant during the time spans much greater than the viscous times .
under some astrophysical conditions , viscously evolving disks formed in binary systems are effectively truncated from outside . in the present work ,
we use the method proposed by @xcite to find the green s functions for the linear viscous evolution equation for the disk with the finite outer radius .
green s functions are obtained for the viscous angular momentum flux @xmath7 and the accretion rate @xmath8 .
two inner boundary conditions at the zero coordinate are considered .
they correspond to the accreting disks and to the disks with no mass sink at the center .
the green s functions allow one to compute @xmath9 and @xmath351 for arbitrary initial surface density profiles that develop into self - similar distributions on sufficiently long time scale .
the analytic formula to calculate the disk evolution with the variable outer mass inflow is derived .
a solution , which can be found with the use of the green s function , has the properties defined by the main equation of the viscous evolution and thus provides the basic description of the transient phenomena in the viscous disk .
the green s functions are found in the form of quickly converging series and can be easily reproduced with standard computer methods .
the c - code written for the case of the accreting disc with the use of the gnu scientific library can be downloadedgalja / lindisk.tgz ] ] .
we present the relations between the rising time , @xmath1-folding time , and disk viscosity . for six bursts in the x - ray novae , which are of fred type ,
we show that the model describes well the peak - normalized light curves before the second maximum .
this favors the mainly viscous nature of their evolution during this period and enables us to obtain an estimate of viscosity @xmath4 , which depends on the outer disk radius .    the models with time - independent viscosity are shown to approximate well the evolution of @xmath3-disks during time intervals comparable to the viscous time .
consequently , estimates of the @xmath3 parameter may be derived for the disks in the binaries with known period and masses .
these estimates rely on the relative half - thickness of the disk at the outer radius .
the estimates of @xmath3 for the five short - period x - ray novae outbursts , groj0422 + 32 in 1992 , a0620 - 00 in 1975 , gs1124 - 68 in 1991 , gs2000 + 25 in 1998 , and xte j1753.5 - 0127 in 2005 , are in line with the values estimated so far for the hot viscous disks  @xcite .
another green s fucntion is found for the disk that has zero accretion rate at the inner radius and acquires angular momentum from the central star . in the steady state
, the mass of the disk can not be very high and the disk has low luminosity .
it radiates most of the rotational power transferred from the central star .
this disk works as a gear transmitting the angular momentum of the central star to the orbital motion .
the ` dead ' stage is ended by an abrupt fall or trickle of the matter on to the star after the centrifugal barrier has moved close enough or beyond the magnetosphere s radius .
the author is grateful to n.  i.  shakura , k.  stempak , j.  l.  varona , m.  perez and the anonymous referee for useful comments .
the work is supported by the russian science foundation ( grant 14 - 12 - 00146 ) .
green s function for the infinite disk , obtained by lp74 , can be written in the dimensionless form with our designations ( c.f .
formula ) as follows : @xmath352 where @xmath353 is the modified bessel function of the first kind , @xmath354 , @xmath355 is some characteristic value of the specific angular momentum , which we use instead of @xmath69 .    to find the physical distribution of @xmath226 in the course of the evolution of the initial narrow ring of matter , initially located at the coordinate @xmath356
, one takes the integral of @xmath357 ( [ eq : fring ] ) with the kernel @xmath358 on the interval @xmath359 $ ] : @xmath360 where @xmath127 is the initial mass of the ring .
the accretion rate at the center can be found from @xmath361 ( see ( [ eq : mdot_from_green ] ) ) .
one finds @xmath362 or @xmath363 where the typical times @xmath364 and @xmath365 are introduced , @xmath366 we find that the accretion rate peaks at the time @xmath367 with the value @xmath368
for infinite disks , see @xcite and @xcite , who have obtained the green s function solution to the viscous evolution of a disk with mass sources / sinks , which are distributed in a fashion along the disk radius .
we have the following problem with the inhomogeneous boundary and initial conditions : @xmath369 let us substitute function @xmath231 by a function @xmath370 using the relation @xmath371 this substitution is used , for example , to solve the problem of finding the temperature distribution in a cylinder , on a surface of which a thermal flux is defined @xcite . for the new @xmath370
the problem is the following : @xmath372 where @xmath373 the corresponding sturm
liouville problem with a free variable @xmath374 is investigated in sect .
[ sec : ev_fin_disk ] : it consists of eq . (
[ eq : sturm - liuv ] ) @xmath375 where @xmath376 or @xmath377 plays a role of an eigenvalue , and the boundary condition @xmath378 at @xmath379 . changing to the free variable @xmath380 , we can write for a non - integer @xmath50 @xmath381\ , , \ ] ] where @xmath382 and @xmath383 depend on the inner - boundary condition . according to the steklov theorem ,
the solution to ( [ eq : problim_inh2 ] ) can be obtained in the form : @xmath384 where coefficients @xmath385 are to be found from the inhomogeneous first - order differential equation over @xmath121 .
let us substitute ( [ eq : steklov ] ) into ( [ eq : sturm - liuv1 ] ) and ( [ eq : problim_inh2 ] ) , multiply by @xmath386 and integrate over @xmath83 from 0 to @xmath138 , using the orthogonality property of the eigenfunctions .
we arrive at @xmath387 with designations @xmath388 where @xmath389 is the norm of eigenfunction @xmath390 .
solution to is as follows @xmath391 let us write down the solution to ( [ eq : problim_inh2 ] ) for the case @xmath392 , that is , @xmath393 : @xmath394 the last expression is equivalent to . in sect .
[ sec : ev_fin_disk ] , we find the coefficients @xmath395 and @xmath82 for the dirac delta function as the initial condition and express @xmath226 using the green s function . in the analogy ,
the solution to the problem with @xmath396 has the view @xmath397 where @xmath220 is given by ( [ eq : green_zerobi ] ) or ( [ eq : green_zeroai ] ) .    substituting ( [ eq : sol_nonzeromdot ] ) into ( [ eq : substit ] ) and using functions @xmath398 and @xmath399 from ( [ eq : problim_inh2 ] )
, we finally arrive at ( [ eq : sol_nonzeromdot_full ] ) .
the coordinate integral in its second term can be taken analytically ; it involves the lommel functions and can be tabulated beforehand for the particular values of @xmath50 and @xmath80 for the sake of the computational efficiency .
the second term of for @xmath94 can be expressed as : @xmath400 \mathrm{d } \tau~,\ ] ] where @xmath401                                                                                , a.  t. 2007 , in astronomical society of the pacific conference series , vol .
367 , massive stars in interactive binaries , ed .
n.  st .- louis & a.  f.  j. moffat , 485 , astronomical society of the pacific conference series