stellar mass loss is of critical importance to the later stages of stellar evolution .
stars evolving up the red giant branch ( rgb ) lose mass via stellar winds .
the remaining mass of a star leaving the rgb tip is the primary determinant of its position on the horizontal branch ( hb ; e.g.  @xcite ) .
subsequent mass loss on the asymptotic giant branch ( agb ) ejects the star s entire hydrogen envelope , creating a post - agb star and ( perhaps ) a planetary nebula ( pn ) .
the timescale and end - point of agb evolution may therefore not be determined by the nuclear burning rate , but by mass loss @xcite .
winds from these stars are the most significant producers of interstellar dust and therefore provided the heavy elements present in population i stars , including the sun and solar system @xcite .
debate exists around many finer points of the mass loss process , particularly its variation with fundamental parameters such as luminosity , temperature and metallicity .
firstly , the amount of mass loss that occurs on the rgb remains largely unmeasured .
rgb mass loss appears to determine the ` second parameter ' ( after metallicity ) that defines the morphology of hbs in globular clusters ( gcs ) @xcite .
secondly , it is not proven where dust formation stars , and whether rgb stars can form significant amounts of dust in their ejecta . if many rgb stars do produce dust
, it would add a large number of dust factories ejecting dust into a galaxy s interstellar medium ( ism ) ( cf .
the missing dust source in many galaxies ; e.g.  @xcite ) .
@xcite suggest that dust production occurs episodically over a large portion of the rgb , while @xcite have presented evidence that dust production is confined to the rgb tip , though gaseous mass loss is still present all along both the rgb and agb @xcite .
dust can also help drive a wind through radiation pressure on grains ( e.g.  @xcite ) , though it is not certain whether winds which lack high - opacity amorphous carbon can be accelerated in this manner @xcite .
dustless stars typically lose mass via magneto - acoustic processes in their chromospheres @xcite .
finally , the composition of circumstellar dust and the chemical processes by which it forms are only partly defined .
of particular interest is the rle that metallic iron grains play in dust formation , and how metallicity affects dust production and composition @xcite .    in this study
, we aim to address these points . to do this , we explore dust production in the galactic globular cluster 47 tucanae ( ngc 104 ) .
globular clusters are unique laboratories with which we can study late - stage evolution of low - mass stars : globular clusters contain stars that are among the oldest in the universe , and most clusters are comprised of a single ( or at least dominating ) stellar population with little internal variance in age or chemical composition .
comparisons between clusters therefore probe differences in mass loss created by metallicity , while comparisons within clusters probe differences created by stellar temperature and luminosity . throughout this paper
, we will refer to similar studies : @xcite , which covers the cluster @xmath0 centauri ( ngc 5139 ; hereafter paper i ) ; @xcite , which covers the cluster ngc 362 ( hereafter paper ii ) ; and a work which accompanies this paper ( submitted apjs ; hereafter paper iii ) .
the cluster 47 tuc itself is among the galaxy s most massive and , co - incidentally , one of the closest to us .
its distance has been estimated photometrically as 41904700 pc @xcite and kinematically as 4000 pc @xcite , with 4500 pc being the median distance determination in the literature .
its mass is 69 @xmath1 m@xmath2 @xcite .
its metallicity is determined to be roughly [ fe / h ] = 0.7 , though various studies have placed it anywhere from [ fe / h ] = 0.89 to 0.66 @xcite .
the cluster suffers from little interstellar reddening ( @xmath3 mag  @xcite ) and is estimated to be over 11 gyr in age @xcite .
the cluster hosts an ionised intra - cluster medium ( icm ; @xcite ) .
this ionised medium is undetected in any other cluster ( m15 has a neutral medium with a dusty component ; @xcite ) . integrating the observed intracluster electron density over the cluster core yields @xmath40.1 m@xmath2 of h ii , representing @xmath55% of the intracluster hydrogen @xcite .
the total stellar mass loss from 47 tuc s stars therefore gives us insight into the fate of gas and dust expelled by stars in this cluster .
paper iii more - precisely examines which stars have mid - infrared excess , indicative of dust , finding again that we can not corroborate the presence of dust around the fainter giants . in this work
, we use the data and results from paper iii to derive mass - loss rates for individual stars and ( by summation ) the cluster as a whole , and investigate the composition of the circumstellar dust produced by analysing archival mid - infrared ( mid - ir ) spectra .
the data for this study come from the accompanying paper iii , which identifies those stars in 47 tuc which have infrared excess , indicative of circumstellar dust .
this list of 22 objects is an exhaustive list of those stars in the cluster which we can be relatively certain are producing dust , given currently available data .
these objects are listed in table [ mdottable ] , which also includes the parameters we determine in later sections .
@l@c@c@c@  c@  c@  c@  c@  c@  c@  c@c@  c@  c@  c@  c@  cl@ star & @xmath6 & @xmath7 & @xmath8 & @xmath9 & @xmath10 & & @xmath11 & @xmath12 & @xmath13 & @xmath14 & & notes +   & ( mag ) & ( km / s ) & ( days ) & ( k ) & ( l@xmath15 & @xmath16 & [ 3.6 ] & [ 8 ] & ( % ) & 10@xmath17m@xmath2yr@xmath18 ) & ( kms@xmath18 & ( k ) & sil . & fe & al@xmath19o@xmath20 & feo &   + v1 & 4.03 & 20 & 221 & 3623 & 4824 & 0.936 & 1.173 & 2.213 & 3.6 & 21.0 & 4.0 & 900 & 12 & 88 & 0 & 0 & + v8 & 1.6 & 16 & 155 & 3578 & 3583 & 0.904 & 1.210 & 2.135 & 5.9 & 14.7 & 4.0 & 900 & 7 & 85 & 7 & 1 & + v2 & 2.78 & 23 & 203 & 3738 & 3031 & 1.103 & & & 4.3 & 12.0 & 3.8 & 900 & 5 & 95 & 0 & 0 & @xmath21 + v3 & 4.15 & 22 & 192 & 3153 & 2975 & 0.954 & 1.198 & & 3.2 & 9.4 & 3.2 & 1000 & 0 & 100 & 0 & 0 & @xmath22 + v4 & 1.5 & 18 & 165 & 3521 & 2603 & 0.932 & 1.167 & 1.902 & 5.6 & 11.7 & 3.6 & 900 & 2 & 95 & 2 & 1 & @xmath4rgb - tip + v26 & 0.3 & & & 3500 & 2541 & 1.073 & 0.807 & 0.995 & 0.9 & 5.9 & 3.8 & 1000 & & 100 & & & @xmath23 + lw10 & 1.2 & & 221 : & 3543 & 2324 & 1.073 & 0.876 & 1.017 & 0.6 & 4.2 & 3.2 & 1000 & 0 & 100 & 0 & 0 & + v21 & 1.0 & 7 & 76 : & 3575 & 2301 & 0.977 & 1.083 & 1.355 & 1.1 & 4.9 & 3.6 & 1000 & 0 & 93 & 5 & 2 & + lw9 & 1.1 & & 74 & 3374 & 2204 & 0.973 & 1.003 & 1.228 & 1.1 & 6.0 & 3.3 & 950 & & 100 & & & + v27 & 0.3 & & 69 & 3374 & 2140 & 1.023 & 0.929 & 1.218 & 1.4 & 3.3 & 3.6 & 1050 & & 100 & & & + lee1424 & & & 65 & 3565 & 2122 & 1.006 & 1.182 & 1.390 & @xmath244 & 4.4 & 4.0 & 1050 & & 100 & & & + a19 & 0.7 & & 60 & 3526 & 2096 & 1.034 & 0.992 & 1.240 & 0.8 & 3.6 & 4.2 & 1100 & & 100 & & & + lw12&0.6&&116 & 3713 & 2079 & 1.082 & 0.967 & 1.093 & @xmath251 & low + x03 & & & & 3816 & 1640 & & 0.958 & 1.156 & 3.2 & 5.1 & 1.8 & 750 & & 100 & & & @xmath26 + lw19 & 0.1 & & 40 & 3738 & 1638 & 1.078 & 0.991 & 1.080 & @xmath251 & low + v20 & 0.8 & & 232 : & 3602 & 1575 & 1.053 & 0.928 & 1.039 & @xmath251 & low&&&&&&&@xmath27 + v22 & 0.5 & & 62 & 3684 & 1528 & & 2.504 & 2.834 & @xmath251 & low & & & & & & & @xmath23 + v6 & 0.6 & 7 & 48 & 3763 & 1374 & 1.001 & 1.000 & 1.167 & @xmath251 & low & & & & & & & + v5 & 0.7 & 8 & 50 & 3741 & 1363 & 1.003 & 1.100 & & @xmath251 & low&&&&&&&@xmath23 + v18 & 0.3 & 5 & 83 : & 3692 & 1297 & 0.989 & 1.004 & & 2.4 & 5.3 & 2.4 & 800 & 33 & 67 & 0 & 0 & @xmath26 + v23 & 0.5 & & 52 & 3775 & 1259 & 1.023 & 0.985 & & @xmath251 & low&&&&&&&@xmath23 + v13 & & 12 & 40 & 3657 & 1029 & 0.990 & 1.106 & & 1.2 & 4.1 & 2.2 & 800 & 0 & 85 & 15 & 5 & @xmath26 +   +      accurate  or even approximate  mass - loss rates are notoriously difficult to measure .
one must assume that the wind has certain properties .
perhaps the most crucial of these is the relationship between wind velocity and distance from the star .
we here assume a constant , spherically - symmetric , purely radiatively - driven wind , starting from a near - stationary , dust - forming layer .
as discussed previously ( paper i ; @xcite ) , this may not be the most appropriate description , but suffices due to lack of contrary evidence : we return to this point in the discussion .    before calculating an estimated mass - loss rate
, we can create a measure of the amount of circumstellar dust surrounding a star , irrespective of most modelling uncertainties . to do this
, we define a quantity @xmath11 , the fraction of the star s luminosity that is reprocessed by circumstellar dust . for low optical depths , this should be directly proportional to the surface area of dust grains , therefore approximately proportional to the mass of circumstellar dust .
the value of @xmath11 can be determined by taking the difference between the observed sed and model photospheric sed in the infrared ( yielding the infrared excess ) and dividing this by the total photospheric luminosity ( the value given in paper iii ) .
we list this quantity in table [ mdottable ] .
errors in the value of @xmath11 are introduced by difficulties in determining the integrated luminosity from photometry alone and from the uncertainty in the photospheric spectrum ( both temperature and luminosity ) .
the former can be approximated as @xmath28 , or roughly 50% of @xmath11 , but is dominated by photometric uncertainty at short ( 2.210-@xmath29 m ) wavelengths .
the latter gives a constant uncertainty of around @xmath30% .
we reduce this latter error by fixing the photospheric luminosity to the @xmath31-band flux for all stars except v26 and x03 , which are fixed to the 3.6-@xmath29 m flux .
the value @xmath11 could not reliably be measured for several stars with low mass - loss rates .
the stars appear to fall into three main categories : those with large ( @xmath45% ) reprocessed flux ( v14 , v8 ) , those with low ( @xmath42% ) reprocessed flux concentrated toward long wavelengths ( v13 and v18 ) , and those stars with low ( @xmath322% ) reprocessed flux concentrated toward shorter wavelengths ( the remainder ) .      given the limitations of our understanding and of the model presented above , we stress that the mass - loss rates and velocities we derive are based only on one model for wind acceleration , but should be comparative both internally and with other works which use these same assumptions .    under these assumptions ,
we calculate mass - loss rates by modelling the sed using the radiative transfer package , dusty @xcite . as inputs
, we use our interpolated marcs model of the stellar spectrum , and a standard mrn grain size distribution @xcite . within our dusty model ,
we vary the dust composition , temperature in the dust - forming region ( the inner edge of the dust envelope ) and optical depth of the wind .
we change these parameters iteratively until a best fit is achieved .
we relate the dusty mass - loss rate ( @xmath33 ) to the actual mass - loss rate ( @xmath34 ) as previously ( paper i ; paper ii ) , using the following formula @xcite : @xmath35 , \label{mdoteq2}\end{aligned}\ ] ] for luminosity @xmath10 , dust - to - gas ratio @xmath36 , grain density @xmath37 , and terminal wind velocity @xmath38 .
we assume that @xmath39 g @xmath40 for silicates , 4 g @xmath40 for aluminium oxides , 5.7 g @xmath40 for feo , and 7 g @xmath40 for metallic iron .
for these three species , we use the optical constants from @xcite ; @xcite ; and @xcite , respectively .
the errors in the above assumptions are substantial , and are discussed in section [ mdoterrorsect ] .
the dust - to - gas ratio , @xmath36 , can be estimated as : @xmath41 } ,    \label{psieqn}\ ] ] which is approximately 1/1000 for 47 tuc .
this assumes that the fraction of condensed metals is identical for all metallicities .
we discuss the accuracy of this assumption in section [ mdoterrorsect ] .
table [ mdottable ] shows our calculated mass - loss rates and inferred wind parameters and dust compositions alongside literature pulsation properties .
the striking fact about the dust composition is the fraction of iron grains present .
the presence of metallic iron in these stars is debatable , as its featureless infrared emission can be mimicked by many other phenomena , including large silicate grains and emission from thick circumstellar molecular layers @xcite .
we consider iron to be the most likely source of emission , though the fractions of metallic iron listed in table [ mdottable ] may not be particularly accurate , and we discuss this further in  [ dustcompsect ] .
many of the mass - loss rates are too low to be sensibly measured : we can merely say that some circumstellar dust probably exists .
those rates that can be measured have their seds shown in figure [ sedsfig ] .
the photometry is taken from the sources cited in paper iii .
mid - ir spectra are included from @xcite ( hereafter vlmo ) , and @xcite ( hereafter lph ) .
these objects are summarised in the following luminosity - ordered list :    = 1em    v1 : this star is clearly identified as the cluster s dominant dust producer by vlmo , lph and @xcite ( hereafter itm ) .
v1 has clear excess at all wavelengths @xmath42 m .
this excess , which we attribute to metallic iron dust , dominates even over its strong silicate feature .
v8 : no silicate feature was seen in this star by vlmo , but a clear silicate and aluminium oxide feature was seen by lph . it has previously been suggested ( lph ; paper i ) that v8 may have variable mass loss .
however , the _ akari _ photometry at 3 , 4 , 7 , 11 , 18 and 23 @xmath29 m shows no variation in the ( dominant ) iron contribution .
this means that if the mass loss is variable , strong variation only occurs in the silicate production .
v8 has a complex ir spectrum showing peaks from silicates , aluminium oxide , and from crystalline minerals at 9.7 , 11.3 , 13.1 , 20 and 32 @xmath29 m .
we discuss the carrier(s ) of these bands in  [ dustcompsect ] .
note that the amplitude of the 10-@xmath29 m silicate peak is substantially over - estimated in our model compared to the 20-@xmath29 m silicate peak .
this implies a lower temperature for the silicate grains than the metallic iron grains , which we also discuss in  [ dustcondsect ] .
v2 : a silicate dust feature was observed by lph , though it was too weak to have been seen by vlmo .
iron dominates the ir flux , however .
the 10- and 20-@xmath29 m silicate peaks are not well modelled , as in v8 ( see also  [ dustcondsect ] ) .
v3 : no silicate emission was seen by either lph or vlmo , but strong , featureless ( metallic iron ) emission is present , which is confirmed by _
akari _ photometry .
v4 : the spectrum in lph shows complex dust mineralogy , akin to that in v8 .
this star is roughly at the rgb - tip , as estimated from the luminosity function .
v26 : blended with a star half its brightness , 2.2@xmath43 to the east , the near - ir , _ akari _ and 24-@xmath29 m photometry are affected .
the amount of ir excess is difficult to determine here , but is likely to be quite small .
no discernible features are present in the spectrum from vlmo .
lw10 : also within the dense cluster core , it appears that the _ akari _ photometry ( at 9 & 11 @xmath29 m ) for this object is affected .
no features are detectable in the spectrum from vlmo .
v21 : observed by lph , the spectrum shows the same features as v8 .
these features , and the silicate contribution , are considerably weaker in v21 .
v11 : no obvious features were seen in the spectrum of vlmo , this star was only covered by _ spitzer _ at 3.6 and 5.8 @xmath29 m .
9-@xmath29 m photometry suggests it may have some excess there , but both the _ spitzer _
8-@xmath29 m photometry and the vlmo spectrum ( which covers 9-@xmath29 m ) show no excess .
the derived mass - loss rate can probably be treated as a rough upper limit .
lw9 & v27 : an excess consistent with circumstellar dust is present in both the _ spitzer _ and _ akari _ photometry
. the chemistry can not be well - constrained by this photometry , but it appears that the flux contribution from components other than iron is comparatively small .
lee1424 : lacks any data redwards of 9 @xmath29 m , making constraining the mass - loss rate , dust temperature and dust chemistry very difficult .
the values given should be treated as order - of - magnitude estimates only .
a19 : as lw9 & v27 , however the _ akari _
7-@xmath29 m data suggests little to no mid - ir excess in this star .
v20 & v22 : not shown in figure [ sedsfig ] , they have no near - ir photometry attributed to them , meaning it is difficult to determine the amount of ir excess present . despite this ,
the dust temperature is relatively - well constrained by the presence of _ spitzer _
24-@xmath29 m photometry .
x03 : located in the cluster core , this star suffers from the same problems as v20 & v22 , but clearly displays a stronger excess , making an approximate determination of the mass - loss rate and wind characteristics possible .
v18 : observed by both lph and vlmo , this star has a clear silicate feature .
this feature was not claimed in vlmo due to poor signal - to - noise , but the detection appears to match that from lph well .
it shows substantial mass loss compared to other stars of similar luminosity . the 10- and 20-@xmath29 m silicate peaks
are not well modelled , as in v8 ( see also  [ dustcondsect ] ) .
v13 : highlighted in both lph and itm , the spectrum of this star shows an excess similar to that of v8 .
the contribution of iron dust is comparatively small , though would appear to be non - zero .
this star is isolated from the cluster core , but has poor - quality photometry due to its comparative faintness .      the accuracy of the wind parameters we find is hard to determine .
the major sources of error are listed below :    1 .   a 30% uncertainty inherent in dusty due to an unspecified gravitational correction .
this uncertainty is exceeded in stars with @xmath44 km s@xmath18 , where the effect of stellar gravity will change the variation of dust particle density with radius from the star ( the stellar mass is unspecified in dusty ) .
2 .   uncertainty in the underlying photospheric spectrum .
as the iron dust contribution is essentially a modified blackbody , one can fit a dust temperature to the 3.6- to 8-@xmath29 m data and extrapolate this fit backwards to @xmath16-band , where the dust contribution is essentially zero .
the photospheric contribution can then be set from the @xmath45 photometry , but the accuracy of this fit depends on the stellar photometry .
( this would equally hold for other possible dust species , e.g.  amorphous carbon . )
uncertainty in the amount of ir excess .
apart from identifying how much flux is generated by the photosphere , missing photometry and photometric errors cause uncertainty in the combined flux of the ( photosphere + dust ) contributions .
4 .   uncertainty in the dust composition .
it is not clear for a number of these stars ( v26 , lw9 , v27 , lee1424 , a19 , x03 and any star marked `` low '' ) what the dust chemistry is .
many of these stars have insufficient coverage in the 923 @xmath29 m region to determine whether they contain any dust components other than metallic iron . lw10 and v11 are not covered by lph : ground - based spectra by vlmo show no obvious spectral features , so we presume a pure metallic iron wind , but the spectra are of low signal - to - noise
. the presence of silicates can significantly increase the implied mass - loss rate due to their low opacity .
any carbonaceous dust present would also be hidden by the modified blackbody of metallic iron .
carbonaceous dust typically has a higher opacity : if it is present , the implied mass - loss rate would be decreased .
uncertainty in grain size distribution and grain shape ( assumed here to be spherical , with radius @xmath46 ) .
the relationship between the grains absorption cross - section to radiation ( @xmath47 ) and their mass ( @xmath48 ) .
a high grain porosity would also imply a decreased dust density and therefore a higher mass - loss rate .
particularly , differences in grain properties between the iron and silicate grains could explain why we model such a high proportion of iron grains in the wind ( see  [ dustcompsect ] ) .
uncertainty in dust temperature .
the dust temperature can only be determined to within @xmath4100 k. as @xmath49 , a 15% error in temperature leads to an error of up to @xmath460% in the flux emitted by the circumstellar dust , depending on the implied error in the radius of the dust shell .
this will be most significant in stars that are losing the least mass , as their dust temperatures are more difficult to determine .
uncertainty in the wind driving mechanism .
we presume here that the winds are sustained by radiation pressure alone .
the very low outflow velocities present in at least three cases ( x03 , v18 and v13 ) would suggest that there is insufficient radiation pressure for it to be the sole driver of mass loss .
additional energy input , for example from pulsation , would increase the wind velocity and ( by virtue of eq .
( [ mdoteq2 ] ) ) decrease the mass - loss rate .
there is a related uncertainty in the dust - to - gas ratio , which is assumed to scale with metallicity .
these uncertainties can not be quantified statistically .
lllll parameter & @xmath50 & @xmath51 & parameter variation used & other parameters affected + grain size@xmath52 & @xmath5310% & @xmath5330% & low : @xmath54 , @xmath55 nm , @xmath56 nm & @xmath57 k , 90% iron , 10% silicates +   &   &   & assumed : @xmath54 , @xmath58 nm , @xmath59 nm & @xmath60 k , 88% iron , 12% silicates +   &   &   & high : @xmath61 , @xmath58 nm , @xmath59 nm & @xmath60 k , 80% iron , 20% silicates + grain density & 50% & 50% & low : 1.63 g @xmath40 & +   &   &   & assumed ( solid ) : 6.52 g @xmath40 & + dust formation & @xmath5310% & @xmath5320% & low : 750 k & @xmath62 , 92% iron , 8% silicates +  temperature &   &   & assumed : 850 k & @xmath63 , 88% iron , 12% silicates +   &   &   & high : 950 k & @xmath64 , 82% iron , 18% silicates + dust : gas ratio & @xmath65 & @xmath66 & low : 1 / 5000 & +   &   &   & assumed : 1 / 1000 & +   &   &   & high : 1 / 200 & + velocity distribution & @xmath4 + 500% & n / a & assumed : 0 @xmath67 v@xmath68 km s@xmath18 & @xmath60 k +   &   &   & high : 10 km s@xmath18 ( constant ) & @xmath69 k + photosphere & @xmath5312% & @xmath532% & assumed @xmath53 5% in luminosity & @xmath70 decreased from 0.45 to 0.42 +   & @xmath5315% & @xmath536% & assumed @xmath53 250 k in temperature & + calculation error@xmath71 & @xmath5330% &   &   &   + total & @xmath72 & @xmath73 &   &   +   +    the easiest way to see how these uncertainties affect a star s determined mass - loss rate is to propagate them through the dusty code . the results for v1 ( table [ mdoterrortable ] suggest that , in absolute terms , our mass - loss rates are little more than an educated guess .
however , since a large number of other studies use the same assumptions , one can use these rates to trace differences between stars , providing these properties are largely invariant between stars .
the parameters listed in table [ mdoterrortable ] have different effects .
the grain size , for example , has little effect on the mass - loss rate but significant effects on other wind properties .
conversely , the dust - to - gas ratio have significant effects on the mass - loss rate and velocity but little effect on the other wind parameters ( temperature , optical depth and minerology ) .
clearly the largest effect on the mass - loss rate is whether the wind has an initial velocity driven by , for example , pulsation .
determining the radial velocity profile of the wind would reduce the error substantially .
while this has not been done accurately for such low - metallicity , dust - producing giants , inferences from chromospheric line modelling ( e.g.  @xcite ; @xcite ; @xcite ) appears to show that wind velocities do not decrease substantially from @xmath410 km s@xmath18 in the stages shortly prior to dust formation .
we point out , however , that the highest mass - loss rates ( and by inference the highest wind velocities ) may be ruled out on evolutionary grounds (  [ dustcompsect ] ) .
we emphasise that the errors derived above are not formal and are only estimates .
these compare relatively well to the adopted error of @xmath45@xmath74 for silicate - dominated winds ( e.g.  @xcite ) , where the grain properties are somewhat better determined due to the long history of their study and constraints from their spectral features .    in relative terms ,
however , we can probably assume that the star - to - star variations are much smaller .
in particular , the driving mechanism and hence the wind s initial velocity should not differ as markedly as the uncertainties we assume .
likewise , the dust - to - gas ratio will only depend on the condensation fraction in the wind , which should not vary as much as a factor of 5@xmath74 between stars at the same metallicity .
this reduces the overall error among stars in the same cluster to a factor of @xmath42@xmath74 and allows us to make direct comparisons between stars .
the wind parameters listed in table [ mdottable ] pose several problems for stellar wind theory , which we discuss in this section .
firstly , the projected terminal velocities we derive are very low , implying radiation pressure from the star is not effective in driving the wind : we must consider whether other processes may contribute to , or dominate , the wind acceleration .
secondly , when one considers the evolutionary timescales of the stars in question , our derived mass - loss rates imply that the stars are losing much more material than is predicted from stellar theory .
finally , we can use these problems to constrain some of the wind parameters listed in table [ mdoterrortable ] .
mass - loss rate , pulsation amplitude and pulsation period are all known to correlate with stellar evolution .
this is demonstrated in figures [ hrdvarsfig ] , [ plfig ] and [ mdotrelfig ] .
these can be used to help solve the above problems .      a major problem when calculating mass - loss rates from metal - poor stars
is that the terminal velocities predicted from radiative transfer codes are very low .
this is a particular worry with dusty , as no correction is input for stellar gravity .
in this situation , it is impossible not to have mass lost from a star when dust is present , whereas the real situation may be that radiation pressure on dust is insufficient for it to overcome the stellar gravity field . assuming the dust layers around these stars
do have some outward velocity , a different method of wind acceleration may be required .
dusty provides a minimum mass for a radiative transfer model .
if a star is less massive than this , the error due to the unspecified correction for stellar gravity will be greater than the 30% inherent in the model .
as this error increases to 100% , it is unclear whether radiation pressure is capable of driving the wind .
this is a particular concern for the less - luminous stars with high mass - loss rates : x03 , v18 and v13 . while expansion velocities of 34 km s@xmath18 are not unreasonable when compared to halo carbon stars @xcite , the lower velocities we derive for those stars losing less mass become increasingly unreliable due to this error . @xcite
predict that stars with metallicities below [ fe / h ] = 1 can not drive winds via radiation pressure , and must rely on acoustic ( pulsation ) driving to eject mass from the star ( see also @xcite ) . at [ fe / h ]
= 0.7 , 47 tuc is near this hypothetical limit , meaning we must consider acceleration by pulsation and other sources of energy .
the velocity amplitude of the photosphere in these variables ranges from 423 km s@xmath18 , with v1 , v2 and v3 showing line doubling in near - ir and/or line emission in the optical due to the propagation of shocks in the atmosphere ( table [ mdottable ] ; @xcite ; @xcite ) .
while this velocity greatly exceeds our calculated wind velocity of @xmath44 km s@xmath18 , the bulk material does not necessarily attain this velocity .
nevertheless , the presence of shocks in the extended atmosphere implies that energy is being transferred to larger radii , and implies that this energy is available for driving mass loss .    a comparison between the presence of circumstellar dust and pulsation reveals some interesting results .
all the stars listed in table [ mdottable ] are variable : lee1424 = acvs002331 - 7222.6 , while x03 shows a variation of @xmath75 mag ( n.  matsunaga , private communication ) .
conversely , 22 of 43 variable stars show at least probable evidence of ir excess .
dust - producing stars therefore appear to be ( in 47 tuc ) exclusively variable , and variable stars tend to be dust producing . in both figures [ hrdvarsfig ] & [ plfig ]
there appears no easy way of discerning rgb - pulsators from agb - pulsators : the outliers in figure [ hrdvarsfig ] are such due to crowding in the cluster core ( paper iii ) , while the pulsation periods in figure [ plfig ] do not show the bimodality that one might expect from having rgb and agb stars at different masses .
however , a bimodality may exist in mass - loss rate with a transition occurring around 25003000 l@xmath2 ( figure [ mdotrelfig ] ; cf .  a similar possible bimodality in the gas - production rate in @xcite ) .
the five stars with elevated mass - loss rates are also those with the highest pulsation velocity amplitudes and among those with the longest pulsation periods . as the presence of pulsation and its amplitude both appear to correlate with mass - loss rate , an obvious conclusion would be that pulsation energy is dissipated in the extended atmosphere and helps drive the wind , though we caution that the strongest variables are also the most - luminous stars with the highest capacity for driving a wind through radiation pressure .    further evidence that pulsation sets the mass - loss rate
can be drawn from the apparent invariance of mass - loss over the narrow range of temperatures we probe ( figure [ mdotrelfig ] ) .
the photospheric temperature sets the radius at which dust can condense and therefore ( if other parameters are invariant ) the density of the wind .
higher - density winds lead to more - efficient dust formation , as the timescale for collision between dust - forming molecules is shorter .
an invariance of mass - loss rate with temperature thus means that the formation of dust has little impact on the mass - loss rates present , providing further evidence that the mass - loss rate is ` hard - coded ' into the wind before the dust - formation radius .
if this is true , the rle of radiation pressure on dust may simply be to affect the outflow velocity of the wind .
if the wind is acoustically driven by pulsations , then the wind velocity and thus mass - loss rate is liable to be several times greater than the values listed in table [ mdottable ] , since these sources of energy are not included in our model .
this poses potential problems , as mass - loss rates this high would lead to the dissipation of the star s atmosphere on much shorter timescales than is inferred from star counts .
this problem may be solved if one invokes different parameters for dust in the wind , which could lead to mass - loss rates much smaller than the listed values ( see  [ dustcompsect ] ) .      in paper iii , we discussed the integrated rgb mass loss for a typical cluster star .
we noted it was likely @xmath320.22 m@xmath2 and possibly , to first order , invariant with metallicity .
this invariance would further suggest that the primary driver of rgb mass loss is not related to dust production .
we had also estimated in paper iii , based on the cluster s luminosity function , that one star reaches the agb tip every 80 kyr .    these figures , however , do not correlate with our derived mass - loss rates .
we find here that the cluster s @xmath422 dust - producing stars are losing a total of @xmath761.2 @xmath77 m@xmath2 yr@xmath18 of dust - containing wind . above the rgb tip ( @xmath762000 l@xmath2 ) , where dust production is ubiquitous , 13 stars are producing @xmath78 m@xmath2 yr@xmath18 .
this implies that a staggering @xmath760.8 m@xmath2 is lost in the @xmath761 myr which a star spends on the agb above the rgb tip .
clearly , this is impossible , as this is the entire mass of the star .
this argues for a mass - loss rate of order 48 times lower than we measure .
if additional wind driving mechanisms are important , this factor will increase further .    mass loss via dust - containing winds between 1000 and 2000 l@xmath2 can be similarly determined . here
, nine out of the 56 stars in this luminosity range are producing @xmath42 @xmath79 m@xmath2 yr@xmath18 , implying 0.16 m@xmath2 is lost over a 4.5 myr timescale .
if one assumes that mass loss from these stars is also lower than measured , and that some or all of these stars are agb stars , this implies that mass loss via non - dusty winds may be of similar or even greater importance in removing the estimated @xmath320.22 m@xmath2 from the star between the middle of the rgb and middle of the agb .    at the rates implied in table [ mdottable ] , the iron in the atmosphere of v1 would be entirely ejected within 1000040000 years , depending on the remaining envelope mass ( 0.080.3 m@xmath2 ) .
v8 , v2 and v3 would similarly disperse the iron in their atmospheres within timescales of @xmath2590000 years .
the padova and dartmouth isochrones used in paper iii predict that these stars ( especially v2 and v3 ) still have several hundred kyr left on the agb , implying that our iron mass - loss rates are ( again ) a factor of several too high .
.maximum elemental abundances in dust [ cols= " < , > , > , > , > , > , > " , ]     in  [ windaccsect ] , we find radiation pressure may not be the dominant driver of winds from metal - poor stars , and is perhaps incapable of driving them in some cases .
the simplest solution is to provide momentum from another source , such as stellar pulsations , in order to increase the wind velocity .
unfortunately , this also has the effect of increasing the mass - loss rate , which is already higher than stellar evolutionary theory can accommodate (  [ dustevolsect ] ) . while we expect another momentum source to be present in the most - metal - poor stars ( simply because radiation pressure on dust is that inefficient ) , the limits implied by the mass - loss rates imply this is not the case in 47 tuc .    a solution to this dilemma
should address both these problems by decreasing the mass - loss rate without significantly decreasing the wind velocity .
examining table [ mdoterrortable ] , the simplest change would be to bring the dust - to - gas ratio closer to unity , which would decrease the mass - loss rate and increase the wind velocity .
alternatively , we may consider measures that would increase the wind velocity without appreciable increasing the mass - loss rate , such as changing the size distribution of grains .
the dust - to - gas ratio , @xmath36 , given in eq .
( [ psieqn ] ) , assumes that the dust fraction scales with metallicity ( cf .
this will hold true if every dust species condenses with the same efficiency .
navely , one may imagine that the fraction of each element condensed would decrease toward low metallicity simply because there is less opportunity for metal - metal reactions to take place in a metal - poor environment . from our conventional understanding , it is difficult to see how the dust - to - gas ratio in the wind could be substantially closer to parity .
we can , however , calculate a maximum amount of dust per unit gas volume , by taking the major dust - forming elements ( o , c , fe , si , mg , ca , al ) plus some other potential dust - forming metals ( ni , ti , cr , mn ) and determining how much dust they can form . for 47 tuc , we assume the abundances of most elements are scaled from the solar values of @xcite by [ z / h ] = @xmath80 dex , following @xcite .
we take the following exceptions , however : [ fe / h ] = 0.7 dex @xcite ; [ he / h ] @xmath76 0.02 dex @xcite ; [ o / fe ] = + 0.2 dex @xcite ; [ mg / fe ] @xmath76 [ si / fe ] @xmath76 [ ti / fe ] = + 0.23 dex and [ ca / fe ] = 0.0 dex @xcite .
table [ psitable ] shows the stages of this calculation , by which we conclude that @xmath81 . here
we make the simplistic assumptions that all dust - forming elements start out in the gas phase , that all carbon forms co , that all silicon forms enstatite ( mgsio@xmath20 ) , that iron and nickel condense entirely as pure metals , and that the leftover metals all form their respective oxides .
chemical reactions must be calculated on the basis of number density , but these must be converted to mass densities in order to calculate a dust - to - gas ( mass ) ratio . for reference ,
the same calculation using solar metallicities arrives at @xmath82 .
this approach neglects several points .
carbon may be enhanced in these stars due to dredge - up 2 m@xmath2 ) , carbon - enhanced stars of @xmath40.8 m@xmath2 exist in globular clusters ( e.g.  @xcite ) .
] , but unless it is sufficiently abundant ( c / o @xmath5 0.75 ) to bind all the remaining free oxygen , this does not affect the end result
. we also neglect water condensation , which could use up free oxygen , but this would serve to decrease @xmath36 .
forsterite ( mg@xmath19sio@xmath83 ) may condense out in preference to enstatite , and calcium and aluminium may form into silicates before enstatite , but these would also decrease @xmath36 where they affect it at all .
the formation of feo and iron - rich silicates would increase @xmath36 , but it seems that only a small percentage of iron condenses into non - metallic forms in these stars .
our stated value of @xmath84 is therefore a relatively robust maximum value .
table [ mdottable ] suggests that the number density of iron grains to other dust species in the winds is very high ( typically 10:1 , where it can be measured ) , but the analysis shown in table [ psitable ] would suggest it is of order 1:2 if all dust is condensed .
this would suggest that the fraction of silicate dust that is condensed is quite low ( @xmath45% ) , however we have not yet accounted for two possibilities : that the temperatures of the grain populations may be limited , or that the size distributions ( specifically the surface - area - to - volume ratios ) of silicates and iron may be different .
differing grain temperatures between dust species may exist when the wind is sparse enough that dust - gas collisions can not equalise the thermal energy .
this would mean that more - opaque grains , like iron , are warmer ; more - transparent grains , like silicates are colder for a given radius , and could explain the difficulties we have in fitting both the metallic - iron continuum , and the 10- and 20-@xmath29 m silicate bands simultaneously (  [ mdotsect ] ) .
more significantly , however , if thermal energy can not be equalised via collisions , this will lead to a significant velocity drift between the gas and dust components of the wind .
this will substantially decrease the projected mass - loss rate , while the dust outflow velocity will also increase accordingly .
if the gas outflow velocity is below the turbulent velocity in the wind , this may allow gas to preferentially fall back , making the star more metal poor .
post - agb stars and planetary nebulae are known to be metal - poor for similar reasons ( e.g.  @xcite ) , but it is not clear that this process is happening while the stars are still on the agb where the outflow velocities are lower and the wind densities at the dust - formation radius are greater .
conversely , if dust grain temperatures are equalised and the formation temperatures of silicates and iron are roughly similar , we may examine the effect of different grain sizes by taking two limiting cases . if we take the limiting case that both silicates and iron grains follow the same continuous distribution of ellipsoids , the silicate condensation fraction is @xmath45% , and we arrive at a value of @xmath85 to 1/3000 if all the iron is condensed
this would increase our derived @xmath34 values by @xmath7660% while decreasing @xmath38 by @xmath7630% .
this is the opposite effect to that needed to reconcile our observations with evolutionary theory , therefore we consider this scenario unlikely . in the converse limit , silicate and iron dust
has condensed efficiently , but the iron grains have average radii 20@xmath74 smaller .
this increases @xmath34 by a moderate ( 3060% ) amount , but increases @xmath38 by much more ( 100200% ) , depending on the exact size distribution .
if one is prepared to accept that iron grains grow preferentially in one direction , e.g.  if magnetic effects are important , needle - like cylinders would be created ( cf .
this would lead to a larger @xmath86 , but in this case without an increased @xmath34 .
however , if identical fractions of silicates and iron are condensed , the length / radius ratio of the iron grains would need to be a staggering @xmath761200 , meaning that the grains must be smaller as well .
for example , a population of grains 5@xmath74 smaller with a length / radius ratio of 10 allows a similar fraction of silicates and iron to be condensed with only a slight increase in @xmath34 .
there is no obvious mechanism that would cause iron to be _ preferentially _ condensed at low metallicities , though magnetic effects and/or the basic problem of forming grains may be important ( an iron grain only needs two iron atoms to bond , whereas a silicate grain requires a series of chemical reactions to form mgsio@xmath20 ) .
higher partial pressures of iron may also assist here , though it is not clear why this would happen in a metal - poor environment .
note that in all cases above , we have assumed grains are homogeneous : our calculation of the wind parameters in  [ mdotsect ] assumes all grains are a uniform mix of the included species , while our discussion above takes the opposite limit : where grains of each species exist separately in pure form . in reality , grains aggregate and form inclusions inside each other , which complicates the above arguments . in general , however , the evolutionary constraints provided in  [ dustevolsect ] would suggest that dust condensation efficiency is high compared to solar - metallicity stars ( suggesting @xmath87 ) , and the iron grains are smaller and/or more - elongated than both the silicate grains and the assumed mrn distribution .
a combination of these effects is likely required to reduce @xmath34 by the factor of 48 required without reducing the outflow velocities , @xmath86 , we derive ( table [ mdottable ] ) .      where we are able to determine it , the composition of circumstellar dust and its variation with luminosity is similar to that we see in other clusters in @xcite ( see also figure [ fehfig ] ) .
the evolution of dust composition in 47 tuc broadly agrees with that seen in other clusters at this metallicity .
metallic iron appears to be ubiquitous around dust - producing stars , while amorphous silicates are produced by most stars without any obvious dependance on luminosity .
the published spectra of several stars ( v4 , v8 , v13 , and weakly in v21 ; see @xcite ) show sharp features around 9.7 , 11.3 , 13.1 , 20 and 32 @xmath29 m .
@xcite attribute the 20-@xmath29 m feature to iron - rich wstite ( mg@xmath88fe@xmath89o ) , where @xmath90 ) .
recent work by niyogi , speck & onaka ( submitted ) suggests that all these features , with the exception of the 13.1 @xmath29 m feature ( which is probably due to an aluminium oxide  @xcite ) , can be explained by iron - rich crystalline silicates . here
, the 20-@xmath29 m feature is still produced by an fe - o bond , but it exists in a large crystalline lattice of fayalite ( fe@xmath19sio@xmath83 ) or ferrosilite ( fesio@xmath20 ) .
iron - rich silicates are not expected to form in circumstellar environments , as their condensation temperatures are below that of both magnesium - rich silicates and metallic iron @xcite .
if iron - bearing crystalline silicates are a significant component of some winds , they must exist alongside metallic iron and amorphous silicate species .
one way this could occur is if cooling of the wind proceeds at a faster rate than dust condensation , which may occur in the tenuous winds of low - luminosity , metal - poor stars .
figure [ fehfig ] shows that globular cluster stars which show a 20-@xmath29 m feature may be confined to luminosities below @xmath44000 l@xmath2 .
it is difficult to say much about their metallicity dependance due to insufficient data , especially given silicate production appears to decline below [ fe / h ] @xmath4 1 .
interestingly , the winds found in such stars ( x03 , v13 and v18 ) appear to be quite cool . in more - luminous stars ,
we find wind temperatures to be higher ( @xmath41000 k ) .
the condensation of metallic iron before silicates remains a mystery : for pressures where iron condenses around @xmath41000 k , magnesium - rich silicates should have already formed closer to the star ( e.g.  @xcite , their figure 2 ) .
significant departures from lte , or unusually - high partial pressures of iron ( caused , perhaps , by magnetic confinement or clumping in the wind ) may offer ways to circumvent this problem , but to invoke them as solutions would be speculative without further data .
by whichever process metallic iron is created , if radiation pressure is effective in wind driving , metallic iron condensation should have a dramatic effect on the wind .
its condensation should cause a rapid acceleration of the wind , dropping the density and preventing other dust species from condensing .
this may be what we see in moderately - luminous stars , such as v3 , lw10 and v11 .    if radiation pressure is not an effective driver , the general absence of other dust species may be explained if silicates form using iron grains as seed nuclei . in this case , denser or more - metal - rich winds effectively condense silicates onto iron grains via collisions , while less - dense or metal - poor winds expand before much condensation occurs , leaving bare iron grains .
this scenario may explain the iron grains found in the interplanetary dust particles known as gems ( glasses with embedded metal and sulphides ; @xcite ) .    in the most - evolved stars ( v1 , v8 , v2 and v4 ) , however , we do see silicate emission making an increasingly - strong presence .
meanwhile , there is a hint that the temperature in the dust formation zone ( as traced by metallic iron condensation ) may be decreasing . as noted in
[ dustcondsect ] , different dust species may also have different temperatures .
more detailed modelling of the winds of these stars , allowing for thermal decoupling of the grain species , is required if we are to speculate on factors which cause this temperature change .
it could , however , be linked to stronger pulsation of the star and opacity shielding within the wind . in a cooler wind , however , silicates once again condense out in increasing preference to metallic iron , which contradicts our tentative finding that the silicate grains are cooler than the iron grains in these stars .
in this work , we have used spectral energy distributions and infrared spectroscopy to establish mass - loss rates and investigate mineralogy around the dust - producing stars in 47 tuc .
we summarise our conclusions as follows :        the total mass - loss rate is estimated to be @xmath41.2 @xmath74 10@xmath91 m@xmath2 yr@xmath18 , mainly in the form of metallic iron , but this is probably over - estimated by a factor of 48 due to unconstrained properties of the dust grains .    it remains unclear whether radiation pressure on dust is capable of driving a wind alone , due to the low projected terminal velocities of the stellar winds .
the existence of circumstellar dust appears to correlate with the presence of stellar pulsation .
we theorise that such pulsation may inject energy into the wind acoustically .
variations in the wind chemistry broadly concur with those found in @xcite , with silicates found preferentially in high - luminosity , high - metallicity winds , while iron is found preferentially in low - luminosity , low - metallicity winds .
iron may condense preferentially at low metallicity .
alternatively , iron grains may be coated by silicates in dense , metal - rich winds , but this coating process may not be effective in winds around lower - luminosity and metallicity stars .    *
acknowledgments : * we are grateful to thomas lebzelter for sharing his original data with us
. this paper uses observations made using the _ spitzer space telescope _ ( operated by jpl , california institute of technology under nasa contract 1407 and supported by nasa through jpl ( contract number 1257184 ) ) ; observations using _ akari _ , a jaxa project with the participation of esa ; and data products from the two microns all sky survey , which is a joint project of the university of massachusetts and ipac / cit , funded by nasa and the nsf .                                                                        , s. 2007 , in astronomical society of the pacific conference series , vol .
378 , why galaxies care about agb stars : their importance as actors and probes , ed
. f.  kerschbaum , c.  charbonnel , & r.  f. wing , 145


actual causation has puzzled philosophers since at least the work by .
one way of phrasing the problem is as follows : suppose we know the causal laws that govern some domain , and that we then observe a story that takes place in this domain ; when should we now say that , in this particular story , one thing actually caused another ?
recent work by has also garnered interest in this topic in the ai community .
their account ( which i will refer to as _ hp _ ) constructs a formal definition that tries to capture this intuition in the context of structural models @xcite . to be more concrete , it defines when some random variables @xmath0 of the structural model having values @xmath1 can be counted as an actual cause for @xmath2 .    in previous work ,
i have tried to show that the knowledge representation properties of pearl s structural models can be improved by borrowing representations and techniques from logic programming . in particular , introduced the probabilistic logic programming language of logic programming with annotated disjunctions , for which riguzzi ( 2008 ; 2010 ) implemented sld and slg based resolution algorithms .
further analysis of this language has lead to a reformulation of its semantics , called _ cp - logic _ , which attempts to clarify its causal aspects and examine its relation to pearl s work @xcite .
a more recent paper @xcite showed that pearl s analysis of interventions and counterfactuals in the context of structural models can be elegantly redone in the context of cp - logic , yielding better results for a number of examples , most notably when cyclic causalities are involved .
the goal of this paper is to examine the notion of actual causation in the context of cp - logic .
section [ sec : mot ] will start with some motivation , by explaining a few of the differences between structural models and cp - logic , and offering some hand - waving arguments for why cp - logic might offer a more appropriate setting for the study of actual causation .
the semantics of cp - logic is briefly recalled in section [ sec : prel ] .
then , section [ sec : hp ] dives into the details of actual causation by discussing the hp definition , while section [ sec : ac ] gives my own account .
the traditional way of testing such a definition is to run through a number of `` tricky '' examples and checking whether the obtained answers are intuitively plausible . in section [ sec : ex ] , i will follow suit . finally , section [ sec : imp ] briefly comments on a naive implementation in prolog that can be downloaded to play with my definitions
a structural model is based on a set of random variables ( rvs ) .
each rv has an associated domain of possible values .
the simplest case are boolean rvs , which have @xmath3 as their domain , and can therefore be thought of as propositional symbols or ground atoms .
boolean rvs suffice for typical examples of actual causation , so i will from now on restrict attention to just these .
a structural model then consists of a set of equations @xmath4 , which define the value of the rv @xmath5 in terms of the rvs @xmath6 by a boolean function @xmath7 .
the rvs that appear in the left - hand side of such an equation are called _ endogenous _ , and the other ones _
exogenous_. typically , these sets of equations are assumed to be acyclic .
their meaning is formalized by the obvious possible world semantics : each assignment of values to the rvs that satisfies all of the equations is a possible world .
an acyclic set of equations has the useful property that an assignment of values to the exogenous rvs uniquely determines a single possible world .
pearl uses structural models to represent causal relations : each equation @xmath8 is taken to mean that the causes for @xmath5 taking value @xmath9 are all assignments @xmath10 for which @xmath11 .
this representation is used by pearl to great effect , studying interventions , counterfactuals , and of course also actual causation .    despite its successes ,
however , there is something peculiar about the structural model representation of causal relations : it does not take into account their dynamic nature .
suppose , for instance , that you make the causal claim that dropping a glass causes it to break .
if i do nt believe you , i might challenge you to prove your claim .
how would you do this ?
presumably , you would first hold out an unbroken glass . then
, you would drop it , so that i could watch it fall , hit the floor , and break . in other words
, you would show me a transition from a state of the world in which the glass is whole to a one in which it is broken .
if you can convince me that it was indeed your dropping the glass that initiated this transition , then you have proven your causal claim .
what this little thought experiment shows is that the idea of a transition from one state of the world to the next is inherently part of the way in which we interpret causal statements .
however , structural models have nothing to do with such transitions . for instance , the causal claim about glasses breaking would just be represented by an equation @xmath12 , which determines two possible worlds : @xmath13 and @xmath14 . in this sense ,
structural models make complete abstraction of the dynamic aspects of causality , until all that remains is a static picture of how the values of different random variables can be defined in terms of each other .
other approaches to causality do not share this static worldview .
for instance , gives an explicitly dynamic account of causation .
he represents causal systems by means of probability trees , in which edges represent transitions between states of the world . for example , in the following picture , the edge going from @xmath15 to @xmath16 represents a transition from a state in which joe hasnt yet taken a swing at the ball to one in which he has and missed : @xmath17\\ & n_1 \ar[ld]_{\text{joe misses ball}}^{0.75 } \ar[rd]^{\text{joe hits ball}}_{0.25}\\ n_2 \ar[d ] & & n_3\ar[d]\\ \ldots & & \ldots } \end{split}\ ] ] the edge @xmath18 represents a transition of the same state @xmath15 to a state where joe has hit the ball .
together , these two edges represent a non - deterministic event , namely that of joe s taking a swing at the ball , which may result in one of these two outcomes .
the edges are labeled with the probabilities of the outcomes : the probability of joe s swing missing is @xmath19 and that of it hitting is @xmath20
.    this paper will use _ cp - logic _ as its formal language , which is essentially just a modular , syntactic representation for such shaferian probability trees . a theory in cp - logic
represents the causal structure of a domain by means of a set of _ causal probabilistic laws _ ( cp - laws , for short ) .
each such cp - law is a blue - print for a class of non - deterministic events .
for instance , the following cp - law : @xmath21 states that , for every player @xmath22 and ball @xmath23 , player @xmath22 s taking at swing at ball @xmath23 causes a non - deterministic event , which has as one possible outcome that @xmath22 hits @xmath23 ( and this happens with probability @xmath20 ) , and as its other possible outcome that @xmath22 misses @xmath23 ( which happens with probability @xmath19 ) .
if @xmath22 and @xmath23 are instantiated to , respectively , a particular player and a particular ball , say joe and the twelfth pitch , the we obtain a description of one particular event ( that may or may not happen , depending on whether joe decides to swing ) : @xmath24 this instantiated cp - law can be seen as a textual representation of picture , provided that , of course , node @xmath15 represents a state of the domain where joe has decided to take a swing at this particular pitch . in this way
, each instantiation of a cp - law describes a piece of probability tree .
as will be explained in more detail later , an entire cp - theory describes a class of probability trees , each of which can be constructed by putting these small pieces together .    unlike structural models
, the formal semantics of cp - logic therefore does provide mathematical objects that represent transitions between states of the world . as argued by ,
these transitions are important for a study of actual causation . indeed ,
when the goal is to figure out what caused what in a given _ story _ , it is obviously convenient to have a language whose formal semantics already offers objects that correspond in a natural way to stories .
a branch of a shaferian probability tree is precisely such an object , because , like a story , it is a description of a sequence of events that change the state of the world .
there is also a second argument in favour of cp - logic .
the goal of actual causation is to explain why things happened .
typically , though , not everything is in need of explanation . a detective solving a murder case , for instance , will be interested in why the victim is dead , but he wo nt care about why he was ever alive in the first place .
the detective s causal model will therefore list causes for dying ( poison , gun shot ,  ) , but not for living ( sexual intercourse , ivf ,  ) . in more technical terms , the detective considers living to be the _ default _ state of a person , and he is only interested in _ deviations _ from this default .
many authors , such as or , have argued that actual causation should be studied under the assumption that each rv has such a default state .
structural models make no distinction between the different values of a rv .
consequently , a rv @xmath25 with values @xmath26 and @xmath27 , and a rv @xmath28 with values @xmath27 and @xmath26 are completely interchangeable . in cp - logic , this is not the case . here , each rv ( i.e. , ground atom ) has * f * as its default value .
this means that the mere existence of the atom @xmath29 implies that the default condition is for @xmath30 to be dead , and that he can only come to life when there is a sufficient cause for this .
by contrast , the detective s theories will contain atoms such as @xmath31 , indicating that living is the default and death is in need of causal explanation . in the probability trees generated by a cp - theory ,
an atom always starts out at its default value , and only deviates from this when it has sufficient cause to do so .
lacking space for a full review of cp - logic , i will only summarize the main ideas and refer to @xcite for details .
the general form of a cp - law is : @xmath32 here , @xmath33 is a first - order formula and the @xmath34 are atoms , such that the tuple of variables @xmath1 contains all free variables in @xmath33 and the @xmath34 .
the @xmath35 are non - zero probabilities with @xmath36 .
such a cp - law expresses that @xmath33 causes some ( implicit ) non - deterministic event , of which each @xmath34 is a possible outcome with probability @xmath35 . if @xmath37 , then at least one of the possible effects @xmath34 must result if the event caused by @xmath33 happens ; otherwise , it is also possible that the event happens without any ( visible ) effect on the state of the world . for the purpose of this paper ,
the propositional fragment of cp - logic suffices , so i will from now on restrict attention to cp - laws in which the tuple of variables @xmath1 is empty .    for a cp - law @xmath38 , we refer to @xmath33 as the _ body _ of @xmath38 , and to the sequence @xmath39 as the _ head _ of @xmath38 .
we denote these objects as @xmath40 and @xmath41 , respectively , and also write @xmath42 for the set of all @xmath34 for which there exists an @xmath35 such that @xmath43 . for cp - laws
that are vacuously caused , @xmath40 may be omitted .
if a cp - law has a deterministic effect , i.e. , it is of the form @xmath44 , it is also written simply as @xmath45 .
a _ cp - theory _ is a finite set of cp - laws .
such a cp - theory describes the non - deterministic evolution of a domain , which is formally represented by a shaferian probability tree .
initially , all rvs of this domain ( i.e. , all ground atoms ) are in their default state .
this means that we can describe the initial state of the domain , which corresponds to the root of the probability tree , by the interpretation that assigns * f * to each of them .
we then extend this root by picking a cp - law @xmath38 whose precondition @xmath40 is satisfied according to this interpretation and creating a child node for each pair @xmath46 in @xmath41 .
the edge to child @xmath47 is labeled with the probability @xmath35 and the corresponding new state of the domain is constructed from the previous state by switching @xmath48 to its deviant state * t*. the cp - law @xmath38 has now happened , and will not happen again .
we repeat this process of adding children to one of the leaf nodes of the current tree , until this is no longer possible , i.e. , until for all leaves @xmath49 of the current tree it is the case that all rules @xmath38 that have not yet happened in @xmath49 have a precondition @xmath40 that is false in @xmath49 .
the resulting trees are called the _ execution models _ of the cp - theory . for a node @xmath50 of the tree ,
i denote by @xmath51 the interpretation that corresponds to the state of the world at that node , and , if @xmath50 is not a leaf , by @xmath52 the cp - law that was used to create the children of this node .
the construction of execution models is quite non - deterministic , in the sense that in any particular node of the tree , there can be many cp - laws that may be used to extend it .
the question is now whether each of these trees actually reflects a sensible way in which a domain described by the cp - theory might evolve .
the answer is a qualified `` yes '' , and depends on precisely how we choose to interpret negation appearing in the body of a cp - law .
consider the following example : @xmath53 we could take the body of this cp - law to mean that this transition may happen in any state where @xmath54 is still at its default state * f * , such as , by definition , the initial state in which suzy has no yet made up her mind about throwing . taking this view ,
every probability tree constructed according to the above principles can indeed be seen as a sensible description of how the domain might evolve . however , this is not very useful .
as argued by , it is more interesting to read negation in a slightly different way , namely , as not just saying that @xmath54 is still at its default value in the current state , but that it can actually _ never _ deviate any more . in other words , according to this reading , the above cp - law will only be applicable _ after _ suzy has decided that she will not refuse the throw .
this idea is formalized in the semantics of cp - logic by means of a construction similar to the gelfond - lifschitz reduct .
we use this to compute , for each state @xmath50 , an overestimate @xmath55 of all atoms that can still be caused in this @xmath50 .
only if an atom @xmath56 does not belong to @xmath57 , do we then say that @xmath58 holds in @xmath50 . if there are no loops containing double negation ( i.e. , some @xmath59 causing @xmath60 and @xmath61 causing @xmath62 ) , then it is the case that , in any branch of a probability tree , each cp - law must either happen at some point , or else become impossible . showed that there is a close connection between the resulting semantics and the well - founded model construction for a logic program
each probability tree @xmath63 defines , in the obvious way , a probability distribution @xmath64 over its leaves . for an execution model @xmath63 of a cp - theory @xmath65
, this distribution @xmath64 induces a probabilistic possible world semantics : the probability @xmath66 of an interpretation @xmath67 is @xmath68 , where the sum is taken over leaves @xmath49 of @xmath63
. showed that each execution model @xmath63 of a cp - theory @xmath65 defines the same possible world semantics @xmath64 .
for instance , the two trees shown on page are execution models of the same theory and , even though they are not isomorphic , they both define the same @xmath64 . in this way
, each cp - theory @xmath65 defines a unique probability distribution , which is denoted as @xmath69 .
the probability of a formula @xmath33 can then be defined as @xmath70 .
the fact that @xmath69 does not depend on the choice of any particular execution model @xmath63 may help to explain why structural models choose to ignore the dynamic aspects of causality in the first place .
indeed , this result shows precisely that , for applications which only care about properties of the final state that the domain will eventually reach , the details of how this final state came about can be safely ignored .
as i attempt to show in this paper , though , actual causation is _ not _ such an application .    like structural models
, cp - logic also makes a distinction between exogenous and endogenous random variables . with @xmath5 the set of all exogenous atoms ,
the semantics of a cp - theory now becomes relative to an interpretation @xmath71 for these atoms .
in particular , an execution model for @xmath65 _ given @xmath71 _ is defined as a execution model that starts not from a root in which _ all _ atoms are * f * , but instead starts with only the endogenous atoms being * f * and the exogenous atoms being interpreted by @xmath71 .
have shown that for each interpretation @xmath71 for the exogenous predicates of a cp - theory @xmath65 , all execution models @xmath63 given @xmath71 define the same probability distribution @xmath64 , which is denoted as @xmath72 .
this section briefly recalls the hp account . their paper starts with this example :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ suppose that two arsonists drop lit matches in different parts of a dry forest , and that both cause trees to start burning , until the entire forest burns down .
both matches are necessary to burn down the forest ; with only one match , the fire would die down . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    it is clear that both arsonists are an actual cause of the forest burning down .
hp reach this conclusion as follows . to represent the causal structure of the example , they use a structural model consisting of a single equation : @xmath73 the particular story under consideration
is then represented by the following assignment of values to the exogenous rvs : @xmath74 .
this of course also uniquely determines the values of the endogenous rvs : @xmath75 .
the hp definition is now reproduced below . in it
, @xmath76 is a structural model with endogenous rvs @xmath77 , @xmath78 an assignment of values to the exogenous rvs , @xmath0 a tuple of endogenous rvs , and @xmath33 a boolean formula in the rvs .
the notation @xmath79 \phi$ ] means that @xmath33 holds in @xmath80 after the intervention of assigning @xmath1 to @xmath0 is performed , i.e. , each @xmath81 has its defining equation removed from @xmath76 and replaced by @xmath82 .
[ hp ] @xmath83 is an actual cause of @xmath33 in @xmath84 if the following three conditions hold .
* @xmath85 ( that is , both @xmath86 and @xmath33 are true in the actual world . ) * there exists a partition @xmath87 of @xmath88 with @xmath89 and some setting @xmath90 of the variables in @xmath91 such that if @xmath92 for all @xmath93 , then both of the following conditions hold : * * @xmath94\lnot\phi$ ] . in words , changing @xmath91 from @xmath95 to @xmath96 changes @xmath33 from true to false . *
* @xmath97\phi$ ] for all subsets @xmath98 of @xmath99 and all subsets @xmath100 of @xmath101 . in words , setting any subset of variables in @xmath99 to their values in @xmath102 should have no effect on @xmath33 , as long as @xmath0 is kept at its current value @xmath1 , even if all the variables in an arbitrary subset of @xmath101 are set to their original values in the context @xmath78 .
* @xmath0 is minimal ; no subset of @xmath0 satisfies conditions ac1 and ac2 .
minimality ensures that only those elements of the conjunction @xmath103 that are essential for changing @xmath33 in ac2(a ) are considered part of a cause .    with @xmath104 , @xmath105 and @xmath106 ,
this definition provides the result that @xmath107 actually caused @xmath108 , since if we change @xmath5 to @xmath109 , while leaving @xmath110 at its original value ( this trivially satisfies ac2(b ) ) , we obtain @xmath111 as required by ac2(a ) .
in other words , in this example , we get actual causation from a simple counterfactual dependency : if it had nt been for @xmath107 , the forest would nt have burned down .
hp also consider a disjunctive variant of this example , where a single match already suffices to burn down the forest ( @xmath112 ) .
this causes the straightforward counterfactual criterion to fail , since stopping only one of the arsonists does not stop the forest burning down .
this motivates the additional machinery of the above definition . by considering the context
in which @xmath113 , that is @xmath110 and @xmath114 , we can re - establish the counterfactual dependency of @xmath115 on @xmath104 .
as we have seen , a question of actual causation can only be asked in the presence of two pieces of information : a causal model of a domain ( the @xmath76 of definition [ hp ] ) and a story that takes place in this domain ( the @xmath78 ) .
my definition will of course assume that the causal model is given in the form of a cp - theory @xmath65 . in the context of cp - logic ,
the most obvious formal counterpart of a `` story '' is a branch of an execution model of @xmath65 .
already , this allows us some more room for nuance than hp , as the following example from shows .    _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ suzy and billy might each decide to throw a rock at a bottle .
if suzy does so , her rock shatters the bottle with probability @xmath116 .
billy s aim is slightly worse and he only hits with probability @xmath117 . _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    this domain corresponds to the following set of cp - laws , where @xmath118 and @xmath119 are exogenous : @xmath120 assuming that suzy and billy both throw , there still exist two different execution models of the theory . representing the states in which the bottle is broken by an empty circle , and those in which it is still whole by a full one
, they look like this : + [ trees ] @h=0.8 cm & & ^0.9|(.3 ) _ 0.1|(.3 ) + & _ 0.8|(.3 ) ^0.2|(.3 ) & & _
0.8|(.3 ) ^0.2|(.3 ) + & & & & @h=0.8 cm & & ^0.8|(.3 ) _ 0.2|(.3 ) + & _ 0.9|(.3 ) ^0.1|(.3 ) & & _
0.9|(.3 ) ^0.1|(.3 ) + & & & & + in the left execution model , suzy s rock reaches the bottle before billy s does , whereas in the right one , it is billy s rock that gets there first . as discussed at the end of section [ sec : prel ] , this difference is irrelevant if we are only interested in the final outcomes that might be reached : the probability of the bottle shattering is @xmath121 in both models .
however , the difference becomes relevant when we want to judge actual causation .
indeed , in the left execution model , it is possible for suzy
s rock to _ actually _ break the bottle even though billy s also would have ( in particular , this happens in the leftmost branch of the tree ) .
according to the execution model on the right , however , this is impossible : here , suzy s rock can only actually break the bottle if billy s rock fails to do so .
goes on to consider the following story :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ suzy and billy both pick up rocks and throw them at a bottle .
rock get there first , shattering the bottle .
since both throws are perfectly accurate , billy s would have shattered the bottle had it not been preempted by suzy s throw . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    this story tells us precisely that we are in the leftmost branch of the left execution model above . hence , suzy s rock should be the actual cause of the bottle breaking , and not billy s . before showing how i reach this conclusion in the context of cp - logic ,
let me first remark that things are more difficult for the hp account .
their paper first tries the following straightforward structural model : @xmath122 here , there is no such thing as one execution in which suzy s rock reaches the bottle first and one in which billy s is first .
hall s story therefore seems to say nothing more than that all five rvs are @xmath123 , and the phrase `` suzy s rock gets there first '' contributes nothing .
of course , because it is precisely this phrase that determines which rock actually broke the bottle , this causal model does not work .
hp fix the problem by introducing two new random variables : @xmath124 ( `` billy s rock hits the ( unbroken ) bottle '' ) and @xmath125 .
the order in which the two rocks actually reach the bottle can then be encoded _ in the structure of the model _ : @xmath126 to me , this does not seem the right way to go .
the order in which the rocks arrive is a purely contingent matter , which belongs to the details of the particular story that is being told , and _ not _ to the general causal structure of the domain . saying that suzy s rock arrives before billy s should not be placed on the same level of causal discourse as the statement that throwing rocks at bottles causes them to break .
this is not just a matter of taste , but also has practical consequences .
if we would want to know whether suzy s rock would still have been the actual cause of the bottle breaking if billy s rock had gotten there first , then  in the hp account  we would not just have to look at a different story in the same domain , but we would have to change the structure of our causal model .
such hand - tailoring of the causal model to the question under consideration is undesirable , and , as i will now show , it is not needed in cp - logic .
my definition too will be heavily based on the intuition of counterfactual dependency from a cause @xmath65 to an effect @xmath127 .
therefore , i first formalize the following criterion :    if all events happen in the way they actually happened _ with the exception that @xmath65 is somehow prevented from occurring _ , then @xmath127 will no longer occur .
( * )    this requires some mathematical machinery .
first , we need to be able to fix the outcome of certain events . for a cp - law @xmath38 of the form @xmath128 , we write @xmath129 to denote the deterministic cp - law @xmath130 .
if we now have a branch @xmath23 that tells us what actually happened , then we can define as follows a theory that fixes the outcome of all events that happened to their actual outcome .
let @xmath131 be a branch of an execution model of a cp - theory @xmath132 .
we define @xmath133 as the union of two disjoint sets @xmath134 and @xmath135 , where @xmath134 contains all cp - laws from @xmath132 that did not happen in branch @xmath23 , i.e. , @xmath136 , and @xmath135 consists of all @xmath137 for which @xmath38 caused @xmath138 in @xmath23 , i.e. , @xmath139    we also need an antonymical transformation , which prevents some @xmath34 from occurring . for an @xmath38 of the same form as above , we write @xmath140 for : @xmath141 to prevent an atom @xmath138 entirely
, it now suffices to apply this transformation to all cp - laws that might cause it .
given a theory @xmath132 , we therefore define @xmath142 as : @xmath143 by combining this transformation with the previous one , we can now construct a theory @xmath144 which corresponds precisely to the counterfactual eventuality that everything happens precisely as it did in branch @xmath23 , with the exception that @xmath65 is somehow prevented from occurring .
i thus formalize the counterfactual criterion ( * ) , by expressing that , according to this new cp - theory @xmath144 , @xmath127 will not occur .
let @xmath131 be a branch of an execution model of a theory @xmath132 . for two atoms @xmath65 and
@xmath127 , such that both @xmath65 and @xmath127 hold in @xmath145 , we say there is a _ counterfactual dependency _ from @xmath65 to @xmath127 if @xmath146 where @xmath147 and , to cover the case where @xmath65 is exogenous , @xmath148 is @xmath149 .    here
, saying that @xmath150 is of course equivalent to @xmath127 being false in each leaf @xmath49 of each execution model of @xmath151 .    this intuition of counterfactual dependency forms the core of the concept of actual causation , but as discussed above , it is in itself not enough .
the additional aspect is the idea of _ relevance_. a causal model might make provisions for a large number of eventualities , many of which may not have been relevant in the actual course of events .
it is typical for judgments of actual causation that truly irrelevant causal mechanisms are ignored , even when they might appear to become relevant in a counterfactual context .
the typical case where this intuition manifests itself is when counterfactual dependencies are masked by _
redundant causation _ : there is some back - up mechanism waiting in the wings , which will ensures that the effect happens anyway , even if we preempt its actual cause .
the example of suzy and billy is a good illustration of this .
the reason why we nevertheless insist that suzy is the actual cause of the bottle shattering is precisely a criterion of relevance : because suzy s rock got to the bottle first , billy s was irrelevant , so we ignore it .    tried to formalize this same intuition by means of the concept of a _ causal beam _ , which is meant to encompass precisely the relevant parts of the causal model .
however , the formal details proved hard to get right , and the refinement that eventually became part of the hp definition seems to be a significant source of complexity , which considerably clouds the otherwise rather simple idea of counterfactual dependency . in the explicitly dynamic context of cp - logic , something much more simple is possible .
let us ask again why intuitions feels that billy s rock is irrelevant if suzy s rock gets to the bottle first and shatters it .
i suggest the blindingly obvious answer : it just got there _
too late_. by the time billy s rock reached the bottle , the damage was already done , the bottle lay in pieces , and there was nothing left to shatter .
in other words , one simply can not cause what is already the case .
my notion of relevance will comprise just this : whatever happened _ after _ the effect is irrelevant , and whatever happened _ on the way to _ the effect is counted as relevant .
of course , this is not yet a complete dichotomy , since it does not rule on the status of those events that did not happen at all . recall that if some cp - law does not happen in a particular branch , this means that , somewhere along the way
, its precondition must have become impossible .
whether an event that did not happen is considered relevant will depend on when its precondition became impossible : if this was _ before _ the effect arose , then it is relevant , otherwise not .
this leads to the following definition .
let @xmath131 be a branch of an execution model of a theory @xmath132 .
let @xmath65 and @xmath127 be two atoms that both hold in the final state @xmath152 of @xmath23 , i.e. , @xmath153 .
@xmath65 is an _ actual cause _ of @xmath127 in branch @xmath23 if @xmath154 with @xmath155 and @xmath156 , where @xmath151 is constructed as follows .
if @xmath157 is the smallest @xmath158 for which @xmath159 , then @xmath160 . in words , @xmath65 is an actual cause of @xmath127 if there is a counterfactual dependency from @xmath65 to @xmath127 , according to the theory @xmath161 that consists of both those events that happened before @xmath127 was caused , and those events that had already become impossible by then .
[ def : actc ]    it is quite easy to check whether this definition is satisfied : you look at the given branch , find the place where @xmath127 first appeared , discard all events that had not yet happened then but still were possible , and check whether the remaining theory exhibits a counterfactual dependency between @xmath65 and @xmath127 or not . to illustrate , consider again the leftmost branch @xmath162 of the left execution model for the billy and suzy example .
the bottle breaks in node @xmath163 , i.e , @xmath164 . therefore , @xmath165 and @xmath166 .
according to @xmath167 , there now is indeed a counterfactual dependency from @xmath118 to @xmath168 , so the first is an actual cause of the second . as this example illustrates , it is important that a branch @xmath169 of an execution model not only records the successive states @xmath170 of the domain , but also the events @xmath171 that caused each of the state transitions .
recall that the hp setting offers no mathematical objects that correspond to a complete story about what happened , so their definition is always just given the final outcome in the form of an assignment of values to the rvs .
in this case , we can not always say with certainty whether some potential cause actually caused an effect or not . indeed ,
if we get only the final interpretation @xmath145 instead of the full branch @xmath172 , then the best we can do is this :    [ partial ] let @xmath132 be a cp - theory and @xmath71 an interpretation for its vocabulary .
let @xmath173 be the set of all branches of all execution models of @xmath132 that end in a state @xmath50 for which @xmath174 .
if @xmath65 actually causes @xmath127 in _ at least one _ branch @xmath175 , we say that @xmath65 is a _
possible _ actual cause for @xmath127 .
if @xmath65 actually causes @xmath127 in _ all _ branches @xmath176 , we say that @xmath65 is a _
certain _ actual cause for @xmath127 .
if , in the bottle breaking example , we are only told that eventually @xmath118 , @xmath119 and @xmath168 all hold , we find ourselves faced with precisely the same problem as hp s first structural model : all that we can say is that both are possible actual causes , but neither is a certain actual cause .
this is typical for redundant causation patterns , and fits well with intuition here : without knowledge about the order in which events happened , we can not say which of the redundant causes actually `` got there first '' .
so far , we have only considered actual causation as it applies to atoms causing atoms . often , it is also interesting to wonder which omissions contributed to an effect ( `` did the doctor s failure to treat the patient cause his death ? '' ) or why some effect was in fact not caused ( `` did the doctor s treatment prevent the patient s death ? '' ) .
extending the framework to also address such questions is easy enough :    * to extend our definition of actual causation to allow also literals @xmath177 to act as effects , we need to specify when such a @xmath177 `` happens '' for the first time , such that we may discard all later events when making counterfactual judgments to determine what caused @xmath177 .
the obvious cut - off point is when @xmath127 no longer belongs to the overestimate @xmath178 . *
to also allow literals @xmath179 to act as causes , we need to define precisely how we will check the counterfactual dependency in this case . to assume that @xmath179 was not the case , we need to assume that @xmath65 has somehow occurred , which we can do formally by just adding a new cp - law `` @xmath180 '' that always causes @xmath65 .
due to space restrictions , formal details are left to the reader .
there is a large literature about actual causation , with many examples , counterexamples , and counter - counterexamples .
while e.g.   have argued that the importance of such small examples should not be exaggerated , it nevertheless remains useful to check that my approach behaves sensibly for them .
due to space restrictions , i will limit myself to those examples that most clearly illustrate the difference between my approach and the hp account .
it is common practice in research on actual causation to formulate examples in terms of neuron diagrams .
a neuron can be in one of two states , one is the default `` off '' state and the other is the deviant `` on '' state in which the neuron `` fires '' or `` is active '' .
different kinds of links between two nodes define how the state of one affects the other . for instance , in the following figure , @xmath127 fires if and only if @xmath181 fires , and @xmath181 fires if at least one of @xmath138 or @xmath65 fires .
[ fill](0,1)0.25 cm [ fill](1,1)0.25 cm [ fill](2,1)0.25 cm [ stroke](0,0)0.25 cm    neuron diagrams typically record not only this causal structure , but also the state of the neurons . in the figure above , nodes that are
`` on '' are represented by full circles and nodes that are `` off '' are shown as empty circles .
so , @xmath138 , @xmath181 and @xmath127 all fire , whereas @xmath65 does not . in the language that we have developed so far , a neuron diagram therefore places us in the partial information setting of definition [ partial ] : we are given a causal model of a domain together with the final state that has been reached , but are not told precisely how this state has come about .
shows a number of counterexamples to hp , and introduces an alternative account , which he formalizes for neuron diagrams only .
one of his counterexamples concerns the following two diagrams :    in both diagrams , the edges from @xmath181 to @xmath182 and from @xmath182 to @xmath127 are _ blocking _ edges : if @xmath181 fires , then @xmath182 will never fire , regardless of its other incoming edges . in the left diagram ,
both @xmath138 and @xmath65 cause @xmath127 : @xmath138 causes it directly and @xmath65 causes it by stopping @xmath183 from preventing @xmath127 . in the right diagram , however , @xmath65 also causes the very `` threat '' to @xmath127 that it prevents . therefore
, hall argues , in this diagram it should not be counted as a cause for @xmath127 .
the hp account correctly handles the left diagram , but fails for the right one , since taking @xmath184 and @xmath185 allows us to create the context @xmath186 in which there is a counterfactual dependency from @xmath65 to @xmath127 .    to see how my definition fares , here are the obvious cp - logic versions . in the first ,
@xmath187 and @xmath183 are all exogenous , while in the second only @xmath138 and @xmath65 are .
+    @xmath188    @xmath189    first , consider the left theory . here
, @xmath127 can only be caused after has already happened and both and have become impossible .
therefore , all these cp - laws are relevant and we end up having to check whether there is a counterfactual dependency from @xmath65 to @xmath127 in the original theory . clearly , this is the case , since no tree that starts from a root in which the exogenous predicates @xmath183 and @xmath138 are * t * and @xmath65 is * f * can produce @xmath127 . in the second theory
, the event may either happen before @xmath127 is caused or after .
this means we either have to check for a counterfactual dependency in the theory @xmath190 or in @xmath191 . in neither theory
we find a counterfactual dependency , so @xmath65 is correctly judged to certainly not be an actual cause of @xmath127 .
hp s problems with this example are caused by the fact that , lacking an explicitly dynamic semantics , they have to resort to interventions to eliminate irrelevant events from consideration . as an undesired side - effect
, they end up allowing the possibility that @xmath183 itself is relevant for judging the impact of @xmath65 on @xmath127 , but the link between @xmath65 and @xmath183 is not .
the following is an example of _ bogus prevention _
@xcite , taken from .
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ assassin is in possession of a lethal poison , but has a last minute change of heart and refrains from putting it in victim s coffee .
bodyguard puts antidote in the coffee , which would have neutralized the poison had there been any .
victim drinks the coffee and survives . _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    here , hp , as well as others such as , erroneously designate the bodyguard s unnecessary antidote as an actual cause for victim s survival .
as i will now show , my account handles this correctly . since the example states that assassin has his @xmath192 before the @xmath193 is administered
, i will not make these exogenous atoms , but instead include them as endogenous atoms that are vacuously caused with some unknown ( and irrelevant ) probability .
@xmath194 the example now tells the story that first event happens , which is then followed by .
however , as soon as happens , both @xmath195 and @xmath196 become impossible , so is considered irrelevant in the actual course of events and will not be part of the theory in which we check for a counterfactual dependency .
hence , preventing @xmath193 in this theory has no effect whatsoever upon the victim s survival , so it is not an actual cause of victim s survival ( but @xmath192 is ) .
note that if the antidote were administered before the assassin s change of heart , then it _ would _ be considered relevant , but still not an actual cause of victim s survival because then would no longer be relevant .
to make the antidote an actual cause of victim s survival , it would have to be administered after the assassin has _ failed _ to have a change of heart .
a prototype implementation can be downloaded from the following url :    http://people.cs.kuleuven.be/~joost.vennekens/actcaus/act.pl[http://people.cs.kuleuven.be/@xmath197joost.vennekens/actcaus/act.pl ]    this small program computes whether an atom is a possible / certain actual cause for an effect in the partial information setting , or an actual cause in the complete information setting .
it was written in swi - prolog , but should also run in sicstus or yap .
currently , it only handles ground theories without disjunction in rule bodies .    in the partial information
setting , this prototype performs a simple backtracking search over all branches that might generate the given observations . obviously , this is not an approach that would scale well for larger examples .
the goal of this prototype , however , is just to allow interested people to experiment with my definition , in order to see whether it corresponds to their intuition . as such ,
it is not meant to handle problems larger than the examples typically considered in the actual causation literature .
future work may investigate better algorithms , e.g. , by means of an integration into riguzzi s @xcite query answering algorithm .
this paper has tried to argue that the hp account of actual causation is flawed for two reasons , both of which stem from their choice of structural models as the formal language to express causal relations .
first , structural models fail to make the distinction between default and deviant values , which has been argued by many authors to play a key role in a correct understanding of actual causation .
second , the static world - view of structural models is ill - suited to handling dynamic concepts , such as the stories that are part of the input to an actual causation problem .
since the hp paper first appeared , it has received a great deal of attention among researchers interested in actual causation , and many counterexample and alternative approaches have been presented .
most of these , such as or , recognize the importance of the deviant / default distinction . the problems caused by the mismatch between the static formalism of structural equations and the dynamic problem of actual causation have achieved less attention , even though they are recently also pointed out by
. nevertheless , also these more recent approaches still use static formalisms such as neuron diagrams or variants of structural models . the main point
i hope to make in this paper is that for problems that , like actual causation , require reasoning about the way in which a domain evolves , it pays to have a language with a formal semantics that contains mathematical objects that correspond to such evolutions .
i have tried to illustrate this by defining a notion of actual causation in the context of cp - logic , a probabilistic logic programming language which can be seen as a modular syntactic representation for shaferian probability trees , which offer precisely the kind of dynamic representation that is perfectly suited for a study of actual causation .
my definition is based on a counterfactual criterion similar to hp s , but is able to leverage the dynamic nature of cp - logic s semantical objects to come up with a very straightforward notion of relevance , namely , it only considers as relevant those events that happened ( or became impossible ) _ before _ the effect first arose .
this is much simpler than the relevance criterion of hp , since i do not have to rely on complex manipulations by means of interventions .
while lacking space for an elaborate review of examples from the literature , i have shown that there are three examples where my definition beats hp : already for simple examples of redundant causation , it offers a more elegant account due to its ability to distinguish the complete and partial information settings ; it is also able to detect fake causes that simply prevent themselves from preventing the effect ; and it also handles bogus prevention .
of course , that is not to say my approach is perfect .
for instance , the railroad switch example from the hp paper can not be handled , because it contains a rv ( @xmath198 ) whose default and deviant values switch in the middle of the story .
i am also offering a prototype implementation of my definition , in the hope that it may help to find further examples where it does not correspond to intuition .
feedback will be appreciated .
inference with logic programs with annotated disjunctions under the well founded semantics . in _ logic programming , 24th international conference , iclp 2008 , udine , italy ,
december 9 - 13 , 2008 , proceedings_. lecture notes in computer science , vol .
springer , 667771 .