in recent years , there has been an increasing interest to link the discussion about complex phenomena in natural science , such as physics or biology , in particular to those in life science , such as sociology , economy or regional planning .
one of the challenges is to reveal cross - links between the dynamic models used in the specific fields , in order to find out pieces for a common theory of self - organization and evolution of complexity .    with respect to economic and urban problems , this perspective seems rather new . in his essay about self - organization in economics , krugman ( 1996b ) states : `` so far this movement has largely passed economics by . ''
however , recent years have also seen numerous attempts for developing a new , evolutionary view of economics ( anderson _ et al .
_ , 1988 , wei  bin , 1991 , arthur _ et al .
_ , 1997 , silverberg , 1997 , schweitzer and silverberg , 1998 )    in a world of fast growing population and intercontinental trade relations , the _ local _ emergence of new economic centers on one hand , and the _ global _ competition of existing economic centers on the other hand has a major impact on the socio - economic and political stability of our future world .
this problem also confronts the natural sciences on their way to a new understanding of complex phenomena .
this paper aims to contribute to this discussion by adapting a model of interactive structure formation , which has already proved its versatility in a variety of applications , to the problem of economic agglomeration .
spatio - temporal pattern formation in urban and economic systems has been investigated for quite a long time @xcite .
there is also a well known german tradition in location theory , associated with names like von thnen , weber , christaller ( 1933 ) and lsch ( 1940 ) .
recent approaches ( sparsely recognized especially by american economists ) tackle the problem of settlement formation , location theory and city size distribution within a stochastic dynamical theory , based on master equations and economically motivated utility functions ( dendrinos and haag , 1984 , weidlich and munz , 1990 a , b , haag _ et al .
_ , 1992 , 1994 , weidlich , 1991 , 1997 ) .
different from these investigations , the current paper focuses on an agent based dynamics of economic concentration .
agent models , which have originally been developed in the artificial life community ( maes , 1990 , steels , 1995 ) recently turned out to be a suitable tool for describing economic interaction ( anderson _ et al .
_ , 1988 , marimon _ et al .
_ , 1990 , holland and miller , 1991 , lane , 1992 , arthur , 1993 , kirman , 1993 , epstein and axtell , 1996 , arthur _ et al .
_ , 1997 ) .
the rational agent model , one of the standard paradigms of neoclassical economic theroy @xcite , is based , among others , on the assumption of the agent s complete knowledge of all possible actions and their outcomes or a known probability distribution over outcomes , and the common knowledge assumption , i.e. that the agent knows that all other agents know exactly what he knows and are equally rational .    in this particular form ,
the rational agent is just one example of a _ complex agent _ , which can be regarded as an autonomous entity with either knowledge based or behavior based rules @xcite , performing complex actions , such as bdi ( belief - desire - interactions ) ( mller _ et al .
_ , 1997 ) .
the complex agent i.e. is capable of specialization , learning , genetic evolution , etc .
however , the freedom to define rules and interactions for the agents , could very soon turn out to be a pitfall , because of the combinatoric explosion of the state space : for 1000 agents with 10 rules , the state space contains @xmath0 possibilities , hence , almost every desirable result could be produced from such a simulation model .
the alternative to the complex agent , promoted in this paper , could be the _
minimalistic agent _ , which acts on the possible simplest set of rules , without deliberative actions . instead of specialization
, the minimalistic agent model is based on a large number of `` identical '' agents , and the focus is mainly on cooperative interaction instead of autonomous action .
a version of the minimalistic agent model , which allows to apply the advanced methods developed in statistical physics and stochastic theory , is based on active brownian particles ( steuernagel _ et al .
_ , 1994 , schimansky - geier _
et al . _ , 1995 , 1997 , schweitzer , 1997 ,
et al . _ , 1998 ) .
generally , active brownian particles are brownian particles with internal degrees of freedom . as a specific action , the active brownian particles ( or active walkers , within a discrete approximation ) are able to generate a self - consistent field which in turn influences their further movement and behavior . this non - linear feedback between the particles and the field generated by themselves results in an _ interactive structure formation _ process on the macroscopic level .
hence , these models have been used to simulate a broad variety of pattern formations in complex systems , ranging from physical to biological and social systems ( lam and pochy , 1993 , schweitzer and schimansky - geier , 1994 , ben - jacob _ et al . _ , 1995 ,
lam , 1995 , schweitzer _
et al . _ , 1997 ,
et al . _ , 1997 ,
schweitzer and steinbrink , 1997 ) . in the following , the active brownian particles are regarded as economic agents with two different internal states and a very simple behavior : employed agents , which are immobile , generate a wage field ( as the result of their work ) , while unemployed agents migrate guided by local gradients in the wage field , i.e. they try to move toward those regions in their vicinity with a high productivity and a high value of the wage field .
further unemployed agents can be employed , dependent on the local economic situation , but employed agents also can be unemployed .    in sect .
2 we present our stochastic approach to the problem of employment and migration . in sect .
3 , the general model is used to derive a special dynamics recently applied by krugman ( 1992 , 1996b ) for economic aggregation . in sect .
4 , we present the economic assumptions involved in our model , while sect .
5 describes the results of computer simulations , which show economic aggregation on two time scales .
6 gives some conclusions and an outlook of possible extension of the model presented .
let us consider a two - dimensional system with @xmath1 economic agents , which are represented by active brownian particles .
these particles are characterized by two variables : their current location , given by the space coordinate @xmath2 , and an internal state , @xmath3 , which could be either one or zero : @xmath4 .
active particles with the internal state @xmath5 are considered employed agents , @xmath6 , active particles with @xmath7 are considered unemployed agents , @xmath8 . with a certain rate , @xmath9 ( hiring rate ) , an unemployed agent becomes employed , while with a rate @xmath10 ( firing rate ) an employed agent becomes unemployed , which can be expressed by the symbolic reaction : @xmath11 employed agents are considered immobile , while unemployed agents are able to migrate .
the movement of a migrant may depend both on erratic circumstances and on deterministic forces which attract him to a certain place . within a stochastic approach
, this movement can be described by the following overdamped langevin equation : @xmath12 @xmath13 describes the local value of a deterministic force , which influences the motion of the agent .
we note here , that the agent is _ not _ subject to long - range forces , which may attract him over large distances , but only to a _ local _ force .
that means the migrant does _ not _ count on global information which guide his movement , but responds only to local information , which will be specified later .
the second term in eq . ( [ langev - red ] ) describes random influences on the movement of the individuals modelled by a stochastic force @xmath14 , which is assumed to be gaussian white noise : @xmath15 @xmath16 is a measure of the strength of the stochastic force . as eq .
( [ langev - red ] ) indicates the unemployed agent will move in a very predictable way if the guiding force is large and the stochastic influence is small , and he will act rather randomly in the opposite case .
the current state of the agent community can be described by the canonical @xmath17-particle distribution function @xmath18 , which gives the probability to find the @xmath17 agents with the internal states @xmath19 in the vicinity of @xmath20 at time @xmath21 .
the change of this probability can be described by a multi - variate master equation @xcite , which considers both changes in the internal states and migration of the agents . here ,
we restrict the discussion to the spatio - temporal densities @xmath22 of unemployed and employed agents , which can be formally derived from the @xmath17-particle distribution function : @xmath23 @xmath24 is the kronecker delta function for discrete variables , which is @xmath25 only for @xmath26 and @xmath27 otherwise , while @xmath28 is dirac s delta function used for continuous variables , and @xmath29 is the system size ( surface area ) . with @xmath4 , we obtain the spatio - temporal density of the employed agents , @xmath30 , and of the unemployed agents , @xmath31 . for simplicity , we omit the index @xmath32 , by defining : @xmath33 we only assume that the total number of agents is constant , while the density of employed and unemployed agents can change in space and time : @xmath34 the density of the employed agents can be changed only by local `` hiring '' and `` firing '' processes , which with respect to eq .
( [ react ] ) , can be described by the reaction equation : @xmath35 the density of the unemployed agents can be changed by two processes , ( i ) migration and ( ii ) hiring of unemployed and firing of employed agents . with respect to eq .
( [ langev - red ] ) , which describes the movement , we can derive the following fokker - planck equation :    @xmath36    the first term of the r.h.s .
( [ fpe - t ] ) describes the change of the local density due to the force @xmath37 , the second term describes the migration of the unemployed agents in terms of a diffusion process with @xmath38 beeing the diffusion coefficient , the third and the fourth term describe local changes of the density due to `` hiring '' or `` firing '' of agents .
by now , we have a complete dynamic model which describes local changes of employment and unemployment , as well as the migration of unemployed agents .
however , so far some of the important features of this dynamic model are not specified , namely ( i ) the deterministic influences on a single migrant , expressed by @xmath37 , ( ii ) the hiring and the firing rates , @xmath39 , which locally change the employment density .
these variables depend of course on local economic conditions , hence we will need additional economic assumptions .
krugman ( 1992 ) ( see also krugman ( 1996b ) for the results ) discusses a dynamic spatial model , where `` workers are assumed to move toward locations that offer them higher real wages '' .
krugman makes no attempt `` to model the moving decisions explicitely '' , but he assumes the following `` law of motion of the economy '' : @xmath40 here , @xmath41 is the `` share of the manufactoring labor force in location @xmath42 '' and @xmath43 is the real wage at location @xmath42 .
@xmath44 is the `` average real wage '' , which is defined as : @xmath45 here , the sum goes over all different regions , where @xmath42 serves as a space coordinate , and it is assumed that `` at any point in time there will be location - by - location full employment '' .
the assumed law of motion of the economy , eq . ( [ krugman ] ) , then means that `` workers move away from locations with below - average real wages and towards sites with above - average real wages '' .    in the present form ,
the assumed law of motion of the economy involves some shortages :    1 .
a constant total number of employed workers is assumed , a change of the total number or unemploment is not discussed .
employed workers are considered to move .
in fact , if a worker wants to move to a place which offers him higher wages , he first has to be a free , that means an unemployed worker , then moves and then has to be reemployed at the new location , again .
this process is completely neglected ( or else , it is assumed that it occurs with infinite velocity ) .
workers at location @xmath42 always exactly know about the average wage in the system .
it is not explained where they get the information about the average wage from .
workers move immediately if their wage is below the average , regardless of the migration distance to the places with higher wages , and with no doubt about their reemployment .    in the following
, krugman s law of motion of the economy shall be derived from the general dynamic model presented in sect .
2 . , in order to elucidate the implicit assumptions which lead to eq .
( [ krugman ] ) .
the derivation is based on four approximations :    * ( i ) * in his paper , krugman ( 1992 ) does not discuss unemployment .
but , with respect to the dynamic model presented in sect . 2 , the unemployed agents exist in an overwhelming large number .
thus , the first assumption to derive krugmans law of motion is , the local change of the density of unemployed agents due to hiring and firing can be simply neglected .    *
( ii ) * further , it is assumed that the spatial distribution of unemployed agents is in a _
quasistationary _ state .
this does not mean that the distribution does not change , but that the distribution relaxes fast into a quasistationary equilibrium , compared to the distribution of the employed agents .    with the assumptions ( i ) and ( ii ) , eq .
( [ fpe - t ] ) for the density of the unemployed agents reduces to : @xmath46 integration of eq .
( [ fpe - st ] ) leads to the known canonical distribution : @xmath47}{{\left\langle \exp \left[\int f(r , t)/d_{n}\ ; dr \right ] \right\rangle } }   \\ { \left\langle \exp \left[\int f(r , t)/d_{n}\ ; dr \right ] \right\rangle } & = & \frac{1}{a }   \int_{a } \exp \left [ \int f(r',t)/d_{n } \ ; dr'\right ] \ ; dr \nonumber\end{aligned}\ ] ] where the expression @xmath48 describes the mean value , and @xmath49 is the mean density of unemployed agents .    *
( iii ) * for a derivation of krugmans equation , we now have to specify @xmath37 in eq .
( [ stat ] ) . here
, it is assumed that a single unemployed agent which migrates due to eq .
( [ langev - red ] ) , responds to the _ total local income _ in a specific way : @xmath50}\ ] ] eq .
( [ respo1 ] ) means that the migrant is guided by _
local gradients _ in the _ total income _ , with @xmath51 being the local income and @xmath52 the local density of employed workers .
we note here again , that the migrant does not count on information about the highest global income , he only `` knows '' about his vicinity . with assumption eq .
( [ respo1 ] ) , we can rewrite eq .
( [ stat ] ) using the discrete notation , which is preferred by krugman : @xmath53 the corresponding equation ( [ react2 ] ) for @xmath54 reads in the discrete notation : @xmath55 since the krugman equation deals with shares instead of densities , we have to divide eq .
( [ react - disc ] ) by @xmath56 , which leads to : @xmath57    * ( iv ) * krugman assumes that the total number of employed workers is constant , so we use this assumption to replace the relation between the hiring and the firing rate in eq .
( [ change - l ] ) .
( [ react - disc ] ) yields : @xmath58 @xmath59 in eq .
( [ change - l ] ) is now replaced by the quasistationary value , eq .
( [ stat - k ] ) .
inserting further eq .
( [ sum ] ) into eq .
( [ change - l ] ) and using the definitions of the average wage @xmath44 , eq .
( [ mean - w ] ) , and @xmath60 , we finally arrive at : @xmath61 eq .
( [ krugman-2 ] ) is identical with krugmans law of motion of the economy , eq . ( [ krugman ] ) , if the prefactor @xmath62 in eq .
( [ krugman ] ) is identified as : @xmath63 .
hence , @xmath62 is a slowly varying parameter , which depends both on the firing rate , which determines the time scale , and on the average wage , @xmath44 , which may change in the course of time .
it is an interesting question whether the assumptions ( i)-(iv ) which lead to krugmans law of motion in the economy , have some practical evidence in an economic context . from a more theoretical perspective , it is noteworthy that krugmans equation , ( [ krugman ] ) , ( [ krugman-2 ] ) has an obvious analogy to a selection equation of the eigen - fisher type , @xmath64 . here , @xmath65 is the fitness of species @xmath66 and @xmath67 is the mean fitness representing the global selection pressure .
it can be proved @xcite that this equation describes a competition process which finally leads to a stable state with only _
one _ surviving species .
the competition process may occur on a very long time scale , but asymptotically a stable coexistence of many different species is impossible .    in the economic context discussed by krugman ( 1992 ) ,
( [ krugman ] ) implies that finally _ all _ workers are located in _
one _ region @xmath42 , where the local real wage @xmath43 is equal to the mean real wage , @xmath44 .
however , in his paper , krugman ( 1992 ) ( see also krugman ( 1996b ) ) discusses computer simulations with 12 locations ( on a torus ) which show the stable coexistence of two ( sometimes three ) centers , roughly evenly spaced across the torus , with exactly the same number of workers ( cf .
fig . 7 in krugman ( 1992 ) and fig . 2.2 . in krugman
( 1996b ) ) .
the paper does not provide information about the time scale of the simulations , and stochastic influences ( apart from a random initial configuration ) are not considered .
krugmans model of spatial aggregation includes of course some more complex economic assumptions , such as consideration of transportation costs , distinction between agricultural and manifactured goods , price index etc .
we are not going to discuss whether the stable coexistence of two centers in krugmans computer simulations might result from those specific economic assumptions in his model or from the lack of fluctuations ( which could have revealed the instability of the ( deterministic ) stationary state ) . due to krugman ( 1992 ) , already a slight variation in the parameters always led to a single center .    in the following ,
we will focus on the more interesting question , whether the general dynamic model of migration and employment introduced in sect .
2 , is able to produce a stable coexistence of different economic centers under the presence of fluctuations .
this may allow us to overcome some of the shortages involved in krugman equation .
the derivation of krugmans equation was based on the assumtion ( ii ) that the time scale for hiring and firing of workers is more determining than the time scale of migration .
if we explicitely consider unemployment in our model , this assumption implies that there are always enough unemployed workers which can be hired on demand . a growing economy , however ,
might be determined just by the opposite limiting case : it is important to attract workers / consumers to a certain area before the output of production can be increased .
hence , the time scale for the dynamics is determined by migration processes .
that means the spatial distribution of the employed agents can be assumed in a _ quasistationary equilibrium _ compared to the spatial distribution of the unemployed agents .
@xmath68 hence , the local density of employed agents can be expressed as a function of the local density of the unemployed agents available , which itself changes due to migration on a slower time scale .
we further assume that the unemployed agent who is able to migrate , eq .
( [ langev - red ] ) , responds to local gradients in the real wages ( per capita ) , instead of gradients in the total income , as assumed for krugmans equation .
the migrant tries to move towards places with a higher wage , but again he only counts on information in the vicinity .
hence , the guiding force @xmath13 is determined as follows : @xmath69 for a further discussion , we need some assumptions about the local distribution of the real wages , @xmath70 .
it is reasonable to assume that the local wages may be a functional of the local density of employed agents , @xmath71 .
some specific assumptions about this dependence will be discussed in the next section . with eq .
( [ change - stat ] ) , we can then rewrite the spatial derivative for the wages as follows : @xmath72 where @xmath73 denotes the functional derivative . using the equations ( [ change - stat ] ) , ( [ respo2 ] ) , ( [ wage - dl ] ) , the fokker - planck equation for the change of the density of the unemployed agents , eq . ( [ fpe - t ] )
can now be rewritten as follows : @xmath74 the r.h.s . of eq .
( [ fpe - eff ] ) now has the form of a usual diffusion equation , with @xmath75 being an _ effective diffusion coefficient _ :
@xmath76 here , @xmath77 is the `` normal '' diffusion coefficient of the unemployed agents , eq .
( [ fpe - t ] ) .
the additional terms reflect that the unbiased diffusion is changed because of the response of the migrants to local differences in the wage distribution . as we see from eq .
( [ d - eff ] ) , there are two contradicting forces determining the effective diffusion coefficient : the normal diffusion , which keeps the unemployed agents moving , and the response to the wage gradient .    if the local wage decreases with the number of employed agents , then , and the effective diffusion increases .
that means unemployed agents migrate away from regions where employment may result in an effective decrease of the marginal output .
however , if , then their wage effectively _ increases _ in regions with a larger number of employed agents , and they are _ attracted _ to these regions . as we see in eq .
( [ d - eff ] ) , for a certain _ positive feedback _ between the local wage and the employment density , the effective diffusion coefficient can be _ locally negative _ and unemployed agents do not leave these areas once they are there . economically speaking , these agents stay there , because they may profit from the local increase in the employment density .
we want to emphasize that this interesting dynamic behavior has been derived _ without _ any explicit economic assumptions .
but for a further discussion of the model , the three remaining functions which are unspecified by now : ( i ) , ( ii ) @xmath9 , ( iii ) @xmath10 , have to be specified , and that is of course where economics comes into play .
in order to determine the economic functions , we refer to a _ perfectly competitive industry _ ( where `` perfect '' means complete or total ) .
in this standard model @xcite , the economic system is composed of many firms , each small relative to the size of the industry .
these firms represent seperate economies sharing common pools of labor and capital .
new competitors can freely enter / exit the market , hence the number of production centers is not limited or fixed .
further , it is assumed that every firm produces the same ( one ) product and every firm uses only _
one _ variable input .
then , the _ maximum profit condition _ tells us that `` firms will add inputs as long as the marginal revenue product of that input exceeds its market price . '' in the case of labor as variable input , the price of labor is the wage , and a profit maximizing firm will hire workers as long a the marginal revenue product exceeds the wage of the added worker .    the marginal revenue product @xmath78 is the additional revenue a firm earns by employing one additional unit of input , where @xmath79 is the price of output and @xmath80 is the marginal product . in a perfectly competitive industry , however , no single firm has any control over prices . instead ,
the price results from the interaction of many suppliers and many demanders . in a perfectly competitive industry , every firm sells its output at the market equilibrium price , which is simply normalized to one , hence the marginal revenue product is determined by @xmath80 .
the marginal product , @xmath80 , can be derived from a _ production function _ - , which describes the relationship between inputs and outputs ( e.i . the technology of production ) .
usually , the production function may include the effect of different inputs , such as capital , public goods or natural resources . here , we concentrate only on one variable input , labor .
thus , @xmath81 is assumed a cobb - douglas production function @xmath82 with a common ( across regions ) exponent @xmath83 , @xmath52 being the local density of employees .
the exponent @xmath83 describes how a firms output depends on the scale of operations .
increasing returns to scale are characterized by a lower average cost with an increasing scale of production , hence @xmath84 in this case . on the other hand , if an increasing scale of production leads to higher average costs ( or lower average output ) , we have the situation of decreasing returns to scale , with @xmath85 . in the following
, we will restrict the discussion to @xmath85 , common to all regions .
the prefactor @xmath86 represents economic details of the level of productivity .
we assume that these influences can be described by two terms : .
@xmath87 should summarize those output dependences on capital , resources etc .
, which are not explicitely discussed here .
the new second term , @xmath88 , considers _ cooperative effects _ which result from interactions among the workers .
since all cooperative effects are _ non - linear _ effects , @xmath88 should be a non - linear functional of @xmath52 : @xmath89 .
hence , the production function depends explicitely only on @xmath52 , now : @xmath90 \;l^{\beta}(r , t)\ ] ] the marginal product ( @xmath80 ) is the additional output by adding one more unit of input , i.e. @xmath91 , if labor is input .
the wage of a potential worker will be the marginal product of labor : @xmath92 if a firm faces a market wage rate of @xmath93 ( which could be bound e.g. by minimum wage laws ) , then , in accordance with the maximum profit condition , a firm will hire workers as long as : @xmath94 to complete our setup , we need to discuss how the prefactor @xmath88 depends on the density of employees , @xmath52 .
here , we assume that the cooperative effects will have an effect only in the intermediate range of @xmath95 . for small production centers ,
the synergetic effect resulting from the mutual stimulation among the workers is too low . on the other hand , for very large production centers
, the advantages of the cooperative effects might be compensated by the disadvantages of the massing of agents .
thus we will assure , that @xmath96 both in the limit @xmath97 and @xmath98 .
these assumtions are concluded in the ansatz : @xmath99 where the utility function @xmath100 describes the mutual stimulation among the workers in powers of @xmath95 : @xmath101 the series will be truncated after the second order .
the constants @xmath102 characterize the effect of cooperation , with @xmath103 .
especially , the case @xmath104 considers saturation effects in the cooperation , e.i .
the advantages of cooperation will be compensated by disadvantages of crowding .
this idea implies that there is an optimal size for taking advantages of the cooperative effect , which is determined by the ratio of @xmath105 and @xmath106 . if one believes that cooperative effects are an always increasing function in @xmath95 , then simply @xmath107 can be assumed
if cooperative effects are neglected , i.e. @xmath108 , then we should obtain the `` normal '' production function , @xmath109 .
the constant @xmath110 is determined by a relation between @xmath111 and @xmath112 which results from the eqs .
( [ prod1 ] ) , ( [ a - u ] ) . within our approach , the constant @xmath111
is not specified .
if we , without restrictions of the general case , choose : @xmath113 then the production function , @xmath114 , eq .
( [ prod1 ] ) , with respect to cooperative effects can be expressed as follows : @xmath115\;l^{\beta}\ ] ] fig .
1 presents the production function , eq .
( [ prod2 ] ) dependent on the density of employees , which can be compared with the `` normal '' production function , @xmath109 .
[ fig1 ]    clearly , we see an increase in the total output due to cooperation effects among the workers .
if we assume @xmath116 , this increase has a remarkable effect only in an intermediate range in @xmath95 . both for @xmath97 and @xmath98 the cooperative effect vanishes .
once the production function is determined , we also have determined the local wage @xmath117 , eq .
( [ wage ] ) , as a function of the density of employees : @xmath118\;\beta \,l^{\beta-1 } \nonumber \\ & + & \frac { \bar{a}}{2 } \;\exp\left(a_{1}\,l+a_{2}\,l^{2}\right )   \,(a_{1}+2a_{2}l)\;l^{\beta } \end{aligned}\ ] ] hence , also the derivative @xmath119 , used for the effective diffusion coefficient , @xmath75 , eq .
( [ d - eff ] ) is determined .
2 shows both functions dependent on the density of employees .
[ fig2 ]    the left part of fig .
2 indicates that , within a certain range of @xmath95 , the derivative @xmath119 could be indeed positive , due to cooperative effects .
with respect to the discussion in sect .
, this means that the effective diffusion coefficient , eq .
( [ d - eff ] ) , can be possibly negative , i.e. unemployed workers will stay in these regions .
the transition rate @xmath9 for hiring unemployed workers is implecitely already given by the conclusion , that firms hire workers as long as the marginal revenue product exceeds the wage of the worker . in accordance with eq .
( [ wage - star ] ) we define : @xmath120 here , @xmath121 determines the time scale of the transitions .
@xmath9 , which is now a function of the local economic situation , is significant larger than the level of random transitions , represented by @xmath121 , only if the maximum profit condition allows hiring .
otherwise , the hiring rate tends to zero .
the firing rate @xmath10 could be simply determined opposite to @xmath9 .
then , from the perspective of the employee , firing is caused by the local situation of the economy , i.e. by _ external reasons _ : workers are fired if @xmath122 .
a more refined description , however , should consider also _ internal reasons _ : a worker can not only _ loose _ his job , he can also _ quit _ his job himself for better opportunities , e.g. because he wants to move to a place where he earns a higher wage .
it is reasonable to assume that the internal reasons depend again on spatial gradients in the wage . due to
( [ langev - red ] ) the unemployed agent migrates while guided by gradients in the wage
. the employed agent at the same location may have the same information .
noteworthy again , this is only a local information about differences in the wage .
if the local gradients in the wage are small , the internal reasons to quit the job vanish , and the firing depends entirely on the ( external ) economic situation .
however , if these differences are large , the employee might quit his job for better chances .
we note that the latter process was already considered in krugmans law of motion in the economy , sect .
3 . different from the assumptions involved in krugmans
( [ krugman ] ) , here the process : employment @xmath123 unemployment @xmath123 migration @xmath123 reemployment is _
explicitely _ modeled . it does not occur with infinite velocity , and there is no guarantee for reemployment .
hence , we define the `` firing rate '' which describes the transition from an employed to an unemployed agent , as follows : @xmath124   + \ ; q\ , \frac{{\displaystyle}\partial \omega(r)}{{\displaystyle}\partial r } \right\}}\ ] ] the additional parameter @xmath125 can be used to weight the influence of spatial gradients on the employee .
eventually , in this section we have determined the variables @xmath37 , @xmath9 and @xmath10 via a production function @xmath126 , which represents certain economic assumptions .
now , we can turn back to the dynamic model described in sect .
2 , which now can be solved by means of computer simulations .
before presenting the results of the computer simulations , the simulation technique should be shortly discussed . the computer program has to deal with three different processes , which have to be discretized in time for the simulation : ( i ) the movement of active brownian particles with @xmath127 due to the overdamped langevin eq .
( [ langev - red ] ) , ( ii ) the transition of the particles due to the rates , @xmath9 , eq .
( [ kplus ] ) , @xmath128 , eq .
( [ kminus ] ) , and ( iii ) the generation of the field @xmath70 .
considering eq .
( [ langev - red ] ) , the new @xmath129-position of a particle @xmath66 with @xmath127 on the two dimensional surface at time @xmath130 is given by : @xmath131 the equation for the @xmath132-position reads accordingly .
@xmath133 is the _ non - constant _ time step , which is calculated below .
@xmath38 is the diffusion coefficient , and _ grnd _ is a gaussian random number with mean equals zero and standard deviation equals unity .    in oder to calculate the spatial gradient of the wage field @xmath134 , eq .
( [ wage2 ] ) , we have to consider its dependence on @xmath52 , which is a local density .
the density of employed agents is calculated assuming that the surface is divided into boxes with the spatial ( discrete ) indices @xmath135 and unit length @xmath136 .
then , the local density is given by the number of agents with @xmath137 inside a box of size @xmath138 : @xmath139 we note that @xmath140 , with @xmath141 being the spatial move of the migrating agent into @xmath129-direction during time step @xmath133 .
that means that the migration process is really simulated as a motion of the agents on a two dimensional plane , rather than a hopping process between boxes .
using the box coordinates @xmath135 , the production function and the wage field can be rewritten as follows : @xmath142 the spatial gradient is then defined as : @xmath143 where the indices @xmath144 ,
@xmath145 , @xmath146 , @xmath147 refer to the left , right , lower and upper boxes adjacent to box @xmath148 .
we further note that for the simulations _ periodic boundary conditions _ are used , therefore the neighboring box is always specified .    using the discretized versions , eq .
( [ disc - l ] ) , ( [ disc - y - w ] ) , ( [ disc - grad ] ) , the transition rates can be reformulated as @xmath149 , @xmath150 accordingly .
they determine the _ average number _ of transitions of a particle in the internal state @xmath3 , located in box @xmath148 , during the next time step , @xmath133 . in a stochastic simulation
however , the actual number of reactions is a _ stochastic _ variable , and we have to assure that the stochastic number of reactions ( i ) does not exceed the actual number of particles available during @xmath133 , and ( ii ) is equal to the average number of reactions in the limit @xmath151 .
this problem can be solved by using the stochastic simulation technique for reactions , which defines the appropriate time step , @xmath133 , as a random variable .
let us assume that we have exactly @xmath152 particles with the internal state @xmath127 in the system at time @xmath153 .
then the probability @xmath154 equals one , and the probability for any other number @xmath155 is zero . with this intial condition , the master equation to change @xmath152 reads : @xmath156 where @xmath17 is the total number of particles , and @xmath157 , @xmath158 are the transition rates for each particle .
the solution of this equation yields : @xmath159 here , @xmath160 is the mean life time of the state @xmath152 . for @xmath161 , the probability @xmath162 to find
still @xmath152 is almost one , but for @xmath163 this probability goes to zero .
the time when the change of @xmath152 occurs , is most likely about the mean life time , @xmath164 . in a stochastic process , however , this time varies , hence , the _ real life time _
@xmath165 is a randomly distributed variable . since we know that @xmath162 has values between @xmath166 $ ] , we find from eq .
( [ m - solv ] ) : @xmath167\big\}\ ] ] ( zero actually has to be excluded . ) .
that means , after the real life time @xmath168 , one of the possible processes which change @xmath152 occurs .
each of these processes has the probability : @xmath169 thus , with a second random number _
rnd_[0,@xmath170 it will be determined which of these possible processes occurs
. it will be the process @xmath125 which satisfies the condition : @xmath171 for the transition probabilities to change @xmath152 , we find in particular : @xmath172 & = & \sum_{i=1}^{n_{0 } }   k^{-}_{i }   t_{m } \nonumber \\   prob[n_{1},t\to n_{1}-1,t+\tau ] & = & \sum_{i=1}^{n_{1 } }    k^{+}_{i } t_{m}\end{aligned}\ ] ] obviously , the sum over these probabilities is one , which means , during the time intervall @xmath168 one of these processes occurs with certainity .
further , using the definition of @xmath168 , we see that @xmath173 yields asymptotically , i.e. after a large number of simulation steps , @xmath174 .    hence , determining the time step as @xmath175 , eq .
( [ tau ] ) , ensures that only one transition occurs during one time step and the number of transitions does not get out of control .
further , in the asympotitc limit , the actual ( stochastic ) number of transitions is equal the average number .
the ( numerical ) disadvantage might be that the time step @xmath176 is not constant , so it has to be recalculated after each cycle before moving the particles with respect to eq .
( [ disc - x ] ) .
this may slow down the speed of the simulations considerably .
let us conclude the procedure to simulate the movement and the transition of the particles :    1 .
calculate the density of particles @xmath177 , eq .
( [ disc - l ] ) 2 .
calculate the production function @xmath178 and the wage field @xmath179 , eq .
( [ disc - y - w ] ) 3 .
calculate the sum over all possible transitions , which determines the mean life time @xmath160 , eq .
( [ m - solv ] ) , of the current system state 4 .
calculate from the mean life time the _ actual _ time step @xmath175 , eq .
( [ tau ] ) by drawing a random number 5 .
move all particles with @xmath127 according to the langevin equation ( [ disc - x ] ) using the time step @xmath168 6 .
calculate which one of the particles undergoes a transition by drawing a second random number , eq .
( [ prob - r ] ) 7 .
update the system time : @xmath180 and continue with 1 .      with the determination of the three functions @xmath37 , @xmath9 and @xmath10 with respect to some economic assumptions ,
we have completed our dynamic model , described in sects . 2 and 4 .
in this section , we want to discuss some features of the dynamics by means of a computer with @xmath181 active brownian particles .
initially , every particle is randomly assigned an internal parameter ( either @xmath27 or @xmath25 ) , and a position on the surface , which has been divided into @xmath182 boxes of unit length @xmath183 .
the diffusion coefficient , which describes the mobility of the migrants , is set to @xmath184 ( in arbitrary units ) , so , a simple brownian particle would approximately need a time of @xmath185 for a mean spatial displacement of @xmath25 ( which is the spatial extension of a box ) .
the minimum wage @xmath93 is set to @xmath186 .
further , @xmath187 and @xmath188 , for the remaining parameters see fig
. 1 .    in the following
we will restrict the discussion to the spatio - temporal evolution of the densities of employed and unemployed agents .
other quantities of interest , such as the spatio - temporal wage distribution , the production function and the local values of the `` hiring '' and `` firing '' rates have been of course calculated for the simulation , but will be discussed in a subsequent paper , which also investigates different parameter sets .
3 shows snapshots of the evolution of the spatial density of employed agents , while fig .
4 shows the corresponding spatial density of unemployed agents .
5 presents the evolution of the total number of employed and unemployed agents , in terms of the total share @xmath189 with respect to eq .
( [ total ] ) .
[ fig3 ]    [ fig4 ]    [ fig5 ]    the simulation shows that the evolution of the densities occurs on two stages :    * ( i ) * during the first stage ( for @xmath190 ) we have a significant higher share of employed agents ( up to 70 percent of the agents community ) .
they are broadly distributed in numerous small economic centers , which have an employment density of about 5 .
this can be understood using figs . 1 , 2 .
the increase of the output due to cooperative effects allows the fast establishment of numerous small firms , which count on the mutual stimulations among their workers and in the beginning offer higher wages ( compared to the model without cooperative effects ) .
this growth strategy , however , is not sufficient for output on larger scales , because the marginal product drastically decreases ( before it may increase again for an above critical employment density ) . as a result the wages may fall ( possibly below the minimum wage rate , @xmath93 ) , which prevents further growth .
so , the first stage is characterized by the coexistence of numerous small economic centers .    *
( ii ) * during the second and much longer stage ( for @xmath191 ) , some of these small centers overcome this economic bottleneck , which allows them to grow further . in the simulations , this crossover occured between @xmath192 and @xmath193 . in the model considered here , we have only one variable input , labor , so the crossover is mainly due to stochastic fluctuations in @xmath52 .    after the bottleneck , the marginal output increases again , and also the wage increases with the density of employment ( cf . figs . 1 , 2 ) .
this in turn effects the migration of the unemployed agents , which , due to eq .
( [ langev - red ] ) , is determined by two forces : the local attraction to areas with higher wages , and the stochastic influences which keep them moving . with @xmath194 , the attraction may exceed the stochastic forces , so the unemployed agents are bound to these regions , once they got there .
this , in turn is important for the further growth of these economic centers , which needs agents to hire . as a result
, we observe on the spatial level the concentration of employed _ and _ unemployed agents in the _ same _ regions .    for their further growth , the economic centers which overcome the bottleneck , _ locally _ attract the labor force at the expense of the former small economic centers .
as a result of the competition process , these small centers , which previously employed about 70 percent of the agents , dissapear . in the new ( larger ) centers
however , only 60 percent of the agents can be employed ( cf .
5 ) , so the competition and the resulting large - scale production effectively results also in an _ increase of unemployment_. but the employment rate of 60 percent seems to be a stable value for the given parameter setup , since it is kept also in the long run , with certain fluctuations .    the concentration process described above , leads to the existence of different _ extended major economic regions _ , shown in figs . 3 , 4 , which each consist of some subregions ( in terms of boxes ) . in the long run ( up to t=100.000 )
we find some remarkable features of the dynamics :    * a _ stable coexistence _ of the major economic regions , which keep a certain distance from each other .
this _ critical distance _ - which is a _ self - organized phenomenon _ - prevents the regions from detracting each other .
in fact , finally there is no force between these regions which would pull off the employed or unemployed agents and guide them over long distances to other regions .
thus , the coexistence is really a stable one . * a _ quasi - stationary non - equilibrium _ within the major economic centers .
as we see , even in the long run , the local densities of employed and unemployed agents do not reach a fixed stationary value . within the major regions ,
there are still exchange processes between the participating boxes , hiring and firing , attraction and migration .
hence , the total share of employed agents continues to fluctuate around the mean value of 60 percent ( cf .
fig . 5 ) .
in this paper , we have proposed a simple dynamic model to descibe migration and economic aggregation within the framework of active brownian particles .
we consider two types of particles ( or agents ) ; employed agents , which are immobile , but ( as the result of their work ) generate a wage field , and the unemployed agents which migrate by responding to local gradients in the wage field .
further a transition between employed and unemployed agents ( hiring and firing ) is considered .
the economic assumptions to be used in the model , are concluded in only three functions : ( a ) the derivative @xmath195 , which describes how the local wage depends on the employment density , ( b ) @xmath9 , ( c ) @xmath10 , which are the hiring and firing rates .
these functions can be determined using a production function , @xmath114 , which describes the output of production dependent on the variable input of labor .
our ansatz for @xmath114 specifically counts on the influence of cooperative interactions between the agents on a certain scale of production .    as the result of our dynamic model
, we find the establishment of distinct economic centers out of a random initial distribution .
the evolution of these centers occurs in two different stages . during the first stage ,
small economic centers are formed based on the positive feedback of mutual stimulation / cooperation among the agents . during the second and much longer stage ,
some of these small centers grow at the expense of others , which leads to the interesting situation that in our model economic growth and decline occur _ at the same time _ , but at different locations , which results in specific _ spatial - temporal _ patterns .
the competition finally leads to the concentration of the labor force in different extended economic regions .
this crossover to large - scale production is accompanied by an increase in the unemployment rate .
although the extended economic regions are in an internal , quasistationary non - equilibrium , we observe the stable coexistence between these regions .
this is an important result , which agrees e.g. with the central place theory of economics ( christaller , 1933 , lsch , 1940 ) .
different from an attempt by krugman ( 1992 ) , who also focused on this problem , we find that ( i ) the coexistence is stable , even at the presence of fluctuations , ( ii ) the centers not necessarily have to have the same number of employed agents , to coexist , ( iii ) the dynamics does not simply converges into an equilibrium state , but the centers exist in an quasistationary non - equilibrium state and still follow a stochastic eigendynamics .
one question , which will be tackled in a forthcoming paper , is about the influence of the different parameters .
we have seen in the simulations , that cooperative effects are able to initiate an increase in the economic output , which leads to the establishment of many small firms .
but the final situation does not include any of those small firms and seems to be independent of this intermediate state .
this is mainly due to the fact , that we considered only _
one _ product .
a more complex production function with _
outputs , however , should also result in the coexistence of small ( innovative ) firms and large scale production .
another important issue is to understand the relation between the response of the agents to local gradients in the wage field , on one hand , and stochastic influences , on the other hand .
most likely , there exist critical parameters which may describe a `` phase transition '' within the agents behavior , which eventually results in different types of spatial - temporal coexistence .
finally , one may argue that the simple model proposed here , is not the simplest one to describe the economic aggregation .
hence , the minimalistic agents discussed in sect .
1 , could be even more minimalistic .
one example is the response of the employed agents to gradients in the wage field , expressed in the parameter @xmath125 which enters the firing rate , @xmath10 . @xmath196
simply means that internal reasons are neglected , and the agent is fired only due to external reasons
. this might be also a reasonable assumption , however , it includes a pitfall . in a growing economy , eventually almost all unemployed agents might be hired at some locations .
but , as long as @xmath197 , employed agents are almost never fired , thus there is a shortage of free agents for further growth .
eventually the dynamics sticks in a dead - lock situation , because agents can not move even if there are locations in their neighborhood which may offer them a higher wage .
therefore , it is reasonable to choose @xmath198 , but it is important to understand the role of @xmath125 , for instance , for the final patterns observed .
so , we conclude that our model may serve as a toy model for simulating the influence of different social and economic assumptions on the spatial - temporal patterns in migration and economic aggregation .
this work has been supported by the deutsche forschungsgemeinschaft ( bonn , germany ) .
i would like to thank j. lobo for discussions on an early draft of this paper , and w. ebeling , l. schimansky - geier and g. silverberg for comments on the final version .
anderson , p.w . , arrow , k.j . , pines , d. ( eds . ) ( 1988 ) .
_ the economy as an evolving complex system _ , reading , ma : addison wesley .
haag , g. ; munz , m. ; pumain , p. ; sanders , l. ; saint - julien , th .
interurban migration and the dynamics of a system of cities : 1 .
the stochastic framework with an application to the french urban system , _ environment and planning a _ * 24 * , 181 - 198 .                                      schimansky - geier , l. ; schweitzer , f. ; mieth , m. ( 1997 ) .
interactive structure formation with brownian particles , _ in : _ f. schweitzer ( ed . ) : _ self - organization of complex structures : from individual to collective dynamics _ , london : gordon and breach , pp .
101 - 118 .
schweitzer , f. ; steinbrink , j. ( 1997 ) .
urban cluster growth : analysis and computer simulation of urban aggregations , in : f. schweitzer ( ed . ) , _ self - organization of complex structures : from individual to collective dynamics _ , london : gordon and breach , pp .
501 - 518 .
silverberg , g. ( 1997 ) .
is there evolution after economics ?
, _ in : _ f. schweitzer ( ed . ) : _ self - organization of complex structures : from individual to collective dynamics _ ,
london : gordon and breach , pp .
415 - 425 .
weidlich , w. ( 1997 ) .
from fast to slow processes in the evolution of urban and regional settlement structures , _ in : _
f. schweitzer ( ed . ) : _ self - organization of complex structures : from individual to collective dynamics _ ,
london : gordon and breach , pp .
475 - 488 .