[ sec : intro ]    the collapsar model @xcite was originally proposed as an explanation for common , `` long - soft '' gamma - ray bursts ( grbs ) with durations near 20 s. this is a typical time scale for the gravitational collapse of the carbon - oxygen core of a massive star and the congruence of these time scales was an argument in favor of the model .
it is well known , however , that forming a disk so quickly requires an unusually large amount of angular momentum to be stored in the core up to the point when the star dies ( e.g. , * ? ? ?
this is presumably one reason why grbs are rare compared with supernovae .
black holes are expected to form in a significant fraction of massive star deaths ( e.g. , * ? ? ?
* ) and it seems possible that , in at least some of these , the black hole will form without generating an outgoing shock . what observable signal would accompany such an event ?
would they all be `` un - novae '' @xcite , stars that simply suddenly disappear from view ? or will some sort of transient display announce the birth of any black hole ? here
we explore the possibility that these black hole systems form an accretion disk , but only after almost all of the star has accreted into the hole .
torques in differentially rotating contracting stars tend to concentrate angular momentum in the star s outer layers .
the surface will also have a large angular momentum in a tidally locked system . on the other hand
, stellar mass loss preferentially removes these same layers , so , just as in the usual collapsar model , reducing mass loss , e.g. , by reducing the metallicity , favors the retention of high angular momentum at the stellar surface .    because the outer layers of the star take longer to fall in - months to years in the case of a red supergiant -
the duration of the transient produced by accretion through a disk in these cases is much longer than for the common grbs studied so far .
they might even be confused with other forms of gamma - ray transients like those arising from super - massive black hole accretion in active galaxies @xcite .    here
we shall consider four varieties of long duration transients categorized by the kind of system where they occur : 1 ) single supergiants with metallicity less than about @xmath4 and low mass loss (   [ sec : bsg ] ) ; 2 ) red or blue supergiants in tidally - locked , or nearly tidally - locked binaries (   [ sec : rsg ] ) ; 3 ) pair - instability induced collapse to a black hole (   [ sec : pair ] ) ; and 4 ) helium stars in tidally locked binaries (   [ sec : helium ] ) . in each case , the duration of the transient that is produced is correlated with the radius of the presupernova star , hence shortest for the helium stars , while the luminosity is inversely correlated with the radius , hence faintest for the red supergiants . as with common grbs , we are unable to predict _ ab initio _ the spectrum of these transients or the efficiency for converting black hole accretion rate into electromagnetic radiation . presumably , as in other forms of transients powered by very rapidly accreting , rotating black holes , accretion energy and black hole spin energy power relativistic bi - polar outflows .
poynting flux , shocks within these outflows , and shocks with the circumstellar medium make x - rays and gamma - rays .
it is common to take an efficiency of @xmath5 to @xmath6 times @xmath7 for the efficiency ( e.g. , * ? ? ?
* ) thus suggesting an integrated power of @xmath8 per solar mass accreted .
interestingly these `` failed supernovae '' @xcite are far more brilliant and powerful than their ( model - wise ) more successful cousins , ordinary supernovae .    the idea that late time accretion in supernovae that make black holes might power long gamma - ray ( or hard x - ray ) transients
was suggested by @xcite .
the transients discussed there , referred to as `` type 2 collapsars '' , were produced because of fall - back in an explosion that launched a relatively strong outgoing shock in a red supergiant progenitor .
two dimensional models explored jet propagation and break out for a variety of fall - back rates , but little attention was given to the angular momentum distribution necessary to form a disk . here we consider a wider variety of progenitors and pay closer attention to the angular momentum , especially in the outer layers .
@xcite also discussed the possibility of very - long duration x - ray transients powered by fall back and accretion through a disk .
he suggested these might look like `` ultra - luminous x - ray sources '' seen in nearby galaxies . these would be much longer in duration ( thousands of years ) and have much lower luminosity than the events considered here , but are an example of the same general idea .    more recently @xcite and @xcite suggested that fallback or a failed explosion in a red supergiant could give a long duration gamma - ray transient that might be confused with a gamma - ray blazar .
quataert and kasen , in particular , suggest this as an alternate interpretation of swift transient sw-1644 + 57 @xcite .
woosley s paper was accepted prior to the discovery of this event .
a variety of other models may lead to long duration bursts .
for example a white dwarf @xcite or helium star merging with a black hole or neutron star @xcite . a neutron star or black hole merging with a a giant star might give a thorne - zytkow object @xcite or a black hole accreting at a rapid rate @xcite .
helium star merger with a compact object also gives a long transient of some sort @xcite .
here we consider only those transients coming from stars that are not actively merging with a companion .
[ sec : black ]    a black hole will form inside a massive star when a sufficient concentration of matter and energy accumulates to form an event horizon .
the simplest possibility is that some portion of a massive star experiences an instability and falls freely , or almost freely , until it reaches such a high density that a trapped surface is formed .
this rarely happens in nature though . only in the most massive stars experiencing a pair instability @xcite or a general relativistic instability @xcite can a trapped surface form before the center reaches nuclear density and bounces .
that is , black hole formation in a massive star is almost always preceded by protoneutron star formation @xcite .    for stars in the range @xmath9 ,
however , the protoneutron star that is formed experiences rapid accretion from the surrounding layers of heavy elements - silicon , oxygen , etc .
if rotation , magnetic fields and neutrinos do not readily produce an outgoing shock , the hot protoneutron star grows beyond a critical mass and collapses to a black hole . the accretion rate is larger for more massive helium cores @xcite , and
black hole formation without an outgoing shock may be a frequent occurrence in stars with main sequence masses over about @xmath10 @xcite . whether _ all _ stars above about @xmath10 produce black holes without outgoing shock waves is uncertain , but the possibility is not precluded by observations or theory
. `` prompt ''
black hole formation is probably more likely in very low metallicity stars that have experienced little mass loss and thus have especially large helium cores , but it could occur for solar metallicity stars , depending upon the uncertain mass loss rate @xcite .
it is also possible that such a weak outgoing shock is generated that most of the star fails to explode and a black hole is made by fall back @xcite .
these were the systems explored by @xcite .
the ones treated here differ in assuming that any outgoing shock is insufficient to eject more than a small fraction of the surface layers of the star .
we shall call such events , where only the surface layers have enough rotation to form a disk and a shock so weak that these layers are not ejected `` type 3 collapsars '' .
this is to distinguish them from ordinary grbs where the disk forms promptly from material deep inside a highly differentially rotating helium core ( type 1 ; * ? ? ?
* ; * ? ? ?
* ) and from systems where longer duration transients result from the delayed fall back of other material also deep inside the star ( type 2 ; * ? ? ?
* ; * ? ? ?
[ sec : stars ]    the preservation of angular momentum in the surface layers of the star can occur either because mass loss is inefficient in a single star or because of tidal interaction with a close companion .
here we consider examples of both sorts .
[ sec : bsg ]    a paucity of metals in a star affects its evolution in several ways - reduced energy generation by the cno cycle in the hydrogen burning shell , reduced opacity , and reduced mass loss . of these effects
, the reduction in mass loss is least well understood , but the lack of atomic lines in the atmosphere on the main sequence and of grains as a red giant seems certain to lead to some reduction .
a common assumption is that the mass loss for massive main sequence stars depends upon some fractional power of the metallicity @xcite . more important and less certain
is the dependence of mass lost as a red or blue supergiant on metallicity .
it is thought that mass loss from cool giants is more dependent upon pulsations and grain formation @xcite .
one might expect therefore a rapid fall off in mass loss below some value necessary for significant grain production .
it has been estimated that red giant mass loss will be significantly less below @xmath4 @xcite .
if stars do not lose mass then they conserve angular momentum .
since the natural course of evolution leads to the contraction and spin up of the inner star , shear instabilities and magnetic torques will concentrate an increasing amount of angular momentum in the outer regions of the star .
whereas the inner part does spin faster due to contraction , it actually loses most of its initial angular momentum by the time it reaches central carbon burning .
@xcite surveyed the evolution of non - rotating zero - metal massive stars from @xmath11 to @xmath12 . in a parallel study currently underway @xcite
, we are surveying the evolution of rotating stars in the same mass range with metallicity @xmath13 , @xmath14 and @xmath15 that of the sun .
the two stars discussed here are extracted from that survey and are typical in the way that they accumulate large amounts of angular momentum in their outer layers before dying as supergiants .
all calculations of stellar evolution presented here used the kepler code @xcite and included angular momentum transport by magnetic torques @xcite .
models v24 and v36 are @xmath16 and @xmath17 main sequence stars with metallicity @xmath18 solar .
both stars rotated rigidly on the main sequence with a moderate speed of about @xmath19 which is about @xmath20 of the keplerian value .
both presupernova stars had hydrogenic envelopes that were , throughout most of their mass , radiative .
model v36 , however , had a low mass surface convection zone that included @xmath21 .
model v36 was consequently a yellow supergiant at death ( @xmath22 ; @xmath23 ; @xmath24 ) .
model v24 was a rather large blue supergiant ( @xmath25 , @xmath26 ; @xmath27 ) .
these two stars ended their lives with cores of helium and heavy elements of @xmath28 and @xmath29 respectively and angular momentum distributions as shown in fig .
[ fig : bsg ] . whereas the lack of sufficient angular momentum within the helium core precludes making a disk around a black hole within that mass
, there is ample rotation in the outer part of the hydrogen envelope to do so .
this could be a very common occurrence . from the survey of @xcite , the total angular momentum of a massive star at birth that has an equatorial rotation speed of about @xmath20 keplerian on the main sequence is @xmath30
if the star does not lose mass , most of this angular momentum becomes concentrated , in the presupernova star , in a nearly rigidly rotating hydrogenic envelope with radius either @xmath31 ( blue supergiant ) or @xmath32 ( red supergiant ) . whereas the density declines with radius in the actual envelope
, one can obtain some interesting scaling relations by assuming constant density . for a moment of inertia , @xmath33 ,
the angular velocity at the surface , @xmath34 , is @xmath35 the specific angular momentum of this rigidly rotating envelope is @xmath36 , or about @xmath37 at the surface of a red supergiant and @xmath38 for a blue supergiant , if no mass is lost .
these can be compared with the angular momentum required to make a disk around a non - rotating black hole , @xmath39 or @xmath40 for a maximally rotating hole ( kerr parameter @xmath41 ) .
all massive stars that do not lose mass have sufficient angular momentum in their outermost layers to make a disk around any black hole formed in their collapse .
of course stars do lose mass and the fact that , for constant @xmath42 , @xmath43 means that the first mass to be lost contains most of the angular momentum . in practice , we find that a red supergiant that loses more than a few percent of its total mass before dying will probably not make a black hole - disk system .
this sensitivity to mass loss is illustrated by the evolution of a @xmath44 model , o25 , with @xmath6 solar metallicity in a calculation that included mass loss at one - half the standard value for that metallicity .
the star lost @xmath45 ( @xmath46 ) of its mass and died as a red supergiant with a radius of @xmath47 .
the value of @xmath48 in its outermost layers was @xmath49 , too little to form a disk around a schwarzschild black hole , but possibly just enough to form a disk around a mildly rotating one .
the point where @xmath48 declined to one - half its surface value in the presupernova star was located @xmath50 interior .
this is borderline case .
stars with less mass loss , i.e. , lower metallicity , and stars that never become red giants could make long gamma - ray transients ( table  [ tab : models ] ) .
[ sec : rsg ]    based upon theoretical models for single massive stars , the red supergiant progenitors of common type iip supernovae are not expected to have rapidly rotating envelopes .
expansion to a red giant slows the surface layers .
torques transfer the angular momentum of the more rapidly rotating inner regions outwards concentrating it in the outermost layers . for stars of sufficient metallicity ,
stellar winds are efficient at removing these layers and , along with them , most of the star s angular momentum .
for example , a @xmath51 solar metallicity star rotating rigidly with an equatorial speed of @xmath19 on the main sequence ends up a @xmath52 supernova progenitor if commonly used estimates for mass loss are employed @xcite .
the surface rotation rate of that red supergiant ( r = @xmath53 ) is @xmath54 and its angular momentum at the surface is @xmath55 .
nowhere in the star does any matter have enough angular momentum to form a disk around a black hole .
if a black hole forms without making an outgoing shock , the star abruptly disappears .
this sort of `` un - nova '' might be a common occurrence in the modern universe @xcite , though see @xcite
.    the outcome can be quite different in a binary system , however .
a stellar merger could lead to a rapidly rotating envelope @xcite , but such mergers probably lead to ejection of a common envelope .
an interesting alternative is detached binary systems with a small enough separation for tidal interaction , but too distant to merge or exchange mass .
the circularization and synchronization time scales for binary stars , one of which is a red giant , have been considered by @xcite , @xcite and @xcite . for a primary with mass ,
@xmath56 , envelope mass @xmath57 , luminosity @xmath58 , and radius @xmath59 , in an orbit having radius @xmath60 , with a companion of mass @xmath61 , the _ circularization _ time , @xmath62 is @xcite @xmath63 where @xmath64 is a number close to unity .
the lifetime of a @xmath51 main sequence star as a red supergiant is approximately @xmath65 , depending upon uncertain convection parameters .
the mass of the envelope that is convective at any one time , however , is only about @xmath66 and the star loses more than this in its lifetime .
a more appropriate time scale might be the time the star spends losing its last @xmath66 , which is about @xmath67 .
the mass of the star when it becomes a supernova is @xmath52 and we assume a companion mass of , e.g. , @xmath68 .
the luminosity of the star as a red supergiant is @xmath69 and its moment of inertia is approximately @xmath70 .
equation  ( [ eq : circ ] ) then implies that during its last @xmath67 , orbits with radius @xmath71 will circularize .
the _ synchronization _ timescale for co - rotation is somewhat different @xcite and , in the present case , is easier to achieve .
@xmath72 most of the angular momentum is in the hydrogen envelope and , from the same @xmath51 model , @xmath73 . with @xmath74 and @xmath75 @xcite ,
@xmath76 which only requires @xmath77 for synchronization in @xmath67 .
in such a system , the period will be shorter than about @xmath78 years ( @xmath79 ) .
systems much closer will merge or experience mass exchange , so we take this to be a typical rotation period for the envelope in a detached binary system that might achieve corotation .
longer period systems will still transfer significant rotation to the envelope , however .
this is about four orders of magnitude more angular momentum in the outer envelope than in the case of the single star .
since the outer several solar masses of the presupernova star is convective , this rotation rate will persist to some depth .
all of that mass would have sufficient angular momentum to form an accretion disk if its helium core collapsed to a black hole . even then though
, the total angular momentum of the envelope , @xmath80 , would be a small fraction of the orbital angular momentum of the binary pair of stars .
the addition of only @xmath5 of the angular momentum estimated here would still result in large quantities of material forming a disk , so the range of binary separations allowed is actually larger , perhaps by a factor @xmath81 , or about two .
since the time scale for material to fall in from @xmath82 is longer than any realistic viscous time scale for the disk , the accretion rate will be given by the collapse timescale of the envelope .
for matter at the base of the hydrogen shell where a disk would first form , the enclosed mass is @xmath83 and the radius about @xmath84 . at the outer edge of the star
the mass is @xmath52 and the radius , @xmath85 .
the corresponding free fall times ( @xmath86 ) are @xmath87 and @xmath88 suggesting that the event would last for months .
accreting @xmath66 in @xmath89 with @xmath6 efficiency for converting accreted mass to outgoing energy would give a jet power of @xmath90 .
conversion of only @xmath6 of the jet power ( @xmath5 of the accreted mass energy ) into gamma - rays would explain the recently discovered transient sw 1644 + 57 @xcite , especially if moderate beaming were involved .
whereas the focus here is on red supergiants because they give the longest transients , tidal locking would also occur in blue supergiants with smaller radii .
the key point is that forced co - rotation with a body in keplerian orbit gives , for a large range of orbital separations , envelope material that is rotating too fast to accrete , without hindrance , into a stellar mass black hole .
since the time scale for accretion is set by the collapse time for the outer part of the star , not the viscous time of the disk , a transient from a tidally locked blue supergiant would be similar in duration within a factor of a few to those studied for single blue supergiants in   [ sec : bsg ] .
blue supergiants , however , will have less mass in the outer layers of the envelope and hence should have less powerful outbursts ( appendix  [ sec : append ] and @xcite ) .
[ sec : pair ]    following helium burning , non - rotating helium cores with masses greater than @xmath91 will collapse directly into black holes with no outgoing shock @xcite .
the collapse is caused by the pair - instability and , above this mass , nuclear burning is unable to reverse the implosion . rotation ,
unless it is very rapid and highly differential , does not change this mass limit greatly since the instability exists in the deep interior where the ratio of centrifugal force to gravity is small .
the collapse of a rotating pair unstable star to a black hole was followed in two dimensions by @xcite .
first , a @xmath92 main sequence star was evolved to the point of instability in the kepler ( 1d ) code .
this star produced a helium core of @xmath93 .
the initial star was assumed to rotate rigidly on the main sequence with a ratio of equatorial speed to keplerian of 20% .
angular momentum transport , especially by shear mixing and eddington sweet circulation , was followed throughout the evolution , but magnetic torques were neglected . as a result , the collapsing helium core rotated sufficiently rapidly to produce a disk outside of a massive black hole once about @xmath94 of the @xmath93 core of helium and heavy elements had collapsed .
@xcite speculated that the accretion of the remaining @xmath95 might produce a very bright , lengthy grb .
we have recalculated the evolution of similar stars , but with magnetic torques @xcite included in the simulation .
models z250a , b , and c are stars with zero metallicity and a mass at the time they formed of @xmath96 .
they differ in the amount of angular momentum each had on the main sequence .
z250a had @xmath97 corresponding to a rotation rate of @xmath98 , or @xmath99 keplerian .
models z250b and z250c had @xmath100 and @xmath101 of angular momentum and rotational speeds of @xmath102 ( @xmath103 keplerian ) and @xmath104 ( @xmath105 keplerian ) .
these made helium cores of @xmath106 , @xmath107 , and @xmath108 respectively , all of which collapsed directly to black holes following helium depletion .
unlike in @xcite however , none of these helium cores had sufficient angular momentum to make a disk around a black hole anywhere in their interior . because mass loss was neglected in these models , the angular momentum lost from the helium core due to magnetic coupling instead ended up concentrated in the outer part of the hydrogen envelope .
[ fig : pair ] shows the evolution of the angular momentum in model z250b .
the evolution of the other two models was similar .
inadequate angular momentum exists in the helium core to form a disk but , for model z250b , the outer @xmath109 can still do so .
as in the tidally locked red supergiant case , the collapse time scale for this material is days to months , so a high energy transient of around @xmath110 , possibly boosted by beaming , could exist for that time .
it should be noted , however , that mass loss is even more uncertain in these models than in the others . due to
rotationally induced mixing , the envelopes of all three models were appreciably enhanced in c , n , and o and thus were red supergiants when they died . neglecting mass loss in such a situation is probably not justified .
also there may be other mechanisms for losing mass in such massive stars besides line and grain driven mass loss @xcite .
[ sec : helium ]    helium stars in close tidally locked binaries have been frequently suggested as progenitors for common grbs @xcite . on the positive side
, such events should occur frequently in nature and , if the core of a massive wolf - rayet star collapses to a black hole while the surface is forced to rotate at a fraction of its keplerian speed , disk formation around a black hole is assured . on the negative side , however , the production of a grb with typical duration @xmath010 s , requires that the _ inner _ part of the helium star rotate rapidly , not just its surface .
the inclusion of magnetic torques in tidally locked stars leads , by itself , to too little rotation to make a disk in that part of the star which might accrete rapidly @xcite . on the other hand , leaving out magnetic torques
makes it difficult to produce the spin rates of ordinary pulsars @xcite , and could lead to an _ overabundance _ of grb progenitors .
this is not to say that close binary systems do not produce _ some _ form of gamma - ray transient though
. the angular momentum in their cores may be adequate to produce a millisecond magnetar , and thus a common grb , just not a ( type 1 ) collapsar .
if a black hole forms , the accretion of the _ outer _ layers of the star may lead to the formation of a disk and grb much longer than typical .
it is this latter possibility that we explore here .
following @xcite , we consider the evolution of helium cores of @xmath111 and @xmath109 whose surface s have become tidally locked with a closely orbiting companion .
@xcite found the shortest possible periods were 2.05 hours ( @xmath112 ) and @xmath113 hours ( @xmath114 ) for @xmath111 and @xmath109 respectively , and they gave examples of systems containing a compact object and a helium core that might approach these limiting values .
an @xmath111 he core in a tidally locked binary with a @xmath115 companion filling its roche lobe had a longer period of @xmath116 hours ( @xmath117 ) . here
we calculate two models for each helium core mass .
model 8a has an enforced rotation rate in its outer @xmath118 of @xmath119 throughout its evolution .
models 8b , 16a , and 16b have rotation rates of @xmath120 , @xmath121 , and @xmath122 respectively .
[ fig : binaryw ] shows the evolution of the angular momentum distribution as a function of mass for these four models . during helium burning ,
magnetic coupling maintains a state of nearly rigid rotation with an angular speed given by the surface boundary condition . for both model 8a and 16a ,
this corresponds to a ratio of rotational speed to keplerian of @xmath123 . for models 8b and 16b
, the ratio is @xmath124 .
roughly the inner @xmath125 ( @xmath111 ) or @xmath126 ( @xmath109 ) of the helium star is convective during helium burning and rigid rotation is enforced in that region .
the outer part of the star is radiative , but magnetic torques there suppress differential rotation . during carbon burning and more advanced stages , the inner core contracts and
rotates more rapidly .
magnetic torques are too inefficient to enforce rigid rotation .
later , during carbon burning , the stars ignite a convective helium burning shell which persists until the star dies .
this convective shell does not extend all the way to the surface though , but does extend beyond the @xmath118 at the surface where corotation is enforced .
thus all stars die with an angular velocity in their helium convective shell given by the surface boundary condition , but have larger rotation rates in their cores , due to contraction during advanced burning stages .
[ fig : binaryj ] shows the corresponding angular momentum in the evolving models .
shear instabilities and magnetic torques transport angular momentum out of the inner core and to the surface . in the absence of a binary companion ,
the surface would have actually spun up , but would have been removed by mass loss . in the present parametrized calculations though , there is no mass loss and the set rotation rate of the surface absorbs any excess angular momentum delivered to it .
consequently the total angular momentum of the star decreases with time . for model 8a at helium ignition , it was @xmath127 , but for the presupernova star it was @xmath128 .
[ fig : binaryjfin ] shows the final distribution of angular momentum in the stars compared with that required to make a disk around a black hole , once the mass interior to the given sample point has collapsed .
four curves are given .
three are for the minimum angular momentum required to form a disk at the last stable orbit of a ) a non - rotating ( schwarzschild ) black hole ; b ) a kerr black hole ; and c ) a black hole with rotation given by the actual angular momentum distribution of the presupernova star .
the fourth curve is the angular momentum distribution in the presupernova star . because the inner parts of the star are not rotating very rapidly , curve a ) and c )
are similar except near the surface . only where the actual angular momentum lies above curve c ) can a disk form .
material with less angular momentum plunges directly into the hole .
none of the models has enough angular momentum to form a disk in its inner core .
thus a type  1 collapsar will not occur .
the angular momentum in the inner @xmath129 of the stars is still appreciable , however .
models 8a and 16a would form @xmath130 ( gravitational mass ) pulsars with periods of @xmath131 and @xmath132 respectively .
models 8b and 16b would form pulsars of @xmath133 and @xmath134 .
model 16b thus qualifies as a `` millisecond magnetar '' candidate and might power a common grb if a pulsar is able to form before accretion turns it into a black hole .
models 8b , 16a and 16b could also make powerful pulsar - powered supernovae .
model 8a and other more slowly rotating systems probably would not .    of greater interest for the present paper , however ,
is the large angular momentum in the _ outer _ layers of three of the models .
these would be relevant if the center of the star collapsed to a black hole rather than a millisecond pulsar .
model 8a has angular momentum sufficient to form a disk in its outer @xmath135 .
two other models ( table  [ tab : models ] ) have similar large rotation in their surface layers .
this is consistent with the results of @xcite who found that their most rapidly rotating @xmath111 and @xmath109 helium cores ended their lives with @xmath48 in outer layers of @xmath136 and @xmath137 .
the radius and gravitational potential of the high angular momentum surface layers allows an order of magnitude estimate of the free fall time scale ( table  [ tab : models ] ) , about @xmath138 .
the tidal interaction actually acts to brake the rotation of the surface layers compared to what they would have had in a star without mass loss and no companion star . to illustrate this , model 8c
was calculated , starting from the rigidly rotating helium burning stage of model 8b ( total angular momentum @xmath139 ) , but with no surface boundary condition and no mass loss .
that is the rotation rate of the surface layers was allowed to adjust to be consistent with whatever angular momentum was transported to them . without mass loss ,
the star was forced to conserve angular momentum overall and the very outer layers ended up rotating very rapidly ( fig .
[ fig : binaryjfin ] ) . the angular speed , rather than being @xmath140 was @xmath141 .
in fact the outer tenth of a solar mass rotated so rapidly that it would be centrifugally ejected . in a realistic calculation , this matter and more underlying matter would probably have been ejected as a disk or a centrifugally boosted wind .
the deep interior still lacked sufficient angular momentum to form a disk around a black hole , but the amount of surface material that could form a disk ( if mass loss did not remove it ) was significant , about 2 @xmath142 , comparable to that in model 8b .
two other models like 8c were calculated that included mass loss appropriate for a wolf - rayet star at solar metallicity and one - tenth that value .
the first ended up with a final mass of @xmath143 and no surface layers with sufficient angular momentum to form a disk .
the second ended up with a mass of @xmath144 and still had enough angular momentum in its outer solar mass to make a disk .
so the major role of the close companion in the calculations done here is to preserve the angular momentum in the outer layers that might otherwise be removed by mass loss .
a companion may also be essential so that the helium core rotates rapidly in the first place , following the loss of the star s hydrogen envelope .
when magnetic torques are included in the calculation , the helium cores of supernova progenitors that are red and blue supergiants rotate too slowly for any portion of that core to make a disk .
a close companion could either preserve angular momentum in the core by removing the envelope very early in the evolution or add angular momentum by tidal locking later on .
[ sec : conclude ]    table  [ tab : models ] and table  [ tab : models1 ] summarize the wide variety of transients that can be produced by type  3 collapsars . in table
[ tab : models ] @xmath145 is the amount of mass capable of forming a centrifugally supported disk around a black hole with a kerr parameter given by accreting all the mass interior to @xmath145 and @xmath59 is the range of stellar radii where that mass exists . the accretion time scale has been estimated using the radius for the center of mass of @xmath145 . the luminosity has been estimated by multiplying the effective accretion rate ( mass / time scale ) by @xmath146 .
this assumes @xmath6 efficiency for making a jet and @xmath6 efficiency for converting jet energy into hard radiation .
most models produce extreme kerr black holes ( @xmath41 , table  [ tab : models1 ] ) , so the radiative efficiencies may be higher than assumed here @xcite .
the numbers in the table are only order of magnitude estimates .
table  [ tab : models1 ] gives the the estimated pulsar period formed in the stellar collapse assuming the angular momentum that exists in the inner @xmath129 of the core , a final gravitational mass of @xmath130 , and a moment of inertia of @xmath147 @xcite . if a black hole forms before the protoneutron star has radiated away most of its gravitational binding energy , the rotation rate will be slower as will the rotational kinetic energy available for explosion . @xmath148 and @xmath149 are the masses of the presupernova star and its helium core respectively and @xmath150 is the lagrangian mass coordinates of the stellar layers that can form a centrifugally supported disk .
the kerr parameter , @xmath151 , is what would develop if the entire presupernova , including its rapidly rotating outer layers , collapsed to a black hole . in case where it exceeds unity
it has been capped at @xmath152 .
the kerr parameter for the helium core alone is @xmath153 .
@xmath154 is the net binding energy , in units of @xmath155 , of the material that could form a disk .
a range of masses can form disks and a large range of time scales characterize the accretion , hence the possible luminosities span many orders of magnitude .
the shorter brighter transients from wolf - rayet stars in tidally locked binaries might comprise an interesting subset of long , but otherwise ordinary grbs , as others have noted .
this channel could be especially important in regions where high metallicity and high mass loss preclude the production of grbs by solitary stars .
the mass of @xmath156ni produced in the collapsar model is sensitive to the temperature , density , time scale , and ejection fraction of the disk and is difficult to estimate .
the lower accretion rate , larger black hole mass , and smaller kerr parameter for the models presented here , however , will give lower temperatures and densities in the disk @xcite .
below about @xmath157 , the nature of the accretion changes and the disk is no longer neutrino dominated .
a thicker disk might qualitatively change the nature of the transient ( macfadyen , private communication ) .
it is thus reasonable to expect that the @xmath156ni production , per solar mass of material accreted , will be different here and perhaps much smaller .
these bursts might not be accompanied by supernovae ( e.g. , * ? ? ?
on the longer end , blue and red supergiants will produce transients of order @xmath158 and @xmath159 respectively .
these could result from either tidally locked systems or more distant systems with low metallicity and mass loss .
these fainter and potentially more common systems may not have been detected yet and certainly should be searched for .
they could easily be confused with flaring gamma - ray blazars @xcite , but would differ in that they only happen once and would be associated with regions of massive star formation , not particularly galactic nuclei .
they too would probably not be accompanied by supernovae and might have a different beaming factor and spectrum from ordinary grbs .
the black holes produced in these supergiant progenitors would frequently be rotating at near the maximum value ( fig .  [
fig : bsg ] and fig .
[ fig : rsg ] ) .
an exciting prospect is that the first generations of stars , those with metallicities from zero to @xmath160 solar might be studied directly using the hard x - rays they emit when they die . whereas predictions are still clouded by the vagaries of mass loss and the explosion mechanism ,
our studies predict that the fraction of massive stars emitting such bursts could be substantial , essentially limited only by the fraction that make black holes promptly .
if none make black holes promptly , then the death of every massive star must make a supernova , though perhaps a faint one .
this too is an interesting possibility worth further exploration .    because of their long duration and potentially high event rate , it is possible that several type  3 collapsars will be visible in distant galaxies at any one time .
@xcite estimate @xmath161 core - collapse supernovae occur each year in a field of view @xmath162 in the redshift range @xmath163 .
this is about @xmath164 supernovae per second spread over the whole sky .
the observation of grbs beyond redshift @xmath165 shows that massive stars were dying at an even earlier time .
black hole formation probably occurs in a significant fraction of these supernovae @xcite , perhaps all those with mass over @xmath10 and metallicity less than @xmath6 solar . in order to estimate , crudely
, an event rate , assume that prompt black hole formation occurs in @xmath6 of core collapses and @xmath5 of these stars either have such low mass loss rates or close companions that their outer layers make a disk .
then @xmath18 of all massive star deaths would make type  3 collapsars . for low metallicity stars
, the latter could be an underestimate .
the universal supernova rate of @xmath164 per second then implies a gamma - ray transient event rate of @xmath166 per second .
if each event lasted @xmath158 , there would @xmath167 events in progress at all times .
most of these would be at high redshift ( low metallicity ) and the spectrum would be softer and duration longer than in the lab frame , so the number of active sources could be higher .
on the other hand , beaming could increase the brightness and decrease the event rate significantly .    in an extreme limiting case where @xmath6 of all massive star deaths made type  3 collapsars with duration @xmath159 with negligible beaming
, one would have @xmath168 active sources at any one time , or about @xmath152 per square degree .
if each source radiated @xmath169 , this could potentially contribute to a `` diffuse '' cosmological background of soft gamma - rays .
the greatest uncertainty in these results is whether the matter with high angular momentum will get ejected , by mass loss either before the explosion or during the explosion , without falling to the center of the star .
as table  [ tab : models1 ] shows , the binding energy of the critical material is very much less than the kinetic energy of a typical supernova ( @xmath170 ) , so a weak explosion or even strong pulsations during silicon burning could eject it .
it seems unlikely that this would happen all of the time in all sorts of models , so perhaps this enters in chiefly as a major uncertainty in the event rates .
we also note that the converse problem , the optical appearance of a `` supernova '' exploding with only @xmath171 of kinetic energy is an interesting one too .
one way or another , it seems unlikely that these stars will all disappear without a trace .
the authors appreciate valuable discussions with christian ott regarding core - collapse in massive stars and with ron taam on tidal interaction in close binaries . this research has been supported by at ucsc by the national science foundation ( ast 0909129 ) and the nasa theory program ( nnx09ak36 g ) .
ah acknowledges support from the doe program for scientific discovery through advanced computing ( scidac ; de - fc02 - 09er41618 ) , by the us department of energy under grant de - fg02 - 87er40328 , by nsf grant ast-1109394 , and by the joint institute for nuclear astrophysics ( jina ; nsf grant phy02 - 16783 ) .
[ sec : append ]    starting from the usual definition of free - fall time , @xmath172 we obtain an order - of - magnitude estimate of the accretion rate from @xmath173 @xmath174 @xmath175 where we used @xmath176 and @xmath177 , and @xmath178 is the mass coordinate and @xmath179 is the radius coordinate of a shell .
note that usually @xmath180 in most parts of the star ; density inversions can exist at the outer edge of convective regions in mixing - length theory especially in rsgs due to recombination and opacity drops . if @xmath181 , as in the case of the outer layers of a bsgs , we obtain @xmath182 , so bsgs may have a more tenuous accretion rate and fainter , softer grbs than rsgs based on the estimate from @xmath183 . on the other hand
the average density of the star is higher as it is more compact , which may be a slight compensation .
density inversions at the edge of rsg envelopes could lead to spikes on accretion rate in this simple approximation .                                                                                of model v24 and
the outer @xmath11 of model v36 will form a disk .
the edge of the helium cores of the two models is apparent in the sharp inflection in the angular momentum at @xmath28 and @xmath29 .
mass loss was included in the calculation , but due to the low metallicity , only @xmath184 and @xmath185 was lost in v24 and v36 , respectively .
note that , if all the surface material accreted here , the black hole would rotate at nearly its maximum allowed value ( i.e. , the red line intersects the green one at the end ) . [
fig : bsg],title="fig:",scaledwidth=47.5% ]   of model v24 and the outer @xmath11 of model v36 will form a disk .
the edge of the helium cores of the two models is apparent in the sharp inflection in the angular momentum at @xmath28 and @xmath29 .
mass loss was included in the calculation , but due to the low metallicity , only @xmath184 and @xmath185 was lost in v24 and v36 , respectively .
note that , if all the surface material accreted here , the black hole would rotate at nearly its maximum allowed value ( i.e. , the red line intersects the green one at the end ) .
[ fig : bsg],title="fig:",scaledwidth=47.5% ]    , the solid dark line shows the actual distribution of angular momentum while the lighter lines show the values required to form a disk around a non - rotating black hole , a kerr black hole , or a black hole with the integrated angular momentum interior to the given mass .
the angular velocity in the envelope is constant with @xmath79 .
note that the black hole formed here would be an extreme kerr hole with @xmath186 .
[ fig : rsg],scaledwidth=90.0% ]                                        ccccccc bsg & v24 & 10.0 & 0.4 - 10 @xmath187 & 20,000 & 5 & 9 + & v36 & 10.2 & 1 - 50 @xmath187 & 60,000 & 1 & 2 + rsg - loz & o25 & 0 . &
1.2 @xmath188 & 10@xmath189 & ... & ... + rsg - bin & s25 & 2.8 & 0.1 - 9 @xmath190 & @xmath191 & 0.01 & 0.02 + pair sn & z250a & 19.6 & 0.9 - 8 @xmath190 & @xmath192 & 0.7 & 1 + & z250b & 20.6 & 0.4 - 6 @xmath190 & 10@xmath193 & 2 & 4 + & z250c & 12.1 & 0.2 - 2 @xmath190 & @xmath194 & 4 & 8 + wr - bin & 8a & 0.18 & 2.5 - 5 @xmath195 & 100 & 20 & 40 + & 8b & 1.19 & 1 - 5 @xmath195 & 100 & 100 & 200 + & 16a & 0.01 & 3 - 4 @xmath195 & 100 & 1 & 2 + & 16b & 0.83 & 2 - 4 @xmath195 & 40 & 200 & 300 + [ tab : models ]    ccccccccc bsg & v24 & 5.4 & 23.96 & 14.0 - 24 & 2.1 & 1 & 8.32 & 0.13 + & v36 & 3.7 & 35.85 & 25.7 - 35.9 & 1.2 & 1 & 14.97 & 0.11 + rsg - loz & o25 & 5.1 & 23.76 & ... & ... & 0.52 & 10.07 & 0.12 + rsg - bin & r25 & 4.6 & 12.09 & 9.3 - 12.1 & 0.0025 & 1 & 9.19 & 0.13 + pair sn & 250a & ... & 250 & 230.4 - 250 & 2.9 & 1 & 142 & 0.05 + & 250b & ... & 250 & 229.4 - 250 & 8.2 & 1 & 166 & 0.09 + & 250c & ... & 250 & 237.9 - 250 & 6.8 & 1 & 222 & 0.13 + wr - bin & 8a & 7.4 & 8.043 & 7.86 - 8.04 & 1.0 & ... & 8.04 & 0.25 + & 8b & 3.2 & 8.043 & 6.85 - 8.04 & 7.8 & ... & 8.04 & 1 + & 16a & 2.5 & 16.087 & 16.08 - 16.09 & 0.08 & ... & 16.08 & 0.18 + & 16b & 1.5 & 16.087 & 15.26 - 16.09 & 12 & ... & 16.08 & 0.48 + [ tab : models1 ]


detecting planetary transits is highly valuable not only because it makes an independent confirmation of a planet from radial - velocity observations but also because the photometric transits provide unambiguous information on planet masses and radii , and thus mean density and interior structure ( e.g. * ? ? ?
* and references therein ) .
various high - precision follow - up studies for transiting planets can also uncover planetary atmospheres and the ( mis-)alignment between the stellar spin and planetary orbital axes via the rossiter - mclaughlin effect ( e.g. * ? ? ?
* ; * ? ? ?
these properties of transiting planets provide valuable hints for planet formation and evolution processes including dynamical interaction between planets .
although transiting planets are as such important , it is difficult to detect them around evolved stars , especially giant stars , due to the large sizes of the host stars
. targets for on - going radial - velocity surveys for planets around giants have typical radii of @xmath110 @xmath6 ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , and then the relative flux variation of such a giant host star caused by a transit of a jupiter - sized planet is only @xmath7 , which is comparable to that for transit of an earth - sized planet across a solar - type star .
it is impossible to detect such a transit from ground - based photometry and thus no transiting planets had been found around giant stars .
our understanding of properties of planets around such stars are thus far behind from those around solar - like stars , although tens of planets have been found around giants by precise radial - velocity surveys ( e.g. * ? ? ?
* and references therein ) .
the _ kepler _ mission was successful in detecting transiting planets with very high photometric precision ( @xmath8 ) since 2009 .
to date it has discovered more than 4000 transiting planet candidates from sub - earth - size to super - jupiter - size in @xmath10.32000 d orbits ( e.g. * ? ? ?
although most of them are orbiting solar - type stars , _ kepler _ s photometric precision is high enough to detect a jupiter - sized planet transiting a giant star .
as expected , _ kepler _ has identified several planet candidates around giant stars with radii larger than @xmath9 @xcite .
the planet candidates have radii comparable to or larger than jupiter s , and it is interesting that many of them have short - period orbits .
since such short - period planets have rarely been found by radial - velocity surveys around evolved stars ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , they will provide us unique opportunities to investigate planet formation and evolution processes around such stars .    here
we report the results of radial - velocity follow - up observations using subaru 8.2 m telescope for two of the _ kepler _ evolved stars , kepler-91 ( koi-2133 ) and koi-1894 .
thanks to the high precision in their radial - velocity measurements , we independently confirmed a jovian planet ( kepler-91b , koi-2133.01 ) previously reported around kepler-91 @xcite and newly found a hint for the existence of additional companion in the system .
we also detected a possible sub - saturn - mass planet around koi-1894 ( koi-1894.01 ) with 23@xmath4 confidence level .
the rest of the paper is organized as follows .
the adopted stellar parameters for the two stars are presented in section [ targets ] .
the observations are described in section [ obs ] and the results of global analysis for the radial velocities and light curves are presented in section [ analysis ] .
section [ summary ] is devoted to discussion and summary .
kepler-91 ( kic 8219268 , koi-2133 ; @xmath10=12.495 ) was identified as a candidate star hosting a transiting short - period jupiter - sized planet ( koi-2133.01 ; @xmath11 , @xmath12 d ) by @xcite .
after that , @xcite reported the confirmation of its planetary nature based on the detailed analysis of orbital brightness modulation seen in the light curve caused by ellipsoidal variation , doppler boosting , and reflected / emitted light from planet ( e.g. * ? ? ?
* ; * ? ? ?
* ) , and then the planet was named kepler-91b with the radius @xmath13 and mass @xmath14 . although @xcite and @xcite claimed a possible non - planertary nature for the system , @xcite and @xcite very recently reconfirmed the planetary nature via radial - velocity measurements with a precision of @xmath1100 m s@xmath2 and @xmath120 m s@xmath2 , respectively , and obtained the planetary mass of 1.09@xmath150.20@xmath16 and 0.73@xmath150.13@xmath16 , respectively .    the stellar parameters ( effective temperature @xmath17 , surface gravity @xmath18 , mass @xmath19 , and radius @xmath20 ) of kepler-91 were reported by @xcite to be @xmath21 k , @xmath22 cgs , @xmath23 , and @xmath24 .
after that , the parameters have been updated by the spectroscopic and asteroseismic analyses in @xcite and @xcite , which are consistent with each other .
thus , we here adopted the values listed in @xcite : @xmath25 k , @xmath26 cgs , @xmath27 , and @xmath28 .
koi-1894 ( kic 11673802 ; @xmath10=13.427 ) was also reported to be a planet - host candidate having a transiting short - period jupiter - sized planet ( koi-1894.01 ; @xmath29 , @xmath30 d ) by @xcite , though radial - velocity follow - up observations for the star have not been reported yet .
the stellar parameters of koi-1894 were derived by @xcite to be @xmath31 k , @xmath32 cgs , @xmath33 , and @xmath34 , and they have been updated by the spectroscopic and asteroseismic analyses in @xcite as @xmath35 k , @xmath36 , and @xmath37 .
here we adopted the updated values in this paper .
recently @xcite reported a non - detection of blended stars for koi-1894 , which could have been physically associated companions and/or responsible for transit false positives if they were within @xmath10.15@xmath382.5@xmath38 separation and with magnitude difference up to @xmath39 , by high - angular - resolution ao imaging .
the parameters for the two stars we adopt here are summarized in table [ tbl : stars ] .
lrr [ tbl : stars ] @xmath17 ( k ) & 4550@xmath1575 & 4992@xmath1575 + @xmath18 ( cgs ) & 2.953@xmath150.007 & 2.87 + @xmath19 ( @xmath40 ) & 1.31@xmath150.10 & 1.410@xmath150.214 + @xmath20 ( @xmath41 ) & 6.30@xmath150.16 & 3.790@xmath150.190
we obtained high - precision radial - velocity data for the stars with the 8.2 m subaru telescope and the high dispersion spectrograph ( hds ; * ? ? ?
* ) in 2013 and 2014 .
we used the setups of stdi2b ( 2013 june 28 - 30 and 2014 july 8 , 10 , 14 - 16 ) and stdi2a ( 2013 december 6 and 2014 july 12 - 13 ) , which simultaneously cover a wavelength region of 35006200@xmath42 and 49007600@xmath42 respectively , the image slicer no.2 ( is#2 ; * ? ? ?
* ) yielding a spectral resolution ( @xmath43 ) of 80000 , and an iodine absorption cell ( i@xmath44 cell ; * ? ? ?
* ) for precise radial - velocity measurements .
we obtained a total of 29 and 18 data points for kepler-91 and koi-1894 with typical signal - to - noise ratio of s / n@xmath453080 pix@xmath2 and 3045 pix@xmath2 , respectively , by an exposure time of 1200 sec depending on weather condition .
the reduction of echelle data ( i.e. bias subtraction , flat - fielding , scattered - light subtraction , and spectrum extraction ) was performed using the iraf software package in the standard manner .
we performed radial - velocity analysis for i@xmath44-superposed stellar spectra ( star+i@xmath44 ) by the method described in @xcite and @xcite , which is based on the method by @xcite and @xcite .
a star+i@xmath44 spectrum is modeled as a product of a high resolution i@xmath44 and a stellar template spectrum convolved with a modeled instrumental profile ( ip ) of the spectrograph .
we obtain the stellar template spectrum by deconvolving a pure stellar spectrum with the ip estimated from an i@xmath44-superposed flat spectrum .
we achieved a radial - velocity precision of 517 m s@xmath2 for kepler-91 and 915 m s@xmath2 for koi-1894 .
the derived radial velocities are listed in table [ tbl : rv - koi2133 ] and table [ tbl : rv - koi1894 ] together with the estimated uncertainties , and are plotted in figure [ fig : rv_koi2133_notrend ] and figure [ fig : rv_koi1894 ] .
lrr [ tbl : rv - koi2133 ] 6473.00247 & 43.20 & 4.61 + 6473.01719 & 33.85 & 4.58 + 6473.03192 & 45.26 & 5.03 + 6473.96075 & 17.35 & 7.36 + 6473.97549 & 8.93 & 7.34 + 6473.99022 & 30.04 & 6.39 + 6474.95764 & @xmath4624.28 & 5.43 + 6474.97235 & @xmath4622.62 & 5.03 + 6474.98708 & @xmath4639.46 & 5.50 + 6633.69753 & @xmath4642.12 & 17.38 + 6847.93547 & 42.30 & 10.56 + 6847.95020 & 37.46 & 13.63 + 6847.96513 & 11.49 & 8.76 + 6847.98006 & 29.77 & 7.96 + 6847.99479 & 17.07 & 6.12 + 6848.00952 & 19.94 & 8.61 + 6848.08121 & 23.25 & 6.59 + 6849.75412 & @xmath4652.12 & 5.23 + 6849.76886 & @xmath4636.10 & 5.06 + 6849.78360 & @xmath4634.47 & 5.48 + 6853.04854 & @xmath4640.23 & 13.48 + 6853.97273 & 5.73 & 5.72 + 6853.98746 & 9.02 & 4.73 + 6854.00219 & 7.66 & 5.17 + 6854.84234 & 2.14 & 4.99 + 6854.85708 & 0.36 & 5.25 + 6854.87181 & @xmath466.69 & 6.03 + 6855.80523 & @xmath4643.37 & 6.63 + 6855.81828 & @xmath4641.91 & 5.28    lrr [ tbl : rv - koi1894 ] 6473.04710 & @xmath469.95 & 9.45 + 6473.06182 & @xmath460.82 & 11.01 + 6473.07656 & @xmath4611.81 & 9.94 + 6474.00523 & @xmath465.31 & 10.97 + 6474.01995 & 1.95 & 10.25 + 6474.03468 & 17.33 & 10.36 + 6474.04941 & 5.73 & 10.49 + 6475.00210 & @xmath463.24 & 10.03 + 6475.01683 & @xmath4621.07 & 9.68 + 6475.03159 & @xmath4613.24 & 8.79 + 6475.04632 & @xmath466.92 & 9.62 + 6848.09625 & @xmath464.14 & 15.02 + 6849.94269 & 0.52 & 10.34 + 6849.95743 & 13.69 & 9.86 + 6849.97216 & 10.09 & 9.04 + 6855.88939 & 0.65 & 9.63 + 6855.90411 & @xmath463.06 & 8.95 + 6855.91885 & 6.48 & 9.09
the variations in the observed radial velocities for both targets show a sign of planetary companions to kepler-91 and koi-1894 . in order to obtain accurate and precise estimates for system parameters of those systems ,
we here present a global analysis that makes use of all the available information from the _ kepler _ photometry and our spectroscopy .
as is well known , a very precise light curve of a star orbited by planet(s ) shows a periodic modulation due to several astrophysical effects : the ellipsoidal variation , doppler boosting , and reflected / emission light from the planet ( e.g. , * ? ? ?
* ; * ? ? ?
the last two effects ( doppler boosting and planetary reflection / emission ) are synchronous with the planet s orbital period @xmath47 , with different peak locations along the planet phase ( @xmath48 for boosting and @xmath49 for the planetary reflection / emission , respectively ) . on the other hand , the ellipsoidal variation , which is caused by a tidal distortion by the planet s gravity , have two flux peaks ( @xmath48 and 0.75 ) within one orbit .
thanks to the different phase dependence of these three effects , a very precise light curve and its modeling enable us to distinguish these three effects , and we can extract physical parameters such as the planet - to - star mass ratio @xmath50 and scaled semi - major axis @xmath51 , where @xmath52 is semi - major axis and @xmath20 is stellar radius . moreover , for transiting systems as in the present cases , incorporating the transit and/or secondary eclipse model into the above phase - curve variation lets us learn more about the planet properties ( e.g. , planet - to - star radius ratio @xmath53 and orbital inclination @xmath54 ) .    to obtain the phase - folded light curves for kepler-91 and koi-1894 , and estimate system parameters , we downloaded all available public light curves ( q0q17 ) from _ kepler _ mast archive .
while only long - cadence data were available for kepler-91 , koi-1894 s light curves involve some short - cadence data . adopting the pdc - sap flux data , for which unphysical artifacts are detrended , we reduced the light curves by the following procedure .
first , after removing planetary transits , we further detrended and normalized the light curve for each quarter by fitting it with a fifth - order polynomial so as to remove the long - term trends that were not removed in the pdc - sap flux .
this process was repeated implementing a @xmath55 clipping to remove outliers .
we then combined the light curves for all quarters and phase - folded them with the ephemerides derived by the official _
kepler _ team using the q1q16 data ; orbital period @xmath56 d and transit center @xmath57 bjd for kepler-91 , and @xmath58 d and @xmath59 bjd for koi-1894 ( from nasa exoplanet archive ) .
we later consider the impact of incorrect ephemeris .
finally , the folded light curves were binned into 300 and 250 phase bins for kepler-91 and koi-1894 , respectively .
these bin numbers were adopted so that each bin approximately covers a cycle span of long - cadence data ( @xmath60 minutes ) . the flux error for each bin @xmath61 , which is the standard deviation of the mean flux ,
was computed based on the dispersion of the flux values within the bin .
the long - cadence and short - cadence flux data were separately folded and binned for koi-1894 , but the binned light curve of the short - cadence data was much noisier than that for long - cadence .
we thus decided to ignore the short - cadence data in the following analysis .
the analysis below is based on the method by @xcite , who employed the evil - mc model @xcite for the phase - curve variation ( i.e. , ellipsoidal variations , doppler boosting , and planetary reflection ) with some revisions ( e.g. , application to an eccentric orbit ) .
this phase - curve model is multiplied by the analytic transit model by @xcite , and the relevant parameters ( e.g. , @xmath51 and @xmath54 ) are simultaneously determined .
in addition , we also model and fit the observed radial velocities @xmath62 .
thus , the @xmath63 statistics in the present case is @xmath64 where @xmath65 and @xmath66 are the @xmath67-th observed radial velocity and its error , and @xmath68 and @xmath69 are the @xmath67-th observed light curve flux and its error , which are the mean flux and its standard deviation in the @xmath67-th bin , respectively ( e.g. , * ? ? ? * ) .
the radial velocity is modeled as @xmath70 where @xmath71 , @xmath72 , @xmath73 , @xmath74 , and @xmath75 are the radial velocity semi - amplitude , true anomaly , orbital eccentricity , argument of periastron , and radial velocity offset of our dataset , respectively . for the light curve model @xmath76
, we refer the readers to @xcite , for details . in the next subsections ,
we describe the fitting procedure for each of the targets .
the binned light curve of kepler-91 in figure [ fig : phase_koi2133 ] shows a clear pattern of phase - curve variation ; the double peaks at @xmath48 and 0.75 are representative of the ellipsoidal variation .
in addition , figure [ fig : phase_koi2133 ] suggests a possible detection of secondary eclipse at @xmath77 , which means the reflected / emitted light from the planet is visible in the folded light curve .
therefore , following @xcite and @xcite , we model the planet s light @xmath78 by the following expression : @xmath79 where @xmath80 is the flux variation amplitude of planet s reflected / emitted light and @xmath81 is the planetary flux offset arising from the homogeneous surface emission .
following @xcite , @xcite , and @xcite , we here introduce the  phase - delay " @xmath82 of the brightest part on the planet surface from the substellar point .
the planetary light @xmath78 is added to the flux model for the beaming and ellipsoidal variation , integrated over the stellar visible hemisphere and wavelength through the _ kepler _ band . in integrating the local flux calculated from the evil - mc model
, we assume that the stellar effective temperature is 4550 k , and adopt the gravity darkening exponent of @xmath83 for equation ( 10 ) of @xcite based on the theoretical calculation by @xcite .
the free parameters relevant to the light curve model are @xmath51 , the transit impact parameter @xmath84 , limb - darkening coefficients @xmath85 and @xmath86 for the quadratic limb - darkening law , @xmath50 ( @xmath87 ) , @xmath81 , @xmath80 , @xmath82 , the overall normalization factor @xmath88 for the folded light curve , @xmath53 , @xmath89 , @xmath90 , and @xmath91 , which represents the small time deviation of the transit center from the ephemeris reported by the _
kepler _ team . in our updated ephemeris ,
the initial transit center becomes @xmath92 .
we fix the orbital period to be @xmath93 d that are derived by the _ kepler _ team based on the q1q16 data ( see section [ method ] ) .    assuming that the likelihood is proportional to @xmath94 in equation ( [ eq : chisq ] ) , we simultaneously model the observed radial velocity and light curve , and compute the posterior distribution for each fitting parameter by implementing markov chain monte carlo ( mcmc ) simulation .
in addition to the above mentioned twelve parameters , we add @xmath75 in equation ( [ eq : rvmodel ] ) to the fitting parameters . note that the radial velocity semi - amplitude @xmath71 in equation ( [ eq : rvmodel ] ) is related to the planet - to - star mass ratio @xmath50 by @xmath95 in computing the posteriors , we do not impose priors on the fitting parameters except for @xmath96 and @xmath97 ; due to the sparse sampling of the long - cadence data ( @xmath60 minutes ) and quality of the binned light curve , the limb - darkening coefficients are poorly constrained in the absence of priors , and thus we decide to put gaussian priors on the limb - darkening coefficients based on the theoretical table by @xcite as @xmath98 and @xmath99 . in our mcmc algorithm , originally developed in @xcite , the step size of each fitting parameter is iteratively scaled so that the overall acceptance ratio falls between 15% and 35% . after running 1,000,000 chains , the best - fit value and uncertainty for each fitting parameter are estimated from the median , and 15.87 and 84.13 percentiles of the marginalized posterior distribution of that parameter .
shortly after we performed the first mcmc trial , we noticed a possible trend or drift in the observed radial velocities .
figure [ fig : rv_koi2133_notrend ] plots the best - fit radial velocity curve as a function of the orbital phase @xmath100 . for clarity ,
the data taken in 2013 are shown in green and those by the 2014 campaign are plotted in red . while we have detected a clear modulation with an amplitude of @xmath101 m s@xmath2 , radial velocities take lower values for the 2014 data , which is evident in the bottom panel indicating the residual of the observed radial velocities from the best - fit model .    in order to confirm the presence of the radial velocity drift , we try to fit the data with an additional parameter : the radial velocity drift @xmath102 . adding the drift term @xmath103 to the right side of equation ( [ eq : rvmodel ] ) , we performed again the mcmc simulation to fit both radial velocity and light curve .
the derived best - fit parameters and their uncertainties are summarized in table [ tbl : best - fit - kepler91 ] .
to compare between the radial velocity models with and without a drift term , we compute bayesian information criteria ( bic ) for the two cases , which is computed as @xmath104 , where @xmath105 is the number of free parameters and @xmath106 is the number of data points . in the absence of a trend , we obtain @xmath107 , and for the case of the radial velocity model with @xmath102 , bic becomes 383 ; @xmath108 is much larger than 10 , meaning that the model with @xmath102 is strongly favored .
therefore , we conclude that a radial velocity drift ( or trend ) is present in our dataset , and report the best - fit parameters for the case with @xmath102 as the final result .
we also tried several periods around @xmath47 and found that the @xmath109 gave slightly better results in the global fitting than @xmath47 did ( @xmath110 ) .
however , the resultant parameters for the two cases are well consistent within 0.2@xmath4 level .
differs by @xmath111 level . ]    figures [ fig : phase_koi2133 ] and [ fig : rv_koi2133 ] plot the folded light curve and radial velocities ( red points ) along with their best - fit models ( blue line ) .
the constant radial velocity offset and trend are both removed from the radial velocity data in figure [ fig : rv_koi2133 ] .
the bottom panel in each figure shows the residuals from the best - fit model . from the posterior distribution of the fitting parameters
, we also estimate the orbital and planetary parameters ( e.g. , the orbital inclination @xmath54 , planet mass @xmath112 and radius @xmath113 ) assuming the stellar properties reported by asteroseismology ( table [ tbl : stars ] ) to be @xmath114 deg , @xmath115 , and @xmath116 .
these estimates are also shown in table [ tbl : best - fit - kepler91 ] .
lrr [ tbl : best - fit - kepler91 ] @xmath51 & 2.253@xmath117 & 2.238@xmath118 + @xmath84 & 0.9041@xmath119 & 0.9062@xmath120 + @xmath85 & 0.737@xmath121 & 0.7370@xmath122 + @xmath86 & 0.493@xmath123 & 0.493@xmath124 + @xmath125 ( 10@xmath126 ) & 4.82@xmath127 & 4.72@xmath128 + @xmath81 ( 10@xmath129 ) & 5.65@xmath130 & 5.63@xmath131 + @xmath80 ( 10@xmath129 ) & 1.41@xmath132 & 1.42@xmath133 + @xmath82 & 0.394@xmath134 & 0.395@xmath135 + @xmath53 & 0.02286@xmath136 & 0.02294@xmath137 + @xmath89 & 0.0280@xmath138 & 0.0277@xmath139 + @xmath90 & @xmath460.043@xmath140 & @xmath460.046@xmath140 + @xmath91 ( 10@xmath141 day ) & @xmath461.09@xmath142 & @xmath461.11@xmath142 + @xmath75 ( m s@xmath2 ) & @xmath4642.5@xmath143 & @xmath4633.2@xmath144 + @xmath102 ( m s@xmath2 day@xmath2 ) & @xmath460.0612@xmath145 &  + @xmath63 & 296 & 396 + bic & 383 & 477 + @xmath54 ( deg ) & 67.37@xmath146 & 67.20@xmath147 + @xmath73 & 0.0519@xmath148 & 0.0535@xmath149 + @xmath74 ( deg ) & @xmath4657.2@xmath150 & @xmath4658.8@xmath151 + @xmath112 ( @xmath16 ) & 0.66@xmath150.06 & 0.65@xmath150.06 + @xmath113 ( @xmath152 ) & 1.40@xmath150.04 & 1.41@xmath150.04      the radial - velocity variations are less visible for koi-1894 owing to the large radial - velocity error compared to the small semi - amplitude .
the phase - folded light curve also shows a very tiny variation , if any , along the orbital phase . to extract possible planetary signals , we simultaneously model the radial velocities and light curve as in the case of kepler-91 . since the observed transit depth and koi-1894 s stellar radius suggest that the radius of the transiting companion ( koi-1894.01 ) is no more than @xmath153 , the reflected / emitted light from the planet
is expected to be very small . according to equations ( 1 )  ( 3 ) in @xcite , the planetary reflection is estimated to be @xmath154 ppm at the location of koi-1894.01 , which is smaller than the expected amplitude of the ellipsoidal variation ( @xmath155 ppm ) .
visual inspection of the binned light curve also suggests the absence of secondary eclipse .
thus for simplicity , we here neglect the planet flux @xmath78 for koi-1894 , and model the folded light curve with the ellipsoidal variation ( including doppler boosting ) and transit only . neglecting @xmath78 also helps to avoid @xmath63 being stuck at a local minimum during the optimization .
fixing the system parameters as @xmath156 k , @xmath157 , the gravity darkening exponent of @xmath158 @xcite , and @xmath159 derived with the q1q16 data ( see section [ method ] ) , we perform mcmc simulations for koi-1894 to infer the posterior distributions for the fitting parameters . in the present case , we have the eleven free parameters : @xmath51 , @xmath84 , @xmath85 , @xmath86 , @xmath50 , @xmath88 , @xmath53 , @xmath89 , @xmath90 , @xmath91 , and @xmath75 .
again , we assume gaussian priors for the limb - darkening parameters as @xmath160 , and @xmath161 from the table by @xcite . due to the weak radial - velocity signal of koi-1894.01 , we found that the global fit to the radial velocities and light curve exhibits a degeneracy between the fitting parameters ( e.g. , @xmath162 and @xmath50 ) .
thus , for koi-1894 , we also impose an additional prior on the host star s density from table [ tbl : stars ] ( @xmath163 ) , assuming a gaussian distribution @xcite .
otherwise , the fitting algorithm is exactly the same as for kepler-91 .
the result of the fit is summarized in table [ tbl : best - fit - koi1894 ] .
the best - fit value for the mass ratio is @xmath164 , implying @xmath165 detection of the planet .
the mass and radius ratios are translated as @xmath166 and @xmath167 assuming the stellar mass and radius in table [ tbl : stars ] .
table [ tbl : best - fit - koi1894 ] also shows the result of our fit in the absence of the gaussian prior on @xmath168 .
as expected , the mass constraint is slightly weaker for this case .
the orbital eccentricity is poorly constrained due to the small planetary signal and small number of radial velocity data , but we can rule out a high eccentricity ( @xmath169 ) with our current datasets .
we also try to fit the data with a circular orbit ( @xmath170 ) , and compare between the eccentric and circular cases .
as shown in table [ tbl : best - fit - koi1894 ] , the circular orbital fit yields a slightly better constraint on the mass ratio ; @xmath171 with the stellar - density prior , leading to a @xmath172 detection .
the mass and radius ratios are translated as @xmath173 and @xmath174 assuming the stellar mass and radius .
based on the @xmath63 values for the best - fit parameters , we compare bic values for the two cases ; bic values of @xmath175 and @xmath176 are obtained for the eccentric and circular cases ( with the stellar - density prior ) , respectively , so that fitting with a circular orbit is favored .
the observed radial velocities and whole light curve are plotted in figures [ fig : rv_koi1894 ] and [ fig : phase_koi1894 ] , with their best - fit models for @xmath170 .
we also show zoomed - in versions of the transit and phase - curve variation in figures [ fig : transit_koi1894 ] and [ fig : ev_koi1894 ] , respectively . in figure
[ fig : ev_koi1894 ] , the binned flux data used for the fit ( 250 bins ) are plotted by the red crosses , and black points with errorbars indicate the flux data binned into 30 bins .
we also tested several periods around @xmath47 and found that about @xmath177 gave the minimum @xmath63 value in the global fitting ( @xmath178 ) , suggesting that the true period might exist around it .
nonetheless , the resultant parameters are well consistent with those for the case of @xmath47 within 0.2@xmath4 level .
differs by @xmath179 level . ]
our global analysis indicates the detection of koi-1894.01 is still marginal , with only @xmath180 level , and the planet mass seems to be mainly constrained by the radial velocity data . in order to see if an independent estimate for the planet mass from the light curve alone gives a comparable result , we fit the folded light curve and estimate system parameters without radial velocity data . as a result
, we obtain @xmath181 , consistent with the above result by the global fit .
we reported the results of high - precision radial - velocity measurements with subaru / hds for the two _ kepler _ evolved stars kepler-91 and koi-1894 , which are candidates to host transiting planets .
based on the simultaneous modeling of the radial - velocity data and _ kepler _ light curves , we independently confirmed the planetary nature of kepler-91b .
the radial velocity semi - amplitude of @xmath182 m s@xmath2 we derived is consistent with what @xcite found ( @xmath183 m s@xmath2 ) with @xmath184 , but incompatible with the result by @xcite who reported @xmath185 m s@xmath2 with @xmath186 , although the reason for this discrepancy is unknown .
we newly detected a drift of @xmath187 m s@xmath2 between the radial velocities taken at @xmath1 1-yr interval thanks to our better measurement precision and longer period of time of observations compared to the previous ones .
the persistent radial velocity drift is suggestive of a possible presence of another companion ( planet ) outside of kepler-91b . due to the lack of data , however , we are not able to pin down the period nor mass of the companion .
intensive high - precision radial - velocity monitoring of the star will uncover the unseen companion and provide hints for formation and evolution of the planetary system .
the estimated parameters for kepler-91b in table [ tbl : best - fit - kepler91 ] are reasonably in good agreement with the previous study by @xcite except for the mass ratio @xmath50 and scaled semi - major axis @xmath51 . this is likely because our radial velocity data show a smaller amplitude and planet mass is estimated to be small . since the flux amplitude due to ellipsoidal variation is approximately proportional to @xmath188 @xcite , a small @xmath50 resulted in the smaller @xmath51 in order to explain the observed flux amplitude .
the westward phase - shift ( @xmath189 ) of the flux maximum in planetary light ( from the substellar point ) is indicative of the inhomogeneous cloud coverage on the planetary surface .
this is consistent with the trend that @xcite found , claiming that westward phase - shifts are preferentially seen for relatively cool close - in planets with equilibrium temperatures of @xmath190 k ; @xcite estimated the equilibrium temperature of kepler-91b to be @xmath191 k , depending on the assumed heat redistribution parameter .
the magnitude of the phase - shift ( @xmath192 ) is , however , unusually large in our best - fit model compared with the previously reported values for other close - in planets @xcite .
the reason is unknown , but the huge coverage of the planet surface illuminated by the host star of the planet surface is always illuminated by the host star due to its huge radius .
] might be responsible .
note that we also detected a dip in the folded light curve ( figure [ fig : phase_koi2133 ] ) around @xmath193 reported in @xcite , and this dip should have more or less affected the fitting result .
lrrrr [ tbl : best - fit - koi1894 ] @xmath51 & 3.79@xmath194 & 3.68@xmath195 & 4.25@xmath196 & 3.24@xmath197 + @xmath84 & 0.634@xmath198 & 0.56@xmath199 & 0.48@xmath200 & 0.68@xmath201 + @xmath85 & 0.724@xmath124 & 0.723@xmath202 & 0.724@xmath123 & 0.7271@xmath203 + @xmath86 & 0.391@xmath204 & 0.3908@xmath205 & 0.392@xmath206 & 0.391@xmath207 + @xmath125 ( 10@xmath126 ) & 1.19@xmath208 & 1.02@xmath209 & 1.29@xmath210 & 0.70@xmath211 + @xmath81 ( 10@xmath129 ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath80 ( 10@xmath129 ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath82 & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath53 & 0.01777@xmath212 & 0.01739@xmath213 & 0.01717@xmath214 & 0.01790@xmath215 + @xmath89 & 0 ( fixed ) & 0.052@xmath216 & 0 ( fixed ) & 0.047@xmath217 + @xmath90 & 0 ( fixed ) & 0.085@xmath218 & 0 ( fixed ) & 0.10@xmath219 + @xmath91 ( 10@xmath141 day ) & @xmath460.25@xmath220 & @xmath460.53@xmath221 & @xmath460.24@xmath222 & @xmath460.57@xmath223 + @xmath75 ( m s@xmath2 ) & @xmath4610.6@xmath224 & @xmath4610.0@xmath225 & @xmath4611.3@xmath226 & @xmath466.8@xmath227 + @xmath63 & 234 & 233 & 235 & 233 + bic & 285 & 295 & 286 & 294 + @xmath54 ( deg ) & 80.4@xmath228 & 80.2@xmath229 & 83.5@xmath230 & 76.2@xmath231 + @xmath73 & 0 ( fixed ) & 0.149@xmath232 & 0 ( fixed ) & 0.186@xmath233 + @xmath74 ( deg ) & 0 ( fixed ) & 61@xmath234 & 0 ( fixed ) & 65@xmath235 + @xmath112 ( @xmath16 ) & 0.18@xmath150.07 & 0.15@xmath236 & 0.19@xmath150.08 & 0.10@xmath237 + @xmath113 ( @xmath152 ) & 0.66@xmath150.04 & 0.64@xmath238 & 0.63@xmath239 & 0.66@xmath240    as for koi-1894 , we did not detect any statistically significant radial - velocity variations with our measurement precision of  915 m s@xmath2 and the number of data points .
we excluded the possibility of a grazing transit by a binary companion for the star and set an upper limit on the mass of koi-1894.01 to be @xmath241 by our non - detection of radial - velocity variations with 3@xmath4 level .
however , we detected possible radial - velocity variations with a semi - amplitude of @xmath115 m s@xmath2 in phase with ellipsoidal variations of the star with 23@xmath4 level .
although we can not say the detection is statistically significant at this stage , it suggests that the koi-1894.01 could be a sub - saturn - mass planet .
figure [ fig : plm - strad ] shows distribution of mass of exoplanets currently known plotted against their host star s radius . as seen in the figure , koi-1894.01 could be one of the lowest mass planets ever discovered around evolved stars together with kepler-56b , a super - neptune - mass planet ( @xmath242 ) detected via ttv ( transit timing variation ) method @xcite .
actually the stellar parameters for koi-1891 are similar to those of kepler-56 ( * ? ? ?
* @xmath243 , @xmath244 ) .
confirmation of koi-1894.01 is highly encouraged in order to uncover such a new population of sub - saturn and super - neptune planets around relatively massive evolved stars , which have rarely been found so far either by radial - velocity surveys or transit ones .
this research is based on data collected at subaru telescope , which is operated by national astronomical observatory of japan ( naoj ) .
we are grateful to all the staff members of subaru for their support during the observations .
this research has made use of the simbad database , operated at cds , strasbourg , france .
is supported by japan society for promotion of science ( jsps ) fellowship for research ( pd:25 - 3183 ) .
is supported by jsps research fellowships for young scientists ( no .
26 - 7182 ) and by the leading graduate course for frontiers of mathematical sciences and physics .
barclay , t. et al .
2014 , , submitted ( arxiv:1408.3149 ) batalha , n.m .
, et al . 2013 , , 204 , 24 butler , r. p. , marcy , g. w. , williams , e. , mccarthy , c. , dosanjh , p. , & vogt , s. s. 1996 , , 108 , 500 claret , a. 1998 , a&as , 131 , 395 claret , a. , & bloemen , s. 2011 , , 529 , a75 da silva , l. , et al .
2006 , , 458 , 609 demory , b .- o .
, et al . 2013 , , 776 , l25 esteves , l.j . , de mooij , e.j.w . , & jayawardhana , r. 2013 , , 772 , 51 esteves , l. j. , de mooij , e. j. w. , & jayawardhana , r. 2014 , , submitted ( arxiv:1407.2245 ) faigler , s. , & mazeh , t. 2011 , mnras , 415 , 3921 faigler , s. , & mazeh , t. 2014 , , submitted ( arxiv:1407.2361 ) hirano , t. , narita , n. , shporer , a. , sato , b. , aoki , w. , & tamura , m. 2011 , , 63 , 531 hirano , t. , et al .
2012 , , 759 , l36 hirano , t. , et al .
2015 , , 799 , 9 huber , d. et al . 2013a , , 767 ,
127 huber , d. et al .
2013b , science , 342 , 331 jackson , b. k. , lewis , n. k. , barnes , j. w. , drake deming , l. , showman , a. p. , & fortney , j. j. 2012 , , 751 , 112 johnson , j.a . , et al .
2007 , , 665 , 785 jones , m.i . ,
jenkins , j.s . , bluhm , p. , rojo , p. , & melo , c.h.f .
2014 , , 566 , 113 kambe , e. , et al .
2002 , , 54 , 865 law , n.m .
2014 , , 791 , 35 lillo - box , j. et al . 2014a , , 562 , 109 lillo - box , j. et al .
2014b , , 568 , 1 liu , y. , sato , b. , takeda , y. , ando , h. , & zhao , g. 2010 , , 62 , 1071 mazeh , t. , nachmani , g. , sokol , g. , faigler , s. , & zucker , s. 2012 , , 541 , a56 noguchi , k. et al .
2002 , , 54 , 855 ohta , y. , taruya , a. , & suto , y. 2009 , , 690 , 1 sato , b. , kambe , e. , takeda , y. , izumiura , h. , & ando , h. 2002 , , 54 , 873 sato , b. , et al .
2008 , , 60 , 539 sato , b. , et al .
2012 , , 64 , 97 sato , b. , et al .
2013 , , 65 , 85 seager , s. & malln - ornelas , g. 2003 , , 585 , 1038 seager , s. & deming , d. 2010 , , 48 , 631 shporer , a. , et al .
2011 , , 142 , 195 sliski , d. h. & kipping , d.m . 2014 , , 788 , 148 tajitsu , a. , aoki , w , & yamamuro , t , 2012 , , 64 , 77 takeda , y. , sato , b. , & murata , d. , 2008 , , 60 , 781 valenti , j.  a. , butler , r.  p. & marcy , g.  w. 1995 , , 107 , 966 .
wang , l. , liu , y. , zhao , g. , & sato , b. 2011 , , 63 , 1035 winn , j. n. 2011 , in exoplanets , ed .
s. seager ( tucson , az : university of arizona press ) , p.5577 zielinski , p. , niedzielski , a. , wolszczan , a. , admow , m. , & nowak , g. 2012 , , 547 , 91