complex systems have attracted much attention among physicists , applied mathematicians , engineers , and social scientists in recent years . in particular ,
agent - based models have become an important part of research on complex adaptive systems @xcite .
for example , self - organized phenomena in an evolving population consisting of agents competing for a limited resource , have potential applications in areas such as engineering , economics , biology , and social sciences @xcite .
the bar - attendance problem proposed by arthur @xcite constitutes an everyday example of such a system in which a population of agents decide whether to go to a popular bar having limited seating capacity .
the agents are informed of the attendance in past weeks , and hence share common information , make decisions based on past experience , interact through their actions , and in turn generate this common information collectively .
these ingredients are key characteristics of complex systems @xcite .
an important step in the recent explosion of research in agent - based models within the physics community , has been the introduction of binary ising - like versions of models of competing populations
examples include the minority game ( mg ) @xcite and the binary - agent - resource ( b - a - r ) game @xcite .    for modest resource levels in which there are more losers than winners , the minority game @xcite represents a simple , yet highly non - trivial , model that captures many of the essential features of such a competing population .
it has been the subject of many theoretical studies @xcite .
the mg considers an odd number @xmath0 of agents . at each timestep
, the agents independently decide between two options ` @xmath3 ' and ` @xmath4 ' .
the winners are those who choose the minority option .
the agents learn from past experience by evaluating the performance of their strategies , where each strategy maps the available global information , i.e. the record of the most recent @xmath5 winning options , to an action .
one important quantity in the mg is the standard deviation @xmath6 of the number of agents making a particular choice .
this quantity reflects the performance of the population as a whole in that a small @xmath6 implies on average more winners per turn , and hence a higher success rate per turn among the agents . in the mg
, @xmath6 exhibits a non - monotonic dependence on the memory size @xmath5 of the agents @xcite .
when @xmath5 is small , there is a significant overlap between the agents strategies .
this crowd effect @xcite leads to a large @xmath6 , implying the number of losers is high .
this is the crowded , or informationally efficient , phase of mg . in the informationally inefficient phase where @xmath5 is large ,
@xmath6 is moderately small and the agents perform better than if they were to decide their actions randomly . in this regime ,
information is left in the resulting bit - string patterns for a single realization of the system . in the inefficient phase
, the mg can be mapped on to disordered spin systems and hence the machinery in statistical physics of disordered systems , most noticeably the replica trick , can be applied @xcite .
however , the replica trick becomes ineffective in the efficient phase .
the crowd - anticrowd theory @xcite gives a physically transparent , quantitative theory of the observed features of the mg in both the efficient and inefficient regimes .
the crowd - anticrowd theory is based on the fact that it is the difference in the numbers of agents playing a given strategy @xmath7 and the corresponding anti - correlated strategy @xmath8 , that dictates the size of the fluctuations and hence performance of the population as a whole .    in agent - based models ,
the non - triviality of the results comes from the actions taken by the agents which are directly related to the decision mechanism .
the decision mechanism depends sensitively on how each strategy performs at the moment of decision . in the efficient phase of the mg ,
no strategy outperforms the others and therefore the relative performance of the strategies oscillates as the game proceeds .
this anti - persistent nature of the strategy performance @xcite is crucial in arriving at a quantitative understanding of the mg s dynamics .
as the systems evolves , it goes from one @xmath5-bit outcome or history bit - string to another .
mathematically , the evolution can be viewed in terms of transitions in the global information ( i.e. history ) space .
the @xmath9 possible history bit - strings for a given value of @xmath5 constitute the nodes in this history space .
as the system evolves , it makes transitions from one node to another .
_ showed that in the efficient ( i.e. low @xmath5 ) phase of the mg , an effective restoring force dominates the strategy - score dynamics yielding a eulerian trail quasi - attractor in history space @xcite .
as @xmath5 increases , a competing bias term  associated with the initial strategy allocation  becomes increasingly important and eventually leads to instability of the eulerian trail quasi - attractor @xcite .
johnson _ et al . _
@xcite subsequently introduced and studied numerically what is known as the binary - agent - resource ( b - a - r ) model , in which the winning group is not necessarily decided by the minority rule .
the b - a - r model features a cutoff parameter @xmath2 , which is referred to as the global resource level @xmath2 ( @xmath10 ) .
the values of @xmath2 and @xmath0 are not known to the agents . at each timestep @xmath11 , each agent decides upon two possible options : whether to access the resource or not .
the winning action is decided by whether the number of agents attempting to access the resource , actually exceeds this resource level .
the mg therefore corresponds to the particular case of @xmath12 in the b - a - r model . in ref.@xcite
, it was found numerically that the population may unwittingly self - segregate itself into groups if @xmath2 deviates sufficiently from @xmath13 . as an example , consider a very high resource level of @xmath14 with each agent holding two strategies : it was found numerically that approximately @xmath15 agents are persistent winners while the rest are persistent losers .    in the present work
we analyze , both numerically and analytically , the transition from modest to high resources in this generic b - a - r model of a competing multi - agent population .
surprisingly , we find that the system exhibits a set of abrupt transitions between distinct yet well - defined states as the resource level @xmath2 increases from @xmath13 to @xmath0 .
in particular , both the highest success rate @xmath16 and the mean success rate @xmath17 among the agents , show abrupt transitions as @xmath2 varies .
in addition , @xmath16 exhibits _ fractional _ values in the plateau regions between each transition .
we show that this behavior can be understood in terms of the system s trajectory in the history space as time evolves . in particular , as @xmath2 increases , the portion of the history space that the system visits becomes increasingly restricted .
we derive analytic expressions for the observed plateaux values and the values of @xmath2 at which the transitions occur .
although the present analysis focuses on a non - networked population , the same elements of ( i ) strategy performance over time , and hence the dynamics of strategy scores , plus ( ii ) the system s trajectory in history space , provide the foundation for a quantitative understanding of a large class of agent - based models , including networked populations @xcite .
we also note that for physicists interested in random walks , the present b - a - r system provides a fascinating laboratory for studying correlated , non - markovian diffusion on a non - trivial network ( i.e. the history space , which corresponds to a de bruijn graph ) .
this non - trivial diffusion is in turn strongly coupled to the non - random temporal patterns arising in the strategy - performance dynamics . finally , we note in passing that the occurrence of abrupt transitions between plateau states , and stable fractions , are known to arise in a multi - electron quantum system as the external magnetic field is increased monotonically ( i.e. fractional quantum hall effect @xcite ) .
however it is very curious to see such ` quantized ' phenomena arise in a classical multi - particle system as a function of a monotonically increasing external control parameter .
the plan of the paper is as follows . in sec .
ii , we define the b - a - r model . in sec .
iii , we present numerical results from extensive simulations .
in particular , we demonstrate the existence of different phases or states at different values of the resource level . in sec .
iv , we discuss how the strategies performance evolves as the system evolves . in sec .
v , we discuss the history space of the b - a - r model and present results for the statistics in the outcome bit - strings at high resource levels . in sec .
vi , we explain the observed numerical features based on the idea that at high resource levels , the system restricts itself to only visit a restricted portion of the history space .
we derive an expression for the highest success rate @xmath16 among the agents and discuss the critical values of resource level at which transitions occur .
we summarize the results in sec .
vii , together with a discussion of how the present approach can be generalized to a wider class of agent - based models .
we consider the binary - agent - resource ( b - a - r ) model @xcite .
the b - a - r model is a binary version of arthur s el farol bar attendance model @xcite , in which a population of agents repeatedly decide whether to go to a bar with limited seating based on the information of the crowd size in recent weeks . in the b - a - r model
, there is a global resource level @xmath2 which is not announced to the agents , where @xmath0 is the total number of agents . at each timestep @xmath11
, each agent decides upon two possible options : whether to access resource @xmath2 ( action ` 1 ' ) or not ( action ` 0 ' ) .
the two global outcomes at each timestep , ` resource over - used ' and ` resource not over - used ' , are denoted by ` @xmath3 ' and ` @xmath4 ' . if the number of agents @xmath18 choosing action @xmath4 exceeds @xmath2 ( i.e. resource over - used and hence global outcome ` @xmath3 ' ) then the @xmath19 abstaining agents win .
by contrast if @xmath20 ( i.e. resource not over - used and hence global outcome ` @xmath4 ' ) then the @xmath18 agents win . in order to investigate the behavior of the system as @xmath2 changes , it is sufficient to study the range @xmath21 .
the results for the range @xmath22 can be obtained from those in the present work by suitably interchanging the role of ` 0 ' and ` 1 ' @xcite . in the special case of @xmath23 ,
the b - a - r model reduces to the minority game .    in the b - a - r model
, each agent shares a common knowledge of the past history of the most recent @xmath5 outcomes , i.e. the winning option in the most recent @xmath5 timesteps .
the full strategy space thus consists of @xmath24 strategies , as in the mg .
initially , each agent randomly picks @xmath25 strategies from the pool of strategies , with repetitions allowed .
the agents use these strategies throughout the game . at each timestep
, each agent uses his momentarily best performing strategy with the highest virtual points .
the virtual points for each strategy indicate the cumulative performance of that strategy : at each timestep , one virtual point ( vp ) is awarded ( deducted ) to ( from ) a strategy that would have predicted the correct ( incorrect ) outcome after all decisions have been made . a random coin - toss is used to break ties between strategies . in the b - a - r model , the population may or may not contain network connections . in the case of a networked population
@xcite each agent has access to additional information from his connected neighbors , such as his neighbors strategies and/or performance . in the present work ,
we focus on the b - a - r model with a non - networked population .    to evaluate the performance of an agent , one ( real ) point is awarded to each winning agent at a given timestep .
a maximum of @xmath2 points per turn can therefore be awarded to the agents .
an agent has a success rate @xmath26 , which is the mean number of points awarded to the agent per turn over a long time window . the mean success rate @xmath27 among the agents is then defined to be the mean number of points awarded per agent per turn , i.e. an average of @xmath26 over the agents .
we are interested in investigating the details of how the success rate changes as the resource level @xmath2 varies in the efficient phase , where the number of strategies ( repetitions counted ) in play is larger than the total number of distinct strategies in the strategy space .
the effects of varying @xmath2 were first reported by johnson _ et al .
these authors studied numerically the dependence of the fluctuations in the number of agents taking a particular option , on the memory size @xmath5 for different values of @xmath2 . for the mg ( i.e. @xmath23 ) in the efficient phase ( i.e. small values of @xmath5 ) the number of agents making a particular choice varies from timestep to timestep , with additional stochasticity introduced via the random tie - breaking process .
the corresponding period depends on the memory length @xmath5 .
the underlying reason is that in the efficient phase for @xmath12 , no strategy is better overall than any other .
hence there is a tendency for the system to restore itself after a finite number of timesteps , thereby preventing a given strategy s vps from running away from the others . as a result ,
the outcome bit - string shows the feature of anti - persistency or double periodicity @xcite . since a maximum of @xmath12 points
can be awarded per turn , the mean success rate @xmath27 over a sufficiently large number of timesteps is bound from above by @xmath28 .    in the b - a - r system with high resource level
, the mean success rate behaves differently . taking the extreme case of @xmath29 ,
the winning action is obviously ` 1 ' ( i.e. access resource ) and in principle every agent could win in every timestep .
the history upon which the agents decide , is persistently @xmath5-bits of ` 1 ' .
however , due to the random initial strategy distribution , some agents may not hold a strategy that predicts the winning option for a history of @xmath5
@xmath4 s .
therefore , there are still losers and @xmath27 is less than @xmath2 .
the number of losers depends on @xmath25 , the number of strategies that each agent holds . for @xmath30 and assuming that the strategies are picked randomly , a mean number of @xmath31 agents in a large population will hold two strategies both predicting the wrong option . the mean success rate is thus given by @xmath32 .
we have carried out extensive numerical simulations on the b - a - r model to investigate the dependence of the success rate on @xmath2 for @xmath21 . unless stated otherwise , we consider systems with @xmath33 agents and @xmath30 .
figure  [ fig : figure1](a ) shows the results of the mean success rate ( dark solid line ) as a function of @xmath2 in a typical run for @xmath34 , together with the range corresponding to one standard deviation about @xmath17 in the success rates among the @xmath0 agents ( dotted lines ) and the spread in the success rates given by the highest and the lowest success rates ( thin solid lines ) in the population . by taking a larger value of @xmath0 than most studies in the literature
, we can analyze the dependence on @xmath2 and @xmath5 in great detail .
in particular , these quantities all exhibit abrupt transitions ( i.e. jumps ) at particular values of @xmath2 . between the jumps , the quantities remain essentially constant and hence form steps or ` plateaux ' .
we refer to these different plateaux as states or phases , since it turns out that the jump occurs when the system makes a transition from one type of state characterizing the outcome bit - string to another . for different runs ,
the results are almost identical . at most , there are tiny shifts in the @xmath2 values at which jumps arise due to ( i ) different initial strategy distributions among the agents in different runs , and ( ii ) different random initial history bit - strings used to start the runs .
these different states are most clearly seen by monitoring the _ highest _ success rate @xmath16 among the agents for given values of @xmath2 and @xmath5 . the most striking feature in fig.[fig : figure1](a )
is that the values of the plateaux in @xmath16 are given by _ simple fractions _ , e.g. @xmath35 , @xmath36 , @xmath37 , @xmath38 , @xmath39 , etc .
this feature strongly suggests that the system goes through different states with different ratios of ` @xmath4 ' and ` @xmath3 ' in the outcome bit - string as @xmath2 varies , as will be discussed in later sections .
figure  [ fig : figure1](b ) shows that the features in the success rates for the simpler case of @xmath40 are similar to those in fig.[fig : figure1](a ) , except that the plateaux in @xmath16 take on fewer values , i.e. @xmath4 , @xmath41 , @xmath39 as @xmath2 decreases .
these values are closely related to the statistics in the outcome bit - string . for large @xmath0 and @xmath40 ,
the outcome bit - string shows a period of 4 bits . for values of @xmath2 with @xmath42
, it turns out that the fraction of the outcome ` @xmath4 ' in a period is exactly @xmath39 . for the range of @xmath2 corresponding to @xmath43 , there are three @xmath4 s in a period of @xmath44 , and so on . for @xmath34 ,
we have also carried out detailed analysis of the outcome bit - string .
for later discussions , we summarize in table [ tab : otsm3s2 ] the values of @xmath16 , the range of @xmath2 corresponding to the observed value of @xmath16 , the ratio of number of occurrence of ` 1'-bits to ` 0'-bits and the period in the outcome bit - string , as obtained numerically from the data shown in fig.[fig : figure1 ] .
hereafter , @xmath16 is used to label the state at a given @xmath2 .
the agents decide based on the best performing strategy that they hold at the moment of their decision .
a strategy s performance is evaluated by its virtual points , which vary as the game proceeds .
it is most illustrative to consider the case of @xmath40 and @xmath30 ( see fig.[fig : figure1](b ) ) since one can readily follow the dynamics for different values of @xmath2 . for @xmath40 , there are only four different strategies in the whole strategy space
. these strategies can be represented by ( @xmath45 ) , ( @xmath46 ) , ( @xmath47 ) , and ( @xmath48 ) , with the first ( second ) index in ( @xmath49 ) giving the action for the history bit - string of ` @xmath3 ' and ` @xmath4 ' , respectively .
the virtual points ( vp ) of the four strategies can then be represented by a @xmath50 matrix vp@xmath51 , with the @xmath49-element giving the vp of the strategy ( @xmath49 ) .
table  [ tab : strgevom1s2 ] shows the time evolution of the vps of the strategies in a few timesteps , the number of agents @xmath52 taking the action ` 1 ' , and the outcome for @xmath40 in a population with _
uniform _ initial distribution of all possible pairs of strategies to the agents .
we illustrate these ideas by considering the case of an initial history of ` 0 ' for given @xmath2 , but the same results are obtained for an initial history of ` 1 ' . for @xmath53 , the system follows the dynamics shown in the left column . in this range of @xmath2 ,
the vps can not run away due to the decision making process of the agents .
the vps restore their values in a few timesteps , as in the mg . for systems of large @xmath0
, the outcome series shows a 4-bit periodic pattern of @xmath54 , with @xmath55 in agreement with the numerical results in fig .
[ fig : figure1](b ) .
this highest success rate is achieved by the agents who hold two identical strategies .
the result also implies that an agent s success rate is determined by the hamming distance between the two strategies that the agent holds @xcite .
one could therefore say that different ` species ' of agent emerge in the population due to the dynamics of the system , with each species characterized by its own hamming distance ` gene ' .    for @xmath56 , two of the four strategies will have runaway vps , with one tending to increase without bound while another tends to decrease without bound ( see central column in table  [ tab : strgevom1s2 ] ) . following the dynamics
, the outcome series shows a 4-bit periodic pattern of @xmath57 , in the limit of large @xmath0 . in this range , @xmath43
as observed numerically and this value is achieved by those agents holding the strategy whose vps increase without bound . for @xmath58 ,
two strategies have their vps increasing ( decreasing ) without bound ( see right column in table  [ tab : strgevom1s2 ] ) .
the outcome series is persistently ` @xmath4 ' , i.e. it becomes effectively period-1 . in this high resource level regime , @xmath32 since three - quarters of the agents
hold at least one strategy which predicts the persistently winning option . in this way , it is possible to follow the dynamics of the system and obtain the number of possible states and the range of @xmath2 for each state .
table  [ tab : wrm1s2 ] summarizes the theoretical results for @xmath40 and @xmath30 , by following the analysis on the dynamics as shown in table  [ tab : strgevom1s2 ] .
the results for the location of the transitions , the values of @xmath27 and @xmath16 are all in good agreement with numerical data ( see fig .
[ fig : figure1](b ) ) .
the important point is that for small @xmath5 , the system undergoes several changes of state with successively higher values of @xmath27 and @xmath16 as @xmath2 increases .
the value of @xmath16 is related to the ratio of occurrences of the two possible outcomes in the outcome series , which in turn is related to the strategies performance and thus the decision - making process for a given value of @xmath2 .    we have carried out similar analysis for the case of @xmath59 and @xmath34 . for @xmath34 ,
the fraction of ` @xmath4 ' in a period is found to take on values in the set @xmath60 , which coincides with the numerical values of @xmath16 obtained by numerical simulations ( see fig .  [
fig : figure1](a ) ) and shown in table  [ tab : otsm3s2 ] .
however , the analysis becomes increasingly complicated for higher values of @xmath5 and/or @xmath25 .
the reason is that there are @xmath61 strategies , and allowing @xmath25 strategies per agent divides the agents into @xmath62 groups according to the sequence in which an agent picks his @xmath25 strategies .
this number increases rapidly with @xmath5 and @xmath25 , and the above microscopic analysis becomes hard to implement .
it turns out the states are closely related to the way in which the system explores the possible histories . in what follows
, we analyze the b - a - r model by an approach that focuses on the transitions between history bit - strings and hence on the path in the history space @xcite .
our approach couples together consideration of the probability of the occurrence of various histories and the ranking in the performance of the strategies @xcite .
the history space consists of all the possible history bit - strings for a given value of @xmath5 . for @xmath34 ,
it includes @xmath63 bit - strings of @xmath3 s and @xmath4 s .
figure  [ fig : figure2](a ) shows the history space for @xmath34 , together with the possible transitions from one history to another .
each history constitutes a node in the history space .
the transitions are marked by arrows together with the outcome necessary for making the transitions .
it will prove convenient to group the possible history bit - strings for a given @xmath5 into columns , in the way shown in fig .
[ fig : figure2](a ) .
each column is labelled by a parameter @xmath64 , which is the number of
s in the @xmath65-bit history ( histories ) concerned . one immediate advantage of this labelling scheme is that the different states characterized by @xmath16 turn out to involve paths in a _ restricted _ portion of the full history space .
for example , the state with @xmath66 is restricted to the @xmath67 portion of the history space , i.e. the @xmath68 history bit - string leads to an outcome of ` 1 ' and hence persistent self - looping at the node @xmath68 in history space .
the states with @xmath69 , @xmath36 , and @xmath37 correspond to different paths in the history space restricted to the @xmath70 and @xmath71 groups of histories , as shown in fig .
[ fig : figure2](b ) .
the states with @xmath72 and @xmath38 have paths extended to include @xmath73 histories .
the state with @xmath74 has paths that cover the whole history space ( @xmath75 ) . in general ,
the deviation of @xmath2 from @xmath13 acts like a driving force in the history space that drifts the system towards an increasingly restrictive portion of the history space bounded by a smaller value of @xmath64 .
one can also view this behavior as the system , in response to the global resource level @xmath2 , effectively avoiding certain nodes in the history space and hence avoiding particular patterns of historical outcomes .
as the game proceeds , the system evolves from one history bit - string to another .
this can be regarded as transitions between different nodes ( i.e. different histories ) in the history space . for @xmath12 in the efficient phase
, it has been shown @xcite that the conditional probability of an outcome of , say , ` @xmath4 ' following a given history is the same for all histories . for @xmath76 ,
the result still holds for states characterized by @xmath74 . note that a history bit - string can only make transitions to history bit - strings that differ by the most recent outcome , e.g. 111 can only be make transitions to either 110 or 111 , and thus many transitions between two nodes in the history space are forbidden .
in addition , these allowed transitions do not in general occur with equal probabilities .
this leads to specific outcome ( and history ) bit - string statistics for a state characterized by @xmath16 .
we have carried out detailed analysis of the statistics of the outcomes following a given history bit - string for @xmath34 , and for each of the possible states over the whole range of @xmath2 .
table  [ tab : tmm3 ] gives the relative numbers of occurrences of each outcome for every history bit - string .
for the state with @xmath77 , for example , the outcomes ` @xmath3 ' and ` @xmath4 ' occur with equal probability for every history bit - string , as in the mg .
for the other states , the results reveal several striking features .
it turns out that @xmath16 is given by the _ relative frequency _ of an outcome of ` @xmath4 ' in the outcome bit - string , which in turn is governed by the resource level @xmath2 . for example
, a ` 0 ' to ` 1 ' ratio of @xmath78 in the outcome bit - strings corresponds to the state with @xmath79 . in table
[ tab : tmm3 ] , we have intentionally grouped the history bit - strings into rows according to the label @xmath64 in fig .  [
fig : figure2 ] .
we immediately notice that for every possible state in the b - a - r model , the relative frequency of each outcome is a property of the _ group _ of histories having the same label @xmath64 rather than the individual history bit - string , i.e. all histories in a group have the same relative fraction of a given outcome .
this observation is important in understanding the dynamics in the history space for different states in that it is no longer necessary to consider each of the @xmath9 history bit - strings in the history space .
instead , it is sufficient to consider the four groups of histories ( for @xmath34 ) as shown in fig .  [
fig : figure2](a ) .
analysis of results for higher values of @xmath5 show the same feature .    for the state characterized by @xmath66 , the outcome bit - string is persistently `
@xmath4 ' and the path in the history space is repeatedly 111@xmath801 .
therefore , the path is restricted to the history labelled by @xmath70 and simply corresponds to an infinite number of loops around the history node 111 .
since there is no ` @xmath3 ' in the outcome bit - string , we will also refer to this state as @xmath81 state .
the system is effectively frozen into one node in the history space . in this case
, there are effectively only two kinds of strategies , which differ by their predictions for the history 111 .
the difference in predictions for the other ( @xmath82 ) history bit - strings become irrelevant . obviously , the ranking in the performance of the two effective groups of strategies is such that the group of strategies that suggest an action ` 1 ' for the history ` 111 ' , outperforms the group that suggests an action ` 0 ' . for a uniform initial distribution of strategies ,
there are @xmath83 agents taking the action 0 and @xmath84 agents taking the action ` 1 ' , since half of the strategies predict 0 and half of them predict 1 . to sustain a winning outcome of ` 1 '
, the criterion is that the resource level @xmath2 should be higher than the number of agents taking the action ` 1 ' .
therefore , we have for the state with @xmath85 that @xmath86 and @xmath87 these results are in agreement with the results obtained by numerical simulations . for @xmath30 , @xmath88 for @xmath58 .
note that equations   and are valid for _ any _ values of @xmath5 .
table  [ tab : tmm3 ] shows that the states with @xmath89 , @xmath36 , @xmath35 have very similar features in terms of the bit - string statistics .
they differ only in the frequency of giving an outcome of @xmath4 following the history of @xmath68 .
note that the @xmath90 and @xmath91 histories do not occur .
the results imply that as the system evolves , the path in history space for these states is restricted to the two groups of histories labelled by @xmath67 and @xmath71 .
the statistics show that the outcome bit - strings for the states with @xmath92 , @xmath36 and @xmath35 exhibit only one @xmath3-bit in a period of @xmath93 , @xmath94 and @xmath95 bits , respectively .
we refer to these states collectively as @xmath96 states , since the portion of allowed history space is bounded by the @xmath71 histories . graphically , the path in history space consists of a few self - loops at the node 111 , i.e. from 111 to 111 , then passing through the @xmath71 group of histories once and back to 111 , as shown in fig .
[ fig : figure2](b ) .
the states with @xmath42 , @xmath38 , @xmath97 involve the other groups of histories and exhibit complicated looping among the histories .
we refer to them collectively as higher ( i.e. @xmath98 ) states .
we now proceed to derive an expression for the observed value of @xmath16 for the @xmath96 states .
recall that each strategy consists of a prediction or action for all of the @xmath5-bit histories .
an important idea is that for states corresponding to paths restricted to a certain portion of the history space , only that part of a strategy corresponding to the histories in question is being used in making decisions .
strategies that only differ in their predictions for the history bit - strings which _ do not _ occur ( i.e. the avoided histories ) are now effectively identical . in the context of the crowd - anticrowd theory @xcite , two previously uncorrelated strategies
could now be correlated when viewed within this restricted history subspace .
as @xmath2 decreases , the system is allowed to explore a larger portion of the history space .
the ranking in the strategies performance becomes more complicated . for the @xmath96 states ,
the paths in history space involve the ( @xmath99 ) histories labelled by @xmath1000 and 1 ( see fig .
[ fig : figure2](b ) ) . for @xmath34 ,
only four out of a total of @xmath101 entries in a strategy corresponding to the histories ` @xmath68 ' , ` @xmath102 ' , ` @xmath103 ' , and ` @xmath104 ' now matter . for general @xmath5 , there are @xmath5 histories with one ` 0'-bit : hence a total of ( @xmath99 ) entries in each strategy now matter . for later discussions , it is useful to first classify all strategies into two groups according to their prediction for the history ` 11  1 ' ( @xmath70 history ) .
these strategies can further be classified according to their @xmath5 predictions for the @xmath71 histories . a strategy having @xmath105 bits predicting 1 and @xmath106 bits predicting 0 for the @xmath5 histories belonging to @xmath71 ,
can be labelled as ( @xmath107 ; @xmath105 , @xmath106 ) where @xmath108 is the prediction for the @xmath70 history .
the first three columns in table  [ tab : strategies10 ] show this classification of strategies for general values of @xmath5 describing the @xmath96 states .    for @xmath109 ,
the outcomes are no longer persistently ` @xmath4 ' , and the system explores both the @xmath70 and @xmath71 groups of histories .
assume that for @xmath2 just below @xmath15 , the outcome must be ` @xmath4 ' for the @xmath71 histories , i.e. the system only visits the @xmath71 histories once in a cycle .
consider the @xmath34 case , for example .
as time evolves , the system exhibits periodic visits to the histories . in each period
, each history in @xmath110 occurs once and the history in @xmath111 occurs @xmath112 times . among these @xmath112 occurrences of the @xmath70 history ,
the outcomes are ` @xmath4 ' for @xmath113 timesteps and ` @xmath3 ' for one timestep ( fig .
[ fig : figure2](b ) ) , since the system must go from ` @xmath68 ' to ` @xmath103 ' after @xmath113 loops in order to sustain the path .
it turns out that paths in the history space for the @xmath89 , @xmath36 , and @xmath35 states ( @xmath34 ) correspond to that shown in fig .
[ fig : figure2](b ) with @xmath114 , @xmath65 , and @xmath44 loops at the @xmath70 history .
since only @xmath99 different histories are involved , the performance of a strategy depends only on the predictions for this subset of histories .
consider the path in fig .
[ fig : figure2](b ) .
the strategy labelled by @xmath115 predicts the correct outcome @xmath116 times for @xmath117 and @xmath118 times for @xmath119 , respectively , in going through the path once . in table
[ tab : strategies10 ] , we list the performance of the strategies labelled by @xmath120 according to the number of successful predictions ( which reflects the vps ) in a closed path ( see fig .  [
fig : figure2](b ) ) consisting of @xmath113 loops at @xmath70 history , i.e. a total of @xmath121 timesteps from @xmath122 to @xmath123 .
it is important to note that there may be overlaps in strategies performance between the @xmath117 and @xmath119 groups of strategies for small values of @xmath113 , i.e. strategies with the label @xmath117 and @xmath119 may win the same number of timesteps in a cycle and hence belong to the same rank in performance of the strategies .
for example for @xmath34 and @xmath114 , strategies labelled by @xmath124 and @xmath125 belong to the same rank in performance .
the number of turns @xmath113 around the @xmath70 history which is consistent with the condition of @xmath126 , is restricted to the range @xmath127 .
this criteria is related to the number @xmath128 of overlapping performances between the @xmath129 and @xmath130 groups of strategies ( see table  [ tab : strategies10 ] ) .
note that for @xmath131 , we have @xmath132 , i.e. strategies with @xmath117 do not have overlapping vps with strategies with @xmath119 .
this implies that the vps of the strategies predicting ` @xmath4 ' for the history in @xmath133 , are always higher than those predicting ` @xmath3 ' .
this further implies that agents will take action ` @xmath4 ' for the @xmath70 history if one of their strategies belongs to the @xmath117 category . for @xmath30
, there will then be @xmath15 agents taking the action ` @xmath4 ' . since @xmath134 ,
the outcome must be ` @xmath3 ' .
therefore , the upper bound for the number of self - loops is @xmath99 , and thus @xmath135 . in other words , for paths with @xmath131
the system must have @xmath58 and the path in history space will be that of an infinite number of loops around @xmath70 history , i.e. the @xmath81 state . for @xmath136 , strategies in @xmath129 will perform worse than or equally well to those in @xmath130 for @xmath137 . as @xmath138 , strategies that predict more @xmath4 s should perform better .
therefore , @xmath136 leads to inconsistency and thus @xmath139 .
thus @xmath127 for the @xmath96 states , with the corresponding @xmath128 of overlapping groups of strategies being in the range @xmath140 .
since each possible allowed value of @xmath113 or @xmath128 gives one @xmath96 state , there are altogether @xmath5 possible @xmath96 states for a given @xmath5 .
the values of @xmath16 for the @xmath96 states can be readily found . for a given value of @xmath113 ,
the best performance among the strategies is to have @xmath141 correct predictions in a path consisting of ( @xmath142 ) timesteps with @xmath113 loops at the @xmath70 history .
therefore , for the @xmath96 states @xmath143 for @xmath34 , we have @xmath144 and hence @xmath145 . there are three @xmath96 states with @xmath89 , @xmath36 , and @xmath35 , exactly as observed in the numerical simulations .
we now derive the critical values of the resource level at which transitions occur from one value of @xmath16 to another for the @xmath96 states .
note that the transition from the @xmath66 state to the @xmath146 state for @xmath34 and @xmath30 , occurs at @xmath147 as predicted in eq.(2 ) .
the condition for transitions from @xmath69 to @xmath148 state is that the value of @xmath2 can no longer support @xmath149 loops at the @xmath133 history before giving an outcome of ` 0 ' for the history 111 . from table  [
tab : strategies10 ] , we note that the performance of the @xmath117 group of strategies becomes increasingly better than the @xmath119 group as @xmath113 increases , and the number of agents taking the action ` 1 ' increases towards @xmath15 .
therefore , the highest number of agents who take the action ` 1 ' _ and _ win will arise at the last turn among the @xmath113 loops where the history 111 is followed by an outcome ` 1 ' . similarly , the lowest number of winning agents will be at the turn when the history 111 is followed by an outcome ` 0 ' , i.e. breaking away from the @xmath113 loops at 111 .    the number of agents choosing the action ` 1 ' given the history 111 , is related to the number @xmath128 ( and hence @xmath113 ) of overlapping performances among the @xmath129 and @xmath130 strategies . from table  [ tab :
strategies10 ] , the number of correct predictions @xmath150 of a strategy @xmath120 in going through a path with @xmath113 loops at the @xmath70 history , is given by @xmath151    for general values of @xmath5 , there are only @xmath152 histories which matter for the @xmath96 states .
therefore there are @xmath153 _ effectively different _ strategies , each of which represents a group of @xmath154 strategies .
the number of effectively different strategies predicting @xmath107 for the @xmath70 history and having @xmath105 predictions of ` 1 ' for the @xmath5 @xmath71 histories , is given by @xmath155 for both @xmath156 and @xmath3 , where @xmath157 is the binomial coefficient .
the performance of the strategies can be ranked by a label @xmath158 , with @xmath159 and @xmath160 representing the best performing group of strategies .
a general situation in which there are @xmath128 overlapping performances between the @xmath129 strategies and @xmath130 strategies , is shown in table  [ tab : vprank10 ] .
for the @xmath96 states , the best performing strategies belong to the @xmath117 group and the worst performing ones belong to the @xmath119 group , with @xmath128 overlapping rankings in between where the allowed range of @xmath128 is @xmath140 .
the ranking @xmath158 of the strategies @xmath120 is related to @xmath105 by the simple relation @xmath161 for a given value of @xmath128 , there are a total of @xmath162 ranks .
therefore , the label @xmath158 is restricted to the range @xmath163 with @xmath164 given by @xmath165 as the number of loops @xmath113 increases , @xmath128 decreases and the strategies spread more widely in terms of performance .
it follows from eqs .   and that the number of effectively different strategies @xmath166 in rank-@xmath158 is given by @xmath167 the fraction of rank-@xmath158 strategies among all the strategies is then given by @xmath168    on the critical turn that determines the minimum value of @xmath2 for sustaining a certain @xmath96 state , only the @xmath117 strategies win . using table  [ tab : vprank10 ] together with eqs .   and
, the number of winning strategies @xmath169 is @xmath170 the corresponding fraction @xmath171 of winning strategies among all strategies of rank - r is given by @xmath172    each agent uses the strategy in his possession which has the best performance record , i.e. the one having the ranking with smaller @xmath158 , in order to make a decision .
assuming a uniform initial distribution of any combination of @xmath25 strategies ( with repetitions allowed ) among the agents , the fraction of agents _ holding _ a rank-@xmath158 strategy as their best performing strategy @xmath173 is @xmath174^s    -\bigg[\sum_{r'=r+1}^{r_{max}}\widetilde{c}(r')\bigg]^s.\end{aligned}\ ] ] with @xmath175 given by eq
..    each @xmath96 state corresponds to a specific value of allowed @xmath113 and hence an allowed @xmath128 . for a given @xmath113 or @xmath128 , the resource level @xmath2 needed to accommodate all the agents that take the action ` 1 ' for the @xmath70 history ,
gives the criterion for the state : @xmath176 note that @xmath164 , @xmath177 , and @xmath173 are all @xmath128-dependent ( see eqs .
, , and ) .
equation gives the lower bounds of @xmath2 for each of the @xmath96 states .
note that the lower bound for a state with a given @xmath128 is also the upper bound for the state with @xmath178 . for @xmath33 , @xmath30 , and @xmath179 , @xmath180 , and @xmath65 ,
gives the lower bounds of @xmath181 , @xmath182 , and @xmath183 for the states characterized by @xmath146 , @xmath36 , and @xmath37 , respectively .
these values are in excellent agreement with those obtained by numerical simulations ( see fig .  [
fig : figure1](a ) and table  [ tab : otsm3s2 ] ) .
we note that our approach of focusing on strategy - performance ranking patterns and the fraction of strategies in each rank , represents a generalization of a similar approach @xcite that has already been successfully applied to the mg to cases in which some of the strategies have runaway vps .
we have studied numerically and analytically the effects of a varying resource level @xmath2 on the success rate of the agents in a competing population within the b - a - r model .
we found that the system passes through different states , characterized either by the mean success rate @xmath27 or by the highest success rate in the population @xmath16 , as @xmath2 decreases from the high resource level limit .
the number of states depends on details of the system such as the memory size @xmath5 and the number of strategies per agent @xmath25
. transitions between these states occur at specific values of the resource level . for small values of @xmath5 , it is possible to explain these states by following the evolution of the performance of the strategies and the decision making dynamics of the system .
more generally , we found that different states correspond to different paths covering a subspace within the whole history space . in the high resource level regime , namely @xmath184 , @xmath66 .
the corresponding path in the history space is one that loops around the history ` 111 ... ' indefinitely . just below the high resource level regime
is a range of @xmath2 that gives @xmath5 states corresponding to the fractions @xmath185 .
this result is in excellent agreement with that obtained by numerical simulations . for these @xmath96 states ,
i.e. the outcome series consists of one bit of  0  in a cycle of @xmath142 bits , the path in history space is restricted to those @xmath5-bit histories with at most one - bit of ` 0 ' and with @xmath113 loops around the ` 111 ... ' history .
this identification of an active portion within the history space implies that only part of each strategy is being used . by considering the performance of the strategies within this restricted portion of the history space ,
the number of loops @xmath113 consistent with the @xmath96 states was found to be @xmath127 .
the number of agents using a strategy that predicts the action ` 1 ' given the history ` 111 ... ' increases as the number of loops @xmath113 increases .
thus a criterion on the resource level for sustaining a state of given @xmath113 can be derived .
the results are again in excellent agreement with numerical results . after passing through the @xmath96 states ,
the system goes into states with more than one ` 0'-bit per cycle in the outcome series as @xmath2 is further reduced .
these @xmath186 states correspond to paths that explore an increasingly larger portion of the history space . while our analysis can also be applied to these @xmath186 states , the dynamics and the results are too complicated to be included here .
a resource level @xmath2 that deviates from @xmath13 acts like a _ driving force _ in the history space . in response to this driving force
, the system effectively adjusts its dynamics to occupy an increasingly restricted portion of the history space as @xmath2 increases . for @xmath12 ,
the system explores the whole history space by passing through trails that are almost eulerian @xcite .
the random initial distribution of strategies and the random initial history bit - string that started the system , provide the seed for the diffusive behavior which develops in the history space as the system evolves .
in fact , results of numerical simulations for @xmath187 show that slightly increasing @xmath2 beyond @xmath13 has the effect of suppressing this random wandering through the history space , and instead locks the system into the eulerian trail . however opposing mechanisms can arise to counteract this driving force , thereby enhancing the diffusive behavior .
for example , this can be achieved by allowing the agents the chance of using a strategy besides the best - performing one @xcite or by allowing some agents to opt out of the system occasionally @xcite .
alternatively , the system can be biased through the initial strategy scores , or by introducing a specially prepared non - random initial allocation of strategies @xcite .
it is this competition between the diffusive and driven behavior that gives the non - trivial global behavior in the b - a - r model and its variations . for this reason ,
the present b - a - r system provides a fascinating laboratory for studying correlated , non - markovian diffusion on a non - trivial network ( i.e. history space ) .
this non - trivial diffusion is in turn strongly coupled to the non - random temporal patterns arising in the vp dynamics .
we note that the present results could also be used to generalize the crowd - anticrowd theory in order to incorporate the effect of restricted history - space dynamics : this would then allow identification of an appropriate set of correlated , uncorrelated and anti - correlated strategies in order to implement the crowd - anticrowd theoretical expressions .    as a side - product ,
our analysis serves to illustrate the sensitivity within multi - agent models of competing populations , to tunable parameters . by tuning an external parameter , which we take as the resource level in the present work
, the system is driven through different paths in the history space which can be regarded as a ` phase space ' of the system .
the feedback mechanism , which is built - in through the decision making process and the evaluation of the performance of the strategies , makes the system highly sensitive to the resource level in terms of which states the system decides to settle in or around .
these features are quite generally found in a wide range of complex systems .
the ideas in the analysis carried out in the present work , while specific to the b - a - r model used , are also applicable to other models of complex systems .    in closing , we remark that besides obtaining analytically the highest success rate @xmath16 and the criteria on the resource level @xmath2 , our treatment can also be extended to obtain the mean success rate @xmath27 for the @xmath96 states .
the analysis is more complicated than that for obtaining @xmath16 .
the procedure is to follow the evolution of the performance of the groups of strategies in _ each timestep _ through a path in the history space .
the number of agents taking a particular action , and hence the number of winning agents , can be found from a strategies performance table like the one shown in table  [ tab : strategies10 ] .
we have carried out the analysis for @xmath27 for @xmath40 , @xmath180 , and @xmath65 , and results are found to be in excellent agreement with numerical results .
our analysis can also be readily extended to consider connected populations in which agents have established links to connected neighbors for collecting additional information @xcite . quite generally , the effect of the links is to modify the number of agents _ using _ a strategy in a particular rank . for connected populations ,
an agent may use a strategy that he does not hold but has access to through his links . for an agent who uses the best performing strategy among his own @xmath25 strategies and those of his connected neighbors
, the success rate behaves in a similar fashion as a function of @xmath2 as that reported here , only that the critical values of resource level at which transitions occur are shifted @xcite .
these results can be understood by incorporating the effects of the linkages into eq .. results for the b - a - r model in a connected population will be reported elsewhere @xcite .
this work was supported in part by the research grants council of the hong kong sar government through grant no .
we acknowledge useful discussions with ho - yin lee , keven k. p. chan , charley s. choe and sean gourley .
.table showing the states characterized by @xmath16 for b - a - r model with @xmath33 agents , @xmath34 and @xmath30 , together with the ratio of number of occurrences of ` 1'-bits to ` 0'-bits in the outcome bit - string and the range of resource level in which the state occurs .
the results are obtained from the numerical data as shown in fig .
[ fig : figure1](a ) .
the last column shows the period observed in the outcome bit - string . [ cols="^,^,^,^",options="header " , ]        ( dark solid line ) as a function of the resource level @xmath2 for a system with @xmath33 agents , @xmath30 strategies per agent , with memory length ( a ) @xmath34 and ( b ) @xmath40 . for each value of @xmath5 ,
data for different values of @xmath2 are taken in a system with the same initial distribution of strategies among the agents .
also shown are the range corresponding to one standard deviation about @xmath27 in the success rates among the @xmath0 agents ( dotted lines ) and the spread in the success rates given by the highest and the lowest success rates ( thin solid lines ) among the agents . ]    .
the nodes correspond to the @xmath9 possible histories .
the transition between nodes are indicated by the arrows , together with the outcome needed for the transitions to occur .
the histories can be grouped into columns labelled by a parameter @xmath64 which gives the number of ` 0 ' bits in the histories .
( b ) for the @xmath96 states at @xmath188 , the system follows a path restricted to histories in the @xmath70 and @xmath71 columns , with @xmath113 loops at the @xmath70 history . ]


two unforeseen developments in recent years have made continued photometry of particularly important : distinctive brightness variations in the near ir accompany its mysterious 5.5-year spectroscopic cycle @xcite and , possibly independent of the cycle @xcite , the star brightened at a surprising rate after 1997 @xcite .
neither of these phenomena has been explained and the data since 1998 are seriously incomplete .
meanwhile , the long - term brightening trend continues ; for general information about @xmath0 car with many references , see @xcite , and @xcite .
photometry of this bright object is difficult for at least two reasons : +    1 .   at visual wavelengths , normal ground - based observations
represent mainly the surrounding `` homunculus '' ejecta - nebula , which appears much brighter than the central star and has structure at all radii from 0.2 to 8 arcseconds .
so far the only available measurements of just the central star have been made with the hubble space telescope ( hst ) .
although the homunculus is primarily a reflection nebula , its apparent brightness measured with respect to the central star has changed greatly . during 1998 - 99 , while ground - based observations showed about a 0.3-magnitude brightening of homunculus plus star , the star itself nearly tripled in apparent brightness @xcite !
this discrepancy presumably involves dust along our line of sight to the star , but , as mentioned above , it is only vaguely understood at this time . 2 .
numerous strong emission lines produced in the stellar wind perturb the results for standard photometric systems . h@xmath1 and h@xmath2 emission , for example ,
have equivalent widths of about 800 and 180 respectively in hst spectra of @xmath0 car . therefore
, broad - band @xmath3 , @xmath4 , @xmath5 , and @xmath6 magnitudes , and most medium - band systems as well , are poorly defined for this object .
photometry around 5500 , e.g. broad - band @xmath7 , is relatively free of strong emission lines , but transformations from instrumental magnitudes to a standard system are imprecise because they involve the other filters ( @xcite , @xcite , @xcite , and references therein ) .    in this paper
we report photometry _ of the central star _
obtained with two hst instruments since 1998 . for reasons noted above , these observations are quite distinct from all ground - based data ; and the star s photometric behavior during 2003 proved to be especially interesting .
first we employ acquisition images produced by the space telescope imaging spectrograph ( stis ) , because these are numerous and internally reliable ; @xmath0 car has been observed many times with this instrument since 1998 , especially during mid-2003 when the most recent `` spectroscopic event '' occurred .
these broad - band red data can not be transformed to any standard photometric system , but they are quite reliable for showing relative fluctuations by the star ( strictly speaking , the stellar wind ) , uncontaminated by light from the surrounding homunculus .
we also report surprising variations in the star s h@xmath1 emission brightness accompanied by measurements of the nearby continuum in the stis slit - spectroscopy data .
we supplement the stis data with similar measurements of a few images obtained in 200203 with hst s advanced camera for surveys ( acs ) .
these represent a different set of wavelength bands .
altogether these stis and acs results probably represent most of the photometry of @xmath0 car that will ever be obtained with the hst .
very few additional stis observations are expected now that the 2003 event has passed .
earlier hst / fos , foc , wfc - pc , and wfpc2 data were not obtained often enough to be suitable for our purposes , and in most cases would be more difficult to compare photometrically , for technical reasons .
we report the stis and acs data now , rather than waiting for two or three more data points in 2004 , because the behavior during 2003 was rather unexpected and therefore deserves to be noted promptly .
in addition , we briefly present a convenient summary of @xmath0 car s photometric record in the aavso archives .
these ` v ' magnitudes represent , of course , the entire homunculus plus star but they are valuable for comparisons with the hst data and also with ground - based photometry using different filter systems , reported by , e.g. , @xcite , @xcite , @xcite , @xcite , @xcite , and @xcite .    in section 2
we present the photometry of the central star based on hst / stis acquisition images and the acs / hrc observations .
this is followed by a section on the ground - based data from the aavso and then a general discussion about how these results relate to the models which have been presented to describe @xmath0 carinae .
normal ground - based imaging and photometry can not resolve the central star and small structures in the homunculus nebula . in this study , we have used the hst / stis acquisition images which have a resolution of about 0.1 arcseconds plus the recent acs / hrc images which have a resolution of 0.05 arcseconds to separate the contribution of central star from the surrounding nebulously .
each set of stis observations of @xmath0 carinae has included a pair of acquisition images .
these images are 100x100 pixel sub - frames ( 5 arcseconds square ) centered on the middle row and column of the ccd @xcite .
each pair of acquisition images includes an initial targeting image and a post - acquisition image .
only the post - acquisition images were used to measure the brightness of the central star because the target is rarely well centered in the initial targeting images , which raises concerns about ghost images and internal reflections that may unpredictably affect them ( see section 2.1.2 ) .
the majority of the stis acquisition images of @xmath0 carinae have been taken using the f25nd3 filter , which is a neutral density filter that covers the wavelengths from 2000 to 11000 .
acquisition images using other filters are ignored since they lack significant temporal baseline or coverage and there is little hope of transforming between filter functions for an object as peculiar as @xmath0 carinae .
the f25nd3 filter includes several prominent features in the spectrum of @xmath0 carinae ( figure 1 ) and changes in them contribute to any measured brightness variations .
we processed the stis acquisition images using the iraf procedure stas.hst_calib.stis.basic2d with version 3.0 of stdas and version 2.13 of calstis .
en lieu of a bias over - scan region , each stis acquisition image has the uniform value of 1510 counts subtracted from each pixel as an approximate bias level correction for the target acquisition process .
a small uncorrected bias level remains in the acquisition images after this initial correction , so we compensate for this tiny effect by modeling it as a function of time and ccd housing temperature .
this residual bias was measured in acquisition images of the standard star after the normal bias correction and dark current were subtracted .
the exposures of agk + 81 266 are short ( 2.1 seconds ) and it is more than 60 degrees above the ecliptic , so zodiacal light and other sources make a negligible contribution to the background @xcite .
we found that the residual bias level in the images of agk + 81 266 ( figure 2 ) is roughly linear with respect to time prior to the failure of the ccd side-1 controller on may 16 , 2001 @xcite .
there may be higher order terms in this relation however they are not significant in this study since they effect the measurements of the central star by less than 0.001 magnitude .
subsequent to the failure of the side-1 controller the stis ccd temperature regulator has been held at a constant voltage permitting the ccd temperature to fluctuate . as a result , after the side-1 failure the residual bias level of the observations depends on both time and ccd housing temperature . using the measurements of the residual bias level in the agk + 81 266 observations , the following relations were obtained for the bias level correction as a function of modified julian date ( mjd ) and ccd housing temperature ( t ) : + mjd @xmath8 52045.5 : @xmath9 mjd @xmath10 52045.5 : @xmath11 when the modeled bias level is subtracted from the agk + 81 266 data the background levels average to @xmath12 counts / pixel / sec ( figure 2 ) .
this correction was applied to the acquisition images for @xmath0 carinae in place of the bias over - scan correction normally implemented in the iraf basic2d procedure .
typically the bias correction amounts to less than 0.005 magnitude for the measurements reported below .
the flux from the central star was measured in each image by summing detector counts within a small circular virtual aperture . in order to avoid pixelization effects ,
the counts were weighted by a function of the distance of a pixel from the aperture center ( @xmath13 ) which decreases smoothly to zero at the aperture radius @xmath5 : + @xmath14 @xmath15 the centering on the star is determined by a peak - up algorithm which locates the position where the sum inside the virtual aperture is maximized . for the stis acquisition images we used @xmath16 pixels which corresponds to a diameter of 0.3 arcseconds and
includes slightly fewer than 30 pixels .
we did not attempt to subtract a background from the measured aperture flux , but this has little effect on the results since the star provides thousands of counts / pixel / second , far greater than any other contribution within the measuring aperture due to scattered and zodiacal light which contribute only about 0.2 counts / pixel / second or less @xcite .
there is no feasible way to account for or remove the contribution from diffuse circumstellar material within 0.1 arcseconds of the star .
however , stis spectroscopic data indicate the level of contamination from this source is much smaller than the other uncertainties .
ghost images exist in the stis optics but do not contribute significantly to our measurements .
when a bright source is centered on the ccd , a ghost image with an peak flux of a few percent relative to the psf maximum appears approximately 0.3 arcseconds ( 6 pixels ) to the right ( increasing column number ) of the peak of the point spread function .
the 0.3 arcsecond diameter aperture used to measure the central star excludes this ghost .
however , any future work that may attempt to measure the flux from the nebula immediately surrounding the central star in these images must confront this problem .
the reduction procedure was tested using 80 pairs of acquisition images of the flux standard agk + 81 266 ( figure 3 ) .
the measured aperture photometry for this standard shows a constant level with an r.m.s scatter of @xmath17 magnitudes .
it should be noted that agk + 81 266 was observed using the f28x50lp `` long - pass '' filter rather than the f25nd3 filter which was used for the @xmath0 carinae observations .
at the same time , the gross spectral energy distribution of agk + 81 266 ( spectral class b5 ) and @xmath0 carinae are very dissimilar .
this is not significant since agk + 81 266 is only used to demonstrate that this method accurately measures a constant brightness over the temporal baseline for an established flux standard .
the brightness of @xmath0 carinae is _ not _ measured relative to agk + 81 266 .
in figure 4 and table 2 we present our results for the stis acquisition images as instrumental magnitudes calibrated for long term sensitivity effects using the photometric constants calculated by the reduction pipeline @xcite and scaled relative to the value on 1999 february 20 ( 1999.14 ) .
these represent an extremely broad - band wavelength sample ( figure 1 ) . for the purpose of comparison with familiar photometric systems , table 1
gives approximate johnson @xmath7 @xcite and cousins @xmath5 and @xmath6 @xcite magnitudes at the reference epoch ( 1999.14 ) .
these were calculated using flux - calibrated stis spectroscopic data @xcite , a flux - calibrated spectrum of vega @xcite , and the appropriate response functions @xcite . for reasons noted in section 1 and uncertainties in the stis slit throughput for an object which is not exactly a point source , entries for @xmath0 carinae in table 1 probably have uncertainties of roughly @xmath18 magnitudes and could possibly be somewhat worse .
it is also possible that the photometric color may trend mildly blue - ward as the object brightens @xcite .
the homunculus is currently of 5th magnitude according to _ ground - based _
johnson @xmath7 photometry ( section 3 below ) , but the central star is fainter than @xmath7 = 7.0 as stated in table 1 .
its apparent continuum is red , while stellar - wind emission lines contribute a significant fraction of its brightness . in the 1999.14 spectrum , h@xmath1 and other emission features
contribute about 16% of the flux in @xmath7 ( 4% from h@xmath1 ) and 38% of the flux in @xmath5 ( 35% from h@xmath1 ) .
meanwhile , there is about a 5% affect on the @xmath6 flux from the hydrogen paschen lines which are in emission . because these emission features are constantly changing ,
it is not possible to formulate a consistent set of transformations to convert the stis magnitudes to more familar johnson @xmath7 and cousins @xmath5 and @xmath6 magnitudes for @xmath0 carinae .
however , the values presented for 1999.14 can serve as a rough zero point .
figure 4 shows that the dramatic brightening reported by @xcite leveled off by 2000 , but may have resumed at least temporarily around the mid-2003 spectroscopic event .
the dip in brightness near the time of the event ( figure 6 ) may be reminiscent of the near - infrared behavior reported by @xcite .
this will be discussed further in section 4 .
eta carinae s broad h@xmath1 emission , produced in the stellar wind at radii of several au , is tremendously strong ( equivalent width usually about 800 ) .
this feature did not vary much relative to the continuum in stis data from 1998 to 2002 , but recently it has faded to a surprising degree as reported below .
since h@xmath1 contributes appreciably to the stis acquisition images , we have investigated the effects on our photometry .
davidson ( 2003 ) has measured the equivalent width of h@xmath1 and estimated the absolute level of the nearby continuum in short - exposure stis slit data .
h@xmath1 was integrated from 6509 to 6650  , and the continuum flux was averaged between 6740 and 6800  ( an interval that is practically uncontaminated by emission lines ) .
these continuum measurements have the advantage of being approximately calibrated in physical units , but on the other hand may be underestimated if the slit was not perfectly centered on the star . fortunately , the h@xmath1 equivalent width measurements are quite insensitive to slit position .
the results are shown in figure 5 .
we estimate that h@xmath1 contributes from 8% to 16% of the stis / f25nd3 counts , varying with time as listed in table 2 .
open triangles in figure 4 represent the instrumental magnitudes that would have been observed if h@xmath1 were not present ; evidently the brightness minimum in early 2003 and the subsequent rapid brightening are not due to this emission line .
the measured brightness fluctuations in the central star appear to be mainly ( though not entirely ) changes in the continuum brightness .    in the six months preceding the 2003.5 event ( december 2002 to june 2003 ) , the equivalent width of h@xmath1 declined by a factor of nearly 2 , while its line profile evolved in an interesting and unexpected way which is beyond the scope of this paper .
the other balmer lines behaved in a similar fashion .
moreover , in july and september 2003 they did not return to the state observed by stis in 1998 at the corresponding point in @xmath0 car s 5.54-year cycle .
a more in - depth account of these effects will be reported in a paper which is now in preparation @xcite .
hst acs / hrs observations of @xmath0 car were obtained for the hst treasury project beginning in october 2002 .
bias corrected , dark subtracted , and flat fielded acs / hrc images were obtained from space telescope science institute via the multi - mission archive ( mast ) .
observations which had overexposed the central star were not used and no drizzle or dithering correction is applied .
the acs / hrc images were taken in four filters which cover near uv and optical wavelengths ( figure 1 ) .
each of these filters is influenced by different types of spectral features .
the uv filters ( hrc / f220w & hrc / f250w ) are heavily influenced by the `` fe ii forest '' @xcite .
the opacity of this `` forest '' is known to dramatically increase during an event @xcite .
the hrc / f330w filter is sensitive to both the balmer continuum and a number of emission features while the hrc / f550 m filter , strategically placed between h@xmath1 and h@xmath2 , almost exclusively measures the continuum .
the relative brightness of the central star in the acs / hrc images is measured with the same weighted 0.3 arcsecond diameter aperture used to measure the stis acquisition images . on
the scale of the acs / hrc this corresponds to @xmath19 pixels and includes about 79 pixels inside the aperture .
the results are given in table 3 and figure 6 .
the magnitudes given are on the stmag system @xcite as calibrated by the photometric keywords calculated by the reduction pipeline @xcite .
note that the hrc / f220w and hrc / f250w curves show a significant dip at the time of the 2003.5 event , due to the absorption as noted above .
interestingly , the visual continuum measured by the hrc / f550 m filter does _ not _ show the decrease in brightness seen at longer wavelengths in the stis data .
while space - based photometry is able to resolve the central star , ground based photometry measures the integrated brightness of both the star and the homunculus nebula . however , ground - based photometry is still valuable because it has a longer temporal baseline and shows many of the same brightness variations , albeit at different amplitudes . also as more publicly funded small telescopes
are closed the data from the aavso and groups like la plata observatory @xcite will become increasingly useful for monitoring objects like @xmath0 car .
one of the few sources for long term monitoring of @xmath0 carinae is the american association of variable star observers ( aavso ) @xcite . from 1968 through mid-2003
there are 6640 observations from 110 observers situated around the globe .
these data are from individual observers who have determined the brightness of @xmath0 carinae by visually comparing it with the brightness of other stars .
ten observers account for 69% of these observations while 36% of the observations are from the two most prolific observers : m. daniel overbeek @xcite and albert f. jones . on average , there are 60 observations per observer but the median is 9 observations per observer .
we averaged the observations in 90 day bins with a median of 32 observations in each bin .
the bottom panel of figure 7 shows binned averages from 1968 through mid 2003 with a typical standard deviation of 0.2 magnitudes , which is consistent with the accepted usual r.m.s .
error of 0.1  0.2 magnitudes for observations of this type .
the top panel of figure 7 shows the photoelectric johnson @xmath7 band photometry ( pepv ) from stan walker @xcite , photoelectric observations from r. winfield jones obtained through the aavso , and the average @xmath20 magnitude from ground based observations reported in figure 2 of @xcite .
the data in both panels follow the same general trend and compare well with the data presented in figure 2 of @xcite .
even the most primitive aavso data track the professional photometry well enough to show irregular behavior .
the maximum noted by @xcite at 1982.11 is visible in all the data , though the errors in the binned aavso data make it harder to recognize .
each dataset also shows an overall brightening trend in @xmath0 carinae over the past thirty years along with a notable jump in brightness in 1998 - 1999 @xcite .    the binned aavso data for the last five years ( figure 8) also follows a trend similar to the brightness variations in the central star measured from the stis acquisition images ( figure 4 ) .
of particular note , both datasets show the brightening of @xmath0 carinae peaked in 2000 , fell off a bit , and then peaked again in early 2002 before brightening through the 2003.5 event .
the data presented here from the hst and ground - based observations confirm the dramatic brightening in both the nebula and central star of @xmath0 carinae which was previously reported by @xcite and @xcite along with a more gradual brightening spanning the past thirty years .
all the available data also show that the apparent brightness of @xmath0 carinae fluctuates with no obvious period .
we also see intervals ( on the order of months ) with short lived increases in apparent brightness similar to those reported by @xcite . the 2003.5 event , as observed in the brightness fluctuations of the central star at visual wavelengths ,
was characterized by a roughly six month rise , then a sharp dip in most filters , followed by a recovery to the highest flux observed for the central star in recent years , which may be the start of another significant brightening episode .
again we emphasize that since 1997 _ the central star has brightened substantially faster than the homunculus _
( compare figures 4 and 7 ) .
the stis , acs , h@xmath1 , and aavso data all have notable trends which began about 200 days before the 2003.54 event .
the visual band photometry from the stis , acs , and aavso began to brighten significantly in january 2003 ( figures 9 & 6 and also @xcite ) which coincided with a steep drop in h@xmath1 luminosity ( figure 5 ) and the beginning of a dip in brightness in the near - uv ( figure 6 ) .
the start time for this pre - event activity is important because most models which invoke a binary companion to trigger the events place the secondary on a very eccentric orbit where it is in close proximity to the primary for only a very brief period of time during each cycle .
@xcite showed that the @xmath21 , @xmath22 , @xmath23 and @xmath24 bands brightened roughly 500 days before the time of the events over the past two decades with an additional sharp rise in the @xmath21 , @xmath22 , and @xmath23 bands about 80 days prior to the 1998.0 event .
the x - ray flux from @xmath0 car ( @xcite and m. corcoran s web page for data on the recent event ) also increased before the 1998.0 and 2003.54 events , preceding them by about 350 and 700 days respectively .
it is not necessary or expected that all of the pre - event activity will temporally coincide since each of these data measure activity which occurs at different radii from the star and a shock front which pushes deeper as the event proceeds will not affect each process simultaneously .
also , as noted below there is compelling evidence that the state and distribution of the material immediately surrounding the central star has evolved over time so while the activity prior to the event should maintain roughly the same sequential order , the exact start times relative to the event could potentially change from cycle to cycle .
the acs / hrc data show that the brief dip in apparent brightness during the 2003.5 event is not achromatic .
the acs data show that the depth of the dip is wavelength dependent ; the dip is deeper at shorter wavelengths and the brightness continually increased in the hrc / f550 m filter through the last data point .
however , @xcite showed that the brief decrease or dip in luminosity occurred at all of the near - ir wavelengths during the 1998.0 event ( figure 9 , reproduced from @xcite ) and @xcite reported similar activity during the 2003.5 event .
this leaves two possibilities which are not mutually exclusive :    1 .
the decrease in brightness occurs later at redder wavelengths and/or 2 .
the decrease in brightness is driven entirely by increased line opacities , decreased emission lines , and balmer and paschen continua formed in the wind , rather than a change in the star s continuum flux .
@xcite show that during the 1998.0 event , the dip in the @xmath21 band occurred about 35 days later than the dip in the x - rays .
they also note that the epoch of the minimum in the dip increases slightly with increasing effective wavelength ( figure 9 ) .
@xcite reported that during the 2003.5 event the dip in the @xmath21 and @xmath24 bands began sometime between mjd 52809.5 ( 2003.467 ) and 52814.5 ( 2003.481 ) , which preceded the observed dip in johnson - cousins @xmath25 by about one to two weeks @xcite .
( note that wavelengths around 2 @xmath26 m observed by @xcite presumably represent free - free emission in the wind at radii substantially larger than the photosphere at visual wavelengths . )
given only the hst data , one might suspect that dip in the hrc / f550 m filter may have occurred between data points and the slight rise observed at mjd 52803.1 ( 2003.45 ) is the rise in luminosity which precedes the dip .
however , the hrc/550 m observation at mjd 52840.1 ( 2003.55 ) shows a _ rise _ in flux roughly coincident with the minimum observed in @xmath25 by @xcite and @xcite .
therefore , it is unlikely that we missed the dip in hrc/550 m due to sampling .
likewise , strong p - cygni absorption is present in the hydrogen lines during the event @xcite .
since the hydrogen emission lines contribute significantly to the total brightness of the central star , anything that changes their profiles will change the brightness measured in a filter which they influence .
the hrc / f550 m filter covers none of the prominent emission lines and it shows no dip concurrent with the dip in the other filters .
therefore , this may indicate that the bolometric luminosity of the central star is largely unaffected by the mechanism which causes the observed decrease in brightness .
@xcite discuss several possible explanations for the variability in the brightness of the central star of @xmath0 carinae .
since the star itself radiates near the eddington limit , there is little room for rapid or dramatic changes in bolometric luminosity .
this is supported by the fact that the surrounding nebula , which acts as a giant calorimeter for the central star , varies in brightness on a significantly smaller scale .
if the bolometric luminosity of the central star stays roughly constant , then significant variations in its brightness may be more appropriately interpreted as changes in the density , temperature , and spatial distribution of intervening circumstellar material .
an expanding shell of ejecta from an event may be able to produce the observed light curves @xcite .
any ejected shell should start hot with high opacities from ionized metals and then later shift these opacities to the red as the shell cools and dust forms . particularly at early stages , this shell may be optically thin and as such would not have a significant impact on the continuum brightness .
the data show that central star of @xmath0 carinae continues to brighten .
we are unsure how long this trend can continue but it appears to be driven by an increase in the apparent continuum brightness which may be caused by changes in the intervening circumstellar material along the line of sight , although how or why the circumstellar material may have been been altered is not completely understood . continued monitoring is critical since coming out of the 2003.5 event there are indications that the central star may undergo another period of sudden dramatic brightening in the near future .
the brief decrease in the star s brightness during the 2003.5 event was more complicated than expected and changes preceding the event in the optical and near - uv started as early as january 2003 .
the wavelength dependence of the dip indicates that it can not be caused by a simple eclipse or occultation .
if the primary star is eclipsed by the presumed secondary we would expect to observe some decrease at all wavelengths including the hrc/550 m band which measures the continuum . on the other hand if an eclipse of the secondary is sufficient to cause a dip of approximately 0.5 mag in the near uv flux of the central star , then we would expect to find some spectroscopic evidence for the secondary in the hst / stis spectra during the 5.5 year cycle which is not the case .
however , the data may support a shell ejection model for the spectroscopic events .
additional data , including spectra gathered by the stis during the 2003.5 event , should further constrain the mechanism powering these periodic phenomena .
this research was conducted as part of the hst treasury project on @xmath0 carinae via grant no .
go-9420 from the space telescope science institute .
we acknowledge kris davidson , roberta m. humphreys , and kazunori ishibashi for directly assisting with the data and measurements and t.r .
gull for preparing most of the hst observing plans .
other participants in the hst treasury project are d. john hillier , augusto damineli , michael corcoran , otmar stahl , kerstin weis , sveneric johansson , fred hamann , h. hartman , nolan walborn , and manuel bautista .
we are also grateful to the aavso international database for the variable star observations contributed by observers worldwide .
we individually thank janet a. mattei and elizabeth o. waagen of the aavso who helped us obtain the data and stan walker and r. winfield jones who provided photoelectric photometry data .
we are very appreciative of m. feast , and p. whitelock for allowing us to reproduce their figures .
and we thank matt gray and qian an for providing the @xmath0 carinae project with computer technical support .
altamore , a. , baratta , g.  b. , cassatella , a. , rossi , l. , & viotti , r.  1986 , new insights in astrophysics .
eight years of uv astronomy with iue , 303 bingham , r.  g.  & cousins , a.  w.  j.  1974 , monthly notes of the astronomical society of south africa , 33 , 15 bohlin , r.  c. , dickinson , m.  e. , & calzetti , d.  2001 , , 122 , 2118 cassatella , a. , giangrande , a. , & viotti , r.  1979 , , 71 , l9 clampin , m. , hartig , g. , baum , s. , kraemer , s. kinney , e. kutina , r. , pitts , r. , & balzano , v./ 1996 , stis instrument science report 96 - 030 , ( baltimore , stsci ) colina , l. , bohlin , r. , & castelli , f.  1996 , instrument science report
cal / scs-008 , ( baltimore , stsci ) cousins , a.  w.  j.  1976 , mem .
soc . , 81 , 25 damineli , a.  1996 , , 460 , l49 davidson , k.  & humphreys , r.  m.  1997 , , 35 , 1 davidson , k. , humphreys , r.  m. , ishibashi , k. , gull , t.  r. , hamuy , m. , berdnikov , l. , & whitelock , p.  1999a , , 7146 , 1 davidson , k.  et al .  1999b , , 118 , 1777 davidson , k. , ishibashi , k. , gull , t.  r. , & humphreys , r.  m.  1999b , asp conf .
ser .  179 : eta carinae at the millennium , 227 davidson , k.  2000 , cosmic explosions : tenth astrophysics conference .   aip conference proceedings , vol .  522 .   college park , maryland , 11 - 13 oct 1999 .   edited by stephen s.  holt and william w.  zhang .
american institute of physics , 2000 .
, p.421 - 430 , 421 davidson , k. et al .
2003 , personal communication davidson , k. et al .
2004 , , in preparation davis , m. , campbell , d. , sticka , r. , faful , b. leidecker , h. , kimbel , r. , & goudfrooij , p.  2001
, stis failure review board final report ( baltimore : stsci ) downes , r. , clampin , m. , shaw , r. , baum , s. , kinney , e. & mcgrath , m.  1997,stis instrument science report 97 - 03b , ( baltimore , stsci ) feast , m. , whitelock , p. , & marang , f.  2001 , , 322 , 741 fernandez lajus , e. , gamen , r. , schwartz , m. , salerno , n. , llinares , c. , farina , c. , amorn , r. , & niemela , v.  2003 , informational bulletin on variable stars , 5477 , 1 gull , t.  r. , davidson , k. , & ishibashi , k.  2000 , cosmic explosions : tenth astrophysics conference .
aip conference proceedings , vol .  522 .
college park , maryland , 11 - 13 oct 1999 .   edited by stephen s.  holt and william w.  zhang .
american institute of physics , 2000 .
, p.439 - 442 , 439 holtzman , j.  a. , burrows , c.  j. , casertano , s. , hester , j.  j. , trauger , j.  t. , watson , a.  m. , & worthey , g.  1995 , , 107 , 1065 ishibashi , k. , davidson , m.  f. , corcoran , k. , drake , s.  a. , swank , j.  h. , & petre , r.  1999 , asp conf .
ser .  179
: eta carinae at the millennium , 266 johnson , h.  l. & morgan , w.  w.  1953 , , 117 , 313 johnson , h.  l.  & morgan , w.  w.  1951 , , 114 , 522 mattei , janet & fraser , brian  2003 , jaavso , 31 , 65 mattei , janet & foster , g.  1998 , iappp communications , 72 , 53 pavlovsky , c. , et al .  2003 , acs instrument handbook , version 4.0 , ( baltimore : stsci ) proffitt , c.  r.  et al .
2002 , the 2002 hst calibration workshop : hubble after the installation of the acs and the nicmos cooling system , proceedings of a workshop held at the space telescope science institute , baltimore , maryland , october 17 and 18 , 2002 .   edited by santiago arribas , anton koekemoer , and brad whitmore .
baltimore , md : space telescope science institute , 2002 .
, p.97 , 97 kim quijano , j. , et al .
2003 , stis instrument handbook , version 7.0 , ( baltimore : stsci ) smith , e. , stys , d. , walborn , n. , & bohlin , r.  2000 , instrument science report stis 2000 - 03 , ( baltimore : stsci ) smith , n. , davidson , k. , gull , t.  r. , ishibashi , k. , & hillier , d.  j.  2003 , , 586 , 432 sterken , c. , de groot , m.  j.  h. , & van genderen , a.  m.  1996 , , 116 , 9 sterken , c. , freyhammer , l. , arentoft , t. , & van genderen , a.  m.  1999 , , 346 , l33 sterken , c. , freyhammer , l.  m. , arentoft , t. , & van genderen , a.  m.  2001 , asp conf .
ser .  233 : p cygni 2000 : 400 years of progress , 39 stys , d.  j. & walborn , n.  r.  2000 , instrument science report stis 2001 - 01r , ( baltimore , stsci ) van genderen , a.  m. , de groot , m.  j.  h. , & the , p.  s.  1994 , , 283 , 89 van genderen , a.  m.  et al .
1995 , , 304 , 415 van genderen , a.  m. , sterken , c. , de groot , m. , & burki , g.  1999 , , 343 , 847 van genderen , a.  m. , sterken , c. , & allen , w.  h.  2003a , , 405 , 1057 van genderen , a.  m. , sterken , c. , allen , w.  h. , & liller , w.  2003b , , 412 , l25 viotti , r. , rossi , l. , cassatella , a. , altamore , a. , & baratta , g.  b.  1989 , , 71 , 983 whitelock , p.  a. , feast , m.  w. , koen , c. , roberts , g. , & carter , b.  s.  1994 , , 270 , 364 whitelock , p.  a. , marang , f. , & crause , l.  2003 , , 8160 , 2 zanella , r. , wolf , b. , & stahl , o.  1984 , , 137 , 79     carinae .
the bottom panel shows the total relative response of each ccd and filter combination used in this study on the same wavelength scale as the top panel . for plotting purposes the curves
are not representative of relative responses between filters .
stis filters are plotted with a dashed line and acs / hrc filters are plotted with a solid line .
the dotted line represents the product of the stis ccd+f25nd3 response curve and the photon flux from the central star .
]                    lllrrrrrrr o4j802qxq&50814.008&1997.999&3.196e-12&13.2&2.35e-12&0.65&0.65&&0.79 + o4j801y4q&50891.414&1998.211&3.720e-12&12.7&3.07e-12&0.49&0.49&&0.64 + o55601hkq&51142.172&1998.898&5.374e-12&14.2&4.47e-12&0.09&0.09&&0.25 + o55602qsq&51230.449&1999.140&5.831e-12&13.7&5.00e-12&-0.00&-0.00&&0.16 + o5f102bkq&51616.484&2000.197&7.476e-12&13.3&&-0.27&-0.28&0.01&-0.13 + o5kz01coq&51623.785&2000.217&7.645e-12&&&-0.29 & & & + o5kz02gbq&51624.285&2000.218&7.580e-12&&&-0.28&-0.26&0.02 & + o5f101yoq&51626.762&2000.225&7.288e-12&&&-0.24 & & & + o5f103a1q&51826.039&2000.769&6.753e-12&&&-0.16&-0.16 & & + o62r01h0q&52016.734&2001.291&6.518e-12&11.9&5.60e-12&-0.12&-0.12&&0.02 + o6ex03a9q&52183.086&2001.747&7.573e-12&13.5&6.52e-12&-0.28&-0.28&0.00&-0.12
+ o6en01c2q&52183.297&2001.747&7.509e-12&&&-0.27 & & & + o62r02fvq&52240.070&2001.903&7.470e-12&13.5&6.10e-12&-0.27&-0.27&&-0.11 + o6ex02tkq&52293.977&2002.051&7.241e-12&14.5&5.47e-12&-0.24&-0.21&0.02&-0.04 + o6ex01wyq&52294.207&2002.051&6.923e-12&&&-0.19 & & & + o6mo02eaq&52459.480&2002.504&6.933e-12&14.1&5.69e-12&-0.19&-0.20&0.01&-0.04 + o6mo01k6q&52459.781&2002.505&7.089e-12&&&-0.21 & & & + o8gm01a1q&52624.043&2002.955&6.840e-12&15.5&5.01e-12&-0.17&-0.17&&0.01 + o8gm12ukq&52682.859&2003.116&6.730e-12&14.9&5.43e-12&-0.16&-0.12&0.03&0.05 + o8gm11zpq&52683.242&2003.117&6.326e-12&&&-0.09 & & & + o8gm21t9q&52727.234&2003.238&6.904e-12&13.5&5.29e-12&-0.18&-0.18&&-0.03 + o8gm41eqq&52764.266&2003.339&7.462e-12&11.7&5.54e-12&-0.27&-0.27&&-0.13 + o8gm33qsq&52776.398&2003.372&7.662e-12&11.2&5.95e-12&-0.30&-0.29&0.00&-0.16 + o8gm32c6q&52778.469&2003.378&7.623e-12&&&-0.29 & & & + o8gm31fnq&52785.797&2003.398&7.617e-12&10.6&6.24e-12&-0.29&-0.29&0.01&-0.17 + o8gm51fuq&52791.035&2003.412&7.752e-12&&&-0.31 & & & + o8gm52kyq&52791.707&2003.414&7.558e-12&&&-0.28 & & & + o8gm63ohq&52812.031&2003.470&8.264e-12&9.0&6.98e-12&-0.38&-0.37&0.00&-0.27 + o8gm61pxq&52812.250&2003.471&8.212e-12&&&-0.37 & & & + o8gm62dwq&52813.703&2003.475&8.222e-12&&&-0.37 & & & + o8ma71jnq&52825.020&2003.506&7.675e-12&8.9&7.05e-12&-0.30&-0.32&0.03&-0.22 + o8ma72liq&52825.355&2003.506&8.043e-12&&&-0.35 & & & + o8ma81mtq&52849.594&2003.573&7.763e-12&8.3&7.05e-12&-0.31&-0.32&0.01&-0.22 + o8ma82b8q&52851.906&2003.579&7.850e-12&&&-0.32 & & & + o8ma91ubq&52903.355&2003.720&8.830e-12&&&-0.45&-0.48&0.03 & + o8ma92ctq&52904.289&2003.723&9.309e-12&&&-0.51 & & & + o8ma83fzq&52960.590&2003.877&9.969e-12&&&-0.58 & & & +    lllrrrr +   + j8gm1aa7q&52561.039&2002.782&1.385e-12&8.547&8.549&0.011 + j8gm1aa8q&52561.039&2002.782&1.402e-12&8.533 & & + j8gm1aalq&52561.047&2002.782&1.366e-12&8.562 & & + j8gm1aamq&52561.047&2002.782&1.375e-12&8.554 & & + j8gm1abmq&52561.098&2002.782&1.426e-12&8.515&8.535&0.015 + j8gm1abnq&52561.098&2002.782&1.404e-12&8.531 & & + j8gm1abyq&52561.105&2002.782&1.394e-12&8.539 & & + j8gm1abxq&52561.105&2002.782&1.374e-12&8.555 & & + j8gm2as1q&52682.547&2003.115&1.562e-12&8.416&8.417&0.007 + j8gm2as5q&52682.551&2003.115&1.574e-12&8.407 & & + j8gm2asaq&52682.555&2003.115&1.555e-12&8.421 & & + j8gm2asoq&52682.566&2003.115&1.549e-12&8.425 & & + j8ma3adhq&52803.055&2003.445&1.316e-12&8.601&8.571&0.018 + j8ma3adlq&52803.086&2003.445&1.372e-12&8.557 & & + j8ma3adqq&52803.090&2003.445&1.373e-12&8.556 & & + j8ma3ae3q&52803.152&2003.446&1.355e-12&8.570 & & + j8ma4aqfq&52840.145&2003.547&8.093e-13&9.130&9.129&0.003 + j8ma4aqj&52840.148&2003.547&8.069e-13&9.133 & & + j8ma4aqnq&52840.156&2003.547&8.109e-13&9.128 & & + j8ma4ar2q&52840.168&2003.547&8.137e-13&9.124 & & + j8ma6ayjq&52957.809&2003.869&9.216e-13&8.989&8.970&0.01 + j8ma6aynq&52957.844&2003.869&9.531e-13&8.952 & & + j8ma6az2q&52957.855&2003.869&9.385e-13&8.969 & & +   +   +   + j8gm1aaaq&52561.039&2002.782&2.935e-12&7.731&7.726&0.011 + j8gm1aabq&52561.043&2002.782&2.995e-12&7.709 & & + j8gm1aaoq&52561.047&2002.782&2.910e-12&7.740 & & + j8gm1aapq&52561.051&2002.782&2.951e-12&7.725 & & + j8gm1abpq&52561.098&2002.782&2.985e-12&7.713&7.722&0.018 + j8gm1abqq&52561.098&2002.782&3.022e-12&7.699 & & + j8gm1ac0q&52561.105&2002.782&2.895e-12&7.746 & & + j8gm1ac1q&52561.105&2002.782&2.937e-12&7.730 & & + j8gm2as2q&52682.547&2003.115&3.181e-12&7.643&7.646&0.006 + j8gm2as6q&52682.551&2003.115&3.198e-12&7.638 & & + j8gm2ascq&52682.559&2003.115&3.174e-12&7.646 & & + j8gm2asqq&52682.570&2003.115&3.148e-12&7.655 & & + j8ma3adiq&52803.055&2003.445&3.058e-12&7.687&7.651&0.021 + j8ma3admq&52803.086&2003.445&3.194e-12&7.639 & & + j8ma3adsq&52803.094&2003.445&3.184e-12&7.642 & & + j8ma3ae5q&52803.156&2003.446&3.208e-12&7.634 & & + j8ma4aqgq&52840.148&2003.547&2.178e-12&8.055&8.063&0.008 + j8ma4aqkq&52840.152&2003.547&2.137e-12&8.076 & & + j8ma4aqpq&52840.156&2003.547&2.171e-12&8.058 & & + j8ma4ar4q&52840.168&2003.547&2.159e-12&8.065 & & + j8ma6aqgq&52957.805&2003.869&2.619e-12&7.855&7.840&0.015 + j8ma6aykq&52957.809&2003.869&2.628e-12&7.851 & & + j8ma6aypq&52957.848&2003.869&2.710e-12&7.817 & & + j8ma6az4q&52957.859&2003.869&2.664e-12&7.836 & & +   +   +   + j8gm1aacq&52561.043&2002.782&3.306e-12&7.602&7.621&0.014 + j8gm1aadq&52561.043&2002.782&3.253e-12&7.619 & & + j8gm1aaqq&52561.051&2002.782&3.243e-12&7.623 & & + j8gm1aarq&52561.051&2002.782&3.189e-12&7.641 & & + j8gm1abrq&52561.098&2002.782&3.378e-12&7.578&7.608&0.021 + j8gm1absq&52561.102&2002.782&3.292e-12&7.606 & & + j8gm1ac2q&52561.105&2002.782&3.273e-12&7.612 & & + j8gm1ac3q&52561.109&2002.782&3.201e-12&7.637 & & + j8gm2as3q&52682.547&2003.115&3.327e-12&7.595&7.598&0.002 + j8gm2as7q&52682.551&2003.115&3.312e-12&7.600 & & + j8gm2asgq&52682.563&2003.115&3.313e-12&7.599 & & + j8ma3adjq&52803.055&2003.445&3.913e-12&7.419&7.379&0.023 + j8ma3adnq&52803.090&2003.445&4.095e-12&7.369 & & + j8ma3adwq&52803.098&2003.445&4.112e-12&7.365 & & + j8ma3ae9q&52803.160&2003.446&4.119e-12&7.363 & & + j8ma4aqhq&52840.148&2003.547&3.671e-12&7.488&7.494&0.003 + j8ma4aqlq&52840.152&2003.547&3.649e-12&7.495 & & + j8ma4aqtq&52840.160&2003.547&3.648e-12&7.495 & & + j8ma4ar8q&52840.176&2003.547&3.642e-12&7.497 & & + j8ma6aqhq&52957.805&2003.869&4.643e-12&7.233&7.222&0.008 + j8ma6aylq&52957.809&2003.869&4.695e-12&7.221 & & + j8ma6aytq&52957.852&2003.869&4.738e-12&7.211 & & + j8ma6az8q&52957.852&2003.869&4.678e-12&7.225 & & +   +   +   + j8gm1aafq&52561.043&2002.782&3.700e-12&7.479&7.485&0.006 + j8gm1aatq&52561.051&2002.782&3.663e-12&7.491 & & + j8gm1abuq&52561.102&2002.782&3.755e-12&7.463&7.475&0.014 + j8gm1ac5q&52561.109&2002.782&3.649e-12&7.494 & & + j8gm1ac4q&52561.109&2002.782&3.745e-12&7.466 & & + j8gm2as4q&52682.551&2003.115&4.129e-12&7.360&7.359&0.004 + j8gm2as8q&52682.555&2003.115&4.162e-12&7.352 & & + j8gm2asjq&52682.566&2003.115&4.134e-12&7.359 & & + j8gm2asxq&52682.578&2003.115&4.118e-12&7.363 & & + j8ma3adkq&52803.059&2003.445&4.963e-12&7.161&7.126&0.020 + j8ma3adoq&52803.090&2003.445&5.156e-12&7.119 & & + j8ma3adzq&52803.102&2003.445&5.168e-12&7.117 & & + j8ma3aecq&52803.164&2003.446&5.206e-12&7.109 & & + j8ma4aqiq&52840.148&2003.547&5.259e-12&7.098&7.099&0.003 + j8ma4aqmq&52840.152&2003.547&5.250e-12&7.100 & & + j8ma4aqwq&52840.164&2003.547&5.234e-12&7.103 & & + j8ma4arbq&52840.180&2003.547&5.272e-12&7.095 & & + j8ma6ayiq&52957.805&2003.869&6.851e-12&6.811&6.807&0.004 + j8ma6aymq&52957.813&2003.869&6.907e-12&6.802 & & + j8ma6aywq&52957.855&2003.869&6.890e-12&6.804 & & + j8ma6azbq&52957.867&2003.869&6.849e-12&6.811 & & +