A study of UK household wealth through empirical analysis and a non-linear Kesten process

We study the wealth distribution of UK households through a detailed analysis of data from wealth surveys and rich lists, and propose a non-linear Kesten process to model the dynamics of household wealth. The main features of our model are that we focus on wealth growth and disregard exchange, and that the rate of return on wealth is increasing with wealth. The linear case with wealth-independent return rate has been well studied, leading to a log-normal wealth distribution in the long time limit which is essentially independent of initial conditions. We find through theoretical analysis and simulations that the non-linearity in our model leads to more realistic power-law tails, and can explain an apparent two-tailed structure in the empirical wealth distribution of the UK and other countries. Other realistic features of our model include an increase in inequality over time, and a stronger dependence on initial conditions compared to linear models.


Introduction
The dynamics of wealth and income inequality is a subject of increasing research interest and public debate, encapsulated by major works such as Piketty's 'Capital in the 21st Century' [1]. The recent COVID-19 pandemic has added to the debate on inequality as some of the very richest, particularly in the tech industry, have gained large quantities of wealth whilst many 'ordinary' households have faced redundancies and reliance on government benefits [2]. Data on standard inequality measures, such as the Gini coefficient as well as wealth or income shares, clearly indicate that inequality has increased since the 1980s in many areas of the world [3]. Potential contributing factors include globalisation, financialisation, decreased taxes, increased tax evasion and avoidance, increased inheritance and domination of the technological sector [1,3,4]. In this paper we summarise these multitude of factors into an idealised growth model for household wealth, dominated by one simple effect: that the wealthier you are, the higher your rate of return (ROR), i.e. the return on wealth you are likely to receive grows superlinearly with wealth. We refer to this type of reinforcement dynamics in our discrete time model as a non-linear Kesten process, which is a generalisation of the work on linear reinforcement initiated by Kesten [5]. The increasing dependence of RORs on wealth has been confirmed in recent studies [6][7][8], and we present further empirical evidence for the UK. Our model uses an agent-based approach, which describes the wealth of individual households as a function of time. The dynamics of individual agents is kept as simple as possible (in our case they evolve independently) and the goal is to predict the collective behaviour of the system via statistical properties of the ensemble of agents. Stochastic agent-based models with multiplicative noise applied to income and wealth dynamics have a long history in economics, with an early major publication in 1953 by Champernowne [9], and since then have been applied extensively and are summarised in several reviews, see for example [10,11]. These models have been used as they exhibit power-law tails, which is a key feature of both income and wealth distributions. Research in the field of Econophysics has focused mostly on exchange of money or wealth, in analogy to energy transfer in models of statistical mechanics (see e.g. [12,13] for an overview). It has been found also in this context, that additive noise leads to Boltzmann-Gibbs type distributions with exponential tails, and heavy tails can result from multiplicative noise or disorder [14]. The focus on pure exchange dynamics has been recognised as unrealistic to model wealth (see [12] page 13), but only very few studies consider both exchange and growth. In [15,16] the authors study growth dynamics of wealth with a global redistribution dynamics, inducing a weak mean-field type interaction between agents. In our model we disregard wealth exchange between households and focus entirely on growth dynamics. This is of course a simplification, but in our view and in line with previous studies mentioned above, growth is clearly the dominant aspect of wealth dynamics for most households, and on average nominal wealth has been growing in an exponential fashion since at least the industrial revolution [17]. Fig 1 shows the tail of the household wealth distribution for the UK from recent wealth and asset survey (WAS) data [18] and rich lists [19,20]. We see here the presence of two power laws in the upper tail: one for the richest in the survey with exponent around 2, and one for the richest in society found in the rich lists with exponent around 1. Such a change in power-law exponent has been observed for other countries [21], and is often argued to be a sampling artefact from survey bias in the data [21,22]. However, due to the particular strength of the effect we believe that the two-tailed structure is a genuine feature of the data. From previous studies [23] linear Kesten processes are known to lead to asymptotic log-normal distributions of wealth. Our non-linear model produces a power-law tail from various generic initial conditions, and in the long run also a two-tailed structure due to a crossover phenomenon resulting from the non-linearity, which we will explain in detail.
We also find that our model has a strong dependence on initial conditions, corresponding to the idea of a low social mobility [24]. It is particularly suitable to describe wealth dynamics since the 1980s, when deregulation of financial markets started to facilitate increasing rates of return for assets typically held by wealthier agents [25], providing increased access to credit and investment opportunities. During the 2007-2008 financial crisis, shortage of available credit temporarily also affected wealth growth for households [26]. But after a relatively short period of adaption and in spite of declining interest rates [27], prices of e.g. housing and financial assets are again increasing at close to pre-crisis levels [28], so the main premise of our model remains valid. While an important macroeconomic question, the mechanisms behind wealth growth are not part of our discussion and we focus on the distribution of wealth among households. Throughout this paper we only model positive wealth, while appreciating that a significant fraction (above 10% [18]) of the UK population has negative wealth, i.e. is in debt. This requires additional modelling and the dynamics we propose do not apply in this case.
Wealth can be defined as assets minus liabilities [29] and is usually measured in a particular currency, GBP in our case. It can be interpreted as the balance sheet of a household, and therefore only assets that can be assigned a monetary value contribute, excluding e.g. health or education of members of a household. We also note that wealth is a stock of value unlike income, which represents a flow of value over time. The WAS categorises wealth into four components: physical, financial, property and pension [30]. In our model we assume that wealth increases on average due to two mechanisms: multiplicative growth due to returns on current wealth, and additive residual savings such as excess salary that is not spent on living costs and other expenses which do not contribute to the balance sheet of the household. In general, wealthier agents can diversify their assets, including riskier strategies with higher average returns [7,8]. Different composition of wealth in different wealth deciles is provided in the WAS [18] and summarised in Section SI.2.4 in S1 Appendix.

Model
We consider independent agents (representing households), whose wealth at discrete time n 2 {0, 1, 2, . . .} (representing years) is denoted by W n > 0. As explained in the introduction, we focus on wealth growth rather than exchange, and model the dynamics of positive wealth only, keeping track of bankruptcy events after which we reset the wealth value of the agent (see Section 4 for details). We assume that the wealth of an agent over the time period n to n + 1 changes via two mechanisms: returns on existing wealth, where R nþ1 2 R denotes the corresponding rate of return (ROR), and residual savings S n+1 � 0, resulting for example from excess earnings which are independent of the current wealth of an agent (see Section 3.3 for details). This leads to the recursion Here the RORs R n and residual savings S n are independent random variables. It is commonly accepted that RORs depend monotonically on wealth [6][7][8], and we assume the

PLOS ONE
following power-law form, where a n 2 R are i.i.d. random variables from some fixed probability distribution, and with small probability can also take negative values. The very simple choice (2) is consistent with empirical data for the UK presented in Section 3.1. We are not claiming that this is the best or most detailed model for RORs, which have been observed in some cases to exhibit an intermediate plateau rather than a strict increase as a function of W n (see e.g. Fig 2 in [8]). But our aim here is to capture the most essential features in a simple model that can also be analysed mathematically, and it is of course possible in simulations to replace (2) by different functions. We find that a non-central t distribution (see Section SI.2.6 in S1 Appendix for details) provides a good match with data for α n , which is discussed in Section 3, Fig 4. Substituting (2) in (1) gives the recursion With γ > 1 we refer to (3) as a non-linear Kesten process. We now summarise theoretical results of (3) for different values of γ. γ = 1. In this case R n = α n and W n+1 = (1 + α n+1 )W n + S n+1 . The stationary version of this linear model has been introduced and studied by Kesten [5], and the non-stationary asymptotic growth case is more recently discussed in [23]. It is easy to see that the asymptotic behaviour of W n is dominated by the exponential e n log j1þa n j , and we present details on the analysis of both cases in Section SI.1.2 in S1 Appendix. In the stationary case with m≔E½logj1 þ R n j� < 0, the model is known to exhibit power-law tails in the limiting distribution, but for wealth dynamics the non-stationary case of asymptotic growth is most relevant, which occurs for μ > 0. Following results in [23], the asymptotics is given by a log-normal distribution such that to leading exponential order where ν 2 ≔ Var[log|1 + R n |] and Z � N ð0; 1Þ is a standard Gaussian. Here the symbol � means that W n ¼ W 0 expðmn þ ffi ffi ffi ffi ffi ffi ffi nn 2 p Z þ oð ffi ffi ffi n p ÞÞ as n ! 1, where o(a n )/a n ! 0 for all sequences ða n : n 2 NÞ with a n ! 1. The rigorous version of this result is subject to further reasonable and mild regularity assumptions on the distributions of parameters (see Theorem 2 (i) in [23]), and the leading order behaviour is independent of the residual savings S n . Since (3) is linear in W n , the model also has a natural scale invariance for the units of wealth (see discsussion in [16]), and the initial condition W 0 enters (4) as a simple multiplicative constant. γ > 1. To our knowledge the non-linear model has not been studied before. Details are given in Section SI.1.3 in S1 Appendix, where we find asymptotic super-exponential growth to leading order, where D is given by a convergent series depending on the distribution of α n and the initial behaviour of the process. Again, we focus on the non-stationary case with W 0 e D > 1. In contrast to the linear case, we see that the asymptotics depend in a strong, non-linear way on the initial conditions and early dynamics of the process. Therefore there is no central limit theorem on the logarithmic scale that leads to (4), and we are not able to predict the asymptotic scaling distribution of W n . But numerical results presented in Section 4 show that the model exhibits power-law tails with realistic shapes on relevant time scales.
For realistic initial conditions and parameters the dynamics follows initially an exponential growth regime, and super-exponential growth sets in when the dominant term in yearly gains in Eq (3) changes from W n to a nþ1 W g n (additive residual savings again do not influence the asymptotic behaviour). This means that the returns from wealth in a single year become of the same order or higher than current wealth, which happens for values around Billionaire return data in Fig 2 below indeed confirm that RORs of around 100% or more can be achieved. From numerical results in Section 4 we see that this crossover leads to a twotailed structure of the distribution of W n similar to what we see in the data in Fig 1, and we think this feature of the model provides a promising explanation for this effect. Since we find in the next section that γ is close to 1, (6) is very sensitive to the value of the random variable α n+1 (which is raised to a large power), leading to a broad crossover region. While this crossover is a realistic feature seen in data from the UK and other countries ( [21], but notably not in the USA, see online Appendix of [21]), the non-linearity also implies that the model is not scale invariant and coefficients will depend on the currency unit.
We further find empirically that α n is mostly positive with a heavy tail, but negative values are possible, see Fig 4 of Section 3.2, and thus W n may become negative. Since our dynamics (1) are not built to describe agents in debt, we replace W n with one of three replacement mechanisms discussed in Section 4.1. We note that bankruptcy events where agents' losses exceed their current wealth are realistic and do occur, but in this paper we focus on modelling the dynamics of agents with positive wealth.
We also note that both, the non-stationary linear and super-linear models, exhibit monopoly, where the wealth fraction of the richest agent in a system of N independent agents tends to 1 as time n ! 1. This behaviour is well known for distributions with heavy tails (see e.g. Table 3.7 in [31]), which include the log-normal distribution in the linear case (4), and is only more pronounced in the super-linear model with heavier tails. We present related numerical results for the Gini coefficient and the top 1% wealth share in simulations, both tending to 1 in the long-time limit. While of course this extreme limit is not realistic currently, inequality measures are well known to increase since the 1980s (see summary in Section SI.2.5 in S1 Appendix). This is consistent with understanding current wealth distributions as transient behaviour of our model, which leads to monopoly if parameters remain unchanged over time. Of course we can only parametrise our model over the current range of wealth values, and in order to get more realistic forecasts for future wealth distributions, we would have to include also the lifetime and inheritance dynamics for agents and the role of external influences (such as war or other catastrophies). The simplified model we present here explains how current wealth distributions can arise naturally from generic initial conditions, and we discuss possible refinements for further study in Section 5.

Data analysis
Before moving on to the simulations of the non-linear Kesten process (3) we undertake some key empirical analysis to parametrise the model. We calculate returns on wealth, R n , and the prefactor, α n , and make statistical fits on these variables. Although residual savings do not evolve with wealth as mentioned above, they are correlated with initial wealth values of an agent as part of their social status or fitness. To infer this dependence, we look at UK income and expenditure data for the year 2016 [32,33].

Statistical properties of returns R n
We rearrange (1) to find the ROR as For wealthy agents, wealth gain is to a large extent dominated by returns on wealth, so that W n+1 − W n � S n+1 and residual savings can typically be ignored. The ROR is then simply given by the wealth growth rate, which we will use to compute R n for billionaires, while we include residual savings to estimate ROR from survey data for other agents.
As mentioned previously, fairly recent work [6][7][8] has suggested an increasing wealth dependence on returns. We also find empirical evidence for this from WAS as summarised in Fig 2, and assume a simple power-law relationship as in (2) which is roughly consistent with the data. According to this we have We fit the power-law exponent γ and the prefactor μ as shown in Fig 2, and also find evidence that returns are independent across time and the variance of returns is proportional to the square of the mean returns as wealth increases (see Fig 3), Such a quadratic scaling relationship of mean and variance is common in multiplicative processes, and consistent with our model assumption (2), as is explained in Section SI.1.1 in S1 Appendix.
Note that the apparent structure in percentile returns data in Fig 2 for individual years does not constitute reliable information in our view, since the variation of the points is artificially decreased due to our numerical procedure as explained in Section SI.2.4 in S1 Appendix. Viewing all years as a combined dataset, we find an increasing wealth dependence of RORs consistent with a simple power-law relationship, which also matches well with data for billionaires. In the next subsection we present a method to estimate a reasonable value of the powerlaw exponent γ so that both, WAS and billionaire return data, can be modelled well with our assumption on returns (2).

Fitting α n
With (3) we have in analogy to (7) As illustrated in Fig 4, we choose the power-law exponent γ = 1.075, such that the return data from the WAS and billionaires can be best explained with a single power law of the form (2). We fit the distribution of the α n (which we assume to be i.i.d.) with a shifted and scaled non-central t-distribution (nct), i.e. we take a n � nctðk; c; l; sÞ:

PLOS ONE
This distribution has four parameters: k > 0 represents the degrees of freedom controlling the heaviness of the tail, c 2 R is the centrality that controls the skewness of the distribution, l 2 R is the shift and s > 0 is the scale, see Section SI.2.6 in S1 Appendix for details.
We find that, while the bulk of the distributions of α n agree well, the billionaire data lead to heavier tails than WAS data. Again, our method of extracting returns from WAS data leads to decreased fluctuations, and therefore we use the parameter values corresponding to billionaire data in simulations in Section 4.

Residual savings S n
We recall that in our model (1) residual savings S n represent all contributions to wealth growth that are independent of the current wealth of an agent. They do not evolve with increasing  [19]. Power law fits according to (8) to the cluster of WAS ROR data combined over all four time periods, leads to μ � 0.003, γ � 1.192 (with both parameters free) and to μ � 0.013 with chosen γ = 1.075 (justified below in Fig 4). We also include γ = 1 for comparison, leading to μ � 0.032, i.e. an average ROR of about 3%. Respective shaded regions are one standard deviation around the power fit means (15) as explained in Section SI.1.1 in S1 Appendix. https://doi.org/10.1371/journal.pone.0272864.g002

PLOS ONE
wealth and only contribute additive noise, which does not influence the long-time behaviour of the dynamics. However, we need to estimate residual savings and their correlation with (initial) wealth to run simulations, and in particular in order to extract empirical RORs from wealth data using (7), which determine the statistics of the crucial parameter α n . [34] presents evidence for recent years in the US, that income and salary are positively correlated with wealth.
We estimate residual savings by equivalised disposable income after expenditure for increasing deciles of median wealth using ONS data sources [32,33]. Equivalised disposable income is household size adjusted income available for spending after tax and deductions, and by expenditure we summarise costs that do not contribute to wealth, such as buying food or paying rent. We fit the dependence on wealth w with a logistic function This is illustrated in Fig 5, where we show data on equivalised disposable income, household expenditure and give the fitted parameter values for (11).
We used (11) as an estimate for additive contributions to wealth growth when calculating percentile returns in Fig 2, see Section SI.2.4 in S1 Appendix and in simulations in Section 4.2 as a function of initial wealth w = W 0 . Note that the logistic fit levels off at κ 1 = 10 6 for large values of w which is an arbitrary cap of 10 6 GBP on wealth independent savings. For most rich  (10) for WAS data percentiles [18] for four time periods along with 2016 billionaire data plotted against wealth W n . We choose γ = 1.075 so that the means of WAS and billionaire data essentially agree (dotted lines). Right: Kernel density of α n+1 for WAS data and 2016 billionaire data as seen in the left Figure. (11) to equivalised disposable income after expenditure, where we choose κ 1 = 10 6 and fit κ 2 = 4.13 � 10 9 and κ 3 = −1.308. ONS data sources used can be found in [32,33]. https://doi.org/10.1371/journal.pone.0272864.g005

PLOS ONE
households, contributions to wealth growth significantly beyond this scale are in the form of wealth returns. It is important to note that none of our results are sensitive to the choice of parameters κ 1 , κ 2 and κ 3 , since savings only really play a role in parameter estimation or simulations on the scales shown in Fig 5.

Simulation results
For all simulations presented in this section we use i.i.d. α n � nct(k, c, l, s) with parameters k ¼ 2:008 ; c ¼ 0:941 ; l ¼ À 0:00156 and s ¼ 0:0112 ; ð12Þ corresponding to data from individual billionaires which represent our best estimate of fluctuations for individual households for γ = 1.075, see Fig 4. We do, however, experiment with changing γ values in which case we multiply the α n by a positive constant to keep the mean at the same level. This is explained further in Section 4.1.

Generic initial conditions without residual savings
To investigate the general properties and dependence on initial conditions of our model over longer time horizons, we consider the following four different initial conditions each with mean 10000: In other words, in I.1 all agents start with initial wealth 10000, in I.2 agents get 5000 plus an exponentially distributed random amount with mean 5000, in I.3 initial wealth is drawn from an exponential with mean 10000 and in I. 4 it is Pareto distributed with scale parameter x m = 5000 and exponent 2.
It is also possible in our simulations for the wealth W n (i) of an agent i to become negative. In this case we choose one of the following replacements for W n (i): R.1. replace with a proportion of the agent's previous positive wealth value pW n−1 (i)>0 such that p is uniformly chosen from (0, 1] R.2. replace with the agent's previous positive wealth value W n−1 (i)>0 R.3. replace with wealth W n (j)>0 of another uniformly chosen agent j We can think of R.1 as the agent losing a random proportion of wealth, R.2 as no change in the agent's wealth and R.3 as the agent being removed from the system and being replaced uniformly with another agent with positive wealth. We note that R.3 is a simple approximation to resampling the agent's wealth from the current wealth distribution. We focus here on simulations with the more realistic compromise mechanism R.1. In Section SI.3.1 in S1 Appendix we will present simulation results for the more extreme replacement mechanisms R.2 and R.3 which lead to similar results, confirming that our model is not very sensitive on the choice of the replacement mechanism.
For each initial distribution we run the simulations iteratively using (3) for N = 10 6 independent agents and zero residual savings S n = 0 with parameters in (12) and replacement mechanism R.1. We choose S n = 0 for convenience in this section, to isolate the effect of the multiplicative dynamics which is dominant in generating the wealth distribution in this model, see Section SI.1.3 in S1 Appendix. Results for empirical tail distributions at times n = 10, 100, 200 and 300 are presented in Fig 6, using the colour code indicated in I.1-I. 4. We also show standard inequality measures (see Section SI.2.5 in S1 Appendix for the definitions), the Gini coefficient g and the top one percent income share s 0.01 for γ = 1.075 up to time n = 300 in the top left and right of Fig 7. We see that all initial conditions eventually lead to monopoly, and for intermediate times power-law tails emerge in the wealth distribution. Due to the crossover (6) to super-exponetial growth, a two-tailed structure emerges for large times and wealth values.
In Fig 8 we show for comparison empirical tails for γ = 1.19 with α n � 0.23 � nct(k, c, l, s), and for γ = 1 with α n � 2.5 � nct(k, c, l, s), so that average ROR values are well approximated with different fits for m ¼ E½a nþ1 � (8) as shown in Fig 2. For γ = 1 we also compute the two inequality measures g and s 0.01 up to n = 400, see bottom left and right of Fig 7 which shows the independence of initial conditions and slower progression towards monopoly. For the higher value of γ = 1.19 we see that the crossover sets in earlier at more realistic wealth values around 10 7 with a two-tailed structure with quite realistic power-law tails (cf. Fig 1). For the linear model with γ = 1 we see no crossover and can fit the distribution for large times well by a log-normal distribution in accordance with (4). In this case there is also no noticeable difference between distributions originating from different initial conditions as we have seen in Fig  7. This is also illustrated in Fig 9, where we also see a clear dependence of final wealth values on initial conditions in the non-linear case with γ > 1.

Realistic initial conditions
In this section we simulate a realistic scenario for the UK, with N = 23 � 10 6 households, initial conditions W 0 extracted from the UK wealth distribution in 2008, and with fixed (non- Comparing Figs 1 to 10 we see that the two-tailed structures differ slightly: While the heavier tail for billionaires with a power-law exponent of about β = 1 is shifting but well preserved, the stability of the lighter power-law tail for millionaires is not well represented in our simulation. This is because we deliberately chose a simple model assuming that average ROR follows a monotone power law with wealth. While this is largely consistent with data, the survey data for RORs show some plateau behaviour for millionaires clearly visible in Fig 2, which has also been suggested for other countries, see Fig 2 of [8]. This may be related to the changing wealth composition of the very rich [35].

Conclusions
The model defined by the iterative Eq (1) represents a generic evolution of household wealth, based on the well motivated assumption that wealth exchange between households does not play an important role. The particular form (3) of a non-linear Kesten process has been motivated by inferring empirically that RORs increase with household wealth, and that this relationship is consistent with a simple power law with exponent γ as in (2), see also Fig 2. We want to stress that the qualitative results and main features of our model do not depend on this particular choice, which we have taken for simplicity and in order to study the effect of the non-linearity with a single parameter. We have seen from theory and simulations that the asymptotic dynamics of the model (3) and the resulting tail of the wealth distribution is dominated by the exponent γ. For the linear case with γ = 1 the RORs do not depend on wealth, and it is known that wealth grows asymptotically with a lognormal distribution (see Section SI.1.2 in S1 Appendix), which does not correspond to power-law tails seen in real data as in Fig 1. As demonstrated by our main results, the non-linear model with γ > 1 exhibits power-law tails   nct(k, c, l, s), γ = 1.075 and right α n � 2.5 � nct(k, c, l, s) and γ = 1. We use initial conditions I. 3 and replacement mechanism R.1. We see a clear dependence on initial conditions for γ > 1, and essentially no dependence for γ = 1.

PLOS ONE
from generic initial conditions, including even perfect equality or light tailed exponential distributions, see Section 4.1. It also leads to a two-tailed structure resulting from a crossover (6) to super-exponential growth for the richest households.
We now summarise the most important theoretical features and differences of the linear (γ = 1) and the non-linear (γ > 1) non-stationary Kesten process (3): • for all γ � 1, including the linear case, the model exhibits monopoly, i.e. for N independent households the wealth fraction of the richest household increases with time and asymptotically approaches 1 (nevertheless, realistic levels of inequality can of course be achieved on intermediate timescales); • the linear model is ergodic, in the sense that the asymptotic exponential growth rate of household wealth does not depend on the initial condition W 0 . The latter only enters as a multiplicative factor and the model is scale invariant, i.e. wealth can be measured in units of W 0 in a dimensionless way; • the non-linear model is not ergodic, i.e. the asymptotic exponential growth rate depends on W 0 and the early dynamics. It is also not scale invariant, and the non-linearity on the right hand side leads to a critical scale (6) where wealth gain per year can exceed current wealth, which is observed in data for the richest households.
Moreover, we would like to stress that our model is phenomenological and not built from first principles, since we simply assume an empirically motivated non-linear relationship between ROR and current wealth. Therefore the model lacks a natural scale invariance and the parameter α n is not universal, but depends on the units of measurement (the currency) and will vary between different countries/economic areas. On the other hand, the non-linearity induces a crossover scale that can be a possible explanation for an apparent two-tailed structure in the data. This is an important aspect of our model which should be investigated further. While not present in data from the USA, the two-tailed structure has been observed [21,22] for several countries which have a less liberal economic system and put more emphasis on social equality. Related political measures such as taxation then lead to a more even wealth distribution and a lighter power-law tail for rich households including millionaires, while the richest in society distribute their wealth globally and can escape such measures, leading to a heavier tail for billionaires.
Other interesting generalisations to make the model more realistic include dynamics for negative wealth, a realistic treatment of bankruptcy events and also household lifetime and fragmentation over longer time periods, or a household dependence of the parameter α n reflecting variations in "fitness" to generate returns from investment. Also, mechanisms of household interaction possibly via a general redistribution or taxation procedure could be included and could lead to interesting effects on the dynamics similar to recent work in [15]. But the aim of this paper was to introduce a simple model, that can explain the main features of wealth distribution and dynamics, and how they can be explained by a non-linear wealth dependent rate of return.