while most observed protoplanetary discs display accretion rates of order @xmath6  m@xmath7  yr@xmath8  @xcite , some objects , in particular fu orionis systems ( e.g. * ? ? ?
* ) , have exhibited episodes of enhanced accretion up to @xmath9  m@xmath7  yr@xmath8  @xcite over observed time - scales of 10 - 100  years .
it has been hypothesised that _ all _ protostellar systems actually undergo such episodic accretion , equating to a small fraction ( the `` duty cycle '' ) of their total lifetime , possibly concentrated in their early stages @xcite .
indeed this would alleviate the apparent discrepancy between the low observed accretion rates and that required to achieve the final mass of the protostars @xcite .
although the mechanism underlying fu orionis outbursts has not been established to date @xcite , a presently popular class of model uses as its only ingredients the two currently favoured drivers of transport in protoplanetary disks , namely the magnetorotational instability ( mri ; @xcite ) and the gravitational instability @xcite .
the `` gravo - magnetic cycle '' , as it has come to be called @xcite , relies on the fact that a disc with a `` normal '' density is typically too cold ( @xmath10  k ) and too dense for its mid - plane to be sufficiently ionised to trigger the mri , giving rise to a _
dead zone_. as the disc is gravitationally stable at such densities , accretion will only be possible in the surface layers where x - ray / uv ( or cosmic ray ) ionisation keeps the mri alive @xcite .
however , a steady - state situation can not be achieved and mass accumulates in the dead zone , until gravitational instability sets in , raising the disc s mid - plane temperature ( @xmath11 ) .
near @xmath12  k , thermal ionisation triggers the mri throughout the entire thickness of the disc , resulting in an accretion burst which lasts @xmath13  years ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
once the previously accumulated mass is dumped on the proto - star , the disc returns to the initial `` low state '' ( or quiescent state ) .
the cycle may then resume , with a duty cycle @xmath14  % at low infall rates ( @xmath15  m@xmath7  yr@xmath8  ) up to @xmath16 at high infall rates ( @xmath17  m@xmath7  yr@xmath8  ) .
the significant changes in the disc s temperature and luminosity during such an accretion burst can chemically process a large fraction of the disc material .
several previous studies have shown the impact of accretion bursts on the chemical evolution of the envelope and disc .
@xcite studied the impact of accretion busts on the co chemistry in the envelopes , in particular the large increase in luminosity , and found that repeated bursts could lead to repeated freeze - out cycles .
these evolutionary co cycles lead to possible observational evidence of accretion outbursts in young star forming cores @xcite .
additionally , @xcite used ` single - point ' models to study the impact of accretion bursts on the local evolution of various chemical species in the envelope , indicating that the outcome depends delicately on the balance between burst and chemical time - scales .
however , what has not been considered in great detail is that this processed material can then be transported to smaller radii in the disc by advection and to larger ones by turbulent diffusion . these effects , which may persist for significantly longer than the duration of the bursts , may lend themselves to observable diagnostics ; whether for extrasolar systems , e.g. with the increased resolution and sensitivity offered by _ alma _ and _ herschel _ , or our own solar system , taking advantage of the rich body of cosmochemical evidence . to this end , in this series we set out to investigate the chemical consequences of accretion bursts .
an astrochemically abundant species of choice in this regard is water @xcite , and in particular its d / h  ratio @xcite .
observations of the interstellar medium and pre - stellar cores suggest that water starts deuterium - rich ( d / h@xmath18 ; @xcite ) . yet
bulk chondrites and comets ( with a further _ in situ _ measurement to come for comet churyumov - gerasimenko @xcite ) in our solar system display lower ( d / h)@xmath0 ratios ( around @xmath19 ) , although some micron - scale d - rich hotspots in clays of the semarkona chondrite do reach ( d / h)@xmath0 @xmath20 @xcite .
this in turn is much higher than the protosolar value ( essentially that of molecular hydrogen ) of @xmath3 @xcite , and reaction kinetics prohibit the conversion of low ( d / h)@xmath0 ratios to higher values in the protoplanetary disc @xcite .
the current interpretation is that the solar system inherited its water from the parent molecular cloud @xcite , and that this water _ partially _ exchanged its deuterium with the nebular gas @xcite .
coupling this deuterium exchange with mixing throughout the disc produced the range of intermediate ( d / h)@xmath0 values we now observe @xcite .
the bottom line is that the ( d / h )  ratio of water is a sensitive tracer of infall , temperature and turbulent mixing conditions .
all previous studies of the ( d / h)@xmath0  ratio assumed a constant turbulence level ( often parametrised in terms of the shakura & sunyaev @xmath21 , @xcite ) , and thus could not reproduce the gravo - magnetic cycle , and had no accretion outbursts .
additionally , only @xcite considered the effects of infall , during which most of the outbursts are expected to take place @xcite , along with the subsequent disc evolution .
@xcite presented steady - state analytical calculations , indicating that in such types of model , a high efficiency of turbulent diffusion - with a schmidt number @xmath22 of 0.2 or smaller - until the snow line and beyond would be required to account for the smallness of the ( d / h)@xmath0 ratio of carbonaceous chondrites and the earth relative to oort cloud comets ; but as yet the realism of such a small schmidt number is unclear @xcite .
however , as isotopic exchange is most efficient for temperatures above @xmath23  k , clearly , outbursts , which may reach temperatures in excess of 1000 k , could significantly impact the ( d / h)@xmath0 profile of protoplanetary disc and possibly account for the low ( d / h)@xmath0 ratios of carbonaceous chondrites .    in this work ,
we examine the effect of accretion outbursts arising from the gravo - magnetic limit cycle on the evolution of the ( d / h )  ratio of water using numerical models presented in section [ sec:1dmodel ] .
since this process has not been studied before , to gain insight we perform numerical experiments where the disc experiences a steady infall rate for the entire simulation ( section [ sec : steady ] ) . secondly , we perform more realistic evolutionary calculations which include a disc building phase for a few hundred - thousand years followed by a stage of isolated evolution for several million years , that are meant to represent the evolution of an actual protoplanetary disc , e.g. the solar nebula ( section [ sec : building ] ) .
the implications are then discussed in section [ sec : discussion ] and the conclusions summarized in section [ sec : conclusion ] .
in order to investigate the effect of accretion bursts on the chemical evolution of the ( d / h)@xmath0 ratio we adopt a simple one - layer approach @xcite . our model is a modified version of that originally presented by @xcite to study the impact of dead - zones in protoplanetary discs and associated accretion bursts , where we additionally include radial transport of heat by radiation and turbulence ( see * ? ? ? * ) .
here we extend the model to include the transport of water due to advection and turbulent diffusion as well as the chemical evolution of the ( d / h)@xmath0 ratio .
the evolution of the surface density @xmath24 for a disc around a forming star , with a dead - zone is given by ( e.g. * ? ? ?
* ) : @xmath25\right\}+\dot{\sigma}^g_{\rm in}\label{eqn : gas}\ ] ] where @xmath26 is the kinematic viscosity due to the mri or self - gravity which we detail in section  [ sec : viscosity ] and @xmath27 is the rate of material falling onto the disc .
@xmath28 is the surface density that is ` mri active ' and is given by : @xmath29 for all our models where we include a dead - zone ( and hence accretion bursts ) we choose @xmath30k and @xmath31 g  @xmath32 . in simulations where we do not include a dead - zone ( and thus no accretion bursts ) we set @xmath33 .
the evolution of the temperature at the mid - plane of the disc is given by @xcite : @xmath34\nonumber\\&&-\left(\frac{2h}{\sigma^g rc_p}\right)\frac{\partial}{\partial r}\left(rf_r\right)\label{eqn : temp}\end{aligned}\ ] ] where @xmath35 is the heat capacity at constant pressure , @xmath36 is the ratio of specific heats and @xmath37 is the scale height of the disc , with @xmath38 the keplerian angular velocity and @xmath39 the isothermal sound speed .
@xmath40 is the ` potential temperature ' of the disc s mid - plane and following @xcite we define it as : @xmath41 noting that equation  [ eqn : temp ] is independent of the choice of scaling for the potential temperature.the net heating rate ( @xmath42 ) is given by : @xmath43 where @xmath44 is the stefan - boltzmann constant and @xmath45 is the vertical optical depth , with @xmath46 the opacity .
the two fluxes @xmath47 & @xmath48 are the radial heat fluxes due to turbulence and radiation respectively ; which are given by : @xmath49 where @xmath50 is the prandtl number ( in most calculations we take @xmath51 , the preferred value suggested by @xcite ) , and @xmath52 with @xmath53 the radiation constant , @xmath54 the speed of light and @xmath55 the mid - plane density . for our model
we adopt the @xcite opacities and equation  [ eqn : temp ] is integrated as described in @xcite .
our heating term ( equation  [ eqn : gamma ] ) does not account for stellar irradiation , which is known to dominate the heating in the outer disc ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
therefore , we do not allow the disc s temperature to drop below what would be obtained for a passively heated disc given by : @xmath56 we set @xmath57 to : @xmath58^{1/4}\ ] ] where the mass scaling is identical to that adopted by @xcite , based on a fit to a t tauri star s birth - line from @xcite .
finally , the temperature is also restricted to be no less than 10  k at very large radius , to account for heating from the surrounding star - forming region . in all calculations
we assume an ideal equation of state with @xmath59 and a mean molecular weight ( @xmath60 ) of 2.35 amu as the majority of the disc material we are interested in remains molecular and neutral .
we treat the viscosity using an ` alpha ' model such that : @xmath61 for mri turbulence we set @xmath62 and for angular momentum transport due to self - gravity we set @xmath63 to @xmath64 where @xmath65 is the toomre parameter @xcite given by : @xmath66 the form of @xmath67 is not important and is merely designed to keep @xmath68 ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) as suggested by simulations of self - gravitating discs ( e.g. * ? ? ?
* ; * ? ? ?
we treat the water component of the disc as a passive scalar species , such that it is advected along with the gas and can diffuse with respect to the gas due to turbulence .
thus , the evolution of the water component is given by ( e.g. * ? ? ?
* ; * ? ? ?
* ) : @xmath69+\dot{\sigma_i}\label{eqn : s_water}\end{aligned}\ ] ] where we evolve equation  [ eqn : s_water ] for the two water species h@xmath70o and hdo .
@xmath71 represents both infall onto the disc and the chemical evolution of the two species .
since the purpose of this work is to understand whether accretion bursts affect the evolution of the ( d / h)@xmath0 ratio rather than make detailed predictions , we do not attempt to model a full chemical network but adopt the simple model proposed by @xcite based on @xcite , where the ( d / h)@xmath0 ratio in water evolves as : @xmath72\label{eqn : water_chem}\ ] ] where @xmath73 is the gas number density at the disc s mid - plane , @xmath74 is the rate constant for the deuteration of water given by @xcite @xmath75 \ ] ] and ( d / h)@xmath76 is the equilibrium d / h ratio at a given temperature , which we take from @xcite assuming the d / h ratio in the background molecular hydrogen dominated disc is @xmath77 .
the latter values were found by @xcite to agree reasonably well with the values obtained from their more detailed chemical network . in all calculations we assume , like @xcite , that the ( d / h)@xmath0 ratio of the in - falling primordial water is @xmath78 .
the infall profile adopted depends on the evolutionary model ( steady infall vs. disc building calculations ) . for the steady infall problem
we use a gaussian infall profile which peaks at a radius @xmath79 and has standard deviation @xmath80 ( e.g. * ? ? ?
for the full evolutionary calculations , where we account for disc building from a uniformly rotating collapsing core , we follow @xcite and adopt the @xcite infall profile : @xmath81,\ ] ] where the instantaneous centrifugal radius ( @xmath82 ) of the infalling material is given by angular momentum balance : @xmath83 where @xmath84 is the core s rotation rate ( considered to be constant in space and time ) , @xmath85 is the temperature of the core and @xmath86 is the mass of the star & disc .
finally , we assume the core follows an isothermal inside - out collapse ( e.g. * ? ? ?
* ) and use a constant infall rate of @xmath87 until the core mass is fully depleted , where @xmath88 is the sound speed for gas at a temperature of @xmath89 in all calculations ( as in * ? ? ?
our numerical code is based on the method detailed in @xcite where advective fluxes are reconstructed using a second order method and a van - leer limiter .
equations  [ eqn : temp ] & [ eqn : gas ] are integrated explicitly using a scheme that is second order in space and first order in time using a staggered mesh and finite volume operators . in order to evolve the transport and chemistry of the water content
we make use of operator splitting : first , equation  [ eqn : s_water ] without the chemical evolution terms is integrated explicitly using a scheme that is second order in space and first order in time .
secondly the chemical evolution is evolved by integrating equation  [ eqn : water_chem ] using an implicit update , such that : @xmath90\end{aligned}\ ] ] such a scheme is stable for equation that are both stiff and non - stiff over a simulation domain and thus appropriate for our problem .
the simulations are initialised with the gas surface density set to the numerical floor of @xmath91 g  @xmath32 and the temperature set to 10  k. we use a non - uniform radial grid , using power - law spacing @xcite . at the inner boundary we use a zero - torque boundary condition for the gas surface density and a free outflow boundary for the water species . for our steady infall calculations
we use 200 cells spaced between 0.1 and 66.8 au and at the outer boundary we use a zero - flux boundary condition for both the gas and water species .
whereas for the disc building calculations we use 400 cells spaced between 0.047  au and 6685  au , at the outer boundary we use a zero shear boundary condition for the gas and a free - outflow boundary condition for the water .    in appendix
a , we benchmark our simplified scheme against the results of the more complicated chemical network used by @xcite and find very good agreement .
therefore , our simple scheme is appropriate for analysing the impact of accretion bursts on the d / h ratio in water .
in order to better understand the effects of accretion bursts we consider a setup with steady infall ; where the star s mass does not grow with time and is fixed to @xmath92  m@xmath7  .
such calculations have been previously used to study the properties of the accretion bursts themselves @xcite as they reach a steadily repeating limit cycle .
initially , only gas falls onto the disc and it reaches a steady - limit cycle after approximately 0.3 - 0.6  myr . after 1  myr
we then allow both water and gas to fall onto the disc and evolve the calculation for a further 3  myr .
we perform large sets of calculations , where in all sets we vary the infall rate from @xmath93 to @xmath94  m@xmath7  yr@xmath8  and the infall radius from 0.32 to 32  au .
we run one set of calculations with no dead - zone and thus no accretion bursts .
additionally , we run five further sets which include accretion bursts : three have a prandtl number of 10 ( the theoretically preferred value of @xcite ) and vary the schmidt number taking @xmath95 , 1.0 & 3.0 ; finally we run two sets with a schmidt number of 1.0 but prandtl numbers of 1 & 100 .      in figure
[ fig : dh_profile1 ] , we show the radial distribution of the d / h ratio in water 0.5 myr after the water has been introduced into the disc ( we plot the profile during quiescence as the d / h ratio shows a transient during the burst ) for a simulation with @xmath96  m@xmath7  yr@xmath8  and @xmath97  au in the six sets of calculations .
profile after 0.5  myr of evolution in a disc with a steady infall rate of @xmath96  m@xmath7  yr@xmath8  occurring at @xmath97  au .
the thin solid line shows the result with no dead zone ( and therefore no accretion bursts ) .
the thick lines show simulations with accretion bursts with schmidt numbers of ( 0.3 - dashed , 1.0 - solid & 3.0 - dot - dashed ) .
the thin dotted lines represent calculations with accretion bursts and a schmidt number of 1.0 but prandtl numbers of 1.0 ( black ) & 100.0 ( grey ) . ]
these profiles clearly show that accretion bursts have a strong impact on the d / h ratio distribution . with accretion
bursts the radius at which the d / h ratio reaches the fully equilibrated value of @xmath98 ( the `` equilibration radius '' ; where the equilibration time - scale is comparable to the transport time - scale c.f .
@xcite ) increases from roughly 1  au to 3  au .
additionally , it decreases the d / h ratio at large radius in the disc from the unprocessed value ( @xmath78 ) with no accretion bursts to @xmath99 with accretion bursts and a schmidt number of 0.3 .
finally , we note that the value of the prandtl number makes little difference .
this is likely due to the fact that while the prandtl number does affect the individual properties of an accretion burst ( e.g. duration ) it does not effect the duty - cycle of the burst significantly @xcite and hence the fraction of disc material processed by the accretion burst when averaged over a time - scale longer than the burst time - scale .    in figures  [ fig : space_time_nb ] & [ fig : space_time_b ] we show the space - time evolution of the temperature ( right - panel ) and ( d / h)@xmath0 ( middle - panel ) , along with the accretion rate evolution ( left - panel ) .
figure  [ fig : space_time_nb ] shows that the ( d / h)@xmath0 ratio quickly reaches a steady state distribution where the equilibration radius occurs at 1  au ( the temperature here is essentially due to viscous dissipation ) , in agreement with the analytic models of @xcite .
however , figure  [ fig : space_time_b ] shows that , in the presence of accretion bursts , the ( d / h)@xmath0 ratio never reaches a steady state ; but like the accretion rate reaches a steadily repeating pattern , where incoming material is accreted into the dead zone region with a ( d / h)@xmath0 @xmath100 , which is then equilibrated in the dead zone region when the disc goes through an accretion burst , which raises the temperature well above the equilibration temperature to @xmath101  k. thus , the dead zone is able to thermally process a larger fraction of the disc material , such that the equilibration radius is considerably further out than predicted by models that do not include accretion bursts . furthermore , as well as being carried by advection along with the gas , this newly equilibrated water in the dead zone region can diffuse out to large radius ( @xmath102  au ) in the disc lowering the ( time - independent ) ( d / h)@xmath0 ratio . in the case with accretion bursts ,
this is more efficient than without because equilibration occurs at a larger radius .              in figures  [ fig : panel_nb ] and [ fig : panel_sc1 ]
we show the ( d / h)@xmath0 ratio measured at different locations in the disc after 0.5 myrs of evolution for the whole range of @xmath103 and @xmath79 examined .
figure  [ fig : panel_nb ] shows the ( d / h)@xmath0 ratio at 1.0 , 3.0 , 10 & 30 au for the simulation set with a schmidt number of 1.0 and no accretion bursts .
this shows that at high accretion rates and small infall radii , most of the in - falling water is equilibrated before it can be transported to large radii resulting in a fully equilibrated disc . at low accretion rates
@xmath104  m@xmath7  yr@xmath8  the disc is too cold to equilibrate the water except at very small radii ( @xmath14  au ) and at intermediate values one gets a profile similar to that shown in figure  [ fig : dh_profile1 ] where the water is equilibrated inside some radius , the equilibration radius , which is an increasing function of mass accretion rate .
essentially , the warmer the disc region onto which the infalling matter arrives , the more equilibrated the water of the disc as a whole will be .    the simulation set with accretion bursts and schmidt number of 1.0
is shown in figure  [ fig : panel_sc1 ] with the ( d / h)@xmath0 ratio measured at 1.0 , 3.0 , 10.0 & 30.0 au after 0.5 myr of evolution .
this shows a similar trend with infall rate and infall radius as previously ; however , with accretion bursts it is much more pronounced .
similarly , the general trend of larger equilibration radius and lower ( d / h)@xmath0 ratio at large radius in the specific case discussed above holds for the majority of the parameter space . comparing these results with those found in figure  [ fig : panel_nb ] indicates that for the majority of the parameter space of interest the terrestrial region will have equilibrated water .
however , in order for the ( d / h)@xmath0 ratio to be reduced from primordial values at large radius requires infall rates @xmath105  m@xmath7  yr@xmath8  , which is expected to occur at early times ( e.g. * ? ? ?
however , the reduced ( d / h)@xmath0 ratio at large radius does appear to sensitively depend on the infall radius , with @xmath106  au resulting in fully equilibrated water at large radius : essentially all the incoming water is thermally processed in the dead - zone before it is transported to large radius .
whereas with @xmath107  au the ( d / h)@xmath0 ratio is only reduced from primordial values by a factor of @xmath108 , where the equilibrated water at large radius comes from the turbulent diffusion of material processed in the dead - zone that has then be transported to larger radius .
s@xmath8 and a schmidt number of unity . ]
we now consider calculations where we build the disc from a uniformly rotating core . initially the star has a mass of 0.2  m@xmath7  and the cloud core contains 0.8  m@xmath7  .
the cloud core then falls onto the disc as described in section  [ sec : infall ] . any material that falls inside the inner boundary
is automatically added to the star , whose mass is subject to an explicit update .
since we have seen above that the prandtl number makes little difference to the results , we only perform calculations with @xmath51 . we then perform a series of simulations where we vary the cloud s rotation rate @xmath84 between @xmath109 & @xmath110  s@xmath8 , the schmidt number between 0.1 & 10 and whether or not the disc experiences accretion bursts .      in figure
[ fig : mass_evolve ] , we show how the individual components ( star , cloud core & disc ) in our disc building calculations evolve during the simulation for a model with @xmath111  s@xmath8 .
the top panel shows a model with no accretion burst and the bottom panel shows a simulation with accretion bursts . in this case
infall takes approximately 0.2  myr .    in the case with
no accretion bursts the star rapidly approaches 1  m@xmath7  after infall ceases and the disc remains below 0.1  m@xmath7  for the majority of the evolution and steadily declines in a power - law fashion .
however , in the case with accretion bursts , the dead zone traps material in the disc .
this results in a more massive disc around 0.1 - 0.2  m@xmath7  for the majority of the evolution , which accretes material onto the star in rapid bursts .
this evolution is similar to that found previously ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) where discs that contain dead zones are generally more massive ( and thus hotter ) than discs that do not .            in figures  [ fig : db_st_nb ] & [ fig : db_st_b ]
we show the evolution of the accretion rate ( left panel ) , centrifugal radius and snow - line radius taken at @xmath112  k ( thick dashed line and thick solid line respectively in middle and right panels ) , ( d / h)@xmath0 ratio ( middle panel ) and disc temperature ( right panel ) as a function of time of the calculations shown in figure  [ fig : mass_evolve ] .
figure  [ fig : db_st_nb ] shows the space - time evolution for a calculation with @xmath111  s@xmath8 and no accretion burst .
the accretion rate evolution shows that once the infall onto the disc ceases the accretion rate drops dramatically ; this coincides with a large decrease in disc temperature , which then slowly declines with time .
the middle panel shows that the ( d / h)@xmath0 ratio has a plateau around 10  au and declines at both large and small radius .
this evolutionary track was seen by @xcite and is also shown in our benchmark calculation discussed in the appendix . at large radius
this disc is primarily composed of material that fell at small radius into the hot disc and then was transported to large radius by decretion ( e.g. * ? ? ?
* ; * ? ? ?
the water that fell at early times into the hot inner disc was equilibrated and then transported to large radius resulting in a ( d / h)@xmath0 ratio close to the equilibrated value at large radius ( although this radius increases with time ) . at small radius
the disc is hot enough to rapidly equilibrate the water , although the equilibration radius moves inwards with time as the disc temperature falls as the disc s mass decreases .
figure  [ fig : db_st_b ] shows the same calculation , but this time with accretion bursts . in this case ,
once infall ceases the accretion rate reaches @xmath15  m@xmath7  yr@xmath8  , which then rapidly burst to @xmath9 m@xmath7  yr@xmath8  occasionally . while the evolution of the temperature and ( d / h)@xmath0 ratio is similar at very large radius to the calculation without accretion bursts with a region of equilibrated water that is transported to large radius by decretion .
however , we see that the presence of a dead - zone and accretion bursts dramatically changes the ( d / h)@xmath0 ratio profile in the inner regions of the disc .
we also see an individual accretion burst has a transient effect that lasts @xmath1130.5  myr .
each accretion burst processes a large fraction of water in the disc , resulting in a lower peak ( d / h)@xmath0 ratio ( as also seen in the steady infall calculations ) and results in a equilibration radius that on average remains constant with time ( at @xmath114  au ) and does not move in with time like in the case without accretion bursts .
the evolution of the ( d / h)@xmath0 profile during an accretion burst is shown in figure  [ fig : pre_post ] .
profile just before ( 0.795  myr , solid ) and after ( 0.819  myr , dashed ) an accretion burst . shown for the simulation with @xmath115 and @xmath111  s@xmath8 where the accretion burst occurs at @xmath1130.8  myr ( see figure  [ fig : db_st_b ] ) . ]    where the solid line shows the ( d / h)@xmath0 profile just prior to an accretion burst and the dashed line shows the ( d / h)@xmath0 profile just after an accretion outburst .
we see that the outburst equilibrates water in the dead zone ( @xmath116  au ) , decreasing the overall ( d / h)@xmath0 ratio of the entire disc .
this newly equilibrated disc material is then transported by advection to smaller radius and to large radius by turbulent diffusion .
this diffusion of newly equilibrated water to large radius results in a lower ( d / h)@xmath0 ratio at large radii in the disc , explaining lower peak ( d / h)@xmath0 ratios compared to burst - free models .          in figure
[ fig : dh_omega ] , we show how the core s rotation rate effects the ( d / h)@xmath0 profile .
we plot the ( d / h)@xmath0 ratio for calculations with a schmidt number of 1 , both with ( dashed ) and without ( solid ) accretion bursts at 1 , 3 , 10 & 30  au after 1  myr of evolution .
these show that there is a critical rotation rate below which the entire disc is equilibrated of @xmath117  s@xmath8 for the calculations without accretion bursts , and @xmath118  s@xmath8 for the calculations with accretion bursts , corresponding to maximum centrifugal radii of @xmath5  au and @xmath119  au , respectively ( see equation ( [ rc ] ) ) . below these critical rotation rates all infalling material lands in regions of the disc where the temperature is above the equilibration temperature .
in the cases with accretion bursts at all rotation rates we find that at radii of @xmath5  au the water is equilibrated , unlike the case without accretion burst where above the critical core rotation rate it returns to near primordial values . at large radii
the trend seen from the steady infall calculations is also repeated here , where accretion bursts results in a lower ( d / h)@xmath0 ratio ; however , the effect is less pronounced .
we note the variability seen in the 1.0 & 3.0 au is indicative of the variability in the ( d / h)@xmath0 ratio with time , in the case with accretion bursts .
ratio as a function of radius for calculations with @xmath120  s@xmath8 for schmidt numbers of 0.1 ( dashed ) , 1.0 ( solid ) & 10.0 ( dot - dashed ) .
the top panel shows calculations without accretion bursts the bottom panel shows calculations with accretion bursts . ]    finally , in figure  [ fig : profile_sc ] we show the impact of the schmidt number on the ( d / h)@xmath0 profile , where we plot our the ( d / h)@xmath0 ratio as a function of radius in our calculations with @xmath121  s@xmath8 after 1  myr of evolution with schmidt numbers of 0.1 ( dashed ) , 1.0 ( solid ) & 10.0 ( dot - dashed ) .
the top panel shows the calculation without accretion bursts and the bottom panel shows the calculations with accretion bursts . at small schmidt numbers
there is enhanced transport of equilibrated water to large radius resulting in a lower ( d / h)@xmath0 ratio at large radius ( as seen by @xcite ) .
conversely , larger schmidt numbers result in sharper transitions from the equilibrated water to the primordial value outward as expected from the reduced efficiency of turbulent diffusion ( e.g. * ? ? ?
* ; * ? ? ?
our calculations have shown that the presence of dead zones in protoplanetary discs , which can give rise to accretion bursts have a significant impact on the evolution of the ( d / h)@xmath0 ratio in water .
this is not surprising since the ( d / h)@xmath0 ratio equilibrates at temperatures @xmath122  k @xcite which lies in - between the quiescent temperature of the dead - zone of @xmath123  k and the dead - zone temperature during an accretion burst of @xmath101  k. hence , un - equilibrated material that resides in the dead zone prior to an accretion burst is brought closer to equilibration by the higher disc temperatures .
our calculations have shown that the cloud core s rotation rate plays a strong role in setting the ( d / h)@xmath0 ratio at large radius in the disc .
thus , given the several order of magnitude spread in observed core rotation rates , with typical values @xmath124
s@xmath8 @xcite then one would expect to observe a large spread in the ( d / h)@xmath0 values measured at large radii in protoplanetary discs .
therefore , herschel ( e.g. * ? ? ?
* ) , alma or tmt measurements of the ( d / h)@xmath0 ratio would allow inferences about the initial angular momentum of the discs .
additionally , if the ( d / h)@xmath0 ratio could be spatially resolved then the properties of accretion bursts could be probed or spatial variations in the ( d / h)@xmath0 ratio then the probe could be used to place constraints on the schmidt number , similar to the method detailed by @xcite for spatial variations at various molecular snow lines .
furthermore , combined observations of the crystalline silicate fraction ( e.g. * ? ? ?
* ) and ( d / h)@xmath0 ratio would allow the outward transport of crystalline materials by decretion during disc building hypothesis ( e.g. * ? ? ? * ) to be tested . indeed discs with high crystallinity fractions at large radius should also possess lower ( d / h)@xmath0 ratios .
thus , future observations of the ( d / h)@xmath0 ratio should be able to probe the angular momentum distribution of the disc during the its building and formation .
additionally , spatially resolved observations of the ( d / h)@xmath0 ratio should be able to determine whether the majority of protoplanetary discs experienced accretion bursts during their early evolution .
to first order , the d / h systematics in the solar system may be subsumed in two salient features @xcite :    * most oort cloud comets show ( d / h)@xmath0 values of @xmath125 ( with some comets showing lower values ) , lower than expected in molecular clouds .
* carbonaceous chondrites display d / h ratios significantly lower than cometary ( around the terrestrial value @xmath126 ) .
both domains of variation are significantly higher than the protosolar value ( @xmath3 ) and significantly lower than the value observed in molecular clouds ( @xmath18 ) .
in our framework , fact ( i ) indicates that infall has occurred on relatively small radii to allow significant equilibration before outward transport , which , as mentioned in the previous subsection , would also account for the presence in comets of crystalline silicates , also deemed to have originated in the inner solar system @xcite . at the same time ,
the relatively large d / h ratio of cometary water compared to the protosolar value suggests that equilibration has not been complete .
judging from figure [ fig : dh_omega ] , that would at face value constrain the rotation rate to be within a factor of several of @xmath127  s@xmath8 ( with a factor of @xmath108 systematic difference between the outburst - free and outburst - laden cases ) .
of course , the infall model is very idealized , so the situation may not be as fine - tuned as it may seem ( whether with or without outbursts , incidentally ) .
it is also possible that comets , whose ( d / h)@xmath0 is somewhat variable , do not mark the _ peak _ ( d / h)@xmath0 value of the disc , and whatever that may be , this latter value may still be consistent with the lower end of those determined for prestellar cores ( e.g. * ? ? ?
the importance of fact ( ii ) may be grasped from examination of figure  [ fig : dh_temp ] which plots ( d / h)@xmath0 against the disc temperature for various disc models , along with the aforementioned solar system data .
the dashed line shows the curve for a simulation without outbursts @xmath128  s@xmath8 and a schmidt number of 1 , for which the maximum d / h ratio roughly matches that displayed by comets at large radii ( @xmath102  au ) .
it is seen that at the snow line ( taken here to correspond to a temperature of 170 k ) , ( d / h)@xmath0 is not very different from this maximum value , that is , in this case diffusion of equilibrated water from the inner solar system is too inefficient to induce a substantial deviation .
however , carbonaceous chondrites , which generally show ample evidence of aqueous alteration and therefore initial presence of ice , should have formed beyond the snow line , and yet their d / h ratio is much lower than this .
this @xcite pointed out as a general problem of standard disc models with @xmath129 .
however , @xcite noted that assuming a schmidt number of about 0.2 ( more efficient turbulent diffusion ) , would alleviate this discrepancy , as illustrated by their analytic profile ( in a quasi - static approximation ) plotted as a red dot - dashed line in the figure ( at higher temperatures ) there . ] . while non - mri and non - ideal turbulence may produce sc values @xmath14
( e.g. * ? ? ?
* ; * ? ? ?
* ) , it is unclear how realistic such low schmidt numbers are , and the required sc values could be yet lower if the d / h values displayed by carbonaceous chondrite clays were compromised by isotopic exchange with organic matter as suggested by @xcite .
if we now allow for bursts , the situation changes : the solid blue line in figure  [ fig : dh_temp ] stands for a simulation with bursts ( @xmath130  s@xmath8 ) also matching the cometary ( d / h)@xmath0 value at large radii , but this time matching the chondritic value at the snow line , despite a schmidt number of 1 .
this is because accretion bursts have been able to process water in the inner disc , affecting the snow line , but the ( d / h)@xmath0 value at large radii has not evolved in pace with it ( see also figure [ fig : pre_post ] ) .
so this may be an alternative way to account for the low d / h ratios of carbonaceous chondrites .
to summarize , in order to explain the distribution of ( d / h)@xmath0 ratios measured in carbonaceous chondrites and comets we need that either : ( i ) the turbulent transport occurs with a rather low schmidt number @xmath131 , and the disc did not experience accretion bursts @xcite ; or ( ii ) turbulent transport occurs with a nominal schmidt number of @xmath5 and the disc experiences several accretion bursts during the initial stages of its evolution .     ratio shown against temperature after 1  myr of evolution
compared to several solar - system constrains : protosolar ( e.g. * ? ? ?
* ) , carbonaceous chondrites ( e.g. * ? ? ?
* ) and the majority of oort cloud comets ( e.g * ? ? ?
* ; * ? ? ?
the snow - line at a temperature of @xmath132  k is marked with the thin vertical dotted line .
the thick solid line shows a simulation with @xmath115 , accretion bursts and @xmath111  s@xmath8 .
the thick dashed line shows a simulation with @xmath115 , no accretion bursts and @xmath128  s@xmath8 .
the thick dot - dashed line corresponds to the best fit solution of @xcite with @xmath133 . ]
while we have shown that accretion bursts can have a significant impact on the evolution of the ( d / h)@xmath0 ratio , our results can not be considered precise predictions .
rather , our simple models have allowed us to explore a large parameter space .
thus , we must discuss the impact our simplifications may have had on the evolution .
in particular , we adopted a rather simple model for the chemical evolution of the water and integrated one rate equation instead of a large chemical network ; however , we have shown our method agrees well with the more detailed calculations of @xcite , although reactions not considered by the latter ( e.g. ion - neutral ) could be important @xcite .
perhaps the biggest simplification is the use of a one - dimensional , one - layer model , where the evolution is controlled by the properties at the mid - plane of the disc .
@xcite performed a series of simulations where they considered the impact of the vertical temperature structure on the evolution of the ( d / h)@xmath0 ratio and found it to be negligible .
however , this may not be true for the layered structure in the dead - zone where the active surface layers that contain little mass , have a temperature different from the more massive dead zone @xcite and there may be turbulent mixing between the active and dead layers , with the possibility that the dead zone itself contains residual levels of turbulence ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
such a 2d structure would be worth taking into account .
we have also assumed throughout the paper that water was tightly coupled to the gas .
however , this may not be true at large distances ( where water is solid ) and late stages of the disc s evolution , especially if water becomes incorporated in growing aggregates of mm to m size ( e.g. * ? ? ?
. it might be _ a priori _ expected that the effect will be to make outward transport of equilibrated water more difficult @xcite but simulations including gas - solid drift , accretion and fragmentation would be worthwhile to address the effect on the ( d / h)@xmath0 ratio .    finally , since we have shown ( in agreement with @xcite ) that the radius and time at which infalling material reaches the disc plays an important role in the late time ( d / h)@xmath0 profile , then more realistic disc formation models should be considered that account for differential rotation and turbulence in the cloud core , as well as the two- ( and possibly three- ) dimensional nature of the infall onto the disc @xcite , along with the associated angular momentum mis - match that can occur between the keplerian disc rotation and the rotation rate of the infalling material .
we have used a 1d model of a layered protoplanetary disc that undergoes accretion bursts due to the gravo - magnetic limit cycle ( a commonly invoked mechanism to explain fu orionis outbursts ) to study their impact on the ( d / h)@xmath0 ratio .
we have performed two types of calculations ; firstly , a set of calculations with a steady infall rate onto the disc where we vary the infall rate , infall radius , schmidt number , prandtl number and whether or not the disc contains a dead zone .
secondly , we perform a series of calculation where we build the disc and star from a uniformly rotating isothermal cloud core , where we vary the rotation rate , schmidt number and whether the disc contains a dead zone . in all these calculations
we find that accretion burst play a major role in equilibrating water with the surrounding hydrogen gas .
thus , the ( d / h)@xmath0 ratio carries an imprint of whether the disc experienced accretion bursts or not .
our main findings are as follows :    1 .   during an accretion burst , water in the dead zone region
is equilibrated , lowering its d / h ratio to @xmath98 , as the dead zone temperature reaches @xmath101  k during the accretion burst , higher than the equilibration temperature .
discs that undergo accretion bursts have , in general , lower ( d / h)@xmath0 ratios as the burst thermally process a larger fraction of the disc material .
equilibrated water from the accretion burst can be carried to large radius by diffusion resulting in lower d / h ratios at 10 - 100s of au .
discs that undergo accretion bursts have larger equilibration radii ( 1 - 5  au ) compared to discs that do not experience accretion bursts . 4 .
the rotation rate of the parent cloud core plays an important role in setting the large scale ( d / h)@xmath0 profile as suggested by @xcite , indicating there should be a significant variation in the d / h ratio in protoplanetary discs .
agreement with the solar - system constraints requires , either : ( 1 ) the solar nebula experienced no accretion bursts and had a low schmidt number @xmath4 ; or ( 2 ) the solar nebula experienced accretion bursts and had a schmidt number close to `` nominal '' ( i.e. unity ) .
given the important role of accretion bursts in the evolution of ( d / h)@xmath0 , we suggest more detailed modelling of their impact on the chemical and physical evolution of protoplanetary discs be performed .
finally , coupling these models with future observations , the ( d / h)@xmath0 ratio and profile could be used as probe for the prevalence of accretion bursts in the earliest stages of protoplanetary disc evolution .
we thank the referee for comments which improved the paper .
we are grateful to phil armitage for useful discussions .
the numerical calculations were performed on the sunnyvale cluster at cita , which is funded by the canada foundation for innovation .
since we have adopted a simplified approach to the chemical evolution in view of the numerous chemical reactions that may govern evolution of the deuterium fraction in water , we need to test whether our method is appropriate for studying the effect of accretion bursts .
therefore , we repeat the ` standard model ' of disc evolution calculation of @xcite who include a much more complete chemical network with 14 species and 48 reactions .
the model parameters are described in table  3 of @xcite ( i.e. @xmath134  s@xmath8 , @xmath135  k & @xmath136 ) .
however , we do not simplify our temperature solver and opacity structure to that used by @xcite , but rather use the more detailed method described in section  [ sec:1dmodel ] .
the results of our calculation are shown every 0.05myr in the left hand panel of figure  [ fig : yang_bench ] and the results from the yang et al.s original calculation are shown in the right - hand panel at 0.1 , 0.2 , 0.3 , 0.5 & 1.0 myr .
we see that the two evolutionary paths agree very well , particularly at late times ( which we are primarily interested in from the point of view of observables ) .
this indicates that our simple method for calculating the evolution ( d / h)@xmath0 ratio produces results comparable with the more detailed chemical method used by @xcite .
this is not surprising since @xcite performed a detailed exploration of the different approaches for the chemical evolution and found they gave similar results .
therefore , we are confident that our method is suitable for exploring the effect of accretion bursts on the d / h ratio in water .
ceccarelli c. , caselli p. , bockelee - morvan d. , mousis o. , pizzarello s. , robert f. , semenov d. , 2014 , to appear in protostars and planets vi , eds
beuther , r. klessen , c. dullemond , th .
henning , university of arizona press ( arxiv:1403.7143 )


unstable particles are a perfect testing ground for theories that try to unify the principle of relativity with quantum mechanics for ( at least ) two reasons .
first , an unstable particle provides a very simple example of a non - trivial interacting quantum system . due to the specific initial condition characteristic to this problem ( there is just one particle in the initial state ) , a rigorous description of the decay is possible in a small hilbert space that contains only states of the particle and its decay products ,
so the solution can be obtained in a closed form .
second , unstable particles played an important role in confirming predictions of einstein s special relativity .
it was demonstrated experimentally that the decay of moving particles slows down in a good agreement with einstein s time dilation formula @xcite .
for a long time it was believed that relativistic quantum mechanics and quantum field theory must reproduce exactly the einstein s time dilation formula for unstable particles @xcite .
however , a detailed quantum relativistic calculation in @xcite has shown that there are small corrections to this formula in the case of particles with definite momentum .
this result was later confirmed in other studies @xcite . although surprising , these findings did not challenge directly the applicability of special relativity to unstable systems .
this is because relativistic transformations connect results of measurements of two inertial observers moving with constant velocities with respect to each other .
then a fair comparison with the time dilation formula requires consideration of unstable states having definite values of velocity for both observers .
such states do not have a definite momentum for , at least , one observer .
their decay laws were considered in recent article @xcite which alleges that the decay accelerates when particle s velocity increases , instead of being slowed down as experiment shows . in the present article
we resolve this controversy by performing a detailed calculation of the decay laws of particles with narrow distributions of velocities observed from moving reference frames .
we found that the result of @xcite is based on an incorrect identification of the subspace of states of the unstable particle in the full hilbert space of the system . in sections [ ss : general ] - [ ss :
instant ] we present a rigorous quantum relativistic framework required for the description of decays . the exact formula for the time dependence of the non - decay probability in a moving inertial frame of reference
is derived in section [ sc : moving - frame ] .
particular cases of this formula relevant to unstable particles with sharply defined momenta or velocities are considered in sections [ sc : mov ] and [ sc : general ] , respectively .
section [ sc : numerical ] is devoted to numerical calculations of the differences between the accurate quantum mechanical result and the standard time dilation formula ( [ eq:1b ] ) .
although these differences are much smaller than the resolution of modern experiments , their presence is in sharp contradiction with special relativity which does not tolerate even small deviations from the einstein s time dilation .
the decay of unstable particles is described mathematically by the _ non - decay probability _ which has the following definition .
suppose that we have a piece of radioactive material with n unstable nuclei prepared simultaneously at time @xmath0 and denote @xmath1 the number of nuclei that remain undecayed at time @xmath2 . then the non - decay probability indicates the rapidity of the observer that measures the non - decay probability , as described below . ]
@xmath3 ( also called the _ non - decay law _ in this paper ) is defined as the fraction of nuclei that survived the decay in the limit of large @xmath4    @xmath5    so , at each time point the piece of radioactive material can be characterized by its _ composition _
@xmath6 , where @xmath7 is the share of undecayed nuclei in the sample and @xmath8 is the share of nuclei that were transformed to the _ decay products_. in this paper , in the spirit of quantum mechanics , we will treat @xmath4 unstable particles as an ensemble of identically prepared systems and consider @xmath3 as a property of a single particle ( nucleus ) , i.e. , the probability of finding this particle in the undecayed state .
in this paper we are concerned with comparison of decay observations made by two observers @xmath9 and @xmath10 that move with respect to each other . without loss of generality
we will assume that observer @xmath10 moves with respect to @xmath9 with velocity @xmath11 along the @xmath12-axis and that at time @xmath0 measured by both observer s clocks the origins of their coordinate systems coincide and all three pairs of coordinate axes are parallel ( @xmath13 , etc . ) .
for discussion purposes we will say that @xmath9 is at rest while @xmath10 is moving .
to simplify formulas we will use the _ rapidity _
parameter @xmath14 instead of velocity @xmath15 .
for example , in this notation the famous relativistic factor takes the form @xmath16 .
if from the point of view of @xmath9 the decay is described by the function @xmath3 , then the non - decay law of the same particle seen by the moving observer @xmath10 will be denoted by @xmath17 .
denotes time measured by the clock belonging to the observer @xmath10 .
] calculation of this function and , in particular , the relationship between @xmath17 and @xmath3 is the major goal of this work .
the purpose of any relativistic theory is to describe the relationships between measurements made by different _
inertial observers _ or from different _ inertial reference frames _ , i.e. , reference frames moving in space with constant velocities along straight lines and without rotation .
the special _ principle of relativity _ tells us that these reference frames are exactly equivalent    * * postulate i. * experiments identically arranged and performed in two different inertial reference frames @xmath9 and @xmath18 always yield the same results .
for any two observers @xmath9 and @xmath18 there is an _ inertial transformation _ that connects @xmath18 to @xmath9 , i.e. , the set of rules that allows one to change from the reference frame @xmath9 to the reference frame @xmath18 . each inertial transformation is a combination of space and time translations , rotations and boosts .
a composition of two inertial transformations is again a valid inertial transformation .
this composition obeys the associativity law .
the inverse of an inertial transformation is also an inertial transformation .
therefore , they form a 10-parameter lie group . the structure of this group is specified in the second postulate    * * postulate ii .
* inertial transformations ( space and time translations , rotations , and boosts ) form the _ poincar group_.    obviously , two different observers obtain different results by measuring observables of the same physical system .
most problems in physics can be understood as translations of descriptions of the physical system between different reference frames .
for example , if we have a full description of the system from the point of view of observer @xmath9 , then the time evolution is obtained by answering the question `` what is the description of the same system from the point of view of observer @xmath19 that is shifted in time with respect to @xmath9 ? ''
it is important to realize that inertial transformations between observers can be divided into two groups : _ kinematical _ and _ dynamical_. kinematical transformations are those whose action on observables does not depend on the interaction .
for example , if two stationary observers look at the unstable particle from different points in space , they would assign the same non - decay probability to the particle .
the same is true for stationary observers having different orientations in space . from this
follows , these forms are not appropriate for the description of non - decay laws of unstable particles . ]    * * postulate iii .
* space translations and rotations are kinematical .
on the other hand , time translations produce non - trivial changes in the system .
the composition of the unstable system looks different for observers @xmath9 and @xmath19 shifted in time with respect to each other .
the exact action of time translations on observables of the physical system should be obtained as a result of solution of dynamical equations which depend on the interaction acting in the system .
hence the following postulate is true for all types of isolated interacting systems .
* * postulate iv .
* time translations are dynamical , i.e. , interaction - dependent .
the above postulates i - iv have overwhelming experimental support .
they are assumed to be valid throughout this paper .
however they are not sufficient for a full description of the unstable particle in a moving reference frame .
such a description requires also knowledge of the nature of boosts .
here we have a choice between two paths forward .
one path is to postulate certain properties of boosts .
this path was taken by einstein .
it leads to special relativity and to the time dilation formula ( [ eq:1b ] ) as explained in the next section . in this paper
we will argue in favor of choosing another path : keep postulates i - iv , add to them well - established postulates of quantum mechanics , and see what are the implications for the transformations of observables ( in particular , the non - decay probability ) with respect to boosts .
this approach is employed starting from section [ ss : general ] throughout the paper .
in addition to postulates i - iv , einstein s special relativity makes two more assumptions regarding the nature of boosts    * * assumption v. * boosts are kinematical .    by postulating the kinematical character of boosts ,
special relativity insists that the internal composition of a compound system does not depend on the velocity of the observer .
for example , if @xmath9 and @xmath10 are two observers moving with respect to each other , then , according to special relativity , both observer will measure the same composition of the unstable system at time @xmath0 .
these statements are often considered as self - evident in discussions of special relativity .
for example , r. polishchuk writes in @xcite _ `` any event that is `` seen '' in one inertial system is `` seen '' in all others .
for example if observer in one system `` sees '' an explosion on a rocket then so do all other observers . '' _ in addition to the universality and interaction - independence of boost transformations , special relativity also postulates the exact transformation laws of physical observables with respect to boosts .
they are referred to as _
lorentz transformations_. the most fundamental are lorentz transformations for space - time coordinates of events .
loses its localization from the point of view of the moving observer @xmath10 @xcite .
so , assumption vi can be applied only in the classical limit . ]    * * assumption vi . *
if from the point of view of observer @xmath9 an event is localized in a space point @xmath20 at time @xmath21 , then from the point of view of observer @xmath10 the same event has space - time coordinates @xmath22 given by lorentz formulas + @xmath23    let us now demonstrate that in special relativity the postulates i - iv and assumptions v - vi are sufficient to unambiguously describe particle decay in different reference frames without invoking any information about the interaction governing the decay .
suppose that from the point of observer @xmath9 the unstable system is prepared in the state with composition @xmath24 at rest in the origin @xmath25 at time @xmath0 .
then observer @xmath9 may associate the space - time point    @xmath26    with the event of preparation .
in accordance with the dynamical character of time translations , the non - decay probability @xmath3 decreases with time . from experiment and quantum mechanical calculations ( see , e.g. , section [ sc : numerical ] ) it is known that the non - decay law has an ( almost ) exponential shape ) here . ]
@xmath27    where @xmath28 is the _ lifetime _ of the unstable particle . at time
@xmath29 the non - decay probability is exactly @xmath30 , so that the composition is @xmath31 .
this `` one lifetime '' event has coordinates    @xmath32    according to the observer @xmath9 .
let us now take the point of view of the moving observer @xmath10 .
due to the assumption v , both @xmath9 and @xmath10 agree that at the `` preparation '' and `` lifetime '' events the non - decay probabilities are 1 and @xmath30 , respectively . however , observer @xmath10 may not agree with @xmath9 about the space - time coordinates of these events . substituting ( [ eq : event1 ] ) and ( [ eq : event2 ] ) in ( [ eq : lorentz - time ] ) - ( [ eq : lorentz - x ] ) we see that from the point of view of @xmath10 , the `` preparation '' event has coordinates @xmath33 , and the `` lifetime '' event has coordinates @xmath34 . therefore , the time elapsed between these two events is @xmath35 times longer than in the reference frame @xmath9 .
this also means that the decay law of the particle is exactly @xmath35 slower from the point of view of the moving observer @xmath10 .
this is reflected in the famous `` time dilation '' formula    @xmath36    this formula was confirmed in numerous experiments @xcite , most accurately for muons accelerated to relativistic speeds in a cyclotron @xcite .
these experimental findings were certainly a triumph of einstein s theory .
however , as we see from the above discussion , in special relativity eq .
( [ eq:1b ] ) can be derived only under two assumptions v and vi , which lack proper justification .
therefore , a question remains whether eq .
( [ eq:1b ] ) is a fundamental exact result or simply an approximation that can be disproved by more accurate measurements ?
our goal in the rest of this paper is to demonstrate that this classical result of special relativity is not exact . in section [ sc : numerical ]
we will calculate corrections to the formula ( [ eq:1b ] ) .
let us now turn to the description of an isolated unstable system from the point of view of relativistic quantum theory . in our approach
we will keep postulates i - iv .
however , we are not going to use assumptions v and vi .
we will consider a model theory with particles @xmath37 , @xmath38 , and @xmath39 , so that particle @xmath37 is massive spinless and unstable , while its decay products @xmath38 and @xmath39 are stable and their masses satisfy the inequality    @xmath40    which makes the decay @xmath41 energetically possible . in order to simplify calculations and avoid being distracted by issues that are not relevant to the problem at hand we neglect the spin of the particle @xmath37 and assume that there is only one decay mode of this particle , i.e. , into two decay products @xmath38 and @xmath39 . for our discussion ,
the nature of the particle @xmath37 is not that important .
for example , this could be a muon or a radioactive nucleus or an atom in an excited state .
when measuring the non - decay law , experimentalists simply count the number of particles and determine their types .
for example , observations of a muon may result in only two outcomes .
one can find either a non - decayed muon or its decay products ( an electron , a neutrino , and an antineutrino ) .
thus the number of particles is a legitimate observable , and for our model system @xmath42 we can introduce hermitian operators for the number of particles @xmath43 , @xmath44 , and @xmath45 . since these observables can be measured simultaneously , we can assume that these three operators commute with each other .
two combinations of their common eigenvalues are allowed in our system :    @xmath46    hence the hilbert space of the unstable system should be represented as a direct sum of two orthogonal subspaces , and @xmath47 are allowed to take any values from zero to infinity
. however , for most unstable particles the interaction between the decay products in the final state can be ignored , and considering the subspace ( [ eq : tensor - product ] ) of the full fock space is a reasonable approximation . ]
@xmath48    where @xmath49 is the subspace of states of the unstable particle @xmath37 which corresponds to the set of eigenvalues ( [ eq : eigenval1 ] ) , and @xmath50 is the orthogonal subspace of the decay products which corresponds to the set of eigenvalues ( [ eq : eigenval2 ] ) .
we can now introduce a hermitian operator @xmath51 ( also known as `` yes - no experiment '' ) that corresponds to the observable `` particle @xmath37 exists '' .
the operator @xmath51 can be fully defined by its eigensubspaces and eigenvalues .
when a measurement performed on the unstable system finds it in a state corresponding to the particle @xmath37 ( i.e. , the state vector is within @xmath52 ) , the value of @xmath51 is 1 . when the decay products @xmath53 are observed ( the state vector lies in @xmath54 ) , the value of @xmath51 is 0 .
apparently , @xmath51 is a projection operator on the subspace @xmath52 . for each normalized state vector @xmath55
the probability of finding the unstable particle @xmath37 is given by the expectation value of the observable @xmath51 @xcite    @xmath56    in relativistic quantum mechanics , the dynamics of the system is described by a unitary representation @xmath57 of the poincar group in the hilbert space @xmath58 @xcite .
it is convenient to express representatives @xmath57 as exponential functions of _ generators _ @xmath59 , and @xmath60 .
then a general unitary operator from the set @xmath61 can be always written as a product    @xmath62    where @xmath63 is the rotation vector , @xmath64 is the velocity of the boost , @xmath65 is the vector of space translation , and @xmath21 is the amount of time translation .
generators are hermitian operators and correspond to certain _ total _
observables of the system .
space translations are generated by the vector of the total linear momentum @xmath66 .
the generator of rotations @xmath67 is the operator of the total angular momentum .
the generator of time translations @xmath68 is the hamiltonian ( total energy ) .
the generator of boosts @xmath69 is called the _
boost operator_. these generators must satisfy the commutation relations of the poincar group lie algebra . in this paper
we will need , in particular , the following commutators    @xmath70 = -i\frac{\hbar}{c^2 } h \delta_{ij }   \label{eq:5.55}\ ] ]    @xmath71 = -i\hbar p_i   \nonumber\ ] ]    where @xmath72 .
they imply the following useful relationships    @xmath73 \label{eq:6.2 } \\
e^ { \frac{i}{\hbar } \mathbf{k } c \vec{\theta } } h e^ { -\frac{i}{\hbar } \mathbf{k } c \vec{\theta } } & = &   h \cosh
\theta + c ( \mathbf{p } \cdot \frac{\vec{\theta}}{\theta } ) \sinh \theta \label{eq:6.3}\end{aligned}\ ] ]    the representation ( [ gen - elem ] ) allows us to relate results of measurements in different reference frames .
let us first take the point of view of the observer @xmath9 and consider a vector @xmath74 the first argument is the rapidity parameter of the reference frame from which the state @xmath75 is observed and the second argument is the time of observation .
this is consistent with the convention adopted for @xmath17 in section [ intro ] . ]
that describes a state in which the unstable particle @xmath37 is found with 100% certainty .
@xmath76    then the time evolution of the state vector @xmath77 in the reference frame @xmath9 is given by    @xmath78    and the non - decay law is given by    @xmath79    from this equation it is clear that the hamiltonian @xmath68 describing the unstable system should not commute with the projection
@xmath51    @xmath80 \neq 0 \label{eq:44b}\end{aligned}\ ] ]    otherwise , the subspace @xmath52 of states of the particle @xmath37 would be invariant with respect to time translations and the particle @xmath37 would be stable ( @xmath81 for all @xmath21 ) .    the moving observer @xmath10 describes the initial state ( at @xmath0 ) by the vector    @xmath82    the time dependence of this state is    @xmath83    then the non - decay law from the point of view of @xmath10 is    @xmath84    where the last equation follows from the property @xmath85 of the operator @xmath51 .
before calculating the non - decay law in the moving reference frame ( [ eq:3b ] ) , let us first consider a simpler case when the interaction responsible for the decay is `` turned off '' .
then the representation of the poincar group @xmath86 acting in @xmath58 is _ non - interacting_. this representation is constructed in accordance with the structure of the hilbert space ( [ eq : tensor - product ] ) as    @xmath87    where @xmath88 are unitary irreducible representations of the poincar group corresponding to particles @xmath37 , @xmath38 , and @xmath39 , respectively .
the generators of the representation ( [ eq : non - int ] ) are denoted by @xmath89 , @xmath90 , @xmath91 , and @xmath92 .
the operator of non - interacting _ mass _    @xmath93    commutes with @xmath89 , @xmath90 , @xmath91 , and @xmath92 . according to ( [ eq : masses ] ) , the operator @xmath94 has a continuous spectrum in the interval @xmath95 and a discrete point @xmath96 embedded in this interval .    from definition ( [ eq : non - int ] )
it is clear that the subspaces @xmath49 and @xmath97 are invariant with respect to @xmath86 .
moreover , the particle number operator @xmath43 and the projection operator @xmath51 commute with the non - interacting generators    @xmath98 = [ t,\mathbf{j}_0 ] = [ t,\mathbf{k}_0 ] =   [ t , h_0]=0 \label{eq : t - comm}\end{aligned}\ ] ]    this implies that the particle @xmath37 is stable with respect to time translations and boosts , as expected    @xmath99    the primary reason for considering the representation @xmath86 is that it allows us to build a convenient basis in the subspace @xmath52 .
let us denote @xmath100 a vector in @xmath52 that corresponds to the particle @xmath37 with zero momentum .
@xmath101    then we find that the vector    @xmath102    describes the particle @xmath37 with definite momentum @xmath103 . indeed , using ( [ eq:6.2 ] ) we obtain    @xmath104 ) | \mathbf{0 } \rangle \\ & = &   \frac{\vec{\theta}}{\theta } m_ac \sinh \theta   e^{-\frac{ic}{\hbar}\mathbf{k}_0   \vec{\theta } }   | \mathbf{0 } \rangle \\ & = & \frac{\vec{\theta}}{\theta }   m_ac \sinh \theta   | \mathbf{p } \rangle \end{aligned}\ ] ]    since particle @xmath37 is spinless by our assumption , the eigenvectors @xmath105 of the total momentum operator @xmath89 form a full basis in the subspace @xmath52 , so that    @xmath106    then any state @xmath107 of the particle @xmath37 can be represented by a linear combination of these basis vectors    @xmath108    where @xmath109 is the wave function in the momentum representation . in order to ensure the normalization @xmath110 ,
the function @xmath111 must satisfy    @xmath112    apparently , vectors @xmath113 themselves are not normalized . if we want to study the non - decay law of a state with definite momentum @xmath114 , we should use a state vector ( which we denote by @xmath115 to distinguish it from @xmath116 ) that has a normalized momentum - space wave function sharply localized near @xmath114 . in order to satisfy eq .
( [ eq : normalization ] ) such a wave function may be formally represented as a square root of the dirac s delta function @xmath117    the action of boosts on the basis vectors @xmath113 and @xmath118 is known from wigner s theory of irreducible unitary representations of the poincar group ( see , e.g. , @xcite )    @xmath119    where    @xmath120    using eq .
( [ eq:129a ] ) and the property    @xmath121    we can find boost transformations for an arbitrary state vector of the form ( [ eq : expansion ] )    @xmath122    these transformations can be viewed as transformations of the corresponding momentum - space wave function , e.g. ,    @xmath123    so that the boost operator can be represented as a differential operator in the momentum space    @xmath124    for further calculations we will need to define the _ newton - wigner position _
operator @xcite in @xmath58    @xmath125    which has the property    @xmath126 & = & i\hbar \delta_{ij } \label{eq : heis}\end{aligned}\ ] ]    that can be verified by direct substitution . according to eq . (
[ eq:7.21 ] )    @xmath127    therefore operator @xmath128 acts as a translation operator in the momentum space . in particular , we can write    @xmath129    where @xmath130 . on the other hand , applying the boost transformation ( [ eq : boost - tr ] ) to the momentum eigenfunction we obtain is the jacobian of the transformation @xmath131 . ]
@xmath132    this suggests that momentum eigenvector ( [ eq : mom - eigen2 ] ) has another useful representation    @xmath133
let us now `` turn on '' the interaction responsible for the decay and discuss the interacting representation @xmath57 of the poincar group in @xmath58 with generators @xmath66 , @xmath67 , @xmath69 , and @xmath68 . according to our postulates iii and iv , the generators of space translations and rotations are interaction - free ,    @xmath134    while the generator of time translations ( the hamiltonian @xmath68 ) contains an interaction - dependent term @xmath135 .
@xmath136    from the commutator ( [ eq:5.55 ] ) it then follows @xcite that the generators of boosts must be interaction - dependent as well    @xmath137    where @xmath138 .
this means that we are working in the so - called dirac s _ instant form of dynamics_. we will further assume that the interacting representation @xmath57 belongs to the bakamjian - thomas form of dynamics @xcite , which is characterized by the property that the interacting operator of mass @xmath139 commutes with the newton - wigner position operator . ]
@xmath140 = 0 \label{eq : nw - comm}\end{aligned}\ ] ]    our next goal is to define the basis of common eigenvectors of commuting operators @xmath89 and @xmath141 in @xmath58 .
, there are other independent operators in the mutually commuting set containing @xmath89 and @xmath141 .
these are , for example , the operators of the square of the total angular momentum @xmath142 and the projection of the total angular momentum on the @xmath12-axis @xmath143 .
therefore a unique characterization of any basis vector requires specification of all corresponding quantum numbers as @xmath144 .
however these quantum numbers are not relevant for our discussion and we omit them . ]
these eigenvectors must satisfy conditions    @xmath145    they are also eigenvectors of the interacting hamiltonian @xmath146    @xmath147    where @xmath148 . that depends on the eigenvalue @xmath149 of the interacting mass operator and @xmath150 in eq .
( [ eq : omega ] ) that depends on the fixed value of mass @xmath96 of the particle @xmath37 . ] in the zero - momentum eigensubspace of the momentum operator @xmath89 we can introduce a basis @xmath151 of eigenvectors of the interacting mass @xmath141    @xmath152    then the basis @xmath153 in the entire hilbert space @xmath58 can be built by formula ( cf .
( [ eq : mom - eigen2 ] ) )    @xmath154 .
these eigenvectors are normalized to delta functions    @xmath155    the actions of inertial transformations on these states are well - known @xcite .
in particular , for boosts along the @xmath12-axis ( cf .
( [ eq:129a ] ) ) and time translations we obtain    @xmath156    where    @xmath157
next we notice that due to eqs .
( [ eq : heis ] ) and ( [ eq : nw - comm ] ) vectors @xmath158 also satisfy eigenvector equations ( [ eq : mom - eigen ] ) - ( [ eq : mass - eigen ] ) , so they must be proportional to the basis vectors @xmath159    @xmath160    where @xmath161 is an unimodular factor . unlike in ( [ eq : alpha ] ) , we can not conclude that @xmath162 , because , generally , the action of @xmath163 on eigenvectors @xmath164 involves multiplication by a unimodular scalar in addition to the shift of momentum .
however , if the interaction is not pathological we can assume that the factor @xmath161 is smooth , i.e. , without rapid oscillations .
this property will be used in derivation of eq .
( [ eq : appr_dec ] ) .
obviously , vector @xmath165 can be expressed as a linear combination of zero - momentum basis vectors @xmath166 , so we can write and @xmath39 , i.e. , discrete eigenvalues of @xmath141 below @xmath167
. then the spectrum of @xmath141 ( similar to the spectrum of @xmath94 ) is continuous in the interval @xmath95 , and integration in ( [ eq : integral ] ) should be performed from @xmath168 to infinity . ]
@xmath169    where function @xmath170 describes the _ mass distribution _ of the unstable particle . for most realistic unstable systems the mass distribution @xmath171 has the breit - wigner form ( see fig .
1 ) ) could be different from the mass @xmath96 that particle @xmath37 has in the absence of decay interaction .
however , we will disregard this possibility here . ]
@xmath172    where parameter @xmath173 is a normalization factor required to ensure that the mass distribution ( [ eq:63 ] ) - ( [ eq:64 ] ) is normalized to unity .
we now use eqs .
( [ eq : alpha ] ) and ( [ eq : integral ] ) to expand the vector @xmath113 in the basis @xmath174    @xmath175    then , from ( [ eq:28a ] ) we obtain a useful formula    @xmath176
suppose that vector @xmath75 in ( [ eq : expansion ] ) describes a state of the unstable particle @xmath37 from the point of view of the observer @xmath9 .
the time dependence of this state seen from the moving reference frame @xmath10 is obtained by applying eqs ( [ eq:44a ] ) , ( [ eq:29a ] ) , ( [ eq:29b ] ) , and ( [ eq : p - expansion ] )    @xmath177    the inner product of this vector with @xmath178 is found by using ( [ eq:35a ] )    @xmath179    introducing new integration variables @xmath180 and taking into account ( [ eq : invar ] ) , this equation can be rewritten as    @xmath181    the non - decay probability in the reference frame @xmath10 is then found by substituting ( [ eq : projector ] ) in eq .
( [ eq:3b ] )    @xmath182    which is an exact formula valid for all values of @xmath183 and @xmath21 .
in the reference frame at rest ( @xmath184 ) , formula ( [ eq:50 ] ) simplifies    @xmath185    let us consider a particular case of this expression when the unstable particle has a well - defined momentum , i.e. , described by the state vector @xmath186 . if the particle @xmath37 is in the state @xmath187 with zero momentum then @xmath188 and from the above equation we obtain    @xmath189    which is the standard formula for the non - decay law of a particle at rest . , this formula can be compared , for example , with eq .
( 3.8 ) in @xcite . ] in particular , if we substitute ( [ eq:63 ] ) - ( [ eq:64 ] ) for the mass distribution @xmath171 we obtain approximate exponential decay ( [ eq : expo ] ) where the lifetime is given by    @xmath190    if the particle @xmath37 has a definite non - zero momentum @xmath191 , then    @xmath192    and    @xmath193    in a number of works @xcite it was noticed that if one interprets the state @xmath118 ( where @xmath194 ) as a state of unstable particle moving with definite speed @xmath195 , then the decay law of the moving particle ( [ eq:35 ] ) can not be connected with the non - decay law of the particle at rest ( [ eq:34 ] ) by einstein s formula ( [ eq:1b ] ) , i.e. ,    @xmath196    this observation prompted authors of @xcite to question the applicability of einstein s special relativity to particle decays .
however , at a closer inspection it appears that this result does not contradict einstein s time dilation formula ( [ eq:1b ] ) directly .
formula ( [ eq:1b ] ) refers to observations made on the same particle from two frames of reference moving with respect to each other .
if from the point of view of observer @xmath9 the particle is described by the state vector @xmath197 which has zero momentum and zero velocity , then from the point of view of @xmath10 this particle is described by the state    @xmath198    which is an eigenstate of the velocity operator[multiblock footnote omitted ]    @xmath199    but is not an eigenstate of the momentum operator @xmath89 .
therefore , strictly speaking , its non - decay law is not described by @xmath200 . in order to compare with einstein s formula ( [ eq:1b ] )
, we need to calculate the decay law in the moving frame of reference @xmath201 .
this is done in section [ sc : general ] .
let us now discuss the non - decay probability @xmath202 at @xmath0 in the moving frame of reference . instead of eq .
( [ eq:50 ] ) , it is more convenient to use the general definition ( [ eq:3d ] ) which expresses @xmath203 as the square of the norm of the projection of @xmath204 on the subspace @xmath52 .
we are going to prove that for a normalized @xmath77 this probability can not be equal to 1 for all non - zero @xmath183 .
suppose that this statement is wrong , so that ( in agreement with assumption v ) for any @xmath205 and any @xmath206 , the vector @xmath207 belongs to @xmath52
. then the subspace @xmath52 is invariant under action of boosts @xmath208 and operator @xmath209 commutes with the projection @xmath51 .
then from the poincar commutator ( [ eq:5.55 ] ) and @xmath210= 0 $ ] it follows by jacobi identity that    @xmath211 & = & \frac{ic^2}{\hbar}[t,[k_z,(p_0)_z ] ] \\ & = &   \frac{ic^2}{\hbar } [ k_z,[t , ( p_0)_z ] ] - \frac{ic^2}{\hbar}[(p_0)_z,[t , k_z ] ] \\ & = & 0\end{aligned}\ ] ]    which contradicts eq .
( [ eq:44b ] ) .
this contradiction implies that the state @xmath207 does not correspond to the particle @xmath37 with 100% probability ; this state must contain contributions from the decay products even at @xmath0    @xmath212    this is the `` decay caused by boost '' .
thus we found that assumption v from section [ ss : relativity ] is not correct .
sometimes this assumption is formulated as @xcite _ `` flavor is the quantum number that distinguishes the different types of quarks and leptons .
it is a lorentz invariant quantity .
for example , an electron is seen as an electron by any observer , never as a muon . '' _ although the above statement about the electron is correct ( because the electron is a stable particle ) , this is not true about the muon which , according to ( [ eq : not - in ] ) can be seen as a single particle by the observer at rest and as a group of three decay products ( an electron , a neutrino , and an antineutrino ) by the moving observer .
unfortunately exact evaluation of eq .
( [ eq:50 ] ) for @xmath213 is not possible , so we need to make approximations .
let us discuss properties of the initial state @xmath77 in more detail .
first , in all realistic cases this state is not an eigenstate of the total momentum operator , so the wave function is not localized at one point in the momentum space ( as assumed in ( [ eq : sqrt ] ) ) but has a spread ( or uncertainty ) of momentum @xmath214 and , correspondingly an uncertainty of position @xmath215 . on the other hand , the state @xmath216 is not an eigenstate of the mass operator @xmath141 , and @xmath77 is characterized by the uncertainty of mass @xmath217 ( see fig .
1 ) that is related to the lifetime of the particle @xmath218 by formula ( [ eq : lifetime ] ) .
it is important to note that in all cases of practical interest the mentioned quantities are related by inequalities    @xmath219    in particular , the latter inequality means that the uncertainty of position is mush less than the distance passed by light during the lifetime of the particle .
for example , in the case of muon @xmath220s and , according to ( [ eq : gg2 ] ) , the spread of the wave function in the position space must be much less than 600 m , which is a reasonable assumption .
therefore , we can safely assume that the factor @xmath171 in ( [ eq:50 ] ) has a sharp peak near the value @xmath221 .
then we can move the value of the smooth function @xmath222 at @xmath223 outside the integral by @xmath149 . )
that the spread of the momentum - space wave function is much less than @xmath224 we could have moved the value of the smooth function @xmath225 outside the integral by @xmath149 .
this would mean that the non - decay law in the moving frame of reference is controlled by the spread of the particle wavefunction in the momentum space rather than by its mass uncertainty which disagrees with experimental observations . ]
@xmath226    now we are going to use formula ( [ eq : appr_dec ] ) to calculate the decay law for an approximate eigenstate of momentum in @xmath52 whose wave function @xmath227 is localized near zero momentum @xmath228 , though inequality ( [ eq : gg ] ) still remains valid .
we denote this state vector by symbol @xmath229 $ ] .
the state @xmath230 $ ] is an approximate eigenstate of the velocity operator ( [ eq : velocity ] ) for all values of @xmath183
. therefore @xmath231}(\theta , t)$ ] in eq . (
[ eq : omega-0-t ] ) can be understood as the non - decay law of a particle with definite speed @xmath232 .
note that the second factor in the integrand in ( [ eq : appr_dec ] ) is a slowly varying function of @xmath191 .
therefore , we can set @xmath233 in eq .
( [ eq : appr_dec ] ) and obtain    @xmath234}(\theta , t )   & \approx &    \vert   \int \limits_{m_b + m_c}^{\infty } dm     \vert \mu(m)\vert^2   e^{-\frac{it}{\hbar}\sqrt{m^2 c^4 + m_a^2 c^4\sinh^2 \theta } } \vert^2   \label{eq : omega-0-t}\end{aligned}\ ] ]    if we approximately ) . ]
identify @xmath235 with the momentum @xmath236 of the particle @xmath37 from the point of view of the moving observer @xmath10 then    @xmath234}(\theta , t )   & \approx &       e^{-\frac{i}{\hbar}\omega_{\mathbf{p } } t}|^2 \label{eq:62}\end{aligned}\ ] ]    so , in this approximation the non - decay law ( [ eq:62 ] ) in the frame of reference @xmath10 moving with the speed @xmath237 takes the same form as the non - decay law ( [ eq:35 ] ) of a particle moving with momentum @xmath238 with respect to the stationary observer @xmath9 .
in this section we are going to perform numerical calculations of the differences between the classical einstein s formula ( [ eq:1b ] )    @xmath234}^{class}(\theta , t ) = \omega_{|\mathbf{0}]}(0 , \frac{t}{\cosh \theta } )   \label{eq : time - dilation}\end{aligned}\ ] ]    and the actual non - decay law ( [ eq : omega-0-t ] ) of a moving particle having either definite momentum or sharply defined speed . in these calculations we assumed that the mass distribution @xmath171 of the unstable system ( see eqs .
( [ eq:63 ] ) - ( [ eq:64 ] ) and fig .
1 ) is centered at the value of mass @xmath239 mev/@xmath240 , the total mass of the decay products was @xmath241 mev/@xmath240 , and the width of the mass distribution was @xmath217= 20 mev/@xmath240 .
these values do not correspond to any real particle , but they are typical for strongly decaying baryons .
it is convenient to measure time in units of the classical lifetime @xmath242 . denoting @xmath243
, we find that einsteinian non - decay laws ( [ eq : time - dilation ] ) for any value of @xmath183 are given by the same universal function @xmath244 .
this function was calculated for values of @xmath245 in the interval from 0 to 6 with the step of 0.1 the calculations were performed by direct numerical integration of eq .
( [ eq:34 ] ) using the _ mathematica _ program shown below    .... gamma = 20    mass = 1000    theta = 0.0    do[print[(1/0.9375349 ) abs[nintegrate [ gamma/(2 pi ) / ( gamma^2/4 + ( x - mass)^2 ) exp [ i t sqrt [ x^2 + mass^2 ( sinh [ theta])^2 ] cosh [ theta ] / gamma ] , { x , 900 , 1010 , 1100 , 300000 } , minrecursion - > 3 , maxrecursion - > 16 , precisiongoal - > 8 , workingprecision - > 18]]^2 ] , { t , 0.0 , 6.0 , 0.1 } ] ....    as expected , function @xmath244 is very close to the exponent @xmath246 . for small values of the argument @xmath245
are too large to be shown on the scale of fig .
next we used eq .
( [ eq : omega-0-t ] ) and the above _ mathematica _ program to calculate the non - decay laws @xmath247}(\theta , \chi )   $ ] in moving reference frames for three values of the parameter @xmath183 ( = ` theta ` ) , namely 0.2 , 1.4 , and 10.0 , that correspond to velocities of 0.197c , 0.885c , and 0.999999995c , respectively .
these calculations qualitatively confirmed the time dilation formula ( [ eq : time - dilation ] ) to the accuracy of better than 0.3% .
however , they also revealed important differences @xmath248}(\theta , \chi )   - \omega ^{class}(\chi ) $ ] which are plotted as thin lines in fig .
2 .    the lifetime of the particle @xmath37 considered in our example ( @xmath249 s ) is too short to be observed experimentally .
so , calculated corrections to the einstein s time dilation law have only illustrative value . however , from these data we can estimate the magnitude of corrections for particles whose time - dependent non - decay laws can be measured in a laboratory , e.g. , for muons ( @xmath250s , @xmath251 , @xmath252 ) . taking into account that the magnitude of corrections is roughly proportional to the ratio @xmath253 @xcite , we can expect that for muons the maximum correction should be about @xmath254 which is much smaller than the accuracy of modern experiments ( about @xmath255 @xcite ) .
so , experimental observation of the predicted corrections requires significant improvements of existing experimental techniques .
note that formulas ( [ eq : omega-0-t ] ) and ( [ eq:62 ] ) for the particle seen from the moving frame are approximate in the sense that they ignore the `` decay caused by boost '' expressed by eq .
( [ eq : not - in2 ] ) . nevertheless , our major approximation ( [ eq : gg ] ) is well justified and it can not explain the discrepancy of our results ( einsteinian time dilation of the decay plus small interaction - dependent corrections ) with the conclusion of ref .
@xcite about the acceleration of the decay of moving particles
. in our view this conclusion does not refer to the experimentally measured non - decay law that is defined as the probability of finding one unstable particle . instead of eq .
( [ eq:3b ] ) , the non - decay probability was defined in @xcite by formula    @xmath256    whose physical meaning remains unclear .
in this article we assumed that the interaction responsible for the particle decay has the bakamjian - thomas instant form of dynamics .
this assumption was used to simplify calculations , but there is no good reason to believe that real interactions have this form .
then , naturally , one may ask a question : `` is there another form of dynamics in which einstein s time dilation formula ( [ eq:1b ] ) is exactly true ? '' the answer is _ no_. in any instant form of dynamics ( including non - bakamjian - thomas instant forms of dynamics ) the boost operators contain interaction terms , so the results ( [ eq : not - in ] ) - ( [ eq : not - in2 ] ) are still valid , and eq .
( [ eq:1b ] ) holds only approximately . in the point form dynamics and rotations @xmath90
are kinematical , while generators of space translations @xmath66 and time translations @xmath68 contain interaction terms . ]
the subspace @xmath52 of the unstable particle is invariant with respect to boosts , so there can be no boost - induced decays ( [ eq : not - in ] ) .
however , due to the interaction - dependence of the total momentum operator @xmath66 , one should expect decays induced by space translations    @xmath257    this prediction is not confirmed by experiments which show that the composition of an unstable particle is not affected by space translations of the observer .
times instead of experimentally observed slowdown @xcite . ]
therefore , the point form of dynamics is not acceptable for the description of decays .
similarly , never observed translation- and/or rotation - induced decays are characteristic for all non - instant forms of dynamics , e.g. , the front form .
therefore only the instant form of dynamics is appropriate for the description of particle decays .
note , however , that the decay slowdown and the absence of translation - induced decays are firmly established only for unstable systems whose lifetimes are long enough to be measured experimentally , i.e. , those governed by electromagnetic and weak interactions . for strongly interacting resonances
these properties are beyond experimental resolution , and the non - instant forms of dynamics can not be ruled out .
in this paper we analyzed the relationships between the non - decay laws in the moving reference frame and in the reference frame at rest .
we used a rigorous quantum relativistic approach that is applicable to any unstable system independent on the nature of interaction governing the decay .
a complete description of dynamics in different reference frames was obtained by using relativistic postulates i - vi and rules of quantum mechanics only .
we found that assumptions v and vi of special relativity are not needed .
moreover , their consequences ( the universal and exact time dilation of the decay of moving particles ) are in contradiction with rigorous calculations .
although the time dilation ( [ eq:1b ] ) of special relativity was qualitatively confirmed by our results , we also found small corrections to this formula that depend on the strength of interaction . in a broader sense
our results indicate that clocks viewed from the moving reference frame do not go exactly @xmath258 slower . the exact amount of time dilation depends on the physical makeup of the clock and on interactions responsible for the operation of the clock .
b. rossi and d. b. hall , _ variation of the rate of decay of mesotron with momentum _
, phys . rev .
* 59 * ( 1941 ) , 223 ; d. s. ayres , a. m. cormack , a. j. greenberg , r. w. kenney , d. o. cladwell , v. b. elings , w. p. hesse , and r. j. morrison , _ measurements of the lifetime of positive and negative pions _ ,
d. * 3 * ( 1971 ) , 1051 ; c. e. roos , j. marraffino , s. reucroft , j. waters , m. s. webster , e. g. h. williams , a. manz , r. settles , and g. wolf , @xmath259 _ lifetimes and longitudinal acceleration _ ,
nature * 286 * , ( 1980 ) , 244 .
j. bailey , k. borer , f. combley , h. drumm , f. kreinen , f. lange , e. picasso , w. von ruden , f. j. m. farley , j. h. field , w. flegel , and p. m. hattersley , _ measurements of relativistic time dilatation for positive and negative muons in a circular orbit _ ,
nature * 268 * ( 1977 ) , 301 ; f. j. m. farley , _ the cern ( g-2 ) measurements _ , z. phys .
c * 56 * ( 1992 ) , s88 .
b. bakamjian and l. h. thomas , _ relativistic particle dynamics .
ii _ , phys . rev .
* 92 * ( 1953 ) , 1300 ; f. coester and w. n. polyzou , _ relativistic quantum mechanics of particles with direct interactions _ , phys . rev .
d * 26 * ( 1982 ) , 1348 .