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: [ sec : s1 ]    the compactness of its core is an important structural characteristic of a presupernova star that affects whether it will explode as a supernova . a shallow density gradient around the iron core , as typically exists in more massive stars , implies a higher accretion rate and  ram pressure  @xcite surrounding the iron core during its collapse that must be overcome by neutrino energy deposition or other energy deposition to turn what is initially an implosion into an explosion . the connection between this structure and the likelihood of explosion has been noted many times @xcite , but little attention has been given to explaining just why the compactness has the values that it does , or why different groups obtain different core structures for models with similar main sequence mass . recently , @xcite quantified this compactness in terms of a parameter , @xmath1    where m=2.5 @xmath0  was chosen as the relevant mass for quantifying the density gradient outside of the iron core , which itself typically has a mass in the range 1.4 to 2.0 @xmath0 . in oconnor et al s definition , the time , t@xmath2 , chosen for evaluating @xmath3 is when the core has collapsed to its maximum ( super - nuclear ) density . as we shall show (   [ sec : s2.3 ] ) , for m = 2.5 @xmath0 , no substantial accuracy is lost if this parameter is evaluated for the `` presupernova model '' defined by when the collapse speed first reaches 1000 km s@xmath4 . the choice of m = 2.5 @xmath0  is justified as being larger than the iron core itself , yet deep enough in the star to sample matter that might accrete , especially in a failed explosion . since @xmath5  is inversely proportional to the radius of a fiducial mass outside the iron core , it is small when the density around the iron core falls off rapidly with radius and greater when the density gradient is shallow . that is , when @xmath5  is large , one need go only a shorter distance to include more mass . the systematics of neutrino - powered supernova explosions were explored as a function of @xmath5  by @xcite and @xcite , both of whom found that a small value of @xmath5  favored explosion . @xcite plotted @xmath5  as a function of main sequence mass for a variety of surveys due to @xcite , @xcite and @xcite , and a subset of that data is shown in fig . [ fig : f1 ] . note especially the distinctly non - monotonic behavior in the vicinity of 25 @xmath0  and 40 @xmath0 .    for a polytrope of constant index and composition supported by a constant fraction ( @xmath6 ) of ideal gas pressure , the compactness defined by equation  [ eq : compact ] will actually _ decrease _ monotonically with increasing mass when evaluated at a constant central temperature and small radius . this is because , for polytropes , the central temperature and density obey a relation ( e.g. * ? ? ? * )    @xmath7    for a given @xmath8 the central density thus declines as @xmath9 and the radius required to enclose a given mass increases with m. as will be discussed in  [ sec : s3.1 ] , equation  [ eq : poly ] works reasonably well for massive stars that are burning hydrogen and helium in their centers . the inner regions of such stars are non - degenerate and can be well represented by polytropes of constant index .    as massive stars evolve , however , they develop nested cores with different compositions and entropies . lighter stars develop more degenerate cores in their late stages , especially during oxygen burning . the large gravitational potential at the edge of these degenerate cores implies a small pressure scale height that results in a steep decline in the local density @xcite . as a result , in heavier stars , the higher entropy actually leads to more extended configurations and larger values of @xmath5 . one does not have to to go out as far in the `` envelope '' of the compact core to encompass 2.5 @xmath0 , while in lighter stars that mass is reached farther out . this accounts for a general tendency of the compactness parameter for presupernova stars to increase with mass , especially for low metallicity stars where mass loss is less important . @xmath10 ) models from @xcite ( thick ) and for the solar metallicity models from @xcite ( dot - dashed ) . see also fig . [ fig : f9 ] of @xcite . note the non - monotonic behavior in both sets of models around 25 @xmath0and 40 @xmath0 . since the mass loss was ignored in the low metallicity models , the mass of the co core and the compactness continues to increase for those models at very high masses . [ fig : f1],scaledwidth=35.0% ]    fig . [ fig : f1 ] also shows that the core compactness has significant non - monotonic behavior above about 20 @xmath0 , and possibly some fine structure below 20 @xmath0 . if this behavior is real and robust , it would have interesting implications for the explosion mechanism(s ) for massive stars , the masses of their compact remnants , and stellar nucleosynthesis . to address the `` real and robust '' issue , it is necessary to understand why this non - monotonic structure exists and its sensitivity to uncertainties and variations in the stellar models . the present study has several parts , some dealing with the systematics of @xmath5  found in new surveys of massive stars , others with the uncertainties one should assign to those results . in   [ sec : s2 ] , three new surveys of massive stellar evolution are presented for stars of solar metallicity , very low metallicity ( 10@xmath11 solar ) , and solar metallicity with suppressed mass loss . these new surveys , though using the same physics as in previous works @xcite , are needed to provide a finer mass grid for examining rapid fluctuations in @xmath5  and to give additional data not archived in previous work . in  [ sec : s3 ] , we address the heart of the matter : why is the compactness a non - monotonic function and , in some places , almost chaotic function of main sequence mass . the timing and location of several carbon and oxygen convective shells are found to play a major role . because there are multiple episodes of shell burning , the final compactness can be quite complex . in  [ sec : s4 ] the sensitivity of these results to uncertain assumptions in the physics used in two stellar evolution codes , kepler and mesa , is explored . these uncertainties , especially the treatment of semiconvection and convective overshoot mixing , lead to large variations in the final co core mass that emerges for a given main sequence mass and account for much of the diversity of published results for presupernova evolution . since these uncertainties cloud the interpretation of the compactness plot , a further study is carried out in   [ sec : s5 ] for bare co stars using both codes . provided the carbon mass fractions at carbon ignition are the same as in the full star models , the same trends seen previously for the compactness parameter are also found in these simpler cores , albeit with an offset due to surface boundary pressure in the full star . thus the non - monotonicity of @xmath5  can be considered robust and the co core mass is a somewhat better indicator of pre - supernova structure ( and to some extent explosion dynamics ) than the main sequence mass . finally ,  [ sec : s6 ] offers some conclusions . [ sec : s2 ]      the models published in @xcite were aimed at sampling supernova nucleosynthesis and light curves for a broad range of non - rotating stars of solar metallicity . the resolution in mass necessary for that purpose was not particularly fine , so that study has been repeated using the same code and physics , but with a finer mass grid . the prescriptions for nuclear reaction rates , mass loss , opacity , convective overshoot mixing , and semiconvection are the same as in the previous survey . for the assumed mass loss rates @xcite , the entire hydrogen envelope is lost for stars with initial masses above 33 @xmath0 . heavier stars thus die as wolf - rayet ( wr ) stars with smaller helium cores than would have existed had mass loss been neglected . for the mass loss rate adopted during the wr - phase @xcite , the core shrinks to the extent that the final compactness becomes small again , thus the compactness of very heavy massive stars when they die depends upon an uncertain prescription for mass loss . is evaluated when the collapse speed anywhere in the iron core first reaches 1000 km s@xmath4 ( see fig . [ fig : f7 ] for other choices ) . notice especially the low values of @xmath5  for stars lighter than 22 @xmath0 , the irregular variations from 18 to 22 @xmath0 , and the subsequent rise to a maximum at 23 - 24 @xmath0  followed by a decline . for heavier stars additional structure , including a second rise , is found at higher masses ( fig . [ fig : f3]).[fig : f2],scaledwidth=35.0% ]    the new survey ( fig .  [ fig : f2 ] ) includes non - rotating stars with solar metallicity and masses in the range 15 to 30 @xmath0  calculated in increments of 0.1 @xmath0using the 1d implicit hydrodynamic code kepler ( table 1 ) . altogether , 151 models were simulated using the standard physics ( for variations see   [ sec : s4 ] ) . nuclear energy generation was calculated , up to oxygen depletion ( central oxygen mass fraction less than 0.04 ) , using a 19 isotope network @xcite . a quasi - statistical equilibrium network and a nuclear statistical equilibrium network were used thereafter . however , a much larger `` adaptive '' network of up to 1500 isotopes was carried in parallel with the main structure calculation @xcite and was used to calculate changes in the electron mole number , @xmath12 . this was particularly important during oxygen burning when the quasi - equilibrium network could not be employed , but substantial weak interactions were already changing the central structure . this is the same approach that was used by @xcite , but a substantial improvement over what was done in @xcite . these models are collectively referred to as the `` s series '' stars - for `` solar '' . in addition to confirming the distinct  bump  previously seen around 23 @xmath0  the new study ( fig . [ fig : f2 ] ) shows other interesting features . fine scale variation of the compactness parameter persists throughout most of the mass range , but is particularly apparent between 19 and 21 @xmath0 . these features are further explored in  [ sec : s3 ] .    for this series of models , @xmath5  is small for main sequence stars lighter than 22 @xmath0 . it then rises rapidly to a peak , but declines again to lower values for main sequence masses between 26 @xmath0  and 30 @xmath0 . above about 30 @xmath0 , results for these solar metallicity stars are clouded by the effects of mass loss reducing the helium core mass , but if mass loss is suppressed , the compactness stays large above 40 @xmath0  with a slight dip around 50 @xmath0  ( fig . [ fig : f1 ] and @xmath13 ) . these features pose a challenge to the conventional notion that the difficulty of exploding a star is a monotonically increasing function of its initial mass and instead imply that models below 22 @xmath0  and between 26 @xmath0  and 30 @xmath0  may be easier to blow - up than other masses . this could alter how we think about the galactic chemical evolution @xcite .    in order to facilitate the comparison with low metallicity stars where the mass loss rate may be negligible , another series of solar metallicity models was calculated for stars above 30 @xmath0  in which mass loss was artificially suppressed . these are the  sh models  covering the mass range of 30 - 60 @xmath0  with varying increments - 31,32,33,35,40, .. ,60 @xmath0(table 1 ) . this survey was truncated at 60 @xmath0  because the effect of the pulsational pair instability became noticeable around 70 - 80 @xmath0 . the 60 @xmath0  model of the sh - series has a helium core mass of 27.5 @xmath0  ( fig . [ fig : f4 ] ) , far greater than that of any of the solar metallicity stars calculated with mass loss .    for all runs , the convective time history was recorded and checkpoints registered at representative points along the evolution including : 1 ) helium depletion - when helium reached 1% by mass in the stellar center ; 2 ) carbon ignition - when the central temperature was @xmath14 k ; 3 ) carbon depletion - when the central carbon mass fraction fell below 1% ; 4 ) oxygen ignition - when the central temperature first reached @xmath15 k ; 5 ) oxygen depletion - when the central oxygen abundance declined to 4% ; 6 ) silicon depletion - when silicon mass fraction reached 1% in the center ; and 7 ) presupernova - when any point in the core collapsed faster than 1000 km s@xmath4 . subsequent discussion will refer to these points as representative of the stellar models at these different times as he.dep . , presn . , etc . [ sec : s2.2 ]    stars with equal masses on the main sequence , but differing initial metallicity can have different presupernova structure for a variety of reasons . most importantly , metallicity affects the mass loss . if the amount of mass lost is very low or zero , the presupernova star , including its helium core , is larger , and that has a dramatic effect on the compactness ( fig . [ fig : f1 ] ) . there are also other less dramatic , but important effects . low metallicity implies a smaller initial helium mass fraction ( and more hydrogen ) . the final helium core mass is sensitive to this and is reduced . nitrogen has an abundance just prior to helium ignition that is directly proportional to the initial metallicity and `` nitrogen burning '' by @xmath16n(@xmath17f(@xmath18o is an important exoergic , convective phase in the star s life that precedes helium burning . low metallicity affects the energy generation during hydrogen shell burning by the cno cycle , and this affects the boundary conditions for the helium core . low metallicity also affects the opacity , and the combined effects of opacity , energy generation , and mass loss determine if the star is a red supergiant or a blue one when it dies . this is especially true for the non - rotating stars studied here . the more compact radiative structure of a blue supergiant envelope places greater surface boundary pressure on the helium core therein and can affect its evolution . to illustrate the systematics of compactness in stars of low metallicity , we include here a previously unpublished set of non - rotating presupernova stars by heger and woosley with metallicity of 10@xmath11 @xmath10 . these used the same physics as the s series , but had an initial composition of 76% hydrogen and 24% helium with only a trace , 10@xmath11 solar , of heavier elements . the masses of these stars were 10 to 95 @xmath0 , though only a subset , 15 - 65 @xmath0 , is considered here . this restricted set , called the u series , contains 86 models with varying mass increments ( 0.2@xmath0between 15 - 25@xmath0 , 0.5@xmath0between 25 - 35 @xmath0 , 1@xmath0between 35 - 45@xmath0  and 5@xmath0between 45 - 65 @xmath0 ; ( table 1 ) . other details of these models not related to their compactness will be published elsewhere . are shown for the s - series ( solar metallicity ; thick line ) and u - series ( 10@xmath11 solar metallicity ; dashed ) stars of mass less than 30 @xmath0 . notice a slight shift of @xmath5  near the first `` bump '' . bottom : same plot for the u - series ( dashed ) and sh - series ( continuous ) models for masses over 30 @xmath0 . [ fig : f3],scaledwidth=48.0% ]    . below 35 @xmath0 , the core masses are not very sensitive to metallicity . above 35 @xmath0 , the helium core is somewhat larger for the higher metallicity stars reflecting a more active hydrogen - burning shell during helium burning . [ fig : f4],scaledwidth=48.0% ]    the compactness of these low metallicity models is shown , along with those from s and sh series , in fig . [ fig : f3 ] . without mass loss , a much larger range of helium and co - core masses is accessible ( figure 4 ) . one sees a continuation of the overall increase in compactness parameter with increasing mass all the way up to 65 @xmath0 , but with new features . the compactness rises rapidly above 30 @xmath0  to a broad peak around 40 @xmath0  and then , following a dip at 50 @xmath0 , resumes its rise . all of the stars above 30 @xmath0 , with no mass loss , will probably be very difficult to explode by any solely neutrino - powered mechanism . mass fraction at the time of core carbon ignition is shown for s ( thick upper curve ) , sh ( dashed ) and u ( thinner lower curve ) stars . note that the carbon mass fraction is smaller in stars that ignore mass loss and that the sh series tracks the u series very well . [ fig : f5],scaledwidth=35.0% ]    below 30 @xmath0 , the compactness curve for the u series shows non - monotonic structure with peaks qualitatively similar to those observed for the s series in figure 2 , but with an offset of about a solar mass . the offset is most pronounced in the range 20 - 25 @xmath0  where the compactness starts to rise in the u series for a lower value of main sequence mass . the main cause of this shift is the lower mass fraction of carbon produced by helium burning in the u stars ( fig . [ fig : f5 ] ) . that value , in turn , affects the critical mass where carbon ceases to burn convectively in the star s center (   [ sec : s3.1 ] ) , shifting it to lower values , which causes the compactness curve to rise earlier (   [ sec : s3 ] ) . thus end their lives as blue supergiants . [ fig : f6 ] , scaledwidth=35.0% ]    this variable nucleosynthesis of carbon for the s and u stars results from differing amounts of helium being mixed into the convective helium core as it grows during the final stages of helium burning . in fact , the carbon mass fractions for s and u stars in this mass range are very similar half way through helium burning , but as the helium mass fraction declines below 10% , the u stars are more effective at bringing additional helium in from their outer regions . each helium convected down turns a carbon into an oxygen . the u stars have somewhat higher helium core luminosities at this point and the structure at the outer edge of the helium convection is also influenced by surface boundary pressure from the hydrogen envelope which is different for red and blue supergiants . below 30 @xmath0 , the u - series stars are blue supergiants while the s - series stars are red ones ( fig .  [ fig : f6 ] ) . the pressure at the edge of the helium core does not decline as steeply in a blue supergiant and the entropy barriers inhibiting convection in the outer helium core are slightly reduced . this makes the growth of the helium convective core easier . convective dredge up in a red supergiant could also possibly reduce the helium core mass , though this effect seems to be small in the present models . with rotation , many of these lighter u - stars would also be red supergiants in the end because of primary nitrogen production at the hydrogen - helium interface , and the effect would be diminished . above 30 @xmath0 the structure of the sh and u stars are very similar because both are red supergiants .    in the past , it has sometimes been assumed that the helium core mass uniquely determines the presupernova structure of a star and hence that structure would be roughly independent of the metallicity for two stars that made the same helium core mass . while this is qualitatively true , the figures here show that this assumption is not very accurate for supernova progenitors below 30 @xmath0 . [ sec : s2.3 ]    with the new surveys , it is possible to address a point of possible concern - the choice of mass ( 2.5 @xmath0 ) and time ( `` presn model '' ) for evaluating the compactness in equation  [ eq : compact ] . the full evolution of a massive star to the point that its iron core collapses and possibly powers an explosion is being characterized by a single number here . what motivates the choice of this particular point in space and time ? @xcite explored the effect of evaluating @xmath19 at different fiducial masses and concluded that 1.75 @xmath0  might be a better discriminant of explosion characteristics rather than the 2.5 @xmath0  chosen here and in @xcite . the time of core bounce rather than initial collapse ( `` presn '' ) also seemed a more relevant time for its evaluation . recently , @xcite also examined the choice of 1.75@xmath0 , since the early neutrino signal is more sensitive to the structure around the neutron star . obviously smaller choices than 1.75 @xmath0  would not be sensible since they often lie _ within _ the collapsing iron core itself , and depend on different physics that occurs after that core has already reached high density , but what about values in between ?    to address these questions , the collapse of the s series was continued until the central density reached @xmath20 g @xmath21 . this is about 100 times greater than the central density of presupernova s - series stars . beyond this point , neutrino trapping and , ultimately , the nuclear equation of state , would influence the dynamics ( fig . [ fig : f7 ] ) . as expected , @xmath19 rises with increasing central density since a smaller radius encloses the same mass . it also rises as the mass chosen for its evaluation decreases , since the average density enclosed by r increases faster than r@xmath22 . running kepler well beyond @xmath23 g @xmath21 to nuclear density gives values of questionable accuracy , but strongly suggests that very little further evolution will occur in the compactness so long as the sampling mass remains greater than 1.75 @xmath0 . as the figure shows , the compactness evaluated at 2.5 @xmath0  changes very little during the collapse from the presupernova star to high central density . at that larger radius , the hydrodynamical response time is longer than the time for the denser part of the core to collapse . ( red ) , and later when the central density reaches @xmath24 g @xmath21 ( blue ) . each pair of curves represent @xmath19 evaluated at the same time in evolution , but at different points in fiducial mass : 1.75 @xmath0(top ) , 2.0 @xmath0  ( middle ) and the default 2.5 @xmath0  ( bottom ) . little evolution occurs during the collapse from about @xmath25 g @xmath21 to @xmath26 g @xmath21 if the compactness is evaluated at 2.5 @xmath0 , but for smaller fiducial masses the contrasts seen at 2.5 @xmath0  are significantly amplified with time . [ fig : f7],scaledwidth=35.0% ]    substantial variation does occur , however , for smaller sampling masses . smaller fiducial masses show a non - linear amplification of the structure in the compactness curve with time . large values get much larger than adjacent smaller ones , suggesting the development of islands of stars that may be hard to blow up . this is due to the tendency of high density regions to collapse faster under the influence of their own gravity and is particularly apparent for stars in the 18 to 22 @xmath0  mass range . this may account for some of the variability in outcome seen by @xcite for supernovae in this mass range . the robustness of @xmath5  evaluated at 2.5 @xmath0 , however , and its strong correlation with the compactness evaluated at smaller masses suggests that we can continue to use our standard choices of time and mass for its evaluation . it should be kept in mind , however , that structures that seem small in some of the plots for presupernova stars may become amplified by the further collapse . [ sec : s3 ]    surveys of stellar evolution find a complex , non - monotonic behavior for the compactness parameter as a function of main sequence mass and metallicity . why is this so ? why does nt the compactness vary smoothly with mass as it would in a polytrope with a single index ?    four factors drive the development of the compactness profile for massive presupernova stars . one obvious effect is the contraction of the core to higher densities in order to burn heavier fuels . this contraction increases the average density inside 2.5 @xmath0  and causes @xmath5  to grow with time . another is the tendency of lighter stars to have lower entropy cores and be more degenerate . degeneracy is responsible for the transition between stars that make planetary nebulae and those that make supernovae around 8 @xmath0 , but the effects of degeneracy on the post - carbon burning evolution continue to be important throughout the entire mass range studied here . third , as has been noted previously @xcite , is the disappearance , around , 20 @xmath0 , of central convective carbon core burning . above 20 @xmath0 , carbon fuses away without contributing a large excess of energy generation over what neutrinos are carrying away in the center of the star . model , a typical supernova mass , ( top : ) from carbon burning until silicon shell ignition , and ( bottom : ) from the end of first oxygen shell burning until 5 minutes before collapse . only the inner 4 @xmath0  is plotted and the time axis is plotted as the logarithm of time until iron core collapse ( @xmath27 ) measured in years . cross hatched regions indicate convection and integer numbers annotate convective carbon burning shell episodes . the convective carbon core burning starts about 1000 years before death and is followed by three episode of carbon shell burning . oxygen burning starts about two years before death and is preceded by a brief stage of neon burning . silicon burning is the last convective episode in the center . at silicon depletion , no oxygen shell is active , only a carbon shell at about 1.9 @xmath0 . shortly afterwards though , the third oxygen shell ignites in the incompletely burned ashes of the second shell . this final vigorous shell merges with the carbon and some neon burning shells shortly before the star dies . the operation of these shells , especially the ones between 1.5 and 2.5 @xmath0 , affects the compactness parameter.[fig : f8 ] , scaledwidth=48.0% ]    fourth and frequently overlooked ( though see barkat 1994 ) are the effects of convective carbon and oxygen shell burning . lighter stars can have three or more carbon burning shells in addition to core carbon burning ( fig . [ fig : f8 ] ) . as these stars evolve , their central regions become increasingly degenerate , especially after carbon depletion in the star s center , and , by oxygen ignition , the concept of a chandrasekhar mass has some approximate meaning , especially for stars lighter than 30 @xmath0 . shells that burn outside the effective chandrasekhar mass , roughly 1.7 @xmath0  depending on thermal corrections , but inside the point where the compactness is measured at 2.5 @xmath0 , will considerably modulate the compactness parameter . [ sec : s3.1 ]    fig . [ fig : f9 ] shows that @xmath5  remains nearly independent of the main sequence mass until after carbon ignites in the core . prior to this time , the fiducial point at 2.5 @xmath0  lies well within a much larger star or helium core that can be characterized by a single polytropic index . the tendency of @xmath28 to decrease as @xmath29 , as given by equation  [ eq : poly ] , is offset by the slight increase of the burning temperature with mass , resulting in a nearly flat curve . as a function of main sequence mass for the s series models up to the point when the carbon mass fraction decreases below 1% in the center of the star ( `` c.dep '' ) . at each successive stage , the star has contracted to higher density and the average curve for @xmath5  rises . below about 20 @xmath0  carbon burns convectively in the core and this reduces @xmath5.[fig : f9 ] , scaledwidth=35.0% ]    by the time carbon has disappeared from the center of the star though , things have started to change . a pronounced dip in @xmath5  develops for the lower mass stars , and the curve shows an abrupt , small rise at 20 @xmath0 . for current code physics and solar metallicity , 20 @xmath0  is the mass below which carbon burns convectively ( exoergically ) at the stellar center as opposed to radiatively ( endoergically ) . convection brings additional fuel into the burning region increasing the effective supply by approximately the ratio of the convective mass to the mass of the burning region . because of the high temperature sensitivity of the carbon fusion reaction , the energy generating region is small , so the enhancement is significant . during this longer time , neutrino losses carry away both energy and entropy , not only from the carbon convective core , but from the hot , overlying helium - rich layers supported by it . this loss of entropy exacerbates the natural tendency of lighter cores to have greater degeneracy and accelerates the development of a compact , white dwarf - like core structure - a dense degenerate core surrounded by a much less dense extended envelope . , averaged over the inner 1 @xmath0  of the star is shown as a function of initial mass at the times of ( top : ) carbon ignition and core carbon depletion , ( bottom : ) oxygen ignition , oxygen depletion and silicon depletion . very negative values of @xmath30 correspond to a non - degenerate gas and large positive values imply degeneracy . the cores of the lower mass models grow much more degenerate due to neutrino losses between carbon depletion and oxygen ignition . the structure around 20 @xmath0  in the curve for oxygen ignition reflects the brief appearance of strong , deeply - sited carbon convective shells in this mass range (   [ sec : s3.2 ] and fig . [ fig : f11 ] ) . these degeneracy curves at oxygen ignition and depletion are nearly mirror images of the @xmath5  curves at these times ( fig . [ fig : f13 ] ) showing the key role played by degeneracy in determining the compactness.[fig : f10],scaledwidth=48.0% ]    as fig . [ fig : f10 ] makes clear , _ central _ carbon burning , by itself , is not the whole story though . most of the increase in central degeneracy occurs after carbon has been exhausted in the star s center . it is during the period between central carbon depletion and oxygen depletion that neutrinos cool the core appreciably and , for masses below about 30 @xmath0 , give it a white dwarf - like structure . during this long cooling - off time , contraction of the inner core is frequently held up by two , or even three vigorous carbon burning shells ( fig . [ fig : f11 ] ) , the first igniting shortly after carbon core depletion typically about 0.5 @xmath0  from the center . as we shall see , the presence and location of these shells is strongly correlated with how carbon burns in the center , so 20 @xmath0  remains a critical mass . it is ultimately the timing and location of these carbon convective shells , however , not just central carbon burning , that account for the structure developed during and after oxygen burning . greatly affects the location and extent of all three convective shells . for slightly larger masses , the first convective shell ( `` 1 '' ) disappears . after that mass carbon burns from the center radiatively out to the base of the former second shell . the location of the third shell is drawn as a dashed line when it that shell exists for less than a year at the end of the life of the star . ( bottom : ) difference between starting times and core collapse for various convective episodes ( oxygen burning - green and carbon burning - blue ) as a function of initial mass . [ fig : f11],scaledwidth=48.0% ]      [ sec : s3.2 ]    fig . [ fig : f12 ] shows the evolution of the compactness parameter for the s series stars from the time carbon is depleted at the center of the star until oxygen is similarly depleted . this is a critical period when the major features of the compactness emerge ( fig . [ fig : f7 ] ) . , as a function of initial mass at the time of central carbon depletion ( lower curve ; see also fig . [ fig : f9 ] ; oxygen ignition ( central temperature equals 1.6 gk ) ; and oxygen depletion ( central oxygen mass fraction equals 0.04 ; top curve ) . non - monotonic behavior caused by carbon burning shells is clearly evident above 18 @xmath0  by the time oxygen ignites , and a broad peak in @xmath5  is emerging around 18 - 24 @xmath0  with a `` notch '' between 19 and 21 @xmath0 . note the strong anti - correlation of @xmath5  with degeneracy below 22 @xmath0  ( fig . [ fig : f10 ] ) especially around 20 @xmath0 . [ fig : f12],scaledwidth=35.0% ]    the typically low values for @xmath5  below 18 @xmath0  result from the operation of the first two convective carbon shells ( fig . [ fig : f11 ] ) . while the star is supported by these shells , its inner regions radiate neutrinos and become cool and degenerate . the rapid variation of @xmath5  in the range 18 to 21 @xmath0  reflects the response of these three shells to the disappearance of the convective core . when that happens , _ the first carbon burning shell disappears _ , and the second and third convective shells move inwards dramatically in response to its loss . these events do not happen simultaneously . the center of the star ceases to be convective at 20.3 @xmath0 , but the first shell has already gone out at 19.9 @xmath0 . the second and third shells  inwards until this mass is reached , and then move outwards for bigger stars . the offsets in these masses and the rapid motion , first in , then out , of shells two and three underlie some of the complex structure in degeneracy and compactness seen around 20 @xmath0 . is plotted in each case and the format is similar to fig . [ fig : f8 ] , but here all stars , except the 20@xmath0  model ( which has a very tiny convective core ) , burn carbon radiatively in their centers . prominent in all figures is a powerful carbon burning shell igniting roughly 10 years before the star dies . note the migration outwards of the base of this shell as the mass increases . the shell that ignites farther out supports the star while the inner core emits neutrinos and becomes increasingly degenerate . the shorter kelvin - helmholtz time for the larger carbon - depleted core in the heavier star also leads to oxygen ignition at a substantially earlier time . the 20@xmath0  and 23 @xmath0  models end up with a large value of @xmath5 , while 20.5@xmath0  and 28 @xmath0  models give a smaller value . see also fig . [ fig : f13 ] . [ fig : f13],scaledwidth=48.0% ]    from 21 to 30 @xmath0 , the compactness is most sensitive to whether , when oxygen core burning ignites , there is , or very recently has been an active carbon convective shell inside of about 2.0 @xmath0  ( fig .  [ fig : f13 ] ) . an active shell supports the overlying star and relieves pressure on the oxygen core , allowing it to evolve as if it had a smaller effective mass . smaller oxygen cores burn longer , emit more neutrinos , cool more effectively , and become more degenerate . the relevant carbon shell can either be the third ( for 18 - 20 @xmath0 ) or the second ( for 21 - 30 @xmath0 ) shell . if that shell is inside 1.4 @xmath0 , the oxygen core lacks sufficient mass to ignite @xcite , so the most important mass range for the convective shell is 1.4 to 2.0 @xmath0 . the location of the boundaries of these convective carbon shells are given in fig . [ fig : f11 ] . the third convective shell is given there as a solid curve only for the mass range where it has a lifetime longer than one year . shorter lived shells do not greatly affect the evolution prior to oxygen depletion . based upon this criterion , the third shell only has a major effect in the range 19 to 21 @xmath0 . for heavier stars , the second shell dominates ( note that , for the sake of continuity , we persist in calling it the `` second '' shell even though the `` first '' convective shell has been replaced by radiative burning ) . this behavior is illustrated by the convective history during carbon and oxygen burning for stars of four masses : 20.0 , 20.5 , 23.0 and 28.0 ( fig . [ fig : f13 ] ) . none of these four stars experience central convective carbon core burning or `` first shell '' burning ( in 20 @xmath0  there is a very tiny convective carbon core ) . the relevant carbon shells are thus `` 2 '' and `` 3 '' , and they are labeled as such in the figure . at oxygen depletion , the 20.5 @xmath0  and 28.0@xmath0  models have small values of @xmath5  ( 0.13 and 0.15 , respectively ) , while the 20.0 and 23.0 have larger ones ( 0.19 and 0.21 ) .    in the 20 @xmath0  model , shells 2 and 3 complete their major burning about a year before oxygen ignites in the stellar center . shell 3 is located at 1.32 @xmath0 , which is too small to allow oxygen to ignite . later , following a brief , inconsequential `` blip '' of convective neon burning , oxygen burning ignites in a core devoid of any nuclear energy sources out to more than 3 @xmath0 . the entropy is higher and the degeneracy is less at both oxygen ignition and depletion ( fig . [ fig : f10 ] ) . oxygen burning ignites late , with less than a year remaining in the star s life . ultimately the star dies with a relatively large compactness parameter , 0.26 . the 20.5 @xmath0  model , though only slightly different in mass , has a very different evolution . the third carbon shell now has its base at 1.62 @xmath0  which is sufficiently large for oxygen burning to ignite before carbon shell burning is over . when oxygen does ignite at the center , the resulting expansion extinguishes the carbon shell , leaving unburned carbon in the region outside 1.62 @xmath0 . some neutrino losses have already occurred in the oxygen before it ignited , but over the course of its burning , which commences much earlier than in the 20.0 @xmath0  model , more cooling occurs and the core becomes very degenerate . by oxygen depletion , the compactness parameter has a value of 0.13 , which declines still further to 0.10 when the star dies .    for 23.0 @xmath0 , the evolution resembles that of the 20.0 @xmath0  model more than it does the 20.5 @xmath0  model . by now , the third carbon shell is igniting so late as to be unimportant . shell 2 , burning at 1.33 @xmath0 , however , again delays oxygen ignition to a late time . the oxygen burns in a core that , at the time has essentially no active burning shell except the helium shell far above it ( off - scale ) . it is easier for the layers above the core to contract since there is a time with no active burning , neither core nor shell . once it ignites , the oxygen core burns rapidly and at low density , producing a less degenerate silicon core . by oxygen depletion , the compactness parameter is 0.21 which rises to 0.42 in the presupernova star , the largest compactness of any model under 30 @xmath0 .    for 28.0 @xmath0  and nearby masses , the evolution resembles that of the 20.5 @xmath0model , but with no third shell . the second shell has now moved out to 1.91 @xmath0  which is enough to allow core oxygen burning to ignite . there is always support by an active burning front , reducing the contraction of the outer core . oxygen burning lasts longer and occurs in a core of smaller effective mass . the compactness at the end of oxygen burning is back down to 0.15 , though this rises to 0.28 in the presupernova star .    beyond 28 @xmath0 , the second carbon shell continues to move out and the effective core mass when it ignites oxygen grows . when the shell crosses 2.5 @xmath0  at 30 @xmath0 , the compactness parameter rapidly increases once more . this marks the end of the epoch of carbon shells . all heavier stars will have large compactness parameters ( fig . [ fig : f3 ] ) with relatively small variations introduced by the _ oxygen _ burning shells . the location and strength of the carbon shells prior to oxygen depletion thus sculpt a highly variable result out of what was a comparatively smooth curve at carbon depletion . in some sense , 20.5 @xmath0  and 28 @xmath0  stars resemble each other more at death than do 20.5 @xmath0  and 20 @xmath0  stars . above 30 @xmath0 , the carbon shells lose their influence . after oxygen finishes burning in the center , powerful convective _ oxygen_-burning shells develop that can also modify the structure either by themselves , or by interacting with existing carbon shells . the continued evolution serves chiefly to amplify the features already present at oxygen ignition ( fig . [ fig : f14 ] ) . having developed a degenerate core , further shell burning pushes matter out at 2.5 @xmath0  and reduces @xmath5 . the location of the first oxygen shell at silicon depletion correlates tightly with @xmath5 . its base is located between 0.8 and 1.2 @xmath0  and is set by the extent of the oxygen convective core . its extent correlates tightly with the location of its base . similarly , the second oxygen shell , when there is one , sets atop the extent of the first . the farther out these shells and the later their ignition , the more extended the core , which results into a shallower density gradient around the iron core later on ( higher @xmath5 ) . the peak in @xmath5  in fig . [ fig : f14 ] , around 23 to 25 @xmath0 ,  is enhanced by the migration outwards in mass of this oxygen burning shell . unlike the carbon burning shells , however , the outwards movement of the oxygen burning shell is not monotonic in mass , since the extent of the oxygen core varies due to the effect of previous carbon shells . as it recedes inwards from 24 to 28 @xmath0 , the @xmath5curve declines , leading to a pronounced dip at 28 to 30 @xmath0 . above 30 @xmath0 , the oxygen burning shells move out rapidly .    ) and silicon depletion . the tendency of existing features to be amplified by the contraction , as seen previously in fig . [ fig : f13 ] , continues . bottom : the location of the outer extent of the oxygen convective core , which is also the base of the first oxygen shell is given as the lower curve labeled `` o.core '' , while the extent of the first oxygen shell is the curve labeled `` 1st o.shell '' . note the correlation between shell locations and the compactness . [ fig : f14],scaledwidth=45.0% ]      [ sec : s3.4 ]    fig . [ fig : f15 ] shows that , while the peak in @xmath5  at 23 @xmath0  continues to grow due to core contraction and the operation of the oxygen shells , the principal qualitative change in @xmath5  during this final stage of the star s life is a `` chopping up '' of the peak that existed at silicon depletion around 20 @xmath0  into finer structures . given that these structures will be amplified in the collapse ( fig . [ fig : f7 ] ) , they may yet affect the outcome of the explosion . , as a function of initial mass at silicon depletion and for the presupernova star . the peak at 24 @xmath0  grows due to core contraction and fine structure is imprinted between 19 and 22 @xmath0  by a strong oxygen burning shells in the last few days of it s life . [ fig : f15],scaledwidth=35.0% ]    this rapid variation for a narrow range of masses results from an almost random decision by the star whether to burn oxygen or silicon in a shell first , the decision depending upon fine details of all that happened before . fig . [ fig : f16 ] shows the history of the convection following central silicon depletion in two stars of the s series with masses 20.1 and 20.2 @xmath0 . note especially the timing of the ignition of silicon shell burning compared with the growing extent of the strong oxygen burning shell with a base at 1.71 @xmath0  in the 20.1 @xmath0  star and 1.59 @xmath0  in the 20.2 @xmath0  star . in both stars this convective oxygen shell is very powerful , and eventually grows to include both the carbon and neon burning convective shells in one large aggregate shell shortly before the star dies . the oxygen shell in the 20.1 @xmath0  star drives a similar linking of convective layers and creates a large positive energy generation rate , but only _ after _ the silicon shell has already burned out and the star is close to collapsing . in the 20.2 @xmath0  star , on the other hand , the strong oxygen shell ignites _ before _ silicon shell burning and operates for a much longer time . ( top ) and 20.2 @xmath0  ( bottom ) stars of s series . each plot shows the evolution inside inner 4.5 @xmath0  during the -2 to -6 in the log of time until collapse . this last few days of evolution roughly corresponds to the time - evolution from the start of silicon ignition in the core until the presupernova stage . shaded regions describe convection . notice in 20.2@xmath0 , where the @xmath5(measured at presupernova ) is low , the silicon ignites int the core earlier and the shell c - o burning is very different . [ fig : f16 ] , scaledwidth=48.0% ]    in the 20.1 @xmath0  star , oxygen shell burning initially generates , only a limited convection zone that slightly alters the structure of the star before silicon shell burning ignites and slows the contraction . in the 20.2 @xmath0  star , the oxygen shell becomes strong first and has already engulfed a large fraction of the co core before the silicon shell ignites . the 20.1 @xmath0  star ends up with a high @xmath5 , the 20.2 @xmath0  star with a low one . the strong oxygen shell in the latter pushes matter out to larger radii and leaves a less extended structure inside . .[fig : f17],scaledwidth=35.0% ]    this sensitivity of the presupernova compactness to the timing of oxygen and silicon shell burning suggests that a correlation might exist between the time from silicon core depletion to death and @xmath5 . this is indeed observed ( fig . [ fig : f17 ] ) . the longer the strong oxygen shell supports the star after the initial iron core has formed from central silicon burning , the smaller @xmath5  for the presupernova star .    these energetic , merged oxygen , neon and carbon burning shells during the last hours of a massive star s life are a robust feature seen in many models of stars above about 15 - 20 @xmath0 , both here and in previous studies @xcite . their study would be an interesting , though perhaps challenging topic for 3d simulation . if such small differences in the central structure of the star at silicon core depletion can lead to major changes in presupernova structure , even for stars differing in mass by only 0.1 @xmath0 , how sensitive might the results be to other `` hidden variables '' , such as the time step criteria and zoning . [ fig : f18 ] shows that , for some masses , a substantial variation in outcome can also be caused by varying these variables . the location of the oxygen shell is a sensitive function of the many convective episodes that went on in the same star during its earlier evolution , and to the physics used to treat that convection .     and b : 20.1 - 20.5 @xmath0 . 1a , b - different zoning : default ( thick ) and @xmath31 of default ( dot - dash ) . 2a , b - time step : default ( thick ) , half ( dot - dash ) and double ( thick - dashed ) . [ fig : f18 ] , scaledwidth=48.0% ]    as previously noted by @xcite , this sensitivity to small changes is especially strong for the stars in the 19 to 21 @xmath0  range . to quantify this behavior better , a subset of the s series was calculated using the same code , compiler , and physics , but slightly altered zoning and time step criteria . in each case , the standard calculations were recomputed using 50% finer resolution ( i.e. , @xmath31 the default zone size ) , and , in separate calculations , the default time step was multiplied by @xmath32 and @xmath33 . even these relatively small changes in operational parameters were frequently sufficient to provoke the star into two different final states . for comparison , we studied two distinct regions of main sequence masses , 17.1 - 17.5 @xmath0- region `` a '' , and 20.1 - 20.5 @xmath0- region `` b '' . as shown in figure 2 , the observed variations of @xmath5  in these two regions are quite different . panels 1a and 2a of fig . [ fig : f18 ] show that the solutions in region a are comparatively robust . decreasing the time step and or changing the zoning does little to alter the final values of @xmath5 . increasing the time step by a factor of @xmath33 does cause some mild variation , @xmath34 . panels 1b and 2b of fig . [ fig : f18 ] , however , show a much more pronounced sensitivity near 20 @xmath0 . the resulting values of @xmath5seem to be bimodal . small variations in mass , zoning , or time step can send the star down one path or the other . typical models in the surveys used 20,000 time steps to reach the presupernova stage and employed approximately 1000 mass shells of variable thickness . it is certainly feasible to double or even quadruple both the resolution and time steps , but the values used here are already conservative . to test the effect of finer time steps , the 21 s - series models from 19 to 21 @xmath0  were rerun with twice as many time steps . while @xmath5  for some individual stars did change significantly , the overall appearance of the pattern , including its rapid , but bounded variation between @xmath5  = 0.1 and 0.25 , did not .      above 30 @xmath0 , if mass loss is neglected , the structure of the core becomes much simpler . the carbon shells have moved far out and no longer affect the solution . the oxygen core is larger and the oxygen shells farther out . because of the larger mass , the degeneracy is reduced and the compactness is always large .    the most notable feature above 30 @xmath0  is a pronounced dip in @xmath5  around 50 @xmath0 . the dip is not large , but is seen in both the u - series and sh - series stars , albeit at slightly different masses ( fig .  [ fig : f3 ] ) . it is also present in the compactness plot for the bare co cores studied later in  [ sec : s5 ] . this behavior can be traced to the presence of a strong , extended convective oxygen burning shell during the post - silicon burning evolution of stars over about 50 @xmath0 . the lighter stars lack this shell ; the heavier ones have it . starting at 50 @xmath0  for the u - series , this shell is present at silicon depletion with a base at 1.8 @xmath0 . moving to heavier masses , the shell grows larger and its base moves outwards , reaching 2.5 @xmath0  at 65 @xmath0 . there is a sharp density decline at the base of the shell and because of this migration outside the fiducial point for measuring @xmath5 , the compactness parameter rises again as the star mass passes about 60 @xmath0 .    whether this shell is present or not depends upon the timing of silicon core ignition and oxygen shell burning . recall the key role played by the carbon shell and oxygen ignition for stars in the range 21 - 30 @xmath0  (   [ sec : s3.2 ] ) . when the carbon shell was situated far enough out , oxygen burning would ignite before carbon shell burning was done with major consequences for the compactness . here , the oxygen shell plays the role of the former carbon shell . if it burns far enough out , the silicon core can ignite earlier . in this case , however , igniting silicon does not blow out the oxygen shell . it persists until the end . and has little effect on @xmath5  in this mass range . the bases of the 1st shells are not plotted for clarity , but they almost always perfectly match with the extents of the cores . notice , how the silicon core size responds as the oxygen core overgrows the silicon shell near 45 - 50 @xmath0 . [ fig : f19],scaledwidth=35.0% ]    fig . [ fig : f19 ] shows the locations of various silicon and oxygen burning episodes as a function of mass for the sh - series models . though it lacks the time dimension of a full convective history plot , the figure shows that the size of the oxygen convective core increases monotonically with mass for these heavy stars . where the oxygen convective shell ignites is pegged to the extent of that convective core . at about 50 @xmath0 , the extent of the oxygen convective core passes the value needed to ignite silicon core burning . an oxygen shell ignites at the edge of that depleted core at about 1.8 @xmath0 . this is sufficiently large that silicon core burning also ignites , and after that , silicon shell burning , but both ignite without putting out the oxygen shell . indeed , the oxygen shell persists until the iron core collapses and forces the compactness parameter down .    above 50 @xmath0 , the compactness parameter continues to increase until about 70 @xmath0 , where the onset of the pulsational pair instability becomes important @xcite . [ sec : s4 ]    studies of massive stellar evolution by different groups yield substantially different presupernova structures and nucleosynthesis . a key quantity is the co core mass . the yield of elements heavier than helium , for an energetic explosion , is the difference between the mass of the collapsed remnant ( neutron star or black hole ) and the co core mass . as we shall see , the compactness is also quite sensitive to this quantity . table 2 shows the disparate values for the co core mass obtained by various groups studying the problem . for example , a 20 @xmath0  star may end up with a co core of anywhere from 2.71 to 6.59 @xmath0  for various choices of convective physics and assumptions about rotation . even for non - rotating stars the range is 2.71 to 4.99 @xmath0 . why are the values so different ? might using the co core mass as a discriminant for the presupernova compactness alleviate some of the sensitivity to uncertain code physics ? the effect of rotation is not surprising . rotation adds new mixing processes ( e.g. * ? ? ? * ) that increase the masses of both the helium and co core . shear mixing in differentially rotating stars can also erase some of the sensitivity to the treatment of semiconvection and convective shell boundaries . the large residual range of co core masses in non - rotating stars reflects these uncertainties in the treatment of convection and warrants discussion . uncertainties in the treatment of convection for non - rotating massive stars can be lumped into three categories : a ) the use of time - dependent mixing length theory ( mlt ) to describe the overall transport ; b ) uncertainties in semiconvection , i.e. , what to do about stellar regions that are unstable to convection by the schwarzschild criterion , but not the ledoux criterion ; and c ) the treatment of convective boundaries ( i.e. , overshoot and undershoot ) . a detailed review is beyond the scope of this paper . most modern studies use time - dependent mlt , if only for the lack of a viable , implemented alternative . semiconvection and convective overshoot can sometimes have similar effects since both allow convective regions on their tops and bottoms to grow smoothly when circumstances warrant it , and prevent the unrealistic fracturing of convective regions into multiple concentric shells . semiconvection additionally slowly mixes regions that , owing to composition barriers , would not have mixed if the ledoux criterion were strictly applied . a relevant case in point is the convective helium core . it has been known for some time that the co core mass is quite sensitive to the efficiency of semiconvection during helium burning @xcite . in situations where semiconvection is very small ( or zero ) , a mathematical instability can develop that leads to the bifurcation of the helium convective core and the production of an unusually small co core . in an algorithm that looks at the composition to determine convective instability ( i.e. , ledoux criterion ) , a single zone that is marginally unstable might temporarily become flagged as stable due to a trivially small increase in the mean atomic weight ( helium plus carbon and oxygen ) in the zone beneath . if this zone is flagged as stable against convection , burning during the next time step will increase the atomic weight farther in the region that is still convectively coupled to the burning at the center of the star without raising it in the outer core . this makes convection between the two shells even more difficult . the helium convective core splits permanently into two pieces . in the outer one , helium never burns to completion and the resulting co core for the star is small . if enough semiconvective mixing ( or rotation ) is included though , the splitting of the the convective core does not happen . in fact , considering the real situation in three dimensions , this numerical instability seems unphysical and should be avoided . overshoot mixing can also help bridge this artificial segregation . it also leads to an increased growth of the helium convective core at the end of helium burning that affects not only the co core mass , but the carbon mass fraction when carbon burning ignites . simulations by meakin and collaborators ( e.g. * ? ? ? * ; * ? ? ? * ) confirm that some type of mixing will always take place at convective boundaries defined by either the ledoux or schwarzschild criterion , although the exact formulation of this mixing remains uncertain . the treatment of convective overshoot and semiconvection in kepler has been discussed previously @xcite . semiconvection is parametrized by a diffusion coefficient that is manufactured from the local convective diffusion coefficient and the radiative one : @xmath35 where @xmath36 is the convective diffusion coefficient the zone would have had based upon the schwarzschild criterion , and @xmath37 is a free parameter that multiplies the radiative diffusion coefficient , @xmath38 . usually @xmath39 so that @xmath38 dominates the transport . the actual value of @xmath37 is not known , or even that such a simple formula captures the essential features of semiconvection . in past surveys using kepler , @xmath37 has usually been taken to be 0.1 , and , unless otherwise specified , that is the value used in this paper . in practice , this generally assigns the semiconvective diffusion coefficient a value of near 10% of the radiative one . overshoot mixing is presently treated very crudely in kepler . single zones at the top and bottom of regions flagged as unstable to ledoux convection are slowly mixed with a diffusion coefficient given by equation  [ eq : semid ] . the convective diffusion coefficient is calculated based upon the schwarzschild criterion : @xmath40 with @xmath41 assumed to be a parametrized fraction , @xmath42 and with @xmath43 having a default of 0.01 . typically the convective diffusion coefficient so calculated is much greater than @xmath44 so that the mixing occurs on a time scale given by this quantity . this prescription lacks a physical basis , but at least allows convective zones the liberty of slowly growing into regions where mild entropy gradients would otherwise disallow mixing .    to explore the sensitivity of the compactness and the co core mass to semiconvection and overshoot mixing , a subset of 16 solar metallicity models in the mass range of 15 to 30 @xmath0  was calculated multiple times using various settings for the both semiconvection and convective overshoot mixing multipliers . the masses studied were the integers from 15 to 30 and a total of 176 new models was calculated ( tables 1 and 2 ) . in these , the semiconvective multiplier , @xmath37 , was varied from 0.001 to 0.1 ( 16 models each with @xmath37 = 0.001 , 0.0025 , 0.005 , 0.0075 , 0.01 , 0.0125 , 0.015 , 0.0175 , 0.02 , 0.05 , and 0.075 in addition to the standard models with @xmath37 = 0.1 ) and are collectively called ss series . additionally , a set of 16 models was run with standard semiconvection ( @xmath37 = 0.1 ) , but zero overshoot mixing ( @xmath45 ) = 0 ( so series ) . all models were generated on the zero age main sequence and run until core collapse . the co core mass does not change significantly after carbon depletion , however , and presupernova values can safely be compared with other studies that stop at an earlier time .    fig . [ fig : f20 ] shows that the compactness exhibits great deal of diversity depending on the efficiency of semiconvection . the mass resolution in these figures is significantly degraded by considering only 16 masses points between 15 and 30 @xmath0 , but the location of the peak in @xmath5 , when one exists , is clear . reduced values of @xmath37 give smaller co cores , more than spanning the space of published values . the carbon mass fraction at carbon ignition is also affected by the choice of convection physics , both in response to the altered co core mass and the mixing of helium at the outer boundary of the helium core that is altered when semiconvection is turned off . consequently the compactness , which is quite sensitive to both the co core mass and the carbon mass fraction , vary significantly .    ) are shown in comparison with the s - series models ( top panel ) with @xmath46 . the s - series is plotted with the same @xmath47 increments as was used in ss series calculations . [ fig : f20],scaledwidth=48.0% ]    for small reductions in @xmath37 , the peaks found for @xmath5  in the standard runs is robust . both the locations and amplitudes do not vary much . though not explored , we expect that small changes upwards of @xmath37 would have given similar results . when @xmath37 declines below 0.02 however , significant changes occur . with smaller @xmath37 , it takes a higher mass main sequence star to give a specific co core mass , so the peak that was at 23 @xmath0  shifts rapidly to higher values , ultimately leaving the grid . by @xmath37 = 0.01 , the peak is no longer clearly discernible in the mass range examined .    continued reduction in @xmath37 ( bottom panel , fig . [ fig : f20 ] ) results in a compactness curve that loses most of its non - monotonic structure and increases steadily with mass . however , such small values imply a co core that is small compared with all previous studies in the literature ( table 2 ) . for @xmath37 = 0.001 , the 15 , 16 and 17 @xmath0  models closely resembled previous models calculated with kepler and other codes for stars with mass near 10 @xmath0  @xcite . oxygen ignition occurs off center in these models and propagates as a convectively bounded flame to the stellar center ( e.g. * ? ? ? a similar result was obtained for @xmath37 = 0.0025 for the 15 and 16 @xmath0  models . such stars are known to have a very steep density gradient surrounding the iron core , which might be a boon to those trying to blow up stars in a wider mass range , but contribute almost nothing to the nucleosynthesis of intermediate mass elements . were the observational limit on the maximum mass supernova to be pegged not far above 18 @xmath0  @xcite , models of this sort would be inconsistent with the theory of stellar nucleosynthesis . a model with 12 @xmath0  and @xmath37 = 0.001 failed even to ignite oxygen burning and would likely become an agb star with a degenerate neon - oxygen core and not a supernova . since this would be in contradiction with a large volume of observational data showing that stars down to about 8 @xmath0 do explode , we must consider values of @xmath37 much below 0.01 as unrealistic .    turning off convective overshoot mixing in kepler also has a significant effect on the co core mass and carbon mass fraction , though not so extreme as that caused by large variations in semiconvection ( table 2 ) . [ fig : f21 ] shows the co core sizes and final presupernova @xmath5  of the so - series models compared with those from s series . the elimination of overshoot reduces the co core sizes ( top - panel ) , and as a consequence , the @xmath5  curve is shifted to higher masses . for example , central carbon burning shifts to being radiative at a higher mass compared to the s - series models with overshooting .    . [ fig : f21],scaledwidth=48.0% ]    these results , which imply that the success or failure of exploding a presupernova model of given mass will vary appreciably depending upon who calculated the model and the description they used for convection physics , are troubling . however , semiconvection and overshoot mixing are real physical processes that will eventually be better understood . some of the very small values for @xmath37 used here may not be physically realistic ( e.g. * ? ? ? * ) and the high threshold mass that they imply for making a supernova is inconsistent with observations . zeroing overshoot mixing is also unrealistic . bifurcation of the helium convective core probably does not occur . real stars rotate , and rotationally induced mixing reduces the sensitivity to semiconvection and overshoot mixing . we believe that our standard choices for @xmath37 and @xmath45 ( 0.1 and 0.01 ) are a good compromise for non - rotating stars since they avoid the unphysical splitting of the convective core , agree reasonably well with the results for rotating models ( table 2 ) , and are not presently at odds with any existing physical prediction . clearly this is an area where more work is needed . see @xcite and references therein for some recent insights .    as a function of initial main sequence mass ( top ) and final co core mass ( bottom ) for the s - series models ( continuous ) , a subset of ss - series models with @xmath48 ( triangles ) and so - series models ( diamonds ) . though the definition of boundary of the co core inside a massive star is ambiguous , its mass can be a better discriminant of @xmath5  than the initial mass in most cases . [ fig : f22],scaledwidth=48.0% ]    our results also suggest that @xmath5 , and presupernova structure in general , will be more similar for two models having the same helium and co core masses than for two that simply started with the same main sequence mass . carrying this to the extreme , the outcome of the explosion of two stars with the same co core mass will be more similar than for those having the same main sequence mass , but different co core masses . fig . [ fig : f22 ] shows this to be approximately true . the offset in peaks seen in fig . [ fig : f21 ] for the stars calculated with and without overshoot mixing is greatly reduced if @xmath5is plotted against co core mass rather than main sequence mass . similarly , the large discrepancy with the s series found for low values of semiconvection ( e.g. , @xmath49 ; fig . [ fig : f20 ] ) is significantly improved . note that the allowed range of co masses in the figure is much more limited for 15 to 30 @xmath0  stars if semiconvection is reduced significantly . just what constitutes the `` co core mass '' inside a massive star is ambiguous however , because there usually is no sharp decline of pressure at its edge (   [ sec : s5 ] ) . here we arbitrarily took the value inside of which the carbon mass fraction rose above 0.20 . differing values of the carbon mass fraction at carbon ignition due to changes in the initial composition or code physics must also be considered . even similar physics often has different implementations in various stellar evolution codes , so to further explore the sensitivity of the compactness curve to the way it is calculated , we repeated part of our study using a completely different code . mesa @xcite is a popular code for the study of stellar evolution which has been successfully applied to a large variety problems in astrophysics and is freely available . mesa has recently been applied to the calculation of presupernova models @xcite . it uses similar nuclear physics to kepler for energy generation during hydrogen through oxygen burning and a very similar equation of state and opacity tables . it can be constrained to use similar mass loss rates . here the red giant mass loss rates of @xcite was used in both codes : @xmath50    however , the two codes have quite different prescriptions for semiconvection and overshooting . given the previous discussion , those differences can be expected to cause changes that might be educational . mesa also uses different nuclear physics during silicon burning and iron - core formation that affect the post - oxygen evolution . in particular , mesa continues to use a 19 isotope network to describe the energy generation and changes in composition after oxygen depletion , while kepler uses a 125 isotope `` quasi - statistical equilibrium '' ( qse ) network @xcite . the latter is more stable and less prone to spurious temperature transients and is also able to follow accurately the appreciable changes in electron mole number that happen in these late phases .    since the major features of the pre - supernova @xmath5(m )  curve , including the existence and location of peaks , have already been imprinted by oxygen depletion , we elected to compare the results of the two codes at that point , not for the pre - supernova star . this way , the post - oxygen burning differences are avoided and also the run time greatly reduced . mesa is a continually evolving code while kepler is more `` mature '' , and for this study we used mesa version 4930 . the stars studied are a subset ( table 1 ) of our s series , that is , non - rotating stars with solar metallicity .    this version of mesa employed ledoux convection . for the mixing length parameter , @xmath51 , the free parameter @xmath52 was taken to be 2.0 , which is within the range found in the literature . in those zones flagged as semiconvective , the @xcite model was used . this calculates the diffusion coefficient as :    @xmath53    where @xmath54 is the specific heat at constant pressure and the efficiency is tuned by the free parameter @xmath55 . we carried out a brief survey varying this parameter between @xmath56-@xmath57 , and found its effect to be rather complicated . for the masses tested , this range of @xmath55 gave a spread of about @xmath58 in the co core mass , @xmath47 difference in final mass , and a one order of magnitude change in radius . in general , for larger values ( @xmath59 ) the mesa results approached that of the pure schwarzschild case . the very low values , on the other hand , gave results that were very similar to one another , but quite different results were obtained using the larger values . the greatest sensitivity for the co core size was found to be between @xmath60 and @xmath61 . in our survey @xmath55 was taken to be @xmath62 as a  midpoint  between these two extremes , which is also well within the range used in literature ( e.g. * ? ? ? we also note that the work to eliminate or greatly constrain this parameter is underway @xcite . overshoot mixing in mesa is accomplished either by fully mixing the zones flagged as overshooting (  step  overshooting ) or by applying an exponential cutoff to the diffusion coefficient calculated for convection in zones that are unstable by the ledoux or schwarzschild criterion . in our survey we have employed the exponential decay formalism , that is :    @xmath63    where @xmath64 is the diffusion coefficient from mlt , @xmath65 is the height , and @xmath66 , the pressure scale height . the efficiency is then controlled by the free parameter @xmath67 ( see @xcite,@xcite for further details ) .    , as a function of initial mass curves at oxygen depletion in the core for each series . [ fig : f23],scaledwidth=48.0% ]    in general , the value of @xmath67 is quite uncertain as well @xcite . a value of 0.016 was used by @xcite for his studies of agb stars . since the stellar context here is quite different , we have performed some additional test runs . the effects of five different choices for @xmath67 ranging from 0 to 0.025 were explored ( table 2 ) . in the end , a value of 0.025 was selected for detailed study since it produced very similar co core sizes to the s - series models calculated using kepler(fig . [ fig : f23 ] ) . a value as small as herwig s had little effect and values much bigger than 0.025 resulted in the bump in compactness around 23 @xmath0  becoming very broadened and shifted to lower masses . though our understanding of overshoot mixing is not adequate to rule this out , there was no compelling region to use larger values that would have given larger co core masses for the heavier stars than previously published by others ( table 2 ) . , evaluated at carbon ignition , carbon depletion , oxygen ignition and oxygen depletion . compare with fig . [ fig : f9 ] and fig . [ fig : f13 ] for the equivalent runs using kepler . [ fig : f24],scaledwidth=35.0% ]    two sets of mesa models were ultimately calculated - one with overshooting ( mov series ) and another without ( m series ) . each set had 16 models in the mass range of @xmath68 with @xmath69 increments , and both used the ledoux criterion ( table 1 ) . as expected , the results from kepler and mesa were nearly identical during hydrogen and helium burning , but started to diverge following helium depletion when the structure became more sensitive to the treatment of semiconvection and overshoot mixing .    as shown in fig . [ fig : f24 ] , by carbon depletion , the discrepancy had become significant , especially around the critical mass for the extinction of central carbon convection , 20 @xmath0 , and continued to grow during oxygen burning . [ fig : f23 ] and table 2 show the comparison at oxygen depletion between the kepler s series and mesa mov series . the quantity @xmath5  for the mov series shows a non - monotonic structure similar to that in the s series , but shifted in mass by about 2@xmath0 . since the co cores are nearly identical for mov- and s - series models , the difference must result from either a different composition at carbon ignition or a different treatment of convection after carbon ignites , or both . the carbon mass fractions at carbon ignition are indeed different in the two studies . mass fractions in the mov series range from 0.34 to 0.24 , declining roughly linearly as the mass rises from 15 to 30 @xmath0 , while the corresponding limits for the s series are 0.21 to 0.18 . the larger carbon abundance probably accounts for most of the shift in the peak of the mov plot . as we shall see (   [ sec : s5 ] ) , however , the treatment of post - carbon burning convection also matters . even bare co cores with identical initial compositions often have significantly different compactness when evolved to oxygen depletion in the two codes . the results from the m series are not shown in fig . [ fig : f23 ] , but table 2 shows that the resulting co cores are much smaller , and the compactness thus differs appreciably , especially for the lighter stars .    despite the shift , it is encouraging that the non - monotonic structure of @xmath5  as a function of mass in fig .  [ fig : f23 ] and , qualitatively , the location of its peak , can be achieved using standard settings for convection physics in two very different codes . the study also highlights the dependence of the outcome on the co core mass . because of the sensitivity of the results to the carbon mass fraction at carbon ignition , outcomes will depend on the rates used for its creation and destruction . three major reaction rates are involved , each of which has some associated uncertainty : 3@xmath72 , @xmath70c(@xmath71o , and @xmath70c+@xmath70c . it is the competition of 3@xmath72 and @xmath70c(@xmath71o that sets the ratio of carbon to oxygen produced by helium burning . increasing the former raises the carbon yield while increasing the latter makes it smaller . as we have shown , small reductions in carbon mass fraction tends to shift the compactness curve downwards , especially the critical mass , 20 @xmath0 , above which carbon burns radiatively . the effects of varying the three uncertain reaction rates on massive stellar evolution , and nucleosynthesis in particular , have been previously explored by @xcite , @xcite , and @xcite for @xmath73 , and @xmath70c(@xmath71o and by @xcite for @xmath70c+@xmath70c , but none of these works explicitly focused on how these rates affect the _ structure _ of presupernova stars . we will not attempt a survey of all possibilities at the present time , but focus here on just the one rate for @xmath70c(@xmath71o .     reaction rate . curves are labeled by a multiplier times the rate of @xcite . a value of 1.2 was used in the rest of the paper . the current 1@xmath74 range is approximately 0.9 to 1.3 . [ fig : f25],scaledwidth=35.0% ]    historically , the large possible range for this reaction rate has been a source , perhaps even _ the _ major source of nuclear physics uncertainty in massive stellar evolution . recent developments ( e.g. * ? ? ? * ) have narrowed the error bar for this rate to less than 20% , which is not so much greater than the uncertainty in 3@xmath72 @xcite . for small errors in the rates , it is the ratio of these rates that matters most , so variation of @xmath70c(@xmath71o can act , qualitatively , as a surrogate for exploring the uncertainty in 3@xmath72 . the standard value used in this paper for @xmath70c(@xmath71o is 1.2 times @xcite who used s(300 kev ) = 146 kev b ; that is , our effective s - factor at 300 kev , near the gamow energy for helium burning , is 175 kev b. @xcite now report laboratory measurements and analysis that gives s(300 kev ) = 161 @xmath75 19 stat - 2 sys + 8 kev b , or a range 0.9 to 1.3 times the buchmann value . a subset of our s - series models was thus recalculated with varying multipliers on the buchmann rate . thirteen different masses of presupernova star were calculated in the mass range of 19 to 25@xmath0  ( 0.5@xmath0  intervals ) for each multipliers of 0.7 , 0.9 , 1.3 , 1.5 and 1.9 times the buchmann rate ( see table 1 ) . these are the sb series ( table 1 ) . [ fig : f25 ] and tables 2 and 3 show the results . decreasing the rate for @xmath76 significantly increases the carbon mass fraction at helium ignition , but affects the mass of the co core very little . the greater ( or smaller ) carbon mass fraction shifts the peak in @xmath5 , formerly seen at 23 @xmath0  for a multiplier of 1.2 to higher ( or lower ) values . for a multiplier of 1.9 , the peak would disappear altogether , and is not plotted for clarity . however , such large values are outside the current experimental range . the current ( 1 @xmath74 ) error bar for the rate translates into a multiplier between 0.9 and 1.3 , but even this smaller range can shift the explodabilty of massive stars enough to have a dramatic potentially observable effect on presupernova masses and nucleosynthesis . a total uncertainty of less than 10% in the combined error of 3@xmath72 and @xmath76 may be needed to pin down the mass range of supernovae that explode by the neutrino transport model to an accuracy of less than 2 @xmath0 . [ sec : s5 ]    the importance of the co core mass motivates a separate study of the compactness of presupernova stars evolved from _ bare _ co stars of constant mass . it is to be emphasized that there is no direct correspondence between the results of evolving an isolated co core and those obtained for a complete massive star with the same co core mass embedded inside a helium core . the lack of a precipitous drop in density and pressure at the edge of a co core inside a massive star results in an evolution that is qualitatively different . the definition of a co core there , even by composition , is ambiguous @xcite . how does one count a partly burned helium shell ( and these are a common case ) ? nevertheless , co cores are simple to evolve and bypass the uncertainties associated with convection during the helium burning phase . if they display the same sort of compactness systematics seen for the full stars , one can have greater confidence in the result .    as a function of mass at central oxygen depletion for bare co stars studied with kepler and mesa . top : ks ( solid ) and ku ( dotted ) ; bottom : ms(solid ) and mu(dotted ) . [ fig : f26],scaledwidth=48.0% ]    a total of 380 co cores were evolved using both kepler and mesa ( table 1 ) . mesa used @xmath77 and @xmath78 and kepler used @xmath79 and @xmath80 . the calculations followed the evolution from central carbon ignition until the presupernova stage for the kepler cores and until the core oxygen depletion stage for mesa models . the core masses ranged from 3 to 25@xmath0  with varying increments : 0.1 @xmath0  increments were used for 3 to 10 @xmath0 , 0.5 @xmath0  was used for 10 to 20 @xmath0 , and 1 @xmath0was used for 20 to 25 @xmath0 . co core masses of 3 to 9 @xmath0  correspond , approximately , to the main sequence masses of 15 to 30@xmath0  studied in the s series . for low metallicity models ( u series ) or other models with no mass loss ( sh series ) , the maximum co core mass ( at 65 @xmath0 ) was 25 @xmath0 , and that sets the upper limit . the physics used to calculate the co core models was identical to that in the full star calculations . the different mass fractions of carbon for the s and u stars ( figure 5 ) , however , necessitated separate surveys . because of convection , the composition of a co core of given metallicity and mass can be taken as roughly constant . a good fit to the abundances in fig .  [ fig : f5 ] is @xmath81 the mass fraction of oxygen was @xmath82 in each case . these 4 sets of calculations are denoted as : ks , ku , ms and mu . , i.e. - ks stands for kepler cores with initial composition from solar metallicity full stars ( s series ) , and mu stands for mesa cores with initial composition from low metallicity full stars ( u series ) . [ fig : f26 ] shows @xmath5  at core oxygen depletion point for the 4 series ( see also table 1 ) . the same sort of non - monotonic variation in @xmath5  seen previously in the full star models is evident in all cases but , as expected , cores with initial composition from u series ( ku and mu ) , have their first peak at a lower mass than the cores with initial composition from solar metallicity stars ( ks and ms ) . this same sort of offset was also seen in the @xmath5  curve for full stars of different metallicities ( fig . [ fig : f3 ] ) .    for presupernova stars derived from bare co cores of a given mass and composition . ( top : ) for the mass range 3 to 12 @xmath0 , corresponding to main sequence masses from about 15 @xmath0  to 30 @xmath0 . note the large bump near 6 @xmath0  corresponding to the similar peak for main sequence masses 22 to 24 @xmath0shown in fig . [ fig : f3 ] . the offset between solar ( ks ; green ) and low metallicity ( ku ; red ) is also similar . ( bottom : ) a similar plot for the heavier co cores . the agreement with fig . [ fig : f3 ] is striking.[fig : f27],scaledwidth=48.0% ]    above co core masses of 11 @xmath0 , the differences between runs using kepler and mesa are minimal . below that , one sees similar structures for ks and ms , for example , but there are differences . the major peak for the ks stars at 7 @xmath0 , is near 5.5@xmath0  for the ms stars and the peak of ku stars is also shifted down about 1 @xmath0  in the mu stars . even after removing the major dependence on co core mass by running bare co cores , there is still some residual difference resulting from the sensitivity of the location of carbon burning shells to convective and semiconvective physics in the two codes . the overall similarity of the patterns is encouraging , however . the compactness is non - monotonic with mass and there are regions of mass that will easier to explode than others . the kepler calculations of co cores were continued until the presupernova stage , and fig . [ fig : f27 ] shows the compactness parameter as a function of core mass for both the ks and ku series at this final stage . a comparison with fig . [ fig : f3 ] is interesting . the peaks at 22 and 23 @xmath0  in the full star models for the u and s series are reflected well in their equivalent mass co cores , 5.8 @xmath0  and 6.2 @xmath0 , respectively . the good agreement continues to heavier masses where the dip at 50 @xmath0  for the u - series appears as a dip at 18.5 @xmath0  in the ku - series . the mass of the co core in the u - model 50 @xmath0  star was indeed 18.5 @xmath0 . similarly , although not so precisely , the co core of the sh - series 55 @xmath0  star was 20.6 @xmath0 , close to the dip at 24 @xmath0  seen in the lower panel of fig . [ fig : f27 ] . the maximum mass co core in the sh series was 27 @xmath0  for the 65 @xmath0  model . for the u series , the corresponding co core mass was 26 @xmath0 . [ sec : s6 ]    while characterizing the complex structure of a presupernova massive star with a single number is a gross over - simplification , considerable insight can be achieved by studying the systematics of the compactness parameter , @xmath5 . the compactness of a given presupernova core turns out to be sensitive to a variety of inputs , including not only the star s initial mass and composition , but the way convection is handled in the code , the nuclear physics employed , and the code used for the calculation . stellar modeling , at least of massive stars , is , in a sense , a statistical science . different groups will almost universally obtain different results for the structure and composition of a non - rotating presupernova star resulting from a given main sequence mass star , even without the complication of mass loss . however , for the same input physics , there should exist a robust pattern of results that can be sampled by running a large set of models . this is a paper about one such pattern , the compactness parameter . our study confirms the non - monotonic nature of @xmath5  noted in previous works , and explores its causes . [ fig : f7 ] shows that the choice of 2.5 @xmath0  as a fiducial point for evaluating the compactness , though arbitrary , is a good indicator of the behavior at other locations deep inside the star , but outside the iron core . whether @xmath5  is evaluated at the  presupernova model ( @xmath83 km s@xmath4 ) , or at core bounce matters little at 2.5 @xmath0 , but for points deeper in , small variations seen in @xmath5  for the presupernova model are amplified during the collapse . small variations in main sequence mass might therefore result in significant differences in outcome for an explosion model @xcite . extensive new surveys of presupernova evolution are reported for both solar metallicity and low metallicity ( 10@xmath11 z@xmath84 ) stars (   [ sec : s2 ] ; table 1 ) . while the emphasis here is on presupernova structure and compactness , this is clearly a rich data set for exploring nucleosynthesis and supernova explosion physics , and this will be done elsewhere . major differences are found in the presupernova compactness for stars of the same initial mass , but different metallicity (   [ sec : s2.2 ] ) . the most obvious differences result from the differing degrees of mass loss expected for the two compositions , but there are subtle secondary effects involving the amount of helium convected into the helium core near central helium depletion . the growth of the helium convection zone results from a different structure at its outer edge which , in turn , is affected by the strength of the hydrogen shell , whether the star is a red or blue supergiant , and how much hydrogenic envelope remains following mass loss ( in the solar metallicity stars ) . even with mass loss suppressed , a substantial offset remains between presupernova stars of solar and low metallicity . generally speaking , the low metallicity stars become difficult to blow up at a lower mass ( fig . [ fig : f3 ] ) and may make more black holes .    the underlying cause for the non - monotonic behavior of @xmath5  with initial mass derives from the interaction of carbon - burning and oxygen - burning shells with the carbon - depleted and , later , oxygen depleted core (   [ sec : s3 ] ) . as has been noted previously , the end of convective carbon core burning at around 20 @xmath0  for standard stellar physics , has a major effect , but not so much due to the entropy decrease during central carbon burning itself as to the effect this has on the ensuing shells . above 20 @xmath0 , both central carbon burning and what was the first convective shell switch to radiative transport . energy generation in excess of pair neutrino losses is unable to support core against further contraction , and the next two carbon convective shells are pulled down . their gradual migration outwards again causes the rise and decline , and rise again of @xmath5  above 20 @xmath0 . above 30 @xmath0 , @xmath5  rises to large values and , in the absence of mass loss , keeps on rising . there is an appreciable dip at about 50 @xmath0 , though , due to the outward migration of the oxygen shell which essentially plays the same role as the carbon shells did for lower mass stars .    broadly speaking , then , the evolution of massive stars separates into four mass intervals characterized by similar values of @xmath5 . for a standard choice of physics and solar metallicity these are stars from 8 to 22 @xmath0 , 22 - 25 @xmath0 , 25 - 30 @xmath0 , and above 30 @xmath0(e.g . , fig . [ fig : f3 ] ) . stars below 22 @xmath0  and between 25 and 30 @xmath0  have low @xmath5  and might more easily explode . between 22 and 25 @xmath0and above 30 @xmath0 , it will be more difficult for the traditional neutrino - transport model to succeed . there can be large variations within those ranges and not all stars lighter than 22 @xmath0  will be easy to explode @xcite , but the distribution of stars that make black holes , for example , or important contributions to nucleosynthesis could be bi - modal . galactic chemical evolution studies using such non - monotonic yields have yet to be done , but may give quite different results , e.g. , for the s - process and oxygen synthesis @xcite . while the pattern seen in @xmath5  _ vs _ mass is robust , the locations of its peaks and valleys are sensitive to the treatment of semiconvection and convective overshoot in the code . these cause variations in the mass of the co core and the carbon mass fraction when carbon burning ignites . reducing semiconvection or convective overshoot makes the co core smaller and this is a leading cause of the large variations seen in the literature for many studies of presupernova ( table  [ tab : cocore ] ) . smaller co cores can mimic the effect of reducing the main sequence mass in a calculation where the convection physics is held constant and thus make cores with steep surrounding density gradients up to a higher mass . however these cores , while easier to blow up will be essentially nucleosynthetically barren . extremely small values for semiconvective efficiency and overshoot mixing can lead to a discrepancy between the lightest stars expected to explode as supernovae based upon observations , about 8 @xmath0 , and that calculated by theory . caution must be exercised when using the co core mass as a discriminant , however , because the outcome is also very sensitive to the carbon mass fraction , and even to the treatment of convection after helium depletion . most of the models in this paper were calculated using the kepler code , but the mesa code was also used to calculate the compactness curve at oxygen deletion where its essential features have already been determined . reasonable agreement is achieved between the two codes provided that some amount of semiconvection and convective overshoot mixing is used in both . the non - monotonic nature of the compactness plot , in particular , is apparent in the results of both surveys . the effects of varying the uncertain reaction rate for @xmath70c(@xmath71o within its error range were also explored and found to have an appreciable effect . a total error of just 10% in the combined error for 3 @xmath72 and @xmath70c(@xmath71o can shift features in the compactness plot by a solar mass or more . to the extent the compactness is related to the difficulty of exploding a given model using neutrinos , the changes for nucleosynthesis and supernova `` cut - off mass '' are significant at that level .    and binding energies for presupernova stars of the s series ( top : ) the binding energy above the iron core shown as a function of initial mass . despite the general increasing baseline trend , the relative variations closely follow that of @xmath5 . ( bottom : ) the final value of @xmath5  is plotted against the binding energy outside of the iron core . [ fig : f28],scaledwidth=48.0% ]    because the co core mass plays such a major role in determining @xmath5 , we explored its systematics using a grid of `` bare '' co - cores . these differ from the co cores of massive stars which are defined by composition changes and do not have sharp declines in density and pressure at their edges . nevertheless , the same non - monotonic compactness curve results . this both confirms the robustness of the @xmath5  relation found in the other surveys and suggests that , for a given presupernova model , its final co core mass , not its main sequence mass might be a more accurate structural indicator . the same behavior was found for co cores studied using mesa . picking the co core mass circumvents some , though certainly not all of the uncertainty introduced by ambiguous treatments of semiconvection and convective overshoot mixing during core helium burning . for wr stars and other stars that have experienced severe mass loss , the presupernova co core mass will be the best indicator of presupernova structure . even for bare co cores , however , the mass fraction of carbon at carbon ignition remains an important `` hidden variable '' that will be sensitive to both the code physics and the metallicity of the star . while this paper has emphasized the relation between compactness of the core and its explodability in terms of the neutrino transport model , the results are actually more general . [ fig : f28 ] shows the binding energy of the matter outside of the iron core of the presupernova star as a function of its compactness . stars with larger values of @xmath5  will be more difficult to explode by _ any _ mechanism , both because of the higher ram pressure immediately after collapse and the larger energy needed to eject the matter outside the iron core . we gratefully acknowledge numerous valuable discussions with alexander heger and bill paxton regarding the kepler and mesa codes and their application to this problem . we also benefited from alex s expert insights into massive stellar evolution on many occasions , and his permission to include unpublished details of the `` u series '' low metallicity models . casey meakin and evan oconnor provided helpful comments on an earlier draft . we thank the anonymous referee for a careful reading of the submitted manuscript and many helpful suggestions . this research has been supported by the national science foundation ( ast 0909129 ) , the nasa theory program ( nnx09ak36 g ) , and the university of california lab fees research program ( 12-lr-237070 ) . kepler & s & zams-@xmath24 g @xmath21 & 15 - 30 & 151 & default + & sh & zams - presn & 30 - 60 & 9 & no mass - loss + & u & zams - presn & 15 - 65 & 86 & @xmath85 , no mass - loss + & ss & zams - presn & 15 - 30 & 176 & weak semiconvection + & so & zams - presn & 15 - 30 & 16 & no overshooting + & sb & zams - presn & 19 - 25 & 65 & varied @xmath76 + & ks & c.ign-presn & 3 - 25 & 95 & co core with s composition + & ku & c.ign-presn & 3 - 25 & 95 & co core with u composition +   + mesa & m & zams-o.dep & 15 - 30 & 16 & default + & mov & zams-o.dep & 15 - 30 & 16 & with overshooting + & ms & c.ign-o.dep & 3 - 25 & 95 & co core with s composition + & mu & c.ign-o.dep & 3 - 25 & 95 & co core with u composition +    ccccccccccc 4.84 & 2.97 & 7.17 & 4.94 & 9.65 & 7.20 & @xmath86 & @xmath87 ( ? ) & c.dep & @xmath88 & @xcite + 4.0 & 2.02 & 6.0 & 3.70 & 8.0 & 5.75 & & o.dep & @xmath89 & @xcite + 4.21 & 2.44 & 6.27 & 4.13 & 8.50 & 6.27 & @xmath90 & @xmath87 & si.dep & @xmath91 & @xcite + 5.68 & 3.76 & 8.65 & 6.59 & 10.0 & 8.63 & @xmath90 & @xmath87 & si.dep & and rotation & @xcite + 3.85 & 2.19 & 5.70&3.54 & 5.18 & 7.73 & @xmath92 & @xmath93 & o.dep & & @xcite + 4.97 & 2.56 & 7.31 & 3.83 & 9.80 & 5.48 & & presn & @xmath94 & @xcite + 5.37 & 3.59 & wr&5.53 & wr&6.63 & & presn & and rotation & @xcite + 4.35 & 3.14 & 6.21 & 4.93 & 8.25 & 6.95 & @xmath95 & @xmath96 & presn & @xmath97 & kepler , s + 4.44 & 3.14 & 7.18 & 5.61 & 9.07 & 7.35 & @xmath95 & @xmath96 & presn & and rotation & kepler + 4.29 & 3.08 & 6.27 & 4.99 & 8.50 & 7.18 & @xmath95 & @xmath96 & presn & @xmath98 & kepler , ss + 4.15 & 2.93 & 6.22 & 4.59 & 8.40 & 6.39 & @xmath95 & @xmath96 & presn & @xmath99 & kepler , ss + 4.09 & 1.96 & 6.00 & 2.95 & 8.30 & 4.40 & @xmath95 & @xmath96 & presn & @xmath100 & kepler , ss + 3.98 & 1.63 & 6.08 & 2.71 & 8.42 & 4.11 & @xmath95 & @xmath96 & presn & @xmath101 & kepler , ss + - & - & 6.22 & 4.96 & 8.28 & 7.00 & @xmath95 & @xmath96 & presn & 1.5@xmath102buchmann & kepler , sb + - & - & 6.21 & 4.90 & 8.26 & 6.98 & @xmath95 & @xmath96 & presn & 1.3@xmath102buchmann & kepler , sb + - & - & 6.17 & 4.96 & 8.23 & 6.99 & @xmath95 & @xmath96 & presn & 0.9@xmath102buchmann & kepler , sb + - & - & 6.16 & 4.94 & 8.21 & 6.99 & @xmath95 & @xmath96 & presn & 0.7@xmath102buchmann & kepler , sb + 4.19 & 2.70 & 6.11 & 4.38 & 8.14 & 6.27 & @xmath95 & @xmath96 & presn & @xmath103 & kepler , so + 3.77 & 1.96 & 5.81 & 3.47 & 7.31 & 3.91 & @xmath95 & @xmath96 & o.dep & @xmath104 & mesa , m + 4.01 & 2.21 & - & - & 7.80 & 5.16 & @xmath95 & @xmath96 & o.dep & @xmath105 & mesa + 4.20 & 2.39 & - & - & 8.06 & 5.47 & @xmath95 & @xmath96 & o.dep & @xmath106 & mesa + 4.41 & 2.61 & - & - & 8.34 & 5.81 & @xmath95 & @xmath96 & o.dep & @xmath107 & mesa + 4.85 & 3.05 & 6.92 & 4.77 & 9.28 & 6.84 & @xmath95 & @xmath96 & o.dep & @xmath108 & mesa , mov + [ tab : cocore ]
Abstract:

Answer: the success or failure of the neutrino - transport mechanism for producing a supernova in an evolved massive star is known to be sensitive not only to the mass of the iron core that collapses , but also to the density gradient in the silicon and oxygen shells surrounding that core . here we study the systematics of a presupernova core s `` compactness '' @xcite as a function of the mass of the star and the physics used in its calculation . fine - meshed surveys of presupernova evolution are calculated for stars from 15 to 65 @xmath0 . the metallicity and the efficiency of semiconvection and overshoot mixing are both varied and bare carbon - oxygen cores are explored as well as full hydrogenic stars . two different codes , kepler and mesa , are used for the study . a complex interplay of carbon and oxygen burning , especially in shells , can cause rapid variations in the compactness for stars of very nearly the same mass . on larger scales , the distribution of compactness with main sequence mass is found to be robustly non - monotonic , implying islands of `` explodability '' , particularly around 8 to 20 @xmath0  and 25 to 30 @xmath0 . the carbon - oxygen ( co ) core mass of a presupernova star is a better , though still ambiguous discriminant of its core structure than the main sequence mass .