Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: gravitational wave ( gw ) observations are expected to provide a new view of relativistic phenomena in the universe . one of the most promising candidates for the direct detection of gws is the merger of compact binaries such as binary neutron stars ( nss ) . the second generation of ground - based gw detectors , such as advanced ligo @xcite , advanced virgo @xcite and kagra @xcite , will reach the sensitivity required to detect gws from the inspiral and coalescence of compact binary systems including binary nss within a few hundred mpc . statistical studies suggest that a few tens of merger events should be observed per year @xcite . electromagnetic counterparts of gw emitters have been recently focused on to maximize a scientific return from the expected detection of gws ( e.g. , * ? ? ? follow - up observations of these electromagnetic counterparts are important to confirm a gw detection and to investigate progenitors and environments . the electromagnetic detection also improves the localization of gw sources because the localization accuracy by photons is much better than that by the ground - based gw detectors @xmath1 deg@xmath2 ( e.g. , * ? ? ? sophisticated simulations have revealed mass ejection associated with the mergers of binary nss by several mechanisms . significant mass is dynamically ejected by gravitational torques and hydrodynamical interactions during the mergers , called dynamical ejecta ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? general relativistic simulations show that these ejecta distribute nearly isotropic compared to newtonian simulations in the cases of binary nss @xcite , while they are anisotropic for ns - black hole ( bh ) mergers @xcite . mass may be also ejected through winds driven by neutrinos @xcite , magnetic fields of and/or amplified by the merged objects @xcite , viscous heating and nuclear recombination @xcite . a traditional electromagnetic counterpart is short - hard gamma - ray bursts ( grbs ; * ? ? ? recent simulations have revealed that a hypermassive ns is formed from the merger of a ns binary ( e.g. , * ? ? ? * ) , which is believed to collapse into a bh at later time . non - collapsed matter and some ejecta falling back to the bh form a torus around the bh ( e.g. , * ? ? ? then , a relativistic jet may be launched from the bh - torus system , which is believed to be the central engine of short - hard grbs . another interesting possibility is a so - called macronova / kilonova , which is thermal emission from ejecta ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? the radiative energy of a macronova is estimated between that of a classical nova and supernova . ejecta can also produce non - thermal emission at later time similarly to supernova remnants @xcite . ejecta may accompany an advanced relativistic part , producing early emission ( @xmath3 hours ; * ? ? ? * ; * ? ? ? emission from macronovae and ns binary merger remnants is almost isotropic and hence different from that of short grbs which depends on the directions of their relativistic jets . moreover , macronovae are closer in time to mergers than emission from merger remnants and do not depend on the properties of circumburst environments . therefore , macronovae are expected to play a crucial role to localize a large sample of gw events @xcite . recently , a macronova candidate following grb 130603b was discovered @xcite . this candidate is widely interpreted as the results of the radioactive decay of @xmath0-process elements produced in the ejecta of a compact binary merger @xcite . we call this scenario _ an @xmath0-process model _ throughout this paper . the ejecta from a merger of binary nss is primarily neutron - rich . then , heavy radioactive elements ( mass number @xmath4 ) are expected to form through neutron - capture onto nuclei ( @xmath0-process nucleosynthesis ) ( e.g. , * ? ? ? although the @xmath0-process nucleosynthesis ends a few hundred millisecond after a merger , synthesized elements release energy due to nuclear fission and beta decays up to @xmath5 days ( e.g. , * ? ? ? * ) . a schematic picture for this model is shown in the left panel of figure [ figure : model ] . if this scenario is correct , the observations also give important insights into the enrichment of @xmath0-process elements in the galaxy evolution ( e.g. , * ? ? ? although the @xmath0-process model explains the observed light curve of the macronova , it is based on the limited observational data and the nuclear heating rate with large uncertainties . required mass of dynamical ejecta to explain the observations is relatively large compared with the simulation results @xcite . in addition , the occurrence of @xmath0-process nucleosynthesis needs the ejecta with low electron fraction ( @xmath6 ) . however , relatively high electron fraction ( @xmath7 ) can be also realized , which has been discussed for neutrino - driven wind ( e.g. , * ? ? ? it is worth considering other possibilities such as the scenarios of an external shock between ejecta and surrounding medium @xcite , a supramassive magnetar @xcite and dust grains @xcite .    in this study , we consider another power source of macronovae , i.e. , energy injection from the activity of the central engine , in addition to the radioactive decay of @xmath0-process elements . this is similar to the early evolution of core - collapse supernovae ( e.g. , * ? ? ? * ; * ? ? ? we call this model _ an engine model _ throughout this paper . there are several motivations to consider that the activity of the central engine contributes to the heating of ejecta . one observational motivation is the extended emission following the prompt emission of short grbs . the origin of extended emission is considered to be the activity of the central engine @xcite because the sharp drop of its light curve is difficult to be reproduced by afterglow emission @xcite . after the merger , a stable ns or a bh is formed . in the case that a bh with a torus ( or disk ) is formed , the energy injection to the ejecta is expected as a form of the jet and/or disk wind ( e.g. , * ? ? ? * ) . in the case that a ns with strong poloidal magnetic field is formed as a result of a merger , the wind of relativistic particles is ejected @xcite . then , the wind collides with the ejecta , and about half of the wind energy converts to the internal energy by the shock - heating . a schematic picture is shown in the right - hand side of figure [ figure : model ] .    the ejecta emission powered by a stable magnetar has already been discussed @xcite . they suggest that the magnetar - powered ejecta emit the brighter optical and x - ray emissions than that of the @xmath0-process model . however , they did not show that the magnetar - powered ejecta explain the detected infrared excess in grb 130603b @xcite . the engine model can provide energy enough to reproduce the detected macronova candidate , grb 130603b . we do not specify the specific heating sources . alternatively , to estimate the luminosity and temperature , we assume that the internal energy @xmath8 erg is injected to the ejecta at the time @xmath9 s after the merger . these values are consistent with typical isotropic energy @xmath10 erg and duration @xmath11 s of the extended emission @xcite . using the velocity of the ejecta @xmath12 , the temperature at @xmath13 is @xmath14^{1/4}$ ] , where @xmath15 is the radiative constant . if we only consider the adiabatic cooling for the cooling process of the ejecta , the evolution of the internal energy @xmath16 and temperature @xmath17 is scaled as @xmath18 and @xmath19 . the luminosity is described as @xmath20 . adopting the ejecta velocity @xmath21 cm s@xmath22 @xcite , the luminosity @xmath23 and the temperature @xmath17 at @xmath24s are @xmath25 and @xmath26 the observations of macronova of grb 130603b give j - band luminosity @xmath27 erg s@xmath22 and the difference between j - band and b - band @xmath28 mag which corresponds to the temperature @xmath29 k at @xmath30 days after grb 130603b in the source rest frame @xcite . therefore , in this estimate , the luminosity and temperature for the engine model is consistent with the observation of the macronova following grb130603b . we model the evolution of luminosity and temperature of a macronova . unlike the previous studies , we treat the model in an analytical manner and formulate a light curve including the early phase ( @xmath31 s ) , which is important for the search of electromagnetic counterparts of gw emitters . we consider shock - heating due to the activity of a central engine as a heating mechanism of the ejecta . for comparison , the @xmath0-process model , which has been discussed in most papers ( e.g. , * ? ? ? * ) , is also formulated . then , we compare the results of our models with observations to constrain the model parameters such as the ejected mass and the velocity of the ejecta . although our models are simplified , it is valuable to make comparison between two heating models . in section [ model ] , we introduce our model assumptions . we describe the analytical models for the evolution of luminosity and temperature in section [ evolution ] . then , we compare our results with observations in section [ discuss ] . implications for the discrimination between two models are also discussed . we summarize our results in section [ summary ] . in appendix a , we summarize the formulae for the observed temperature and bolometric luminosity . significant mass of material @xmath32 is ejected during a binary merger . we model ejecta by following the results of the general relativistic simulations of ns - ns mergers in @xcite . the simulations show that ejecta expand in a nearly homologous manner ( see also * ? ? ? the morphology of the ejecta is quasi - spherical in the case of a merger of binary nss . according to these results , we assume an isotropic and homologous expansion for the ejecta . then , the velocity of ejecta @xmath12 is @xmath33 where the radius @xmath0 originates the central engine and the time @xmath34 is measured from the time when a compact binary merges . note that in the case of a merger of ns - bh binary , the ejected mass expands with significant anisotropy @xcite . we do not consider such anisotropic ejecta in this work . @xcite found that the profile of ejecta obtained from simulations by @xcite can be well fitted by a power - law function @xmath35 . the power - law index of snapshot density @xmath36 is more or less independent on the dynamics of mergers , which is in the range of @xmath374 for @xmath38 , where @xmath39 and @xmath40 are the velocities of the outer and inner edges of the ejecta , respectively . we choose the middle of this range @xmath41 in this study . we also fix the maximum velocity @xmath42 from simulation results @xcite . the maximum velocity @xmath39 is comparable with the escape velocity of the system . the minimum velocity @xmath43 is mainly determined by complicated dynamics at the initial stage of the merger @xmath44 . for the mass density profile at the front of the ejecta , we assume the discrete boundary and the mass density @xmath45 at the region @xmath46 . although this profile may be far from the actual one , ) . such a profile is difficult to calculate precisely with current numerical calculations . ] our main aim is to compare two models for energy sources , so that our conclusions are not affected . in section [ outer ] , we discuss the dependence on the mass density profile at the outer region of the ejecta for the observed light curve . here , we only consider the evolution after the initial stage of the merger ( @xmath47 ) and treat @xmath43 as a model parameter . because of homologous expansion , the density decreases as @xmath48 . then , the density profile is described by @xmath49 where @xmath50 and @xmath51 are normalization factors . the factor @xmath52 is related to the total mass of the ejecta @xmath53 as following , @xmath54,\end{aligned}\ ] ] where we use @xmath55 from equation ( [ sec2:v ] ) . we also introduce the radius of ejecta outer edge @xmath56 and their inner edge @xmath57      the inner part of the ejecta is optically thick , and therefore the propagation of radiation in the ejecta can be regarded as a diffusion process . photons can diffusively escape from the region which satisfies that the diffusion time , @xmath58 , is smaller than the dynamical time @xmath34 , @xmath59 the medium in this region is called to be effectively thin @xcite . for convenience , we introduce a diffusion radius @xmath60 which is the radius satisfying the condition @xmath61 . furthermore , we divide ejecta into two regions called the effectively thin ( @xmath62 ) and effectively thick ( @xmath63 ) regions . near the diffusion radius , the optical depth is @xmath64 . we consider random walk for photons so that the mean number of scatterings to propagate for the distance @xmath65 is @xmath66 , where @xmath67 is the mean free path for a photon . hence , the diffusion time @xmath58 for the propagation distance @xmath65 is @xmath68 in the right hand of equation ( [ sec2 - 1:t_diff ] ) , we use @xmath69 . we calculate the diffusion radius @xmath70 from the condition @xmath71 . since the mass density profile of the ejecta is described by a decreasing power - law function ( equation [ sec2 - 1:rho ] ) , the diffusion time @xmath58 is negligible in an outer part . thus , in order to calculate the diffusion radius @xmath70 , it is a good approximation to only consider scatterings near @xmath70 ( @xmath72 ) . however , in the early phase , the distance from the outer edge of the ejecta @xmath73 to the diffusion radius @xmath70 is smaller than the diffusion radius @xmath74 . therefore , we should take the propagation distance as @xmath75 we call the first the thin - diffusion phase and the second the thick - diffusion phase throughout this paper . we schematically show these two phases in figure [ figure : diffuse ] . note that in the thin - diffusion phase , since the size of the effectively thin region is much smaller than the size of the ejecta ( @xmath76 ) , the calculation of the radiative transfer using monte carlo technique ( e.g. , * ? ? ? * ; * ? ? ? * ) requires a large number of realizations to follow the temporal evolution , which do not seem to have been considered properly so far .    to obtain the diffusion radius , we need to calculate the optical depth @xmath77 of photons which propagate a distance @xmath65 . using equations ( [ sec2 - 1:rho ] ) and ( [ sec2 - 1:deltar ] ) , the optical depth @xmath77 is described as , @xmath78^{-1 } \nonumber \\ & & \times   { \displaystyle \left[\left(\frac{r_{\rm diff}}{v_{\min}t}\right)^{1-\beta}-\left(\frac{v_{\max}}{v_{\min}}\right)^{1-\beta}\right]},\end{aligned}\ ] ] in the thin - diffusion phase , and @xmath79^{-1 } \nonumber \\ & & \times   { \displaystyle \left(\frac{r_{\rm diff}}{v_{\min}t}\right)^{1-\beta}(1 - 2^{1-\beta})},\end{aligned}\ ] ] in the thick - diffusion phase , where @xmath80 is the opacity of the ejecta . for simplicity , we use a grey approximation and a spatially uniform value of the opacity @xmath80 . from the results of @xcite ( see also * ? ? ? * ) which consider the contribution from all @xmath0-process elements to the opacity of merger ejecta , the evolution of the bolometric luminosity can be approximately described by the constant value of the opacity , @xmath81 cm@xmath2g@xmath22 . following their results , we use this value for the opacity of the ejecta . note that the exact value of the opacity of the ejecta has some uncertainties in the production efficiency of @xmath0-process elements and its spatial distribution . moreover , if the ejecta temperature is low enough for dust formation ( @xmath83 k ) , the opacity significantly increases @xcite . from these reasons , we consider the dependence on @xmath80 in section [ evolution ] . our model is based on the formulation of the light curves of supernovae ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , but there are several differences . in the case of type ii supernovae , the opacity is significantly reduced due to hydrogen recombination ( e.g. , * ? ? ? however , since the ionization potentials of the lanthanides included in the @xmath0-process elements are generally lower than that of hydrogen and the iron group , the opacity remains high at relatively low temperature @xcite . therefore , we do not consider the recombination effects for the opacity .    as far as we know , the supernova studies ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) have not taken into account the thin - diffusion phase , which is necessary for treating the thickness of the diffusion length appropriately and estimating the physical quantities by the values at the outer edge of the ejecta in the analytical formulae . this phase may be also important for the case of supernovae . some supernova studies consider the planar phase @xcite in which the evolution of the ejecta is approximately planar as long as its radius do not double . in the case of the ns - ns merger , since the initial length scale of the merger system is small @xmath84 cm and the velocity of the merger ejecta is subrelativistic , the planar phase is irrelevant for the observations . one of the two heating mechanisms we consider is nuclear heating by @xmath0-process elements . since the beta decay products of @xmath0-process elements produced in ns binary mergers naturally heat ejecta , this mechanism is considered to power the emission of a macronova ( e.g. , * ? ? ? the nuclear heating rate is calculated in several works @xcite . the derived heating rates per unit mass @xmath85 are described by the following formula @xmath86 in this study , we use @xmath87 and @xmath88 obtained by @xcite . the value of @xmath89 has been obtained by simulations under some simplified assumptions with only limited parameter regions . thus , we should note that the value of @xmath90 has uncertainties . the injected internal energy by the nuclear decay is @xmath91 in the region @xmath63 . on the other hand , the injected energy in this region is decreased by adiabatic cooling . the time - evolution of internal energy due to the adiabatic cooling is proportional to @xmath92 . comparing the two temporal evolution , the index of the adiabatic cooling is smaller than that of the increase of internal energy due to the nuclear decay for @xmath93 . since we use @xmath87 , we neglect the injected internal energy in the region @xmath63 . unlike the @xmath0-process model , energy injection occurs only within the time @xmath13 in the engine model . we only consider adiabatic cooling as a cooling process of ejecta after @xmath13 , and therefore , the temperature distribution at time @xmath34 is , @xmath94 where the index @xmath95 is a parameter for a snapshot distribution and @xmath96 is a normalization factor described later . the time dependence of @xmath92 is the effect of adiabatic expansion . the normalized value @xmath96 is determined by using the relation of total injected internal energy @xmath97 as @xmath98,\end{aligned}\ ] ] where we use @xmath99 . for the temperature index @xmath100 , the innermost region of ejecta has dominant internal energy . as will be shown in section [ evolution ] , since the luminosity and temperature always depend on the product of @xmath97 and @xmath13 , we treat @xmath101 as a parameter . thus , the engine model has two parameters , @xmath95 and @xmath101 instead of @xmath90 and @xmath102 in the @xmath0-process model . energy injection is not always a single event and the shock does not always get through the whole ejecta . it is considered that the activity of the central engine accompanies violent time variability . in this case , multiple shocks propagate into the ejecta . some of the shock may not catch up with the outer edge of the ejecta . current general relativistic simulations can not calculate the evolution of ejecta for such a long time after merger ( @xmath103 s ) , so that the index @xmath95 of temperature distribution is highly uncertain . therefore , we treat the temperature index @xmath95 as a parameter .    unlike the case of core - collapse supernova @xcite , it is difficult to determine the temperature distribution of heated ejecta by the activity of a central engine . in the case that the activity of a central engine injects the energy into the ejecta , the radiation - dominated shock ( where the internal energy behind the shock is dominated by radiation ) is formed in the ejecta . the ejecta are heated during the propagation of the shock . this situation is similar to the initial phase of core - collapse supernovae ( e.g. , * ? ? ? * ; * ? ? ? in the cases of core - collapse supernovae , the kinetic energy of ejecta before the shock heating is much smaller than the injected internal energy . in such ejecta , the relation between velocity and mass density was obtained by @xcite ( in the non - relativistic case for the velocity of the ejecta ) . using sakurai s ( 1960 ) solution and the equipartition between the kinetic energy after the shock heating and the internal energy @xcite , the distribution of the temperature distribution is derived . however , in the case of compact binary mergers , the merger ejecta have a large velocity ( @xmath104 ) before the shock heating @xcite . then , injected internal energy is not always larger than the kinetic energy of ejecta so that it is not clear whether we can use the equipartition to estimate the distribution of internal energy or not . the kinetic energy of the ejecta @xmath105 is described as @xmath106}{(5-\beta)\left[1-\left(\frac{v_{\max}}{v_{\min}}\right)^{3-\beta}\right]}}.\end{aligned}\ ] ] note that if the injected internal energy @xmath97 is larger than the kinetic energy of the ejecta , it is expected that some of the internal energy converts to the kinetic energy of the ejecta . as a result , the internal energy and the kinetic energy are equal as in the case of core - collapse supernovae . then , the mass density distribution and the maximum velocity of the ejecta derived from simulations may be changed because the injection time may be long @xmath107s compared to that calculated by simulations @xmath108s @xcite . for simplicity , we only consider the case @xmath109 . clccc + symbol & & fiducial model & minimum mass model & hot interior model + @xmath53 & ejecta mass & @xmath110 & @xmath111 & @xmath112 + @xmath43 & minimum velocity & @xmath113 & @xmath114 & @xmath115 + @xmath39 & maximum velocity & @xmath116 & @xmath116 & @xmath116 + @xmath36 & index of the density profile & 3.5 & 3.5 & 3.5 + @xmath80 & opacity & 10 cm@xmath2 g@xmath22 & 30 cm@xmath2 g@xmath22 & 10 cm@xmath2 g@xmath22 + @xmath90 & nuclear heating rate at 1 day & @xmath117 erg s@xmath22 g@xmath22 & @xmath118 & @xmath117 erg s@xmath22 g@xmath22 + @xmath102 & index of nuclear heating rate & 1.3 & @xmath118 & 1.3 + @xmath97 & internal energy at @xmath13 & @xmath119 erg & @xmath120 erg & @xmath121 erg + @xmath13 & injection time & @xmath122 s & @xmath122 s & @xmath122 s + @xmath95 & index of the temperature profile & 1.6 & 1.1 & 2.7 +    [ tab : parameter ]    in this section , we present the evolution of the observed temperature and luminosity of a macronova using our model introduced in the previous section . in sections [ thin ]  [ transparent ] , we focus on the parameter dependence of the evolution using some approximations . in section [ fiducial ] , we calculate the temperature and luminosity using the fiducial model with parameters summarized in the first column of table 1 .    to calculate the luminosity and temperature , we assume that the emission is well described by the blackbody radiation ( e.g. , * ? ? ? * ) . for simplicity , we assume that the observed temperature equals to the temperature at the diffusion radius @xmath70 . we also assume that the temperature is not so different from the diffusion radius @xmath70 to @xmath123 so that in the thick - diffusion phase ( @xmath124 ) , we only consider the emission from @xmath70 to @xmath123 to calculate the observed luminosity for both the @xmath0-process and engine models . in some studies ( e.g. , * ? ? ? * ) , the observed temperature is approximated by the temperature at the radius of the photosphere @xmath125 where the optical depth is unity . since the velocity of the ejecta is near the light speed , the optical depths at the diffusion radius @xmath70 and its twice @xmath123 are @xmath126 . therefore , our assumed temperature approximately equals to the temperature at the photosphere .    in section [ diffusion ] , we introduced two phases , the thin- and thick - diffusion phases ( figure [ figure : diffuse ] ) , depending on the size of the region where photons make the diffusion in the ejecta @xmath65 . we also introduce another phase @xmath127 , the transparent phase , in which photons can diffuse out from the entire of the ejecta . thus , we divide the evolution into these three phases for the values of the diffusion radius @xmath70 as described below .      the size of the effectively thin region gets larger with time . at the early phase of a macronova , the diffusion radius @xmath70 , which is the inner radius of the effectively thin region , is near the outer edge of the ejecta @xmath73 . in this early phase , we take the propagation distance @xmath65 of a photon as @xmath128 . since we assume that the density is a homologous function of the velocity @xmath35 , the density can be approximated as @xmath129 in the region @xmath130 . using the escaping condition for the diffusing photons @xmath131 , equation ( [ sec2 - 1:t_diff ] ) and approximation on the optical depth @xmath132 , the propagation distance @xmath65 can be estimated as @xmath133 in the discussion of parameter dependence ( sections [ thin ]  [ transparent ] ) , we only consider the dominant term . for example , we neglect the second term in the right - hand side of equation ( [ sec2 - 1:m_ej ] ) to derive the parameter equation ( [ sec2 - 2 - 1:deltar ] ) because the index of the mass density is @xmath134 in our model . in section [ fiducial ] , we include the subdominant terms to calculate the light curves numerically . first we consider the @xmath0-process model . the evolution of temperature @xmath135 is obtained by the internal energy density @xmath136 at the radius @xmath137 . using equations ( [ sec2 - 1:rho ] ) and ( [ sec2 - 1:m_ej ] ) , the parameter dependence of the density is @xmath138 . the observed temperature is @xmath139 for @xmath87 , the observed temperature evolves as @xmath140 . this is because in the thin - diffusion phase the ejecta is effectively a single expanding shell with @xmath141 and the injected energy @xmath142 is almost constant so that the observed temperature approximately follow adiabatic cooling @xmath19 . note that in this phase the observed temperature does not depend on the opacity . the bolometric luminosity @xmath143 for the radioactivity is described as the product of the mass within the thickness @xmath65 in equation ( [ sec2 - 2 - 1:deltar ] ) and the nuclear heating rate @xmath89 in equation ( [ sec2 - 1:dotepsilon ] ) so that @xmath144 for @xmath87 , the evolution of the bolometric luminosity is @xmath145 . next we consider the engine model . we should take into account the freedom of the temperature index @xmath95 in the temperature distribution ( equation [ sec2 - 1:t ] ) . since we only consider a dominant term in the right - hand side of equation ( [ sec2 - 1:e_int ] ) ( the first term for @xmath146 or the second term for @xmath147 ) in this subsection , the parameter dependence of the temperature @xmath96 is described as @xmath148 substituting @xmath149 into equation ( [ sec2 - 1:t ] ) , the observed temperature is described as @xmath150 since the observed temperature @xmath135 approximately equals to the temperature at the outer edge of the ejecta , @xmath151 , the evolution of the observed temperature and the luminosity are also determined by the adiabatic cooling . the bolometric luminosity in the effectively thin region is equal to the total radiation created by thermal emission in this region @xcite . using equation ( [ sec2 - 2 - 1:t_obs , sh ] ) , the bolometric luminosity is described as @xmath152 the time evolution of bolometric luminosity for the engine model does not depend on the temperature index @xmath95 in the thin - diffusion phase . comparing the engine model with the @xmath0-process model in the thin - diffusion phase , the bolometric luminosity and the observed temperature decrease faster in the engine model than those in the @xmath0-process model . these time - dependence do not depend on the indices of the density and temperature . note that the light curve may depend on the detailed profile of the front of the ejecta in this thin - diffusion phase . the profile of the ejecta front is difficult to calculate by the numerical simulation due to its low density , and hence has large uncertain @xcite . we discuss its dependence in section [ discuss ] . we consider diffusion to evaluate the diffusion radius in the thick - diffusion phase . we take the propagation distance @xmath72 after the time when the difference between the radius of the outer edge of the ejecta @xmath73 and the diffusion radius @xmath70 is larger than the diffusion radius , @xmath153 , since the optical depth of the outer part is negligible for the density profile in equation ( [ sec2 - 1:rho ] ) . in this thick - diffusion phase , mass density significantly deviates from @xmath154 . substituting equations ( [ sec2 - 1:t_diff ] ) and ( [ sec2 - 1:tau2 ] ) into @xmath61 , the diffusion radius @xmath70 is calculated as @xmath155^{\frac{1}{\beta-2 } } \nonumber \\ & \propto&\kappa^{\frac{1}{\beta-2}}m_{\rm ej}^{\frac{1}{\beta-2}}v_{\min}^{\frac{\beta-3}{\beta-2}}t^{\frac{\beta-4}{\beta-2}},\end{aligned}\ ] ] where we use the relations @xmath72 and @xmath156 . the latter is obtained from the assumption of the homologous expansion . regarding the optical depth @xmath77 , the second term is neglected in the right - hand side of equation ( [ sec2 - 1:tau ] ) to focus only on the dominant term to study parameter dependence in sections [ thin ]  [ transparent ] . for @xmath41 , the diffusion radius decreases with time ( @xmath157 ) . then , emission from the region with relatively high mass density can be observed progressively in this phase ( @xmath158 ) . we introduce the transition time @xmath159 between thin- and thick - diffusion phases , which satisfies the relation @xmath160 . substituting equation ( [ sec2 - 2 - 2:r_diff ] ) and @xmath161 into the relation @xmath160 , we can obtain the transition time @xmath159 as @xmath162 where @xmath163 cm@xmath2 g@xmath22 , @xmath164 , @xmath165 and @xmath166 . as seen above , the transition time @xmath159 is typically several days . this timescale is expected to allow the follow - up observations @xcite . thus , we should consider both phases to predict something useful for follow - up observations . if we fix @xmath80 and @xmath39 and use @xmath41 , equation ( [ sec2 - 2 - 2:t_times ] ) gives @xmath167 . if we increase the total mass of the ejecta @xmath53 and the velocity at the inner edge of the ejecta @xmath43 , the mass density of the ejecta @xmath168 and hence the optical depth are increased . as a result , the transition time @xmath159 becomes large . first , we consider the @xmath0-process model . here , we introduce the velocity @xmath169 based on the homologous relation . using the velocity @xmath170 and equation ( [ sec2 - 2 - 2:r_diff ] ) for the mass density ( equation [ sec2 - 1:rho ] ) , the evolution of temperature is @xmath171 for @xmath87 and @xmath41 , the evolution is described by @xmath172 so that the observed temperature increases with time . the bolometric luminosity is described as the product of the mass between @xmath70 and @xmath123 and the nuclear heating rate @xmath89 . using equations ( [ sec2 - 1:dotepsilon ] ) and ( [ sec2 - 2 - 2:r_diff ] ) , we obtain the bolometric luminosity as , @xmath173 for @xmath87 and @xmath41 , the evolution of the bolometric luminosity is @xmath174 . next , we consider the engine model . using equations ( [ sec2 - 1:t ] ) and ( [ sec2 - 2 - 2:r_diff ] ) , the evolution of the observed temperature is described as @xmath175 for @xmath41 , the value @xmath176 is the boundary whether the observed temperature increases with time ( @xmath100 ) or not ( @xmath146 ) . the evolution of the luminosity equals to the total radiation created by thermal emission in the sphere with radius @xmath70 . using the relation @xmath177 and equations ( [ sec2 - 2 - 2:r_diff ] ) and ( [ sec2 - 2 - 2:t_obs , sh ] ) , we obtain @xmath178 if we take @xmath41 and @xmath179 , the evolution of the bolometric luminosity is @xmath180 . this is almost the same dependence as in the @xmath0-process model . note that even if the inner part of the ejecta has the larger internal energy ( @xmath100 ) , the bolometric luminosity does not always increase with time . using the relation @xmath181 , from the adiabatic cooling , the evolution of bolometric luminosity for a given mass shell with @xmath12 is @xmath182 , where @xmath183 is the total internal energy for the mass shell with a given expanding velocity @xmath12 . since @xmath184 and @xmath185 , the bolometric luminosity increases with time for the value of the temperature index @xmath186 . once the diffusion radius reaches the inner edge of the ejecta @xmath187 , all photons emitted from the ejecta can diffuse out within dynamical timescale . if energy is not injected into the ejecta in this transparent phase , the internal energy in the ejecta runs out immediately . the transition time from the thick - diffusion phase to the transparent phase @xmath188 is described as @xmath189 where we use the diffusion radius @xmath190 . first we consider the @xmath0-process model . the observed temperature equals to the temperature at the inner edge of the ejecta , @xmath191^{1/4}$ ] . using equation ( [ sec2 - 1:rho ] ) , we obtain @xmath192 since the energy is continuously injected due to the nuclear heating in the @xmath0-process model , the bolometric luminosity from the entire ejecta is described as @xmath193 . however , the outer part of the ejecta emits photons with lower temperature and/or x - rays and @xmath194-rays produced directly in radioactive decays . although such emission contributes to the bolometric luminosity , we here focus only on the optical and infrared emissions . in the thick - diffusion phase , the observed emission comes from the region between @xmath195 and @xmath196 . in the transparent phase , we assume that the time evolution of the diffusion radius @xmath70 is the same as the thick - diffusion phase until @xmath197 and the observed luminosity comes from the region from @xmath198 to @xmath123 for simplicity . then , the bolometric luminosity is described as @xmath199 although it appears that this time evolution directly reflects the nuclear decay rate , when we calculate the mass between @xmath198 and @xmath123 the evolution of the upper limit of the integration @xmath123 makes the decrease of the luminosity faster than @xmath200 ( see a dashed line in the middle panel of figure [ figure : evolution ] ) . in addition , the evolution of @xmath70 depends on the index @xmath36 ( see equation [ sec2 - 2 - 2:r_diff ] ) , so that the mass between @xmath198 and @xmath123 also depends on the index @xmath36 . next we consider the engine model . we assume that the internal energy is exhausted when the diffusion radius reaches @xmath197 . for the observed temperature @xmath135 , we assume the relation @xmath201 and use equation ( [ sec2 - 1:t ] ) , @xmath202 the bolometric luminosity is described as @xmath203 since the internal energy at the innermost region almost equals to the total internal energy @xmath204 and determines the bolometric luminosity @xmath205 for the temperature index @xmath100 , the bolometric luminosity does not depend on the mass @xmath53 and velocities @xmath39 and @xmath43 . this luminosity always corresponds to the maximum luminosity for @xmath147 , so that we can impose the lower limit on the parameter @xmath101 . -process model and the engine model , respectively . for comparison , we also plot the bolometric luminosity from the whole ejecta for the @xmath0-process model after the transparent phase ( @xmath206 in equation [ sec2 - 2 - 3:t_tr ] ) as a blue long - dashed line in the middle panel.,width=264 ]        we show the temporal evolution of the diffusion radius @xmath70 , the bolometric luminosity @xmath143 and the observed temperature @xmath135 in figure [ figure : evolution ] under the fiducial parameter set . the parameters are summarized in the first column of table 1 . here , we do not use approximations @xmath207 and @xmath208 at the thin - diffusion phase as in section [ thin ] . instead , the diffusion radius @xmath70 is calculated from equations ( [ sec2 - 1:t_diff ] )  ( [ sec2 - 1:tau ] ) without approximations . using the obtained diffusion radius @xmath70 and the relation @xmath169 , we calculate the observed temperatures in the thin- and thick - diffusion phases , @xmath209^{1/4}$ ] ( equation [ sec2 - 2 - 2:t_obs , nuc ] ) , and @xmath210 ( equation [ sec2 - 2 - 2:t_obs , sh ] ) for the @xmath0-process and the engine models , respectively . in the transparent phase , the temperature in equations ( [ sec2 - 2 - 3:t_obs , nuc ] ) and ( [ sec2 - 2 - 3:t_obs , sh ] ) are evaluated with @xmath211 . equations on observed temperature and bolometric luminosity for both models are summarized in appendix . the set of parameters we choose here explains the observed optical and infrared light curves of grb 130603b ( see next section ) . the vertical dash - dotted lines in figure [ figure : evolution ] show the time @xmath212 ( equation [ sec2 - 2 - 2:t_times ] ) . the diffusion radius is plotted only up to the transition time @xmath213 ( equation [ sec2 - 2 - 3:t_tr ] ) .    in the thick - diffusion phase , the diffusion radius ( @xmath214 for @xmath41 ) moves inward in the ejecta ( @xmath215 ) . since the observed luminosity and temperature are determined at the diffusion radius @xmath70 , the time evolution of luminosity and temperature strongly depends on the indices of the profile , @xmath36 and @xmath95 . for the @xmath0-process model , the bolometric luminosity decreases with time ( @xmath216 ) in the thick - diffusion phase ( see equation [ sec2 - 2 - 2:l_bol , nuc ] ) , which is more rapid than that in the thin - diffusion phase ( @xmath217 , see equation [ sec2 - 2 - 1:l_bol , nuc ] ) . since the index of the mass density @xmath41 is close to 3 , in which the mass of each shell with a certain size @xmath218 is the same value in logarithmic scale , the mass between the diffusion radius @xmath70 and its doubled value @xmath123 does not significantly change with time . the luminosity is mainly determined by that mass , so that the evolution of the luminosity is slow compared with the evolution of nuclear heating rate ( @xmath200 ) in the thick - diffusion phase . on the other hand , bolometric luminosity and observed temperature increase with time in the engine model with the parameter set of the fiducial model . these mainly reflect the profile of the temperature distribution ( @xmath219 ) . in fact , using equation ( [ sec2 - 2 - 2:l_bol , sh ] ) , the index of the time @xmath34 for the bolometric luminosity is @xmath220 for the engine model . after the transition time @xmath221 , the luminosity and temperature are almost determined by the quantities at the inner edge of the ejecta . then , the evolution of the luminosity and temperature does not significantly depend on the indices of profile @xmath36 and @xmath95 as in the case of the thin - diffusion phase ( except for the case @xmath146 of the engine model , equation [ sec2 - 2 - 3:l_bol , sh ] ) . since our used profile of mass density has an artificially steep cut - off at the inner edge of the ejecta ( figure [ figure : diffuse ] ) , bolometric luminosity in both models rapidly declines after the time @xmath221 . in the bottom panel of figure [ figure : evolution ] , the observed temperature in both models has a steep cutoff at @xmath197 . for comparison , we also consider the time evolution of bolometric luminosity from the whole ejecta @xmath222 in the @xmath0-process model . time evolution is shown in the middle panel of figure [ figure : evolution ] as a blue long - dashed line . this luminosity evolution ( @xmath223 ) is significantly slower than that of the engine model in the transparent phase . in section [ implication ] , we discuss the implication for discriminating the @xmath0-process model and the engine model using these temporal behaviors . -process model ) . since two regions are overlapped , the color looks like purple for the @xmath0-process model . these regions are only a schematic view . we fix @xmath42 , @xmath41 , @xmath224 erg s@xmath22 g@xmath22 and @xmath87 . for the opacity @xmath80 , we use the range @xmath225 cm@xmath2 g@xmath22 . for the engine model , we treat @xmath101 and @xmath95 as free parameters to fit the light curve . the circle , square and triangle denote the case of the fiducial model , minimum mass model and hot interior model , respectively ( see table 1).,width=377 ]          we compare the results with the optical and infrared observations of short grb 130603b in figure [ figure : observation ] . the fiducial parameter set in table 1 is adopted . the @xmath0-process model and the engine model result in similar light curves at the optical and infrared bands . both of them satisfy the observational data of grb 130603b . note that the detection point at f606w band at @xmath226 s is consistent with the afterglow of grb 130603b modeled as a smoothly broken power law ( blue dashed line , * ? ? ? we regard this detected value as an upper limit for the luminosity of emission from the ejecta . the detection point at f160w band at @xmath84 s exceeds the extrapolation of the afterglow emission ( red dashed line , * ? ? ? * ) , so that we regard this detected emission as a thermal radiation from the ejecta . the range of the model parameters @xmath43 and @xmath53 to satisfy the constraints obtained from the observation of grb 130603b is shown in figure [ figure : parameter ] as colored areas ( red area for the @xmath0-process model and blue area for the engine model ) . note that the red area has a completely overlap with the blue area . we fix the other model parameters @xmath42 , @xmath41 , @xmath224 erg s@xmath22 g@xmath22 and @xmath87 as in the fiducial model . we take into account the uncertain range of the opacity , @xmath225 cm@xmath2 g@xmath22 to constrain the parameters , @xmath43 and @xmath53 . in the engine model , @xmath95 and @xmath101 are additionally treated as free parameters to derive the allowed area in figure [ figure : parameter ] .      in the @xmath0-process model , the luminosity becomes smaller for smaller ejecta mass @xmath53 . the small ejected mass @xmath227 can not reproduce the infrared excess of grb 130603b ( figure [ figure : observation ] ) . the required ejecta mass is relatively large compared to the mass indicated by recent numerical simulations for a merger of binary nss ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? note that in @xcite , @xmath228 is required to explain the observed infrared excess , which is a factor @xmath32 smaller than our results . their theoretical light curves are based on the study of @xcite . in @xcite , a broken power - law mass density profile with the index @xmath229 for the inner layer and @xmath230 for the outer layer of ejecta is adopted , in which the mass of the ejecta is efficiently concentrated at the transition point of the density index . therefore , the luminosity of ejecta is evaluated as the heating rate multiplied by the total ejecta mass at the moment when the diffusion radius reaches the transition point . on the other hand , the index of our mass density profile of the ejecta is @xmath231 , which is indicated by general relativistic simulations by @xcite . this profile is quite different from the profile adopted in @xcite ; the index is close to 3 , in which the mass in each logarithmic radius is constant . then , at the time @xmath232 ( equation [ sec2 - 2 - 3:t_tr ] ) , the mass contributing to the luminosity is about @xmath360% of the total ejecta mass in the case @xmath233 . this profile predicts luminosity dimmer than that other studies . in fact , @xcite tried to explain the observed infrared excess using the mass profile which is almost the same with ours . in the case of a binary ns merger with ejecta mass @xmath30.02 @xmath234 , even if they use a larger nuclear heating rate ( larger by a factor of 2 ) , their predicted luminosity is slightly smaller than the observed infrared excess ( in the left panel of their figure 3 ) . this result is consistent with our model , i.e. , our model requires larger mass than most of previous studies .    in the engine model , the injected internal energy which determines the luminosity does not depend on the ejecta mass ( except for the limit in equation [ sec2 - 1:e_int ] ) . however , the luminosity declines rapidly after the transition time @xmath235 which depends on the ejecta mass as in equation ( [ sec2 - 2 - 3:t_tr ] ) . the condition @xmath236 s in the observer frame is required to reproduce the excess observed from grb 130603b in the near - infrared band . this condition gives the lower limit for the ejecta mass in the engine model , @xmath237 with the opacity @xmath238 cm@xmath2 g@xmath22 . note that the observed upper limit on the infrared luminosity at @xmath239 s in the observer frame ( figure [ figure : observation ] ) , which corresponds to @xmath240 s , gives the upper limit on the ejecta mass for both models . however , this limit is not important for the range @xmath241 in the range of the opacity @xmath225 cm@xmath2 g@xmath22 . the smaller minimum velocity @xmath43 gives the smaller bolometric luminosity at certain time in the @xmath0-process model ( see equations [ sec2 - 2 - 1:l_bol , nuc ] and [ sec2 - 2 - 2:l_bol , nuc ] ) . the small minimum velocity enlarges the size of ejecta ( when we fix the maximum velocity @xmath39 ) . then , the diffusion time @xmath58 of photons emitted from the inner region of the ejecta becomes long for the small velocity @xmath43 ( equation [ sec2 - 2 - 3:t_tr ] ) . the mass between @xmath70 and @xmath123 ( or @xmath73 ) increases toward inner region of the ejecta ( as long as @xmath134 ) so that the mass is reduced for the small minimum velocity @xmath40 at certain time . in fact , the dependence of the mass on the minimum velocity is @xmath242 . as a result , smaller minimum velocity gives smaller luminosity to reproduce the observed infrared excess of grb 130603b . moreover , smaller minimum velocity gives larger temperature @xmath135 at certain time ( equations [ sec2 - 2 - 2:t_obs , nuc ] and [ sec2 - 2 - 3:t_obs , nuc ] ) because mass density at a shell with small velocity is large . the difference between the detected luminosity at f160w band and the upper limit on the luminosity at f606w band at @xmath84 s in the observer frame gives the upper limit on the observed temperature ( @xmath243 k ) . to satisfy the observed upper limit on the temperature from grb 130603b , a lower limit of @xmath244 is obtained for @xmath245 . the smaller minimum velocity @xmath43 gives higher temperature @xmath135 in the engine model ( equations [ sec2 - 2 - 2:t_obs , sh ] and [ sec2 - 2 - 3:t_obs , sh ] ) . the observational limit for the temperature at @xmath84 s in the observer frame indicates that the range of the minimum velocity @xmath40 is limited in the engine model ( @xmath246 for @xmath245 ) . we discuss the dependence on the value of @xmath80 . as mentioned in section 2.2 , we use the temperature - independent opacity @xmath80 with the grey approximation . in general , the @xmath0-process line opacity depends on frequency and changes with temperature and ionization state of the ejecta @xcite . the indicated grey opacity is @xmath225 cm@xmath2 g@xmath22 .    in the case of the @xmath0-process model , the luminosity significantly depends on opacity @xmath80 . the larger opacity causes larger diffusion time @xmath58 , so that larger time is required to observe the inner region of the ejecta for given ejecta mass @xmath53 and minimum velocity @xmath43 . in fact , two transition times @xmath159 and @xmath188 are proportional to @xmath247 ( equations [ sec2 - 2 - 2:t_times ] and [ sec2 - 2 - 3:t_tr ] ) . then , the mass around the diffusion radius @xmath70 is small at certain time , so that the luminosity is reduced . as a result , in order to explain the infrared excess observed in grb 130603b , larger mass @xmath53 is required for the larger value of opacity @xmath80 . for the opacity @xmath248 cm@xmath2 g@xmath22 , total ejecta mass @xmath249 is required to reproduce the observed excess , which is much larger than the simulation results of mergers of binary nss ( e.g. , * ? ? ? * ) . on the other hand , the transition time @xmath188 is smaller for the smaller value of the opacity . then , the luminosity significantly increases at @xmath226 s. for the opacity @xmath250 cm@xmath2 g@xmath22 , there is no parameter set which gives smaller luminosity than the detection at f606w band ( @xmath226 s in the observer frame ) and the luminosity comparable to the observed excess at f160w band simultaneously in the @xmath0-process model .    in the case of the engine model , a larger value of opacity @xmath80 reduces the lower limit for the mass @xmath53 to explain the observed excess . for certain temperature and luminosity , the opacity @xmath80 and the ejecta mass @xmath53 always degenerate in the form @xmath251 ( see equations [ sec2 - 2 - 1:t_obs , sh ] , [ sec2 - 2 - 1:l_bol , sh ] , [ sec2 - 2 - 2:t_obs , sh ] , [ sec2 - 2 - 2:l_bol , sh ] , [ sec2 - 2 - 3:t_obs , sh ] and [ sec2 - 2 - 3:l_bol , sh ] ) . this dependence comes from the optical depth ( equation [ sec2 - 1:tau ] ) because the internal energy in the ejecta does not depend on the opacity and the ejecta mass , contrary to the @xmath0-process model . we present a parameter set to give the minimum ejecta mass @xmath53 in table 1 as the minimum mass parameter set . we also plot the value of @xmath53 and @xmath43 of this model in figure [ figure : parameter ] as a square . this ejecta mass is naturally realized in general relativistic simulations ( e.g. , * ? ? ? although the larger value of the opacity @xmath80 reduces the lower limit for the ejecta mass @xmath53 , the kinetic energy is @xmath252 erg ( equation [ sec2 - 1:e_kin ] ) which is close to the initial injected energy @xmath253 erg for the minimum mass parameter set . the lower ejecta mass @xmath53 reduces the kinetic energy of the ejecta , @xmath254 in equation ( [ sec2 - 1:e_kin ] ) , so that the required energy @xmath97 may exceed the kinetic energy of the ejecta for the larger opacity . for the small value of the opacity , larger mass and smaller minimum velocity is required to satisfy the condition @xmath236 s ( @xmath255 in equation [ sec2 - 2 - 3:t_tr ] ) . for the opacity @xmath256 cm@xmath2 g@xmath22 , the condition corresponds to @xmath257 . the observational constraint for the temperature also requires a large value of the minimum velocity @xmath43 . then , there is no solution to explain the observed excess within the parameter range shown in figure [ figure : parameter ] for the opacity @xmath258 cm@xmath2 g@xmath22 . therefore , for the small opacity @xmath259 cm@xmath2 g@xmath22 the engine model can not explain the observed excess .      since the engine model has additional free parameters , @xmath95 and @xmath101 , the allowed region of the parameters is larger than that of the @xmath0-process model . we can impose the lower limit on the parameter @xmath101 by regarding the infrared luminosity @xmath27 erg s@xmath22 at @xmath260 day as bolometric luminosity in the source rest frame with equation ( [ sec2 - 2 - 3:l_bol , sh ] ) . the derived limit is @xmath261 . to satisfy the optical upper limit at @xmath226 s and the detected luminosity at @xmath84 s in the observer frame ( figure [ figure : observation ] ) , we find the lower limit on the index of the temperature profile @xmath262 . for a smaller value of the index @xmath95 , emission from the ejecta with relatively high temperature can be observed at time @xmath226 s , so that luminosity at f606w band is larger than the observed upper limit of grb 130603b . in addition , the smaller value of @xmath95 decreases relative internal energy in the inner edge of the ejecta . to reproduce the luminosity at time @xmath263 s in the observer frame when observed emission comes from the inner ejecta , the smaller @xmath95 requires the larger initial internal energy @xmath97 which exceeds the kinetic energy of the ejecta @xmath105 in some cases . several deep optical observations of short grbs give stringent upper limits on the luminosity of macronovae @xcite . we compare the results with two deep optical observations of short grbs , grb 050509b and grb 080905a . for the fiducial parameter set , the luminosity exceeds the observational upper limits on these two observations . in the engine model , we can reduce the luminosity in the early phase @xmath264 s without reducing the luminosity in the late phase @xmath84 s by utilizing the steep temperature profile ( large @xmath95 ) .    here , we introduce the hot interior parameter set with larger value of index @xmath95 than that of the fiducial parameter set . since emission from the inner part of the ejecta is observed at the later time , the luminosity at the early phase decreases and avoids the observational limits if most of the internal energy is injected to the inner part of the ejecta . we show the light curve of the hot interior parameter set in figure [ figure : observation2 ] . we choose the parameters as @xmath265 and @xmath266 ( the right column of table 1 ) . from figure [ figure : observation2 ] , the light curves are consistent with all three observations using the same model parameters . a possible scenario for the hot interior parameter set is that the shock produced by the activity of the central engine may not be able to catch up with the outer part of the ejecta because the velocity of the ejecta is close to the light speed ( @xmath42 ) . then , only the inner part of the ejecta will be heated . for comparison , we also show the light curves in the @xmath0-process model with the parameter set of the hot interior in figure [ figure : observation2 ] as dashed lines . the luminosity of the @xmath0-process model exceeds the observed upper limits in two observations , grb 050509b and grb 080905a ( left and middle panels of figure [ figure : observation2 ] ) if we choose the parameter set of the hot interior model . we are not able to find any parameter set in the @xmath0-process model , which simultaneously satisfies the observed limits of the three observations . note that we do not argue that the @xmath0-process model is excluded from these results because we need to take into account the variations of the model parameters for each event . note that the extended emission was not detected in three short grbs . however , it is not unreasonable to miss the extended emission of these bursts . one possibility is a selection effect . observationally , the fraction of short grbs with extended emission is significantly larger at softer energy bands : @xmath325% in the swift bat samples ( @xmath26715 kev ; * ? ? ? * ) and @xmath37% in the batse samples ( @xmath26720 kev ; * ? ? ? this suggests that observations with a low energy threshold may dramatically increase short grbs with extended emission @xcite . this will be further tested by future soft x - ray survey facilities such as wide - field maxi ( 0.7 - 10 kev ) @xcite . the three referred short grbs were detected by swift bat and therefore swift bat could not detect extended emission by chance . alternatively , the outflow following the main short grb jet could not breakout the ejecta . @xcite and @xcite investigated the propagation of jets in merger ejecta . they found the cases that relativistic jets can penetrate merger ejecta and produce the prompt emission of short grbs , but in the late energy injection cases , outflow fails to breakout the ejecta . therefore , some extended emission may not be observed , although the central engine works actively . -process model on the shape of the front of the ejecta . we plot the case of grb 130603b . thick - dashed lines denote the fiducial model . thin - solid , thin - dashed , and thin - dot - dashed lines correspond to the exponential profiles of mass density ( equation [ sec3:rho ] ) with @xmath268 and @xmath269 , respectively . we also plot the observational data from @xcite.,width=302 ]    in the thin - diffusion phase , the light curve strongly depends on the density profile of the ejecta surface . the density profile is determined by the complex merger dynamics @xcite , so that the density profile of the ejecta can not be analytically derived as mentioned in section 2.3 . since the outer part of the density profile is difficult to calculate precisely , little attention has paid on the mass profile at the outer region in current numerical simulations . in order to investigate the dependence of the light curve on the mass profiles in the thin - diffusion phase , we consider other forms of the mass profile and compare the light curve with that of equation ( [ sec2 - 1:rho ] ) . we adopt an exponential profile @xmath270 we introduce an additional free parameter @xmath271 and the ejecta expand @xmath272 . the calculations are the same except for using @xmath273 . we fix the mass with velocity larger than @xmath274 and calculate three models for @xmath275 . for the other parameters , we adopt from the fiducial parameter set in table 1 . we show the results in the @xmath0-process model in figure [ figure : exp_comp ] . since we fix the mass with velocity larger than @xmath274 , the bolometric luminosities at the time @xmath276 are almost the same values . in the thin - diffusion phase ( @xmath277 ) , both the luminosity and the temperature are smaller than the fiducial model . this is because the density at the front of the ejecta is reduced in this mass profile . since the maximum velocity effectively becomes large and the adiabatic cooling becomes efficient , these effects for luminosity and temperature should be also seen in the engine model . we conclude that the luminosities in the thin - diffusion phase ( @xmath277 ) have uncertainties at least with @xmath278 mag , which originates from the uncertainty of the outermost mass profile . note that the emission from the ejecta with the mass profile discussed here reduces the tension between the light curve in the @xmath0-process model and the upper limits of the deep optical observations ( grb 050509 and grb 080905a ) as discussed in section [ other ] . especially , the optical luminosity in the case @xmath279 ( thin - dot - dashed line ) significantly decreases after @xmath280 s. therefore , the ejecta with relatively shallow mass distribution at the front of the ejecta is able to explain the current optical follow - up observations in the @xmath0-process model .      in the fiducial parameter set , the light curves for two models in the optical and infrared bands are similar ( figure [ figure : observation ] ) . in figure [ figure : parameter ] , the allowed parameter region to explain the observation of grb 130603b for both models are also overlapped . recall that emission from ejecta is described as blackbody radiation for the two models , whose spectrum is narrow in bands . therefore , no excess at other wavelengths is expected and also the prediction is consistent with radio observations @xcite . therefore , it is difficult to discriminate two models from the currently available observational data .    in the @xmath0-process model , the light curve with mass profile @xmath35 ( @xmath281 ) and a parameter set which explains the infrared excess detected from grb 130603b can not explain the upper limits obtained from the deep optical observations of some short grbs . therefore , if both stringent optical upper limits at @xmath226 s and bright infrared emission at @xmath84 s are simultaneously obtained from a single event ( with a difference larger than two magnitudes @xmath282 mag ) , the @xmath0-process model is significantly restricted . for the engine model , these observations give a constraint for the temperature distribution in the ejecta , which may give new insights into the activity of the central engine .    as shown in the middle panel of figure [ figure : evolution ] , the bolometric luminosity in the @xmath0-process model from the whole ejecta ( blue long - dashed line ) , including low temperature and/or x - rays and @xmath194-rays produced directly in radioactive decays @xcite , declines more gradually than that for the engine model . this is because there is no energy injection after the time @xmath283 for the engine model . then , the luminosity significantly decreases when photons at the inner edge of the ejecta begin to diffuse out ( see the middle panel of figure [ figure : evolution ] and figure [ figure : observation ] ) . the luminosity from the whole ejecta can be described as @xmath193 in the transparent phase . the index of time @xmath34 is determined by the nuclear heating rate , @xmath284 . therefore , the two models are distinguishable by observing the temporal evolution of bolometric luminosity from the whole ejecta in this phase . we calculated the light curves of macronovae by developing analytical models . we modeled the ejecta based on the results of numerical simulations for a merger of binary nss . in addition to the nuclear decay of @xmath0-process elements ( the @xmath0-process model which is often discussed ) , we considered another heating mechanism for the ejecta , the engine - driven shock ( engine model ) . we compared the results with the optical and infrared observations of the first macronova candidate associated with grb 130603b , and showed that both models can explain the observations . in order to reproduce the observed light curve , the @xmath0-process model requires relatively large ejecta mass @xmath285 which is mainly determined by the observed infrared luminosity @xmath27 erg s@xmath22 at @xmath84 s. in the engine model , the internal energy of ejecta , which mainly determines the observed luminosity , does not depend on the ejecta mass . then , unless the entire of the ejecta is effectively thin ( the diffusion time is smaller than the dynamical time , @xmath286 , at the inner edge of the ejecta ) , s in the engine model , the luminosity rapidly decreases after the time @xmath287 when photons can diffuse out from the inner edge of the ejecta ( section [ transparent ] ) . ] the required ejecta mass is @xmath237 , which is comparable to the recent numerical simulation results . the initial internal energy @xmath97 and the injection time @xmath13 are required as @xmath288 , which is consistent with the observed extended emission of short grbs , @xmath10 erg and @xmath11 s. the required minimum velocity is about @xmath289 for both models , which is mainly determined by the constraint for the observed temperature @xmath29 k at @xmath84 s. for the range of the opacity @xmath250 cm@xmath2 g@xmath22 , it is difficult for both models to explain the observations of macronova associated with grb 130603b by the ejecta mass less than @xmath241 . if macronovae are identical , the upper limits on the luminosity obtained in the deep optical observations of other short grbs give stringent constraints on the @xmath0-process model . on the other hand , the engine model satisfies these constraints if the temperature profile is centrally concentrated in the ejecta ( large @xmath95 ) . thus , if the difference between the optical magnitude at @xmath226 s and the infrared magnitude at @xmath84 s is larger than @xmath290 mag in a single event , the @xmath0-process model is difficult to explain the observations unless the front of the ejecta has much shallow mass distribution . another difference in the light curves between two models is the bolometric luminosity at the transparent phase when dynamical time is smaller than the diffusion time at the inner edge of the ejecta @xmath198 . although the optical and infrared luminosities rapidly decrease in the transparent phase , the bolometric luminosity from the whole ejecta , including lower frequency than near - infrared band and/or x - rays and @xmath194-rays produced directly in radioactive decays , is determined by the energy injection rate of nuclear decay , @xmath291 . for the engine model , the bolometric luminosity decreases rapidly in this phase ( faster than @xmath292 ) . therefore , we expect that the light curve of the bolometric luminosity from the whole ejecta can distinguish between two heating mechanisms . our results show that early light curves depend on the density profile of the outermost edge of the ejecta . it is necessary to develop a method to calculate the low - density region of the ejecta in either the analytical or numerical ways in order to precisely predict the early light curves of macronovae . we thank k. asano , k. kashiyama , k. kiuchi , h. nagakura , t. nakamura , y. sekiguchi , m. shibata for fruitful discussions . this work is supported by kakenhi 24103006 ( s.k . , k.i . ) , 24.9375 ( h.t . ) , 24000004 , 26247042 , 26287051 ( k.i . ) . we summarize the formula for the observed temperature and bolometric luminosity . the detailed derivation of equations in this section is described in sections 2 and 3 . since we assume that the observed luminosity and temperature approximate the luminosity and temperature at the diffusion radius @xmath70 ( section 3 ) , we need to calculate the diffusion radius . for the dynamics of the ejecta , we assume an isotropic and homologous expansion . then , the velocity of ejecta @xmath12 is described by equation ( [ sec2:v ] ) @xmath293 where the radius @xmath0 originates the central engine and the time @xmath34 is measured from the time when a compact binary merges . as in section [ fiducial ] , we calculate the diffusion radius @xmath70 from the condition that the diffusion time equals the dynamical time , @xmath71 . the diffusion time is described by equation ( [ sec2 - 1:t_diff ] ) as @xmath294 where @xmath295 is the speed of the light , @xmath77 is the optical depth described by equation ( [ sec2 - 1:tau ] ) as @xmath296 and @xmath65 is the width of the diffusion region described by the equation ( [ sec2 - 1:deltar ] ) as @xmath297 in the calculation of the optical depth @xmath77 , we use the spatially uniform value of the optical depth @xmath80 with grey approximation and the ejecta mass density @xmath298 described by equation ( [ sec2 - 1:rho ] ) as , @xmath299 where @xmath50 and @xmath51 are normalized factors , and @xmath43 is the velocity at the inner edge of the ejecta . the radius @xmath73 is the outer edge of the ejecta , described by equation ( [ sec2 - 1:r_out ] ) as @xmath300 where the velocity @xmath39 is at the outer edge of the ejecta . the radius at the inner edge of the ejecta @xmath198 is described by equation ( [ sec2 - 1:r_in ] ) as @xmath301 the normalization factor @xmath52 in the profile of the mass density is determined by the ejecta mass @xmath53 ( in equation [ sec2 - 1:m_ej ] ) as @xmath302    we consider two heating sources of the ejecta . in the @xmath0-process model , the internal energy of the ejecta is determined by the nuclear heating rate of the @xmath0-process element described in equation ( [ sec2 - 1:dotepsilon ] ) as @xmath303 for the engine model , we assume the temperature profile of the ejecta @xmath304 as a result of the activity of the central engine , described in equation ( [ sec2 - 1:t ] ) as @xmath305 where @xmath96 is the normalization factor . this factor is determined by injected internal energy @xmath97 at the time @xmath13 , which is described in equation ( [ sec2 - 1:e_int ] ) as @xmath306    observed temperatures in the @xmath0-process model and the engine model are given by @xmath307^{1/4 } } \\ ~~~~~~~~~~~~~~~(r_{\rm diff}>r_{\rm in } ) \\    \\ { \displaystyle \left[\frac{\dot{\epsilon}t\rho(t , v_{\min})}{a}\right]^{1/4 } } \\ ~~~~~~~~~~~~~~~(0.5r_{\rm in } < r_{\rm diff } \le r_{\rm in } ) \\ \\ { \displaystyle 0 } \\ ~~~~~~~~~~~~~~~(r_{\rm diff } \le 0.5r_{\rm in } ) , \\ \end{array } \right.\end{aligned}\ ] ] and by @xmath308 respectively . note that we do not use the approximation @xmath309 and @xmath310 in the thin - diffusion case . the bolometric luminosities for the @xmath0-process model and the engine model are given by @xmath311 and @xmath312 respectively . an example of the calculated result is shown in figure [ figure : evolution ] . we present the numerical values with the parameter dependence for later use . unlike equations ( [ sec2 - 2 - 2:t_times ] ) and ( [ sec2 - 2 - 3:t_tr ] ) , we include the contribution from subdominant terms to the numerical values when we integrate equations . some of the subdominant terms include the ratio @xmath313 . hereafter , the value @xmath314 in subdominant terms are fixed and are not included in the parameter dependence . we introduce the normalized quantities @xmath164 , @xmath315 , @xmath166 , @xmath316 cm@xmath2 g@xmath22 , @xmath317 erg and @xmath318 s. for other parameters , we fix the index of the mass density profile @xmath41 and the parameters of the nuclear heating rate @xmath224 erg s@xmath22 g@xmath22 and @xmath87 . we also introduce the normalized time @xmath319 s and @xmath320 s. the values of observed temperature and bolometric luminosity in the thin - diffusion phase are @xmath321 and @xmath322 respectively . the transition time from the thin - diffusion phase to the thick - diffusion phase @xmath159 is @xmath323cv_{\max}}\left(\frac{v_{\max}}{v_{\min}}\right)^{3-\beta } } \nonumber \\   \nonumber \\ & \sim&4.53\times10 ^ 5~\kappa_{10}^{0.5}m_{\rm ej,0.1}^{0.5}v_{\min,0.1}^{0.25}v_{\max,0.4}^{-0.75}~{\rm s}.\end{aligned}\ ] ] the values of the observed temperature and bolometric luminosity in the thick - diffusion phase are @xmath324 and @xmath325 respectively . note that since the diffusion radius @xmath70 can not be analytically described in the thin - diffusion phase , we use the approximations @xmath309 and @xmath310 in equations ( [ app : t_obs , thin ] ) and ( [ app : l_bol , thin ] ) . these approximations make discontinuity at the transition time @xmath159 . the ratios of the temperature in the thick - diffusion phase to the temperature in the thin - diffusion phase for the @xmath0-process @xmath326 and the engine model @xmath327 at the time @xmath159 are @xmath328 and @xmath329 respectively . the ratios of the luminosity in the thick - diffusion phase to the luminosity in the thin - diffusion phase for the @xmath0-process model @xmath330 and the engine model @xmath331 at the time @xmath159 are @xmath332 and @xmath333 respectively . the transition time from the thick - diffusion phase to the transparent phase @xmath188 is @xmath334cv_{\min } } } \nonumber \\ \nonumber \\   & \sim&7.62\times10 ^ 5~\kappa_{10}^{0.5}m_{\rm ej,0.1}^{0.5}v_{\min,0.1}^{-0.5}~{\rm s}.\end{aligned}\ ] ] we introduce another transition time @xmath335 when the upper limit of the integral for the luminosity @xmath123 reaches the inner edge of the ejecta @xmath198 , @xmath336 the values of the observed temperature and bolometric luminosity in the transparent phase ( @xmath337 ) are @xmath338 and @xmath339 } \\ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(r-{\rm process } ) \\    \\ { \displaystyle 1.33\times10^{41}~~{\rm erg~s}^{-1 } }   \\ ~~{\displaystyle \times e_{\rm int0,51}t_{\rm inj,2}t_6^{-2}\left[1-\left(\frac{t}{t_{\rm tr2}}\right)^{1.333}\right ] } \\ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~({\rm engine},~\xi=1 ) , \\    \\ { \displaystyle 1.00\times10^{41}~~{\rm erg~s}^{-1 } }   \\ ~~{\displaystyle \times e_{\rm int0,51}t_{\rm inj,2}t_6^{-2}\left[1-\left(\frac{t}{t_{\rm tr2}}\right)^{6.667}\right ] } \\ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~({\rm engine},~\xi=2 ) , \\    \\ { \displaystyle 1.00\times10^{41}~~{\rm erg~s}^{-1 } }   \\ ~~{\displaystyle \times e_{\rm int0,51}t_{\rm inj,2}t_6^{-2}\left[1-\left(\frac{t}{t_{\rm tr2}}\right)^{12}\right ] } \\ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~({\rm engine},~\xi=3 ) , \\ \end{array } \right.\end{aligned}\ ] ] respectively .
Abstract:

Answer: a macronova ( kilonova ) was discovered with a short gamma - ray burst , grb 130603b , which is widely believed to be powered by the radioactivity of @xmath0-process elements synthesized in the ejecta of a neutron star binary merger . as an alternative , we propose that macronovae are energized by the central engine , i.e. , a black hole or neutron star , and the injected energy is emitted after the adiabatic expansion of ejecta . this engine model is motivated by extended emission of short grbs . in order to compare the theoretical models with observations , we develop analytical formulae for the light curves of macronovae . the engine model allows a wider parameter range , especially smaller ejecta mass , and better fit to observations than the @xmath0-process model . future observations of electromagnetic counterparts of gravitational waves should distinguish energy sources and constrain the activity of central engine and the @xmath0-process nucleosynthesis .