a decade of extrasolar planet discoveries has shown that the process of planet formation is more complex than originally anticipated .
it leads to a remarkable diversity of planetary configurations , ranging from migrating hot jupiters to eccentric giant planets , as well as our own `` circular '' solar system . following these discoveries ,
progress has been made in understanding planet formation , but the theory is still incomplete ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the leading scenario for the formation of giant planets is the core accretion mechanism .
icy planetesimals located beyond the snow line in their host disk collide repeatedly to grow a core with a modest gaseous atmosphere .
if this core succeeds in reaching a critical mass of a few tens of earth masses , runaway accretion of a massive gaseous envelope proceeds and leads , ultimately , to the formation of a gaseous giant planet @xcite .
some evidence supporting this scenario has emerged in recent years , in the form of a metallicity trend for stars hosting planets ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , the discovery of a high density hot jupiter @xcite and that of a surprisingly low - mass micro - lensing planet @xcite .
however , a long - standing difficulty for the core accretion scenario , which has not yet been fully elucidated , is the fact that the timescale required to build - up critical core masses and thus large gaseous envelopes ( @xmath0-@xmath1  yr ) is comparable to the lifetimes of proto - planetary disks ( e.g. * ? ? ?
* ; * ? ? ?
orbital migration adds a layer of complication to theories of planet formation . as a result of gravitational interactions with their gaseous disk @xcite , the orbits of planets in the terrestrial mass range
are predicted to decay on timescales ( @xmath2  yr ) short compared to disk lifetimes ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
migration is slower for planets of much smaller or much larger masses : in the first case because the torque causing migration is quadratic in planet mass , and in the second case because the planet opens a gap and then migrates on the disk s accretion timescale , which can be comparable to its lifetime . while it is possible or probable that many terrestrial planets form by agglomeration of smaller bodies after the gas is gone , this is not an option for the solid cores of jovian planets since , in the core accretion scenario , the cores must form before the gaseous envelopes can be accreted .
the prevalence of jovian planets with orbital periods of only a few days deepens the mystery as it suggests that these planets did migrate but stopped short of merging with their stars at orbital radii where even the accretion timescale would seem to have been very short ( e.g. * ? ? ?
* ) .    analytic calculations and most hydro - dynamical simulations of migration usually assume a disk that is laminar apart from the waves and shocks excited by the planet itself ( e.g. * ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
? * ; * ? ? ?
* ; * ? ? ?
but the effective viscosity of disks probably involves turbulence .
@xcite , @xcite , and ( * ? ? ?
* hereafter n2005 ) have found in 3d simulations of magneto - rotational turbulence that the instantaneous torque exerted on a planet in the terrestrial mass range is subject to fluctuations many times its mean value , apparently caused by turbulent density fluctuations in the planet s vicinity .
in fact , no obvious secular decay manifests itself in the orbits of planets with mass @xmath3 , although because the simulations are limited to @xmath4 planetary orbits
 compare @xmath5 for jupiter during the lifetime of the proto - solar nebula  and the predicted decay in semimajor axis is @xmath6 over this period , the expected secular trend might be difficult to discern amid the fluctuations .
the possibility that planetary torques are stochastic raises several issues :    1 .
is the time - averaged torque exerted on a planet of a given angular momentum affected ?
n2005 finds that turbulence excites planetary eccentricity , @xmath7 . in linear theory
, eccentricity allows the planet to interact with the disk at additional resonances that do not couple to circular orbits @xcite .
also , at second and higher orders in @xmath7 , the radial positions of the original resonances shift if one compares eccentric and circular orbits at fixed angular momentum ; they do not shift at fixed orbital energy since keplerian energies and mean motions vary only with semi - major axis .
thus for eccentric planets , one must distinguish the rate of change of the planetary angular momentum ( torque ) from the rate of change of orbital energy .
smoothing of the planetary potential along epicycles tends to reduce the strength of the highest-@xmath8 lindblad resonances that dominate the torque on circular orbits .
all of these perturbative effects are expected to become significant when @xmath9 , where @xmath10 is the disk thickness .
planets of terrestrial mass achieved such eccentricities in n2005 s simulations . on the other hand , at about the same eccentricity the epicyclic velocities become supersonic , which will probably raise significant shocks , a non - perturbative effect .
thus both the sign and magnitude of the changes in the mean torque due to turbulently excited eccentricity are unclear .
we are not aware of any systematic investigation of this issue , nor have we attempted one . throughout this paper ,
we treat the planetary orbit as circular .
2 .   even for circular orbits , stochastic changes in semi - major axis on timescales
long compared to the period but short compared to the nominal migration time may significantly alter the mean orbital lifetime .
this is particularly likely when the local mean torque varies strongly with radius , as it does in the alpha disks studied by @xcite . because of rapid changes in opacity in certain temperature regimes , the disks exhibit sharp peaks in the local migration timescale . under laminar conditions ,
orbital lifetimes can be dominated by the time required to drift across these peaks .
turbulent fluctuations in the torque , however , could allow planets to `` jump over '' the peaks , thus shortening the lifetime .
on the other hand , diffusion to large radii where the local migration time is long might prolong the orbital lifetime .
3 .   finally , orbital diffusion should allow the occasional planet to survive much longer than the mean lifetime defined by its birthplace .
if sufficiently many planets are formed , this might in principle help to reconcile theoretical predictions of rapid _ mean _ migration with the observations at least of the longer - period planets ( though it seems unlikely to explain hot jupiters ) .
the methods adopted in this paper allow us to explore the second and third issues above , but not the first . in
[ sec : fp ] , we develop a simple advection - diffusion or fokker - planck equation for the probability distribution function of planetary orbital angular momentum in the presence of torque fluctuations .
parametrizations of the advection and diffusion coefficients are calibrated against semianalytic formalisms and the simulations of n2005 , respectively .
3 explores analytic steady - state and time - dependent solutions to this equation in steady disks whose properties are power laws in radius .  4 and  5 present numerical , time - dependent solutions in such a disk ( the minimum mass solar nebula ) and in alpha disks with realistic opacities .
a summary of our main results and conclusions is given in  6 .
planets migrate as a result of torques exerted by the disk . because of turbulence ,
the torque is a stochastic function with a mean part @xmath11 and a fluctuating part @xmath12 . for a given planetary mass ( @xmath13 ) and a given time - averaged surface density and thickness of the disk [ @xmath14 respectively ] , the mean torque depends only on the semimajor axis , @xmath15 , of the planetary orbit , which is assumed to be circular .
the diffusion equation is most conveniently derived in terms of the orbital angular momentum , where @xmath16 is the mass of the central star . ]
@xmath17 , rather than @xmath15 , as the spatial independent variable .
the fluctuating part of the torque depends on time as well as @xmath18 : @xmath19 . by construction , @xmath20 .
the fluctuating torque is taken to be a temporally stationary stochastic process with a finite amplitude and correlation time , so that the integral @xmath21 exists .
@xmath22 will play the role of diffusion coefficient .
the correlation time can then be defined by @xmath23 ^ 2}\,.\ ] ] temporal stationarity implies that @xmath22 , @xmath24 ^ 2}$ ] , and @xmath25 are independent of @xmath26 .
the desired advection - diffusion equation derives from a fokker - planck equation , and its formal derivation is completely standard ( _ e.g. _ * ? ? ?
* ) , but it is useful to repeat the derivation here in order to emphasize the underlying assumptions and to ask whether they are justified in this case .
two important assumptions have already been introduced : the finiteness of eqs .
( [ ddef ] ) & ( [ tcorr ] ) .
let @xmath27 be the probability that a planet of initial angular momentum @xmath18 suffers a change @xmath28 during time interval @xmath29 .
moments of the change are @xmath30 of course , except for @xmath31 , these depend implicitly on @xmath18 and @xmath29 .    next , let @xmath32 be the probability that the planet has angular momentum @xmath18 at @xmath33 given some initial condition at @xmath34 .
the probability at @xmath35 is @xmath36\,d\delta j \nonumber\\ & = & \sum\limits_{n=0}^\infty\frac{(-1)^n}{n ! } \frac{\partial^n}{\partial j^n}\left\{\vphantom{\frac{a}{a } } \overline{(\delta j)^n}f(t , j)\right\}\,.\end{aligned}\ ] ] in the fokker - planck approximation , the changes are presumed small compared to @xmath18 itself but cumulative , so that one can ignore the terms @xmath37 in the sum above .
more will be said about the justification for this below .    clearly @xmath31 , @xmath38 , and @xmath39 ^ 2}\\[1ex ] & = & \bar\gamma(j)^2(\delta
t)^2~+~\int\limits_0^{\delta t}\,d t_m \int\limits_{-2 t_m}^{2 t_m}\,d \tau~   \overline{\delta\gamma ( t_m-\,\tau/2,j)\delta\gamma ( t_m+\,\tau/2,j)}.\end{aligned}\ ] ] the next important assumption is that @xmath25 is short compared to the timescale on which @xmath18 changes by order itself  let us write @xmath40 for the latter timescale  so that we may take @xmath29 in the range @xmath41 . in that case , the double integral above @xmath42 . at this point , @xmath43\delta t ~+~o\left[(\delta t)^2\right],\ ] ] where the @xmath44 $ ] part includes the @xmath45 contribution to @xmath46 . neglecting these terms on the grounds that @xmath47 , one has the desired advection - diffusion equation : @xmath48 f_j & \equiv & \bar\gamma(j ) f(t , j)~-~\frac{\partial } { \partial j } \left[\vphantom{\frac{1}{2 } } d(j)f(t , j)\right]\,.\end{aligned}\ ] ]    an interesting feature of the above result is that if the diffusion coefficient @xmath22 varies with @xmath18 , then it influences the _
mean _ migration rate .
that is to say , one can rewrite the flux in eq .
( [ diffusion ] ) as @xmath49 so that @xmath50 contributes to the advection term .
we have tried to abstract a value and scaling for @xmath22 from the mhd simulations of n2005 ; the results of @xcite were difficult to translate into our framework .
nelson reports that the typical timescale of fluctuations in @xmath51 is of order half an orbital period .
hence we take @xmath52    a direct estimate of the correlation function of the data suggests , however , that correlations may persist to much longer than an orbital period , and it is conceivable that the integral ( [ ddef ] ) for @xmath22 does not even exist ( nelson 2005 , private communication ) .
if the correlation time does not exist , then the methods of this paper are inapplicable to migration .
such long - term correlations imply a `` memory '' in the turbulence , presumably involving persistent structures in addition to the planet of interest : for example , long - lived vortices ( e.g. * ? ? ?
* ; * ? ? ?
* ) , or of course other planets .
pending further numerical evidence , we adopt as a working hypothesis that @xmath25 _ does _ exist but with uncertain magnitude .    to complete the diffusion coefficient
, we require a parametrization of the variance of the fluctuating torque in terms of time - averaged disk properties .
this entails some guesswork , as the simulations of n2005 and @xcite have explored only a limited range of disk models and planetary radii .
a natural scale for the gravitational force exerted on the planet by the local gas is @xmath53 , where @xmath54 is the surface density of the disk .
this is the force that the planet would feel if suspended just above the disk .
the corresponding natural scale for the torque is @xmath55 .
therefore , we postulate that the variance of the torque fluctuations is @xmath56 ^ 2}\approx \left(\mathcal{c}_d\times 2\pi g\bar\sigma r_p m_p\right)^2,\ ] ] where @xmath57 is a dimensionless coefficient depending upon the strength of the turbulence .
since the parametrization ( [ amptorque ] ) underlies all of our results , it is worth more discussion .
the turbulence is accompanied by density perturbations of r.m.s .
amplitude @xmath58 on scale @xmath59 .
an individual perturbation of this scale exerts a gravitational force per unit mass @xmath60 at distance @xmath61 if @xmath62 , and @xmath63 if @xmath64 . presuming that the forces from different perturbations add in quadrature , the mean - square force is dominated by the closest such perturbation , at @xmath65 , since the number of @xmath59-scale perturbations within distance @xmath61 increases more slowly with @xmath61 than the mean - square contribution from individual perturbations decreases : @xmath66 for @xmath67 and @xmath68 for @xmath69 , versus @xmath70 .
thus , the r.m.s
. force due to scale @xmath59 is @xmath71 for @xmath62 and @xmath72 for @xmath64 .
the contributions to @xmath22 from @xmath62 are probably unimportant because of the weighting by @xmath59 and because smaller scales are likely to have shorter correlation times .
if the density power spectrum were as `` blue '' as white noise , @xmath73 , and if the correlation time @xmath74 , then there would be equal contributions from all scales @xmath75 , with the result that the total contribution from these scales would be larger than the contribution from @xmath76 by a logarithmic factor .
in reality , the density power spectrum is surely red , since local magnetorotational ( henceforth mri ) simulations find that most of the turbulent magnetic and kinetic energy resides at scales @xmath77 @xcite .
thus , the dominant contribution to the force comes from scales @xmath77 . since @xmath78 in this regime , the total r.m.s .
force fluctuation is proportional to the total r.m.s .
surface - density fluctuation , @xmath79 .
it remains to discuss how @xmath80 and hence @xmath81 may depend upon the strength of the turbulence .
density fluctuations may arise from fluctuations in gas pressure or in entropy .
local mri simulations indicate that @xmath82 @xcite , as might be expected from equipartition , since the internal energy per unit mass associated with a density fluctuation is @xmath83 . in a statistically steady accretion disk , the local thermal time @xmath84 , so if the correlation time of mechanical fluctuations is @xmath85 , entropy fluctuations due to dissipation are expected to be at most @xmath86 .
entropy fluctuations may also arise by advection of material in the presence of a background entropy gradient , @xmath87 .
the simulations cited above indicate that the turbulent reynolds stress scales with the magnetic stress so @xmath88 , and we assume that the largest displacements have frequencies @xmath89 , so the r.m.s .
displacement @xmath90 in any direction .
vertical and horizontal entropy gradients are likely to have scale lengths @xmath91 and @xmath92 , respectively , which would seem to make vertical displacements more important .
however , the azimuthal force , and hence torque , exerted by a mass element changes only at second order as it is displaced vertically through @xmath93 , so vertical displacements contribute at @xmath94 .
radial displacements contribute at first order , @xmath95 .
the vertically uniform disk of n2005 has a radial entropy gradient , since the sound speed is constant but the density is not ; however , since n2005 obtains @xmath96 ( in agreement with previous work ) , it happens that @xmath97 in these models , so that we can not distinguish between @xmath98 and @xmath99
.    it would be useful to calculate torque fluctuations in a variety of disk models with varying entropy gradients , thicknesses , and ( if possible ) strengths of turbulence .
lacking such information , we provisionally assume @xmath100 , where @xmath101 is some function ( probably linear ) .
this leads to the parametrization ( [ amptorque ] ) if @xmath102 is a universal constant .
as will be seen , diffusion then dominates mean migration , _
i.e. _ inward drift at the rate predicted for a laminar disk , in the outer parts of plausible disk models .
if @xmath79 is also proportional to @xmath103 , diffusion becomes even more dominant in the outer parts of a flaring disk .
n2005 reports that the r.m.s .
fluctuating torque per unit planetary mass is @xmath104 in code units . in the same units , @xmath105 ( nelson 2005 , private communication ) ;
the latter is independent of radius because @xmath106 .
therefore , @xmath107 is taken as the reference value of the dimensionless factor in equation ( [ amptorque ] ) .    following eq .
( [ tcorr ] ) , the diffusion coefficient @xmath22 is the product of the correlation time @xmath25 ( [ tcorrsim ] ) and the torque variance ( [ amptorque ] ) .
both factors are uncertain , but within the framework of this advection - diffusion approach , they combine into a single parameter .
we hope that future numerical simulations will refine our estimates of this parameter and , more importantly , test the validity of the scalings ( [ tcorrsim ] ) and ( [ amptorque ] ) .    for comparison , @xcite s expression for the mean lindblad torque density
, @xmath108 leads to the scaling @xmath109 for consistency with our calibration of torque fluctuations , we also calibrate this relation against nelson s simulations of laminar disks , in which the type i torque per unit mass is @xmath110 in code units . the constant of proportionality in eq .
( [ eq : meantorque ] ) is then @xmath111 . using eq .
( [ eq : meantorque ] ) implies that the mean torque is always negative .
this is fine for plausible power - law disks but might be misleading for alpha disks , where the surface density and temperature can be highly structured , allowing in principle for positive ( or at least much reduced ) @xmath11 at some radii .
thus when modeling alpha disks we use eq .
( [ eq : ward ] ) directly , as described in @xcite .
power - law disk models such as the minimum - mass solar nebula ( hereafter mmsn ; @xcite ) allow the advection - diffusion equation to be solved analytically , at least in steady state .
some aspects of time - dependent evolution can also be determined analytically for this class of disk models . throughout this section ,
all disk models are considered to be infinite in radius as well as steady . while unrealistic , this is technically convenient , and we are still able to identify results that are likely to be sensitive to the details of the outer boundary conditions : for example , the mean orbital lifetime . for reference ,
the mmsn model adopted here has mean surface density @xmath112 and aspect ratio @xmath113 .
a steady - state distribution of planets , @xmath114 , satisfies @xmath115 where the right - hand side is a source term representing the rate of formation of planets at radius @xmath15 with angular momentum @xmath17 .
analytic solutions to eq .
( [ eq : steadystate ] ) exist when mean torque and the diffusion coefficient are both power - laws in @xmath18 : @xmath116 and the source term is a dirac delta function @xmath117 where @xmath118 is the rate of planet formation at @xmath119 .
we see immediately that the flux of planets interior @xmath120 and exterior @xmath121 to the source are both constant and related to each other by @xmath122 first we consider the limiting behavior as @xmath123 .
the local timescale for secular inward migration is given by @xmath124 , and the local diffusion timescale is given by @xmath125 .
the ratio , @xmath126 indicates that the mean torque becomes asymptotically negligible if @xmath127 .
( this condition holds for the minimum - mass solar nebula , in which @xmath128 and @xmath129 . ) for such disks the solution to eq .
( [ eq : steadystate ] ) at large @xmath18 exterior to the source becomes @xmath130 , where @xmath131 is a constant of integration .
we insist that @xmath132 since there should not be an incoming flux of planets from infinity .
thus the numerator eventually becomes negative .
the diffusion coefficient @xmath22 is unquestionably positive , so @xmath101 becomes negative at sufficiently large @xmath18 .
but this is absurd , since @xmath133 is a probability .
thus we can not have both a positive flux and a positive @xmath101 at large @xmath18 in steady - state , and it follows that @xmath134 and @xmath135 .
all the flux generated at the source is accreted through the inner edge of the disk .
we now discuss eq .
( [ eq : steadystate ] ) in full , emphasizing behavior at small and large @xmath18 .
let @xmath136 so that @xmath137 $ ] is an integrating factor for eq .
( [ eq : steadystate ] ) . then @xmath138.\ ] ] as before , the upper ( lower ) sign on @xmath139 and @xmath140 applies to @xmath141 ( @xmath142 ) .
when @xmath143 and @xmath22 are power laws , @xmath144 can be evaluated explicitly , but we will not write it out .
normally the mean torque @xmath145 , so that @xmath146 . then we can choose the constant of integration in eq .
( [ eq : mu ] ) so that @xmath147 for @xmath148 and @xmath149 for @xmath150 .
all disks of interest in this paper are dominated by the mean torque at small radii and by diffusion at large @xmath18 .
hence @xmath151 $ ] becomes large at small @xmath18 , and so for a well - behaved solution , the contents of the square brackets in eq .
( [ eq : fss ] ) must tend to zero as @xmath152 .
this defines a relationship between @xmath153 and @xmath154 . for @xmath141 , @xmath155 , so the integral increases at least as fast as @xmath156 . as before ,
if the disk is indefinitely extended then we must have @xmath157 to prevent @xmath158 from becoming negative at large @xmath18 .
the main conclusion of this discussion is that _ as long as the mean torque is nonpositive in the outer disk , the flux of planets diffusing to infinite radius vanishes .
_ we have reached this conclusion by considering steady solutions but will assume that it applies generally .
the time - dependent problem is more difficult , but we can determine some general characteristics of the solution by analytic means , most importantly that _ the probability of surviving without accretion in an infinite steady disk , though declining monotonically with time , can have a long power law tail_. for the special case of the mmsn , where @xmath129 and @xmath159 , this probability declines only as @xmath160 .
the asymptotic form of the solution at late times is self - similar except near @xmath161 and is described explicitly in dimensionless units by equation ( [ eq : ssim ] ) .
the evolution of a planet created with @xmath162 at @xmath34 can be represented by solving eq .
( [ diffusion ] ) with a source term @xmath163 .
the laplace transform @xmath164 reduces the time - dependent problem to @xmath165- \frac{\partial}{\partial j}\left[\bar\gamma(j)\hat f(s , j)\right ] -s\hat f(s , j)=-\delta(j - j_s).\ ] ] unfortunately , although eq .
( [ eq : odes ] ) is just an ordinary differential equation in @xmath18 , we can not solve it completely even for power law disks described by eq .
( [ eq : powers ] ) .
the difficulty is in part that physically interesting cases have @xmath166 , which is precisely the condition for eq .
( [ eq : odes ] ) to have an irregular singular point at @xmath161 .
however , a formal solution can be obtained in powers of @xmath167 , @xmath168 replacing @xmath169 by its taylor series in the laplace transform ( [ eq : laplace ] ) shows that @xmath170 is the @xmath171 moment of @xmath32 with respect to @xmath26 , @xmath172 more interesting than the moments of @xmath101 are the moments of the corresponding flux : @xmath173 the probability that the planet survives longer than @xmath26 before accretion , @xmath174 , is related to the flux at the origin by @xmath175 .
therefore , @xmath176 is the @xmath171 moment of the lifetime , and in particular , @xmath177 is the mean lifetime .
thus we are motivated to find the @xmath170 .
substituting eq .
( [ eq : expans ] ) into eq .
( [ eq : odes ] ) leads to @xmath178''-[\bar\gamma(j)f_0(j)]'&=&-\delta(j - j_s ) , \\
\label{eq : recur }    [ d(j)f_n(j)]''-[\bar\gamma(j)f_n(j)]'&=&-n f_{n-1}(j ) \qquad n>0,\end{aligned}\ ] ] where the primes indicate derivatives with respect to @xmath18 .
the solution to eq .
( [ eq : f0 ] ) , which is formally identical to the steady - state case ( [ eq : steadystate ] ) , is a green s function for eq .
( [ eq : recur ] ) .
it can be constructed from homogeneous solutions @xmath179 and
@xmath180 with fluxes @xmath181 and @xmath182 , respectively : @xmath183'~-\bar\gamma(j)y_0(j ) & = & 0,\\    \label{eq : y1 } [ d(j)y_1(j)]'~-\bar\gamma(j)y_1(j ) & = & 1.\end{aligned}\ ] ] the green s function is @xmath184 the formal solution to the system ( [ eq : f0])-([eq : recur ] ) is then @xmath185 the existence of the moments ( [ eq : moment ] ) therefore depends upon the convergence of the integrals in ( [ eq : ssol ] ) .
an important special case is the mmsn , where @xmath186 and @xmath187 .
it is convenient to choose @xmath188 as the unit of angular momentum and @xmath189 as the unit of time so that @xmath190 and @xmath191 .
the homogeneous solutions ( [ eq : y0 ] ) become @xmath192 & \approx & \cases {     j^2-(3!!)j^4+(5!!)j^6-\ldots & $ j\to0^+$ \cr     1-\sqrt{\frac{\pi}{2}}\sum\limits_{n=0}^\infty\frac{1}{(2n)!!}j^{-2n-1 } ~+\sum\limits_{n=1}^\infty\frac{1}{(2n-1)!!}j^{-2n } & $ j\to\infty$. } \ ] ] the first expansion is only asymptotic , in keeping with the fact that @xmath161 is an irregular singular point of eq .
( [ eq : y1 ] ) , whereas the second converges for all @xmath193 .
armed with these results , we can now determine the convergence of the integrals ( [ eq : ssol ] ) . because @xmath194 as @xmath195 , it follows from ( [ eq : green ] ) that the integrand for @xmath196 is proportional to @xmath197 as @xmath198 .
thus the integral over @xmath199 is logarithmically divergent .
consequently _ the mean lifetime defined by eq .
( [ eq : tmom ] ) is not finite_. since the divergence is only logarithmic , one may guess ( correctly ) that @xmath200 at late times in the infinite mmsn .
the long tail in the distribution of lifetimes results from diffusion to large radii where torques are weak and planets linger .
it is natural to seek a self - similar solution to eq .
( [ diffusion ] ) that reflects this behavior .
unfortunately @xmath201 defines a preferred scale that prevents self - similarity . omitting @xmath11 from eq .
( [ diffusion ] ) on the grounds that this term is unimportant at @xmath202 removes the obstruction : @xmath203,\ ] ] in the present dimensionless units with @xmath204 as in the mmsn .
the tilde distinguishes the solution to the approximate equation ( [ eq : nogamma ] ) from the exact solution @xmath32 of eq .
( [ diffusion ] ) .
we require @xmath205 as @xmath195 .
on the other hand , @xmath206 should be positive and finite at @xmath161 for the following reason . at late times such that @xmath207 , we expect @xmath32 to match onto a multiple of the constant - flux steady - state solution defined by eq .
( [ eq : y1mmsn ] ) : @xmath208 .
we may replace @xmath101 with @xmath209 in this relation if @xmath202 , and in that region @xmath210 .
( [ eq : nogamma ] ) can be solved with these boundary conditions by separation of variables . and @xmath211 , which produces the bessel equation of order one on the righthand side .
hankel transforms can then be used . ] for arbitrary positive initial conditions , @xmath212 as promised , this has self - similar form : @xmath18 enters only in the combination @xmath213 [ or @xmath214 in dimensionful units ] .
it can be verified that eq .
( [ eq : ssim ] ) is an exact solution of eq .
( [ eq : nogamma ] ) , and its integral over all @xmath18 is @xmath215 , in agreement with our expectation for @xmath174 at late times .    for other values of the indices
@xmath102 & @xmath216 satisfying @xmath217 , we have the following situation .
again , @xmath11 is negligible at large @xmath18 , so eqs .
( [ eq : y0 ] ) & ( [ eq : green ] ) yield @xmath218 at @xmath219 . therefore the integral ( [ eq : ssol ] ) for the first moment @xmath196 diverges if @xmath220 , and consequently , the mean lifetime converges only if @xmath221 . in other words ,
the mean lifetime exists only if @xmath222 increases more slowly than @xmath223 .
the following technical result helps to interpret differences between the late - time behaviors of the planetary distribution in the mmsn and in alpha disks .
the large-@xmath26 behavior of @xmath32 is dominated by the singularities of @xmath224 at largest @xmath225 .
these correspond to eigenvalues @xmath226 obtained by separating variables in eq .
( [ diffusion ] ) with assumed time dependence @xmath227 . if the eigenvalues are real and discrete , @xmath228 , then @xmath229 at late times , where @xmath230 is the eigenfunction corresponding to @xmath59 .
( we assume that these modes are complete . ) on the other hand , if the eigenvalue spectrum is continuous and extends to arbitrarily small positive @xmath59 , then the late - time decay is slower than exponential , perhaps a power law . defining @xmath231 and setting @xmath232 and @xmath233 , with @xmath234 as in eq .
( [ eq : mu ] ) , one can re - cast the separated version of eq .
( [ diffusion ] ) in self - adjoint sturm - liouville form @xmath235~+\lambda w(j ) g_\lambda(j)=0.\ ] ] in all disks of interest to us , @xmath236 as @xmath152 and the integral of @xmath101 over @xmath18 is finite , whence @xmath237 and @xmath238
. a theorem of @xcite then asserts that if the limit @xmath239^{1/2 } \left[\int_{j}^\infty \frac{dx}{p(x)}\right]^{1/2}\ ] ] vanishes , then the spectrum is discrete and positive .
[ the theorem also requires @xmath240 and @xmath241 to be positive and continuous , which is always true for us , and it encompasses `` quasilinear '' generalizations of eq .
( [ eq : sl ] ) . ]
it is easy to show that in disks where @xmath242 tends to zero as @xmath152 and to infinity as @xmath195 , the limit ( [ eq : dk ] ) is equivalent to @xmath243^{1/2}.\ ] ] in the mmsn [ @xmath204 ] , this limit diverges , so a continuous spectrum extending to @xmath245 is allowed and indeed must exist , or else we would not have the self - similar behavior discussed above . in the alpha disks studied numerically in
[ sec : numerical ] , on the other hand , it appears that @xmath246 as @xmath195 [ fig .
( [ fig:4dprop ] ) ] , so the limit vanishes and a discrete spectrum is implied .
this is clearly true of power - law disks with @xmath127 and @xmath247 , so the survival probability @xmath174 should decay exponentially in these disks .
we conjecture that the continuum exists and @xmath174 follows a power law for all @xmath248 , because further consideration of eqs .
( [ eq : ssol ] ) indicates that the @xmath171 moment of the lifetime exists only if @xmath249 .
the diffusion equation ( [ diffusion ] ) with a source term reads : @xmath250 we find numerical solutions with an implementation of the generic relaxation code described in @xcite .
the finite - difference scheme is fully implicit in time . given a distribution @xmath32 at time @xmath26 , we only need to solve the one - dimensional boundary - value problem in @xmath18 to find the solution at time @xmath251 .
we thus relax successively from one solution to another over a variable timestep , @xmath252 .
controls over the size of the timestep ensure that , at each grid point , the change in each variable over @xmath252 is within an acceptable tolerance . the system is described by two first - order differential equations equivalent to eq .
( [ fpsource ] ) , supplemented with a third differential equation that defines the numerical domain . for simplicity , we do not implement an adaptive mesh but use instead a fixed grid on the numerical domain , chosen to be uniformly spaced in @xmath253 .
we take as variables the set @xmath254 .
numerical stability of the solutions is achieved with centered - differencing schemes for @xmath101 and @xmath18 , and an upwind - differencing scheme for @xmath255 .
three boundary conditions are required .
the numerical domain extends from the inner edge of the disk , at @xmath256 , out to @xmath257 . at the inner boundary we trivially set the value of @xmath18 and are free to fix @xmath258 . at the outer boundary we demand the flux
in steady - state , as discussed in  [ subsec : steady ] , this is a boundary condition that must be satisfied if diffusion dominates at large radii . to avoid potential complications near the disk edge ,
we extend our numerical domain to sufficiently large radii that it does not influence our results .
we have tested our numerical results for power - law models against the analytic steady - state solutions described in  [ subsec : steady ] .
the numerical solutions converge satisfactorily for a modest grid spacing ( typically 1000 grid points over five decades in radius ) .
additional tests at ten times this resolution show clear numerical convergence of the steady - state results .
agreement with the late - time behavior expected from the time - dependent self - similar solutions described in  [ subsec : timedep ] will also be demonstrated below ( see fig .  [
fig : mmsnsv ] ) .
we study the combined effects of migration and diffusion on the orbital evolution of low - mass proto - planets for two separate classes of disk models : the minimum mass solar nebula ( mmsn ) and t tauri alpha disks .
the mmsn model adopted here has a mean surface density @xmath112 and an aspect ratio @xmath259 , which are normalized to @xmath260 and 0.1 at @xmath261 , respectively .
the t tauri alpha disk model adopted is identical to the irradiated disk described by ( * ? ? ?
* see their fig .  1 ) , accreting steadily at a rate @xmath262 .
we take as our reference model the alpha disk with viscosity parameter @xmath263 . for comparison
, we also consider a more massive and cooler alpha disk , with @xmath264 .
the inner edge of the disks is located at @xmath265 and @xmath266 in the mmsn and the alpha disks , respectively .
the mass of the primary is @xmath267 in the mmsn and @xmath268 in the alpha disks .
these slight model differences have very little influence on our results .
note that , while the disk models used are assumed to be steady , we are studying the time - dependent orbital evolution of planets embedded in these disks with equation ( [ diffusion ] ) . as a consequence
, our method can only apply to proto - planets of small enough mass that they exert a negligible feedback on the disk structure ( broadly speaking , @xmath3 ; e.g. * ? ? ?
turbulence is assumed to be present everywhere in our disk models .
the magnitude of the corresponding diffusion coefficient is uncertain through both the amplitude and the correlation time of the turbulent torque fluctuations .
we quantify these uncertainties with a single parameter , @xmath269 .
if we define @xmath270 then @xmath271 ^ 2}=\epsilon\,(2\pi)^3(0.046)^2\frac{\bar{\sigma}^2j^7}{g^2m_*^4m_p^5}.\ ] ] the calibration performed in
[ subsec : calib ] suggests that @xmath272 is a reasonable value , which we have assumed to be independent of radius .
however , we also investigate the effects of varying this parameter by several orders of magnitude .
it is useful to discuss a few general scalings specific to our disk models before we discuss the details of time - dependent numerical solutions . to quantify the relative importance of diffusion _ vs. _ migration in the models ,
we define a local diffusion time @xmath273 ( for @xmath274 ) and a local migration time , @xmath275 according to eq .
( [ eq : meantorque ] ) , so that @xmath276 in the mmsn . the ratio of these timescales has a particularly simple dependence upon disk and planetary properties : @xmath277    one expects the evolution of planets located in regions of a disk with @xmath278 to be dominated by migration and vice - versa .
while the diffusion time is independent of planet mass , the migration time is inversely proportional to @xmath13 .
therefore , lower mass proto - planets will be more sensitive to the effects of turbulent torque fluctuations .
in addition , the two timescales differ in their dependence on the mean surface density , @xmath279 , which tends to make planets in denser disks more susceptible to the effects of diffusion .
as we shall now see , a general consequence of these scalings is that one expects the effects of diffusion to dominate only in the outer regions of proto - planetary disks ( at least for the classes of disk models considered here ) .
figure  [ fig:4dprop ] shows profiles of various important quantities in the mmsn and our two alpha disk models .
panel i shows profiles of mean surface density , @xmath280 . both the slope and normalization of @xmath280 profiles matter for the relative importance of migration _ vs. _ diffusion .
panels iii and iv show profiles of local migration and diffusion times in the various disk models .
while @xmath281 increases linearly with @xmath18 ( @xmath282 ) in the mmsn , it decreases with increasing radius in the alpha disk models .
this is a consequence of the steeper @xmath280 profile in the mmsn .
the profiles of @xmath283 , derived from eq .
( [ eq : meantorque ] ) for the mmsn and from the detailed calculations of @xcite for the alpha disk models , show local migration times which decrease toward small radii in both classes of models .
finally , panel ii shows the ratio of timescales , @xmath284 , which is important to characterize the behavior of planets embedded in these disks .
as mentioned earlier , in all three models , this ratio grows above unity ( indicated by the horizontal dash - dotted line in panel ii ) only at radii @xmath285 if @xmath286 , indicating that diffusion should dominate over migration only in the outer regions of these disks ( and this despite the increasing value of @xmath281 with radius in the mmsn ) .
although very useful , the profiles shown in figure  [ fig:4dprop ] do not capture all of the complexity of the migration - diffusion problem .
first , all the quantities plotted in panels ii - iii - iv were derived for a planet mass @xmath287 and a torque fluctuation normalization @xmath288 .
the profiles are easily rescaled for different values of these parameters ( as indicated by labels in fig .
[ fig:4dprop ] ) , but it is clear that the exact location of the transition from migration- to diffusion - dominated evolution depends on the proto - planetary disk model adopted , the value of the planet mass , and the normalization @xmath269 of the diffusion coefficient . nevertheless , the dominance of diffusion over migration in the outer regions of these disks is robust because @xmath289 tends to increase with radius in disk models of interest [ eq .
( [ eq : migvsdiff ] ) ] .
secondly , as shown in panel iii , in a disk model with realistic opacities , sudden and non - monotonic variations in the value of @xmath283 at radii @xmath290 could complicate the orbital evolution of planets substantially , by stalling their migration as they reach the disk inner regions .
this is best explored with global numerical solutions of the time - dependent fokker - planck equation , as seen in
[ subsec : burst ] .
a third complication arises from the radial dependence of the diffusion driven by turbulent torque fluctuations . as mentioned in
[ sec : fp ] , the radial variation of the diffusion coefficient leads to an additional contribution to the advected flux of planets . since @xmath291 for @xmath274 , and since @xmath292 in the mmsn and @xmath293 in the alpha disk models , this additional term tends to advect planets inward in all the disk models considered .
we can define a local advection timescale associated with this new contribution , @xmath294 , and compare it to the values of @xmath283 and @xmath281 defined above .
we find that this advection timescale is everywhere exactly equal to @xmath281 in the mmsn , while it is shorter than @xmath281 everywhere in our alpha disk models ( as shown explicitly for the @xmath263 case by the short - dashed line in panel iv of fig .
[ fig:4dprop ] ) . consequently , even in disk regions where diffusion nominally dominates over migration in our models ( as measured by the ratio @xmath295 ) , this additional contribution due to radially inhomogeneous diffusion should still cause a rather significant effective advection of planets towards the star .
we have studied two different types of time - dependent models for the diffusive migration of planets in proto - planetary disks , corresponding approximately to a continuous and a burst scenario for planet formation . in both cases , for the sake of clarity , the source of planets was localized , taking the form of a sharply peaked gaussian approximating the delta source function of eq .
( [ eq : sourceterm ] ) .
we have verified that our results do not depend sensitively on the detailed shape of the gaussian source function adopted .
as we shall see , the predictions for these rather different scenarios for planet formation are broadly consistent with each other .
figure  [ fig:4steady ] shows steady - state distributions obtained by numerically solving eq .
( [ fpsource ] ) with a source term centered at @xmath296 , steadily producing earth - mass planets at a rate @xmath297 .
the topmost curve in each panel is the steady - state distribution .
the enclosed curves , on the other hand , are snapshots of the decaying distributions at selected time intervals after the source is shut off .
we have verified that the steady - state distributions for the mmsn in panels i and iii are fully consistent with the analytic solution , eqs .
( [ eq : mu])-([eq : fss ] ) , for power - law disks .
a fiducial normalization of @xmath288 for torque fluctuations is assumed in all cases except in panel iii , where @xmath298 has been scaled down to extend the migration - dominated regime out to @xmath299 in the mmsn .
thus , at @xmath300 , the disks in panels i and iv are diffusion - dominated and the disks in panels ii and iii are migration - dominated .
the distinction between diffusion - dominated and migration - dominated evolutions is most clearly illustrated by comparing panels  i and  iii . while the initial steady - state peaked distribution is eroded and advected for the most part in panel iii , with little diffusion at large radii , there is significantly more diffusion at large radii in panel i. that distinction is not as clear in the alpha disk models , however .
note that , in both panels  ii and  iv , the initial steady - state distribution is strongly peaked not at the source radius , but at radii where migration stalls ( as seen from panel iii in fig .
[ fig:4dprop ] ) .
in addition , even though diffusion is more important in the model with @xmath264 in panel iv , the effect appears to be minor in smoothing out only slightly the same sharp distribution features as in panel ii .
this is partly explained by the strong contribution of radially inhomogeneous diffusion to inward advection of planets , which is always important in our alpha disk models ( as explained in the previous subsection ) .
even though diffusion is equally important at @xmath296 in the models of panel i and iv ( see @xmath301 in fig .  [ fig:4dprop ] ) , it is more inhomogeneous in the alpha disk model and thus promotes a significantly more efficient inward advection of planets .    our models with a burst of planet formation
confirm these general trends . to study this formation scenario ,
we solve the diffusion equation ( [ diffusion ] ) without source term and impose a specific initial condition for the distribution of planets at @xmath34 instead .
the initial distribution , @xmath302 , is modeled as a narrow gaussian centered on @xmath303 .
figure  [ fig:4spike ] shows the evolution of such distributions of earth - mass planets initially located at @xmath296 in the same four disk models as in fig .
[ fig:4steady ] .
as before , a comparison between the two mmsn models shows that the effects of diffusion are much more important in the model with @xmath304 ( panel i ) than in the one with @xmath298 ( panel iii ) .
in particular , a sub - population of planets diffusing at large radii is clearly visible in panel i. on the contrary , while diffusion still acts to widen the initial distribution in panel iii , the dominant effect is a global advection of the distribution of planets .
finally , in agreement with our previous results , figure  [ fig:4spike ] shows that the most significant feature of alpha disk models is the effect of rapid variations in the migration rates of planets , which causes them to stall at specific locations in the disk and results in time - varying peaks in the distributions of panels ii and iv . even though diffusion is comparatively as important in panels i and iv , at @xmath300 , its effects appear to be much less important in the alpha disk model of panel iv .
we attribute this difference , once again , to the more inhomogeneous diffusion in alpha disk models , which contributes substantially to the inward advection of planets into regions of the disk where diffusion becomes gradually less important and the effects of stalling radii are felt more strongly .
two important aspects of the orbital evolution of planets which are not well captured by figure  [ fig:4spike ] are the properties of the surviving minority of planets at late times and large radii , and the representative planetary lifetimes in burst - formation models
.    the fraction of planets that diffuse to orbital radii greater than their formation radius , @xmath305 , can be characterized by the time - dependent quantity @xmath306 where @xmath307 is the angular momentum at a fiducial radius @xmath308>r_s$ ] .
hereafter , the integral over the initial distribution appearing in the denominator of eq .
( [ eq : chia ] ) will be normalized to unity .
the maximum values of @xmath309 and @xmath310 , reached typically just after the majority of planets were accreted by the star , are listed in table [ tb : chia ] for the four models of figure  [ fig:4spike ] .
not surprisingly , models with a strong role for diffusion tend to have a fairly large maximum fraction of planets diffusing beyond @xmath311 ( @xmath312 ) and @xmath313 ( @xmath314 ) , while models with little diffusion have much reduced maximum fractions of planets at large radii ( @xmath315 beyond @xmath311 , @xmath316 beyond @xmath313 ) .
lccc disk model & torque fluctuations & @xmath317 & @xmath318 + mmsn & @xmath288 & @xmath319 & @xmath320 + alpha - disk ( @xmath263 ) & @xmath288 & @xmath321 & @xmath322 + mmsn & @xmath323 & @xmath324 & @xmath325 + alpha - disk ( @xmath264 ) & @xmath288 & @xmath326 & @xmath327 +    we compute three additional time - dependent quantities to help us better characterize the properties of surviving planets in our models .
the first is the survival probability , computed as @xmath328 the second is the median orbital radius of surviving planets , @xmath329 , defined by @xmath330 finally , since the mean lifetime is often ill - defined (  [ subsec : timedep ] ) , we calculate the median lifetime @xmath331 defined by @xmath332 .
figure  [ fig : mmsnsv ] shows how the survival probability and the median radius of surviving planets evolve with time in burst formation models for earth - mass planets in the mmsn .
three separate models with burst formation at @xmath333 , @xmath334 and @xmath313 are shown for comparison ( solid lines ) .
note that the middle solid line corresponds to the exact same model as the one shown in panel i of figure  [ fig:4spike ] . for reference ,
vertical dotted lines bracket typical lifetimes for gaseous proto - planetary disks ( @xmath335@xmath1  yrs ; e.g. * ? ? ?
the median lifetime of the population of planets , which is such that @xmath336 , is indicated by a horizontal dash - dotted line in the upper panel .
the lifetimes differ substantially in the three models shown , from much shorter to somewhat longer than the typical disk lifetimes .
the self - similar @xmath337 behavior predicted in
[ subsec : timedep ] is indeed achieved at late times in the three models .
as shown in the lower panel of figure  [ fig : mmsnsv ] , the median radius of surviving planets does not evolve much at early times and remains close to the formation radius , as the distribution of planets mostly diffuses out to both small and large radii .
once a significant fraction of planets have been accreted by the central star , however , the surviving planets are found at a median radius increasing almost quadratically with time independently of where they were formed .
this is just what one would expect from the self - similar solution ( [ eq : ssim ] ) : @xmath338 .
figure  [ fig : alpsv ] shows the evolution of the same two quantities in burst formation models for earth - mass planets in the alpha disk ( with @xmath263 ) this time .
six separate models with burst formation at @xmath333 , @xmath339 , @xmath334 , @xmath340 , @xmath341 and @xmath342 are shown for comparison ( solid lines ) .
note that the fourth solid line from the top corresponds to the exact same model as the one shown in panel ii of figure  [ fig:4spike ] .
the results are qualitatively consistent with the mmsn ones , with a few important differences .
there is no self - similar behavior at late times in this case and the more dominant role of advection in the alpha disk model tends to make the median radius of surviving planets move inward initially , until the majority of them are lost to the central star . among the surviving minority
, however , @xmath343 asymptotes to @xmath344 at late times , independently of the formation radius , @xmath305 .
this behavior is explicable if , as predicted by the theorem cited in
[ subsubsec : eig ] , the eigenmodes of eq .
( [ diffusion ] ) are discrete : then @xmath32 converges upon the most slowly decaying mode , and @xmath329 tends to the median radius of this eigenmode .
figures  [ fig:4dprop ] and  [ fig : alpsv ] suggest that this radius is close to that at which @xmath345 [ more accurately , @xmath346 , where @xmath347 is minimized , as might be expected since planets should linger longest where the combined effects of diffusion and migration are weakest . in the mmsn ,
both timescales increase without limit as @xmath195 , and the spectrum of eigenmodes is apparently continuous .
it is also worth emphasizing that the largest peak value of @xmath348 at sub - au distances in alpha disk models ( see panel iii in fig .  [ fig:4dprop ] ) effectively limits the migration time of earth - mass planets into their star to values @xmath349  yrs ( for @xmath263 ) , which exceeds typical disk lifetimes .
based on this general understanding of diffusive migration in our burst - formation models , let us now focus more systematically on several important model dependencies .
figure  [ fig : medlrs ] illustrates how the median lifetimes of earth - mass planets formed at different radii in the mmsn and the alpha disk are affected by changes in the normalization of torque fluctuations , @xmath269 . for the nominal @xmath304 value adopted ( indicated by a vertical dash - dotted line in fig .
[ fig : medlrs ] ) , diffusion mostly affects the orbital evolution of planets formed at the largest radii , by reducing somewhat their median lifetime .
as @xmath269 is increased above this nominal value , the trend becomes clearer .
median lifetimes for all formation radii are considerably reduced at larger values of @xmath269 , more so in the mmsn than in the alpha disk .
the long - dashed line in figure  [ fig : medlrs ] results when we neglect the mean torque @xmath11 in the governing equation ( [ diffusion ] ) .
the convergence of the solid curves for the alpha - disk model to a value independent of @xmath305 is related to the value of the diffusion time , @xmath222 , at the disk inner boundary . for sufficiently large @xmath269 , @xmath350 everywhere on the grid . but @xmath222 increases toward small radius in the alpha disks ( fig .
[ fig:4dprop ] ) .
so the median lifetime in the large-@xmath269 regime is comparable to @xmath222 at the inner boundary , @xmath351 ; we have verified that @xmath352 decreases as @xmath351 is increased .
nevertheless , in all cases , the effect of diffusion is to reduce the median lifetime .
figure  [ fig : medlmp ] illustrates how the median lifetimes of planets of various masses , all formed at @xmath296 in the mmsn and in the alpha disk , are affected by changes in the normalization of torque fluctuations , @xmath269 . for the nominal @xmath304 value
, diffusion mostly affects the orbital evolution of the lower mass planets , and reduces somewhat their median lifetimes .
this results from the larger values of @xmath348 at lower planet mass , whereas @xmath222 is independent of planet mass .
as @xmath269 is increased , diffusion progressively dominates and the distinction among planets of various masses gradually disappears . for large enough values of @xmath269 , the median lifetime of all planets with masses
@xmath353 formed at @xmath354 is reduced below typical disk lifetimes . as in the previous figure , however , the lifetimes for @xmath355 in the alpha disk are dominated by the diffusion time at the inner boundary .
as anticipated by @xcite , @xcite , and @xcite , our fokker - planck treatment confirms that diffusion driven by turbulent torque fluctuations in proto - planetary disks can greatly influence the orbital evolution of embedded low - mass proto - planets .
contrary to some expectations , diffusion does not promote planetary survival but instead systematically reduces the lifetime of most planets in the disk .
however , it does help a small fraction of planets diffuse to large radii where they can survive for extended periods .
( we consider that a planet `` survives '' if its orbit remains larger than the inner edge of the disk . )
our results point to a potentially important role for non - deterministic effects in planet formation scenarios . in models with a significant level of diffusion , most proto - planets in the earth - mass range could end up being accreted by their stars , while only a small fraction of them would survive by diffusing to large radii . in scenarios where all proto - planetary disks are efficient at forming planets ,
this could be interpreted as leading to only a fraction of all potential host stars being successful in keeping a system of surviving planets .
even systems starting with very similar initial conditions may end up with very different planetary orbital configurations once their gaseous disks disappear .
neglecting the effects of direct gravitational interactions between proto - planets , which are likely to happen from differential rates of migration and diffusion , is one of many model limitations in our work .
we have already mentioned that our study is restricted to proto - planets of low enough masses that they do not affect in any way the structure of their host disk .
this excludes the late stages of giant planet formation , which are obviously of considerable interest .
in particular , it has been suggested that migration and diffusion , each independently , could accelerate the rate of growth of cores until they reach the critical mass at which runaway envelope accretion occurs ( e.g. * ? ? ? * ; * ? ? ?
* ; * ? ? ?
studying these possibilities within the context of global models such as ours would be interesting .    even within the strict regime of applicability of our models ,
significant uncertainties exist .
let us mention a few .
we have focused our work on strictly circular orbits .
we have assumed that the underlying proto - planetary disk structure does not evolve with time , and that the amplitude @xmath81 of surface - density fluctuations is constant with radius .
we have assumed that torque fluctuations can simply be added to a laminar mean torque which is not modified by the turbulence .
the laminar torque itself is calculated from a simple extension of the two - dimensional lindblad torque theory , which does not account for the role of corotation torques .    despite all these uncertainties
, we believe that our results should be broadly valid , provided that turbulence and torque fluctuations in proto - planetary disks satisfy the various stochasticity criteria enunciated in  [ sec : fp ] . understanding better the nature of turbulence in proto - planetary disks and in particular indications for long - term correlations in mhd simulations may turn out to be crucial in this context .
we would like to issue a plea to the simulation community to study these questions , and in particular , whether the dimensionless amplitude and correlation time of density and torque fluctuations depend only upon the strength of the turbulence ( _ via _
@xmath102 ) or upon the thickness @xmath103 as well (  [ subsec : calib ] ) . the answers will matter both quantitatively ( in determining the relative importance of migration _ vs. _ diffusion in global models such as ours ) and qualitatively ( in justifying a fokker - planck approach or in falsifying it ) .
we have seen that details of the protoplanetary disk structure ( mmsn _ vs. _ alpha disks ) do influence our results at the quantitative level .
in that respect , the possibility that a magnetically layered structure with quiescent dead zones exists in weakly - ionized regions of proto - planetary disks @xcite may be important .
indeed , the results of @xcite suggest that the nature of turbulence in dead zones differs qualitatively from that in mhd turbulent regions , perhaps leading to comparatively reduced torque fluctuations ( see also * ? ? ? * for other suggested effects of dead zones in planet formation scenarios ) .
the general tendency for a fraction of planets to diffuse to large radii in proto - planetary disks may have consequences of some observational relevance .
we note , for instance , that @xcite have looked into scenarios for outward planet migration in an attempt to explain observed structures in dust disks as resulting from perturbations by planets located at radii beyond those at which they are thought to be able to form . rather than invoking outward migration ,
this may be expected in scenarios with a prominent role for diffusion , as illustrated by our various models .
we thank richard nelson for help with understanding his simulations .
this material is based on work supported in part by the national aeronautics and space administration under grant nag5 - 11664 issued through the office of space science .
alibert , y. , mordasini , c. , benz , w. & winisdoerffer , c. 2005a , , 434 , 343 alibert , y. , mousis , o. , mordasini , c. & benz , w. 2005b , , 626 , l57 armitage , p.  j. & rice , w.  k.  m. ( 2006 ) , to appear in a decade of extrasolar planets around normal stars ( astro - ph/0507492 ) balbus , s. a. & hawley , j. f. 1998 , rev .
phys . , 70 , 1 beaulieu , j .-
2006 , nature , 439 , 437 bender , c.  m. & orszag , s.  a. 1978 , advanced mathematical methods for scientists and engineers ( new york : mcgraw - hill ) brandenburg , a. , nordlund , a. , stein , r.  f. , & torkelsson , u. 1995 , apj , 446 , 741 dangelo , g. , kley , w. & henning , t. 2003 , apj , 586 , 540 drbek , p. & kufner , a. , 2005 , proc .
soc . , 134 ( 1 ) , 235 .
fleming , t. & stone , j. m. 2003 , apj , 585 , 908 fromang , s. & nelson , r.  p. 2005
, , 364 , l81 fromang , s. , terquem , c. & balbus , s. a. 2002 , mnras , 329 , 18 gammie , c. f. 1996 , apj , 457 , 355 gilliland , r.  l. et al . 2000 , , 545 , l47 glassgold , a. e. , najita , j. & igea , j. 1997 , apj , 480 , 344 [ erratum : 485 , 920 ] goldreich , p. & tremaine , s. 1980 , , 241 , 425 goldreich , p. & sari , r. 2003 , , 585 , 1024 haisch , k.  e. , lada , e.  a. & lada , c.  j. 2001 , , 553 , l153 hameury , j .- m . , menou , k. , dubus , g. , lasota , j .-
, & hur , j .- m .
1998 , , 298 , 1048 hawley , j.  f. , gammie , c.  f. , & balbus , s.  a. 1995 , , 440 , 742 hayashi , c. 1981 , prog . theor .
, 70 , 35 ilgner , m. & nelson , r.  p. 2006
, , 445 , 205 johansen , a. & klahr , h. 2005 , , 634 , 1353 kley , w. , dangelo , g. & henning , t. 2001 , apj , 547 , 457 korycansky , d. g. & pollack , j. b. 1993 , icarus , 102 , 150 laughlin , g. , steinacker , a. , & adams , f.  c. 2004 , , 608 , 489 lin , d. n. c. , bodenheimer , p. & richardson , d. c. 1996 , nature , 380 , 606 lin , d. n. c. & papaloizou , j. 1986 , apj , 309 , 846 lufkin , g. et al .
2004 , mnras , 347 , 421 marcy , g. w. , cochran , w. d. & mayor , m. 2000 , in protostars and planets iv ( tucson : university of arizona press ) , eds .
mannings , v. , boss , a.p . , russell , s. s. ) , p. 1285
matsumura , s. & pudritz , r. e. 2003 , , 598 , 645 menou , k. & goodman , j. 2004 , , 606 , 520 mizuno , h. 1980 , prog . of theor .
64 , 544 nelson , r.  p. 2005
, , 443 , 1067 ( n2005 ) nelson , a. f. & benz , w. 2003a , apj , 589 , 556 nelson , a. f. & benz , w. 2003b , apj , 589 , 578 nelson , r.  p. & papaloizou , j.  c.  b. 2004 , , 350 , 849 papaloizou , j. c. b. & terquem , c. 1999 , apj , 521 , 823 papaloizou , j.  c.  b. & terquem , c. 2006 , rep .
69 , 119 pollack , j. b. et al .
1996 , icarus , 124 , 62 rice , w.  k.  m. & armitage , p.  j. 2003 , , 598 , l55 safronov , v. s. 1969 , evoliutsiia doplanetnogo oblaka i obrazovanie zemli i planet ( moscow : nauka ) ( english transl
. evolution of the protoplanetary cloud and formation of the earth and planets [ nasa tech .
f-677 ] [ jerusalem : israel program sci .
[ 1972 ] ) sano , t. , inutsuka , s. , turner , n. , & stone , j.m .
2003 , , 605 , 321 .
santos n. c. , israelian g. , mayor m. 2004 , , 415 , 1153 sato , b. et al .
2005 , , 633 , 465 schaefer , c. , speith , r. , hipp , m. , kley , w. 2004 , , 418 , 325 van kampen , n. g. 1981 , stochastic processes in physics and chemistry ( amsterdam : north - holland ) veras , d. & armitage , p. 2004
, , 347 , 613 ward , w. r. 1997a , icarus , 126 , 261 ward , w. r. 1997b , apj , 482 , l211 weldrake , d  .t .
f. , sackett , p.  d. , bridges , t.  j. & freeman , k.  c. 2005 , , 620 , 1043 wuchterl , g. , guillot , t. & lissauer , j. j. 2000 , in protostars and planets iv ( tucson : university of arizona press ) eds .
mannings , v. , boss , a.p . ,
russell , s. s. , p. 1081