Tax Evasion Dynamics via Ising Model Spin S = 1

In this work, we study the problem of tax evasion on Erdös Rènyi random graphs. Here, we consider that the agents may be in three different states, namely honest tax payers, tax evaders, and undecided that are individuals in an intermediate class among honests and evaders. Every individual can change his/her state following an Ising model dynamics with spin S = 1. In addition, we consider the punishment rules of the Zaklan econophysics model, for which there is a probability a p of an audit each agent is subject to in every period and a length of time k detected tax evaders remain honest. The dynamic of temporal evolution of the Zaklan model was studied initially via the equilibrium Ising model with two opinions (−1 and +1), and recently via a non-equilibrium three-state kinetic agent-based model on a fully-connected population. Here, through Monte Carlo simulations, we study the problem of the tax evasion fluctuations using an Ising model with spin S = 1 (−1, 0, and +1) on Erdös Rènyi random graphs in the dynamic of the temporal evolution of the Zaklan model. Then, we found that the Ising model with spin S = 1 is as efficient as the Ising model and non-equilibrium three-state kinetic agent-based model in controlling the tax evasion fluctuations. This control is even better when we use strong punishment values k even for low audit probabilities a p .


Introduction
In several branches of science, it is common to formulate problems using structures in networks (graphs). This approach considers the forming parts as well as the interactions of these parts and is able to explain the emergence of effects arising from emergent behavior. In economic sciences, they are present in a system of interaction between consumers and producers, a network of suppliers for the synthesis of high technology and high added value products. Research in complex networks benefits from the increase of processing power, which allows the analysis of real network data. Initial research of real networks was guided by the modelsŕandoms from Erdös & Rènyi (1959). An Erdös Rènyi (ER) random graph is a set of N vertices (sites or nodes) connected by c links (bonds) (Erdös et al., 1959(Erdös et al., , 1960(Erdös et al., , 1961. The probability p that a given pair of sites is connected by a In social systems, humans tend to adopt herding behavior and follow the crowd because they feel more comfortable when their decisions are supported by the like decisions of other people. In financial markets, this behavior is reinforced by so-called noise traders agents, who follow trends and over-react to both good and bad news when they buy and sell. In contrast, some agents find that following the global minority brings the best return, i.e., they buy when noise traders depress prices and sell when noise traders push prices up. These contrarian traders are also called fundamentalists, sophisticated traders or α-investors (Lux, 1996(Lux, , 1999(Lux, , 2000Bornholdt, 2001;Kaizoji et al., 2002;Takaishi, 2005;Long et al., 1990;Day et al., 1990;Mukherjee & Chatterjee, 2016). They base their decisions not on market euphoria but on rational expectation, and they push prices toward fundamental values. Another classification of agents in terms of two basic financial market strategies has also been investigated using realistic models. These models considered the presence of two subgroups containing individuals who are optimistic or pessimistic about the future development of the market (Lux et al., 1999(Lux et al., , 2000. Does Econophysics make sense? (Stauffer, 1999) Econophysics is an approach to quantitative economy using ideas, models, conceptual and computational methods of statistical physics (Ausloos, 2013). The name econophysics, a hybrid of economyand physics, was coined to describe applications of methods of statistical physics to economy in general. In practice, majority of the research concerned the finances. In such a way, physicists entered officially and scientifically the field of financial engineering. On top of similar statistical methods used by financial mathematicians, physicists concentrated on the analysis of experimental data using tools borrowed from the analysis of real complex systems. Bloomquist (2006), Föllmer (1974, Andreoni (1998), Lederman (2003), Slemrod (2007), Wintrobe & Gërxhani (2004) through an analysis of the social and economic behavior of a community of people fulfilled indicate that tax evasion in a community is a major cause of concern for governments and through em-  Gächter (2006), Frey & Torgler (2006) have provided that the group members or neighborhood of tax evaders are important in deciding whether or not to pay taxes. In 2008 physicist D. Stauffer andeconomists Zaklan et al. (2008a, 2008b) and Westerhoff (2008) developed an economics model to study the problem of tax evasion dynamics on a people community using the equilibrium Ising model with spin S = 1/2 (two states +1 and −1) on a square lattice (SL) and Monte-Carlo simulations with the Glauber algorithms. On SL the Ising model presents a phase transition well-defined second-order phase transition at critical temperature 2.269 c T =  . The problem of tax evasion fluctuations can be better described as a non-equilibrium than by an equilibrium system. So, different than Zaklan et al. (2008aZaklan et al. ( , 2008b, Oliveira (1992), Lima (2010) proposed the use of the non-equilibrium Majority-Vote Model (MVM) to study the tax evasion in any complex network (Lima, 2012a). He points out that people do not live in a social equilibrium, and any social noise can drive to a government or market chaos. Then, it is reasonable to expect that a non-equilibrium model (MVM) explains better events of non-equilibrium. In order to show that the problem of tax evasion fluctuations is better described by a non-equilibrium model than by an equilibrium model, Lima (2012bLima ( , 2012cLima ( , 2015Lima ( , 2016) studied the Zaklan model (ZM) on Apollonian, Opinion-Dependent, Solomon, and Small-World Networks. In all these cases the tax evasion problem was analyzed using the two-state version of MVM, where the honest agent was rated +1 and the evaded value −1.
In 2014, Crokidakis (2014) studied the problem of tax evasion via ZM on a fully-connected population. In his work, the agents may be in three different states (−1, 0, and +1), namely honest taxpayers (+1), tax evaders (−1) and undecided (0). The undecided agents are individuals in an intermediate class among honests and evaders. Each agent can change their state following a kinetic exchange opinion dynamics, where the agents interact by pairs with negative (probability p) and positive (1 − p) affinity parameters ij µ , representing agreement/ disagreement between pairs of agents. Using Monte Carlo simulations, Lima (Lima, 2012d) has studied the Ising model with spin S = 1/2 and 1 on ER random graphs, with z neighbors for each spin. In the case with spin S = 1/2, the ER graphs present a spontaneous magnetization in the universality class of mean field theory, where on ER random graphs the model presents a spontaneous magnetization at which is the percolation threshold. For ER graphs with spin S = 1, he finds a first-order phase transition for z = 4 and 9 neighbors.
Here, we study the behavior of the tax evasion on an agents community of honest citizens, tax evaders, and undecided where the agents are positioned on sites of ER random graphs, but now using an Ising model with spin S = 1 for the temporal evolution of ZM. The article is organized as follows: We present, in the next section, the model and some of the Monte Carlo simulation details. In Section 3, the results are presented, and in the final section, some concluding remarks are discussed.

ZM via Equilibrium Dynamics of Ising Model Spin S = 1
We use the ZM via ferromagnetic spin S = 1 Ising model to study the tax evasion fluctuations on a community of homogeneous agents located on ER random graphs. In every time period, each ER random graphs site is inhabited by an agent (individuals) with opinion i σ , who can either be an honest taxpayer It is assumed that initially, everybody is honest. Each period the agents can rethink their behavior and have the opportunity to become the opposite type of agent they were in the previous period or undecided. In the same way, an undecided individual may become an honest or evader in the next time step. Then, each agent's social on ER random graphs may either prefer tax evasion, reject it or remain indifferent. The evolution in time of these systems is given by a single spin-flip like dynamics with a probability i P given by ( ) where T is the temperature, B k is the Boltzmann constant, and i E is the energy of the configuration obtained from the Hamiltonian where the sum runs over all neighbor pairs of sites and the spin S = 1 variables i σ assume values ±1 and 0. In the above equation, J is the exchange coupling. The spin S = 1 case is well known in the literature (Onsager, 1944;Kaufmann, 1949;Baxter, 1982).
The ZM presents an enforcement mechanism that consists of two components: a probability of an efficient audit a p ; and a number k of periods. Then once the tax evasion is detected, the tax evaders, where N is the total number and honest N is the honest number of agents. The

Results and Discussion
Following the ZM, we start with a population of 4 10 N = individuals on ER random graphs, where initially everyone is considered honest Then we apply the rules of the ZM as the punishment period k and the audit probabilities a p .
In Figure 1 (1) and (2) lead the system to a disordered state with an equal fraction 1/3 (30%) for ±1 and 0 states. Thus, from Figure 1 that if the audits are efficient ( 90% a p = ) the tax evasion can be considerably  reduced to ≈18% for 10 k = and to ≈5.0% for 50 k = . This behavior of fluctuations of tax evasion is identical to that reported in ZM on regular lattices and networks (Zaklan, 2008a(Zaklan, , 2008bLima, 2010Lima, , 2012. For the cases where 5% a p = , considered more realistic by the literature, the punishment is more efficient when the penalty duration is high as 50 k = . In this case, the tax evasion can be reduced for values around 30%. Notice that when the value of T decrease, the fraction of tax evaders also decreases, this is, for high T, the fraction of opinions −1 is greater than in the cases of lower values of T. In Figure 3(a) and Figure 3 Zaklan et al. (2008a) using Monte Carlo Simulation and the Ising model in twodimensional with two states (−1 and +1) (equilibrium model) performed the first numerical attempt to model the fluctuations tax evasion problem. Then Ising model was used as a tool for the dynamic evolution of the Zaklan model. However, the Ising model cannot be used in some geometries, such as Apollonius and Barabasian networks due to the lack of phase transition. Lima (2010) overcame this obstacle by proposing to use the majority vote model (non-equilibrium model) as a temporal evolution of the Zaklan model. Both proposals were successful in their particular topologies.

Conclusion
Here, we studied the behavior of the tax evasion via the Ising model with spin S = 1 (equilibrium model) that presents three states −1, 0, and +1. As we saw above, −1 represents evaders, 0 undecided, and +1 the honests. We have also verified that the fraction of undecided agents favors the reduction of tax evasion. This fact, together with the control given by the public policies may lead to low levels of evasion.
The behavior of the Ising model with spin S = 1 in relation to controlling the tax evasion fluctuation is very similar to those found by the Ising and majority vote models. Then, its model is another useful tool for studying the tax evasion problem in a community. We still verify that our results are in agreement with the results of Crokidakis (2014) for the non-equilibrium model with three different states (−1, 0, and +1) disagreement between pairs of agents. Therefore, we found the plausible result that tax evasion is diminished by stronger punishment k and audit probability a p .