Why multilingual, and how to keep it-An evolutionary dynamics perspective.

While many languages are in danger of extinction worldwide, multilingualism is being adopted for communication among different language groups, and is playing a unique role in preserving language and cultural diversities. How multilingualism is developed and maintained therefore becomes an important interdisciplinary research subject for understanding complex social changes of modern-day societies. In this paper, a mixed population of multilingual speakers and bilingual speakers in particular is considered, with multilingual defined broadly as zero, limited, or full uses of multiple languages or dialects, and an evolutionary dynamic model for its development and evolution is proposed. The model consists of two different parts, formulated as two different evolutionary games, respectively. The first part accounts for the selection of languages based on the competition for population and social or economic preferences. The second part relates to circumstances when the selection of languages is altered, for better or worse, by forces other than competition such as public policies, education, or family influences. By combining competition with intervention, the paper shows how multilingualism may evolve under these two different sources of influences. It shows in particular that by choosing appropriate interventional strategies, the stable co-existence of languages, especially in multilingual forms, is possible, and extinction can be prevented. This is in contrast with major predictions from previous studies that the co-existence of languages is unstable in general, and one language will eventually dominate while all others will become extinct.

Introduction

As the world becomes increasingly globalized, more pan class="Species">people become multilingual across the continents. It is reported that now more than half of the world’s population speak at least two languages [1]. While many languages are in danger of extinction, multilingualism is being adopted as a common way of communication among different language groups, and is playing a unique role in preserving language and cultural diversities [2]. Therefore, how multilingualism is developed and maintained becomes an important research subject in linguistics, social studies, and beyond [3]. The promotion and protection of multilingualism has also been a hot topic long discussed in public [4].

Languages compete and spread among their speakers, as genes are inherited and passed down to biological populations, where some are selected while others become extinct [5]. Genes may be carried over in mixed forms. So are languages by multilingual speakers. Genes flow in and out of a population and are subject to constant mutations among their variants. Speakers also migrate and can as well be converted from speaking one language to another. In this paper, a mixed population of multilingual speakers and bilingual speakers in particular is considered, with multilingual defined broadly as zero, limited, or full uses of multiple languages or dialects, and an evolutionary dynamic model for its evolution is proposed, similar to that for genetic evolution [6].

Unlike genetic evolution though, the uses of languages are not only dependent of competition, but also subject to various societal interventions, common in social or cultural evolution. The proposed model consists of two different parts accordingly, formulated as two different evolutionary games, respectively. The first part accounts for the selection of languages based on the competition for population and social or economic preferences. It assumes that the speakers are able to adjust their ways of using single or multiple languages, although in reality, it may take time and extra learning efforts. The second part relates to circumstances when the selection of languages is altered, for better or worse, by forces other than competition such as public policies, education, or family influences. Under these circumstances, the number of speakers of one language may be increased or decreased due to the impact of certain public policies or the migration of the speakers in or out of the population, etc.

Here are a few examples where both language competition and societal intervention may play a role for the development and maintaining of multilingualism: In twenty-five European countries surveyed, 56 percent of the population speak at least two languages, and 28 percent speak at least three languages [1]. In the United States, over 20 percent of the US population are bilingual [7], and as high as 70 percent of immigrant families speak a language other than English at home [8, 9]. In China, almost all different provinces have their own local Chinese dialects, yet by 2015, over 73 percent of the population have learned Mandarin Chinese, a once northern Chinese dialect later promoted as the official Chinese by the government [10]. However, Cantonese is better preserved in Canton than all other Chinese dialects, as the locals have kept many newspapers, radio broadcastings, movies and TV series in Cantonese [11]. Hebrew, an extinct ancient language, revived in 19th century and later became one of the official languages of Israel, with 9 million speakers worldwide today [12].

Much work has been done on modeling language competition, although not specifically for the evolution of multilingualism. A well known model was proposed by Abrams and Strogatz in 2003 for the study of pan class="Disease">language death [13]. The model was later extended to more general and complex cases by several other groups [14-17]. The models along this line focus mainly on language competition for population and social or economic preferences. In general, the models predict that one language will eventually dominate the population while all others become extinct, and the co-existence of languages is unstable and hard to sustain [13, 16, 17]. While successfully applied to some language populations, these models have not explicitly distinguished language competition from possible societal interventions that may reverse the course of language changes. By combining language competition with possible societal interventions, this paper shows how multilingualism may evolve under these two different sources of influences. It shows in particular that by choosing appropriate interventional strategies, the stable co-existence of languages, especially in multilingual forms, is possible, and extinction can be prevented, as seen in many multilingual communities.

The paper has three core sections titled Evolutionary Models, Dynamic Analysis, and Dynamic Simulation. The evolutionary models for populations with competition or intervention are discussed in Evolutionary Models. In Dynamic Analysis, the dynamic behaviors of populations under different conditions are analyzed. The stability conditions of equilibrium states are also justified for different interventional strategies. In Dynamic Simulation, a two-dimensional dynamic simulation scheme based on the proposed evolutionary model is described. Two sets of simulation results are presented showing potential geographical impacts on the evolution of multilingual populations. To keep it simple and easy to follow, the paper is focused more on bilingual populations, with models for general multilingual populations provided as an appendix in the end. To make the paper more accessible, all formal definitions, theorems, and especially proofs are also organized into the appendix sections. Nonetheless, these appendices are considered to be necessary contents as well for supporting the results and conclusions presented in the paper.

Evolutionary models

To keep the description simple and easy to follow, consider a bilingual population, the simplest yet most common form of multilingualism. Assume it is large and well mixed, i.e., every individual speaker can interact with all others in the population. If an individual speaker uses two languages A and B with frequencies x and x, respectively, this individual is called an (x, x)-speaker, where 0 ≤ x, x ≤ 1, and x + x = 1. Thus, a (1, 0)-speaker is an A speaker, a (0, 1)-speaker is a B speaker, and an (x, x)-speaker with x, x ≠ 0 is a speaker of both A and B, with frequency x for A and x for B.

Likewise, if language A and B are used with frequencies y and y in average in the whole population, this population is called a (y, y)-population, where 0 ≤ y, y ≤ 1, and y + y = 1. Thus, a (1, 0)-population is an A population, a (0, 1)-population is a B population, and a (y, y)-population with y, y ≠ 0 is a bilingual population of A and B, with average frequency y for A and y for B. Here, if the population size is n and speaker i uses A and B with frequencies and , respectively, then and .

Evolution with competition

Let P(y) and P(y) be the payoff functions for A and B speakers in a (y, y)-population, respectively, defined in terms of population percentages and social or economic preferences, with P increasing in y and P in y, meaning that the larger the population percentage of a language, the more benefit the language provides. Then, the payoff function for a general (x, x)-speaker in a (y, y)-population can be defined in terms of the average use of A and B by this speaker: For example, the payoff of a (0.5, 0.5)-speaker, i.e., a half-A and half-B speaker, in a (0.5, 0.5)-population, i.e., a half-A and half-B population, should be 0.5P(0.5) + 0.5P(0.5).

Considering (x, x) as the strategy of an individual and (y, y) the strategy of the population, an evolutionary game (Appendix A) can be defined, where every individual tries to maximize his/her payoff. The latter can be achieved when an optimal strategy is found for every individual. The strategy for the whole population then becomes as well, and a Nash equilibrium is reached. A Nash equilibrium of this game is thus a strategy such that

Assume for a (y, y)-population that y and y vary over time as they reach their equilibrium. Then, a system of replicator equations (Appendix A) can also be defined for the population: meaning that the changing rate of average frequency y (or y) of language A (or B) depends on the payoff P (or P) compared with the payoff P (or P): if it is higher, y (or y) increases; otherwise, it decreases. Based on evolutionary game theory, a Nash equilibrium of the evolutionary game as conditioned in (2) is always a fixed point of the system of replicator equations given in (3) [18-20] (Appendix A).

Evolution with intervention

Now consider the situation where some societal influences other than simple competition intervene in the use of languages. Assume that the societal influences or interventions are implemented to counter the arbitrary increase or decrease of either language in a bilingual population. Let and be the payoff functions for A and B speakers in a (y, y)-population, respectively, defined in terms of societal influences and interventions, with decreasing in y and in y, meaning that the smaller the population percentage of a language, the more incentive or less penalty the speakers of the language receive. Here, incentive means more benefit from the society such as more work opportunities or public supports, etc., and penalty means less advantage in the society such as less access to public education or discouragement from the public, etc. Then, the payoff function for a general (x, x)-speaker in a (y, y)-population can be defined in terms of the average use of A and B by this speaker:

It follows that another evolutionary game can be defined, with the Nash equilibrium being a strategy such that and also a corresponding system of replicator equations:

Note that the game with competition in (2) has three possible equilibrium strategies. They can be obtained as the fixed points of the system of replicator equations in (3): It is not so hard to prove that the strategies in (a) and (b) are evolutionarily stable while the one in (c) is not (see Appendix C). This means that under competition, one of the languages is expected to die eventually while the other takes over the whole population. The co-existence of the two languages is not sustainable.

On the other hand, the game with intervention in (5) also has three possible equilibrium strategies. They can be obtained as the fixed points of the system of replicator equations in (6): However in this case, it can be proved that the strategy in (c) is evolutionarily stable while the ones in (a) and (b) are not (see Appendix C). This means that with intervention, it is possible to keep the population of either language not too large or too small. A bilingual population can be maintained.

Evolution with competition and intervention

Now the two types of games described above can be combined to obtain a complete model for the evolution of a bilingual population. A simple combination can be obtained by merging the systems in (3) and (6) into the following system: where λ ∈ [0, 1] is a scalar. If λ = 1, the system is reduced to the one in (3) with competition only. If λ = 0, it is reduced to the one in (6) with intervention only. If λ ∈ (0, 1), the system is in some sense a convex combination of competition and intervention, with more weight on competition when λ > 0.5 while more on intervention when λ < 0.5. For convenience, however, a simpler system, with an equal weight given to competition and intervention, will be considered in the following discussion: In any case, with a combined system, the impacts from language competition and influences of societal intervention are both included for the uses of the languages. The first part of the system, based on language competition, pushes the population to monolingual, while the second part, counted for societal interventions, prevents the population from losing either language. The combination of the two is then expected to provide a more comprehensive account on the changes of the languages and the evolution of the bilingualism.

Note that the system in (10) can be reformulated in a more compact form: where It is then a system of replicator equations with payoff functions and , and corresponds to an evolutionary game, with the Nash equilibrium being a strategy such that where is the payoff function for the game, and for an (x, x)-speaker in a (y, y)-population

In general, P and P can be any increasing functions and and be any decreasing functions. However, based on previous studies, they can be defined as the following: where c and are scaling constants, α, , s, s, , are all parameters, . The parameters α, determine the order of dependency of the payoffs on the population percentages. Since 1 < α ≤ 2, the payoffs from P and P increase with increasing population percentages. On the other hand, since , the payoffs from and decrease with increasing population percentages. The parameters s, s are used to define the payoffs from competition. They are indicators of social or economic impacts on the payoffs. The larger these values, the more benefit for the corresponding language groups. The parameters are used to define the payoffs from interventions. They are rates for language reversing due to interventions. The larger these values, the faster the reversing rates.

The definitions of the functions in (16) and (17) are actually based on previous work on language competition and, in particular, the work by Abrams and Strogatz 2003 [13]. In fact, with these functions, the systems in (3) and (6) are both equivalent to an Abrams-Strogatz system as shown below, differing only in the ranges of the α values and the parameters c, s, s versus , , : where P is a function, P(y, s) = cy s for some constant c. Abrams and Strogatz 2003 surveyed a large number of regions in UK and South America for language competition and estimated the α value around 1.31 ± 0.25 [13]. In general, the system can be adopted for modeling language competition for 1 < α ≤ 2 as well as language reversing for 0 ≤ α < 1, as discussed above.

Dynamic analysis

Continue the assumption of bilingual populations and consider their possible dynamic behaviors based on the models described in previous sections. Assume that the payoffs are defined by the functions in (16) and (17).

Dynamics with competition

For competition only, an evolutionary game and a corresponding system of replicator equations are given in (2) and (3). With P and P given in (16), the game in (2) is actually a so-called potential game [18, 20] (Appendix B), i.e., there is a potential function f such that Indeed, f(y, y) = (P(y, s) + P(y, s))/α, where P(y, s) = cy s, 1 < α ≤ 2. Therefore, an equilibrium strategy for the game must be a KKT point of the following potential maximization problem: Further, if the equilibrium strategy is evolutionarily stable, it must be a strict local maximizer of the potential maximization problem, and vice versa [18, 20] (Appendix B).

For a given strategy (y, y), the Hessian H of f is given as follows. According to the theory on constrained optimization [21, 22], as a KKT point of the potential maximization problem in (20), if an equilibrium strategy for the game in (2) is a local maximizer, the Hessian projected on the null space of the active constraints at this strategy is necessarily negative semidefinite; conversely, if the projected Hessian is negative definite, the equilibrium strategy must be a strict local maximizer of the potential maximization problem [21, 22]. Thus, the evolutionary stability of the equilibrium strategy can be justified by the negative definiteness of the projected Hessian at the strategy (Appendix B).

The game in (2) has three possible equilibrium strategies corresponding to three fixed points of (3), as given in (7). It is easy to see that strategies (a) and (b) are evolutionarily stable (Appendix C). In either case, one pan class="Disease">language dies while the other takes over the whole population. If s > s, the chance for (a) will be greater than (b), and if s > s, the chance for (b) will be greater than (a). For (c), it is easy to verify that the projected Hessian at this strategy is always positive definite (Appendix C). It follows that the strategy can never be a local maximizer of the potential maximization problem in (20), and can never be evolutionarily stable.

The above behaviors of the game can be further demonstrated in Fig 1, where P and P are plotted against y, and the changing directions of y are pointed with arrows. The graph in Fig 1(a) shows the behavior of a population with α = 2, when P and P are simple linear functions. The curves of P and P are displayed over y ∈ [0, 1] with s = 0.75 and s = 0.25, when language A dominants B. There is an equilibrium strategy in (0, 1) such that , with and . It is unstable as y decreases to 0 if it starts from below while increases to 1 if it starts from above . Note that is relatively small. Therefore, y increases to 1 even if y starts small.

Fig 1

Dynamic behaviors with competition.

Payoff functions P and P are plotted against y and the changing directions of y are pointed with arrows. In (a), α = 2, s = 0.75, s = 0.25, , and . In (b), α = 3/2, s = 0.25, s = 0.75, , and .

The graph in Fig 1(b) shows the behavior of a population with α = 3/2, when P and P are square root functions. The curves of P and P are displayed over y ∈ [0, 1] with s = 0.25 and s = 0.75, when language B dominants A. There is an equilibrium strategy in (0, 1) such that , with and . It is unstable as y decreases to 0 if it starts from below while increases to 1 if it starts from above . Note that in this case is very large. Therefore, y decreases to 0 even if y starts large.

Dynamics with intervention

For intervention only, an evolutionary game and a corresponding system of replicator equations are given in (5) and (6). With and given in (17), the game in (5) is also a potential game, i.e., there is a potential function such that Indeed, for , , and for , , where . Therefore, an equilibrium strategy for the game must be a KKT point of the following potential maximization problem: Further, if the equilibrium strategy is evolutionarily stable, it must be a strict local maximizer of the potential maximization problem, and vice versa [18, 20] (Appendix B).

For a given strategy (y, y), the Hessian of is, Again, as a KKT point of the potential maximization problem in (23), if an equilibrium strategy for the game in (5) is a local maximizer, the Hessian projected on the null space of the active constraints at this strategy is necessarily negative semidefinite; conversely, if the projected Hessian is negative definite, the equilibrium strategy must be a strict local maximizer of the potential maximization problem [21, 22]. Thus, the evolutionary stability of the equilibrium strategy can be justified by the negative definiteness of the projected Hessian at the strategy (Appendix B).

The game in (5) has three possible equilibrium strategies corresponding to three fixed points of (6), as given in (8). It is easy to see that strategies (a) and (b) are unstable (Appendix C). For (c), it is easy to verify that the projected Hessian at this strategy is always negative definite (Appendix C). It follows that this strategy must be a strict local maximizer of the potential maximization problem in (23), and must be evolutionarily stable.

The above behaviors of the game can be further demonstrated in Fig 2, where and are plotted against y, and the changing directions of y are pointed with arrows. The graphs in Fig 2(a) shows the behavior of a population with , when and are reciprocal functions. The curves of and are displayed over y ∈ [0, 1] with and , when language A is more promoted than B. There is an equilibrium strategy in (0, 1) such that , with and . It is stable as y increases to if it starts from below while decreases to if it starts from above .

Fig 2

Dynamic behaviors with intervention.

Payoff functions and are plotted against y and the changing directions of y are pointed with arrows. In (a), , , , , and . In (b), , , , , and .

The graph in Fig 2(b) shows the behavior of a population with , when and are inverse square root functions. The curves of and are displayed over y ∈ [0, 1] with and , when language B is more promoted than A. There is an equilibrium strategy in (0, 1) such that , with and . It is stable as y increases to if it starts from below while decreases to if it starts from above .

Dynamics with competition and intervention

Now consider the dynamic behaviors of a bilingual population under the influences of both inter-language competition and external interventions, as modeled by the system in (10) and the corresponding game in (14). The payoff functions and are simply combinations of those for competition and intervention as given in (12) and (13). With these two payoff functions, the game in (14) is again a potential game, and therefore, there must be a function such that Indeed, , with f and defined in previous sections. It follows that an equilibrium strategy for the game must be a KKT point of the potential maximization problem: Further, if the equilibrium strategy is evolutionarily stable, it must be a strict local maximizer of the potential maximization problem, and vice versa [18, 20] (Appendix B).

For a given strategy (y, y), the Hessian of is Again, as discussed in the previous two sections, if an equilibrium strategy for the game in (14) is a local maximizer of the potential maximization problem in (26), the Hessian projected on the null space of the active constraints at this strategy is necessarily negative semidefinite; conversely, if the projected Hessian is negative definite, the equilibrium strategy must be a strict local maximizer of the potential maximization problem [21, 22]. Thus, the evolutionary stability of the equilibrium strategy can be justified by the negative definiteness of the projected Hessian at the strategy (Appendix B).

For convenience, set for the following analysis. Let , , be an equilibrium strategy for the evolutionary game in (14). Then, , i.e., If , then This is a necessary and sufficient condition for , , to be an equilibrium strategy for the evolutionary game in (14). For a specific population, for example, for α = 3/2 and , it can be simplified to and for α = 3/2 and , to

Note that at , , the potential maximization problem in (26) has only one active constraint y + y = 1, and the Hessian projected on the null space of this constraint will be If , if and only if This is a sufficient condition for an equilibrium strategy , , to be evolutionarily stable. For a specific population, for example, for α = 3/2 and , the condition can be simplified to: and for α = 3/2 and , to

An equilibrium strategy , , corresponds to a bilingual population with both language A and B co-existing. In previous sections, competition-only and intervention-only populations have been discussed. For the case with competition-only, such a strategy is never stable, and therefore, the two languages can never co-exist, even in bilingual forms. For the case with intervention-only, such a strategy is always stable, and the two languages can co-exist. By combining the two, it is hoped that such an equilibrium strategy still exists and is also stable. Indeed, such an equilibrium can be achieved by choosing appropriate interventional strategies:

Strategy 1: If , i.e., , the stability condition in (33) can be satisfied easily by choosing appropriate reversing rates and : Since , and , and therefore, the stability condition in (33) is satisfied if and are sufficiently large, say and , where 1 ≤ t ≤ min{1/s, 1/s}. With such a choice of and , one can prove that is also unique (see Appendix D for verification and Fig 3(a) for demonstration).

Fig 3

Dynamic behaviors with competition and intervention.

Payoff functions and are plotted against y and the changing directions of y are pointed with arrows. In (a), α = 3/2 and , and , and and . In (b), α = 3/2 and , and , and and .

Strategy 2: If , i.e., , also set and for any 1 ≤ t ≤ min{1/s, 1/s}. Then, , , is an equilibrium strategy for all three type of games, the competition-only game, the intervention-only game, and the combination of the two, i.e., the game with both competition and intervention. Mathematically, if and only if and . In addition, the stability condition in (33) is satisfied for . One can then prove that is evolutionarily stable and also unique (see Appendix D for verification and Fig 3(b) for demonstration).

Strategy 3: In general, given a desired equilibrium strategy , , not necessarily optimal for either the competition-only or intervention-only game, it can be made a true equilibrium strategy for the combined game by choosing appropriate interventional strategies: For , set Then, and . It follows that , and becomes an equilibrium strategy for the game. In addition, let , and assume that is selected such that or . Then, the condition in (33) is satisfied at , and the strategy is also evolutionarily stable (see Appendix D for verification and Fig 4 for demonstration).

Fig 4

Dynamic behaviors with competition and intervention.

Payoff functions and are plotted against y and the changing directions of y are pointed with arrows. In (a), α = 3/2 and , s = 0.6 and s = 0.4, and , and and . In (b), α = 3/2 and , s = 0.4 and s = 0.6, and , and and .

Note that in the last case, when is close to 1, the functions and are likely to have more than one intersections over the interval (0, 1) for y, as demonstrated in the examples in Fig 4. This implies that there may be more than one equilibrium strategies, among which is the desired one. In this case, the interval (0, 1) is divided by the values of these strategies, and the value of the selected strategy is converged only within a small subinterval around it. When is much smaller than 1, however, the functions and are more separated, and usually there is only a single intersection, corresponding to a unique equilibrium strategy.

The dynamic behaviors of the combined game in (14) can be demonstrated through examples shown in Figs 3 and 4. The graph in Fig 3(a) shows the behavior of a population with α = 3/2 and . The curves of and are displayed over y ∈ [0, 1] with and , when language A has social or economic advantages over B. If without intervention, the game would have an equilibrium strategy with and , but it is unstable. However, with intervention, the equilibrium strategy becomes with and . It is stable as y increases to if it starts from below while decreases to if it starts from above .

The graphs in Fig 3(b) shows the behavior of a population with α = 3/2 and . The curves of and are displayed over y ∈ [0, 1] with and . If without intervention, the game would have an equilibrium strategy with and , but it is unstable. However, with intervention, the equilibrium strategy remains to be with and , but becomes stable as y increases to if it starts from below while decreases to if it starts from above .

Fig 4 shows examples for how to obtain desired equilibrium strategies by choosing appropriate interventional strategies. The graph in Fig 4(a) is an example with α = 3/2 and . The curves of and are displayed over y ∈ [0, 1] with s = 0.6 and s = 0.4, when language A has social or economic advantages over B. By setting and , a desired equilibrium strategy with and is obtained. It is stable as y increases to if it starts from below while decreases to if it starts from above .

The graph in Fig 4(b) is an example with α = 3/2 and . The curves of and are displayed over y ∈ [0, 1] with s = 0.4 and s = 0.6, when language B has social or economic advantages over A. By setting and , a desired equilibrium strategy with and is obtained. It is stable as y increases to if it starts from below while decreases to if it starts from above . Note that in this example, there are three equilibrium strategies in (0, 1), as marked by half circles. The desired strategy is one of them. This happens when is close or equal to 1.

Dynamic simulation

Computer simulation can always be used to validate theories when physical experiments are not accessible, and to reveal additional insights into system behaviors which may otherwise be hard to see. This section describes further dynamic behaviors of bilingual populations through computer simulation with the combined model given in (10) and (14). To carry out the simulation, the population is assumed to be distributed in a two-dimentional space, and the changes of the bilingual level of the population will then be tracked across time and space.

More specifically, to carry out the simulation, a 2D torus-shaped lattice of n × n cells is constructed first, with each cell assumed to be occupied by an individual speaker. A random individual can then be selected repeatedly from the lattice, and a game is played for the individual against the population of the lattice. Let (x, x) be the current strategy for the individual, and (y, y) the strategy for the population. Let and be the payoffs for A and B speakers, respectively. Then, the payoff for the individual, , is computed. If p > π, x is increased by setting x = y if x < y. On the other hand, if p < π, x is reduced by setting x = y if x > y.

Initially, all individuals are assigned with a random strategy. The game is played n × n times for the population to complete a generation. The game is repeated for 100 generations to make sure the population reaches its equilibrium. In general, the game can be played in a neighborhood of each selected individual. Let the neighborhood be an m × m sub-lattice, with the selected individual located at the center. Then, the game can be carried out for each selected individual only against the population in its neighborhood of this size, with the population strategy (y, y) computed from the population in the neighborhood. Such a game may in fact be more realistic, as pan class="Species">people usually interact only with a small group of others around them.

Recall that the population is assumed to be large and well mixed in the game models discussed in previous sections. This means that every individual should be able to interact with all others in the population. With such an assumption, the dynamic behaviors and especially the equilibrium states of the population can be computed by directly solving a corresponding system of replicator equations. The 2D simulation described here can be considered as a discrete solution of the system of replicator equations if the individual in each cell of the lattice is allowed to interact with the individuals in all other cells. However, the 2D simulation is more flexible, allowing restrictions on interactions. It can therefore be used to see how a population behaves when each individual interacts only with a selected group of other individuals, which a continuous model would not be able to reveal.

Fig 5 shows the results from two sets of simulations in two columns of graphs, respectively. The results in the first column are for a population with α = 2 and , and and . Those in the second column are for a population with α = 3/2 and , and and . In both cases, the population is distributed on a 75 × 75 lattice. For each case, the game is played three times with neighborhood size equal to 75 × 75, 25 × 25, and 5 × 5, respectively. The final distribution of the individual frequency in the population is displayed for each play in the corresponding order.

Fig 5

Dynamic simulation results.

The distributions of the A speaking frequencies in the 2D lattice at equilibrium are displayed in graphs from top to bottom, with corresponding neighborhood sizes equal to 75 × 75, 25 × 25, and 5 × 5. Each of the graphs is a 75 × 75 2D lattice. The x-axis and y-axis of the graph represent the 75 units of the lattice in the horizontal and vertical directions, respectively. In column (a), α = 2 and , and , and and . In column (b), α = 3/2 and , and , and and .

In the first column of graphs, the one on the top shows the result from the game with the neighborhood size equal to 75 × 75, when each individual interacts with all others in the whole population. The equilibrium frequency of the population in this case is approximately equal to 0.75, which agrees with the direct prediction from the continuous model described in previous sections. In addition, the distribution of the individual frequency in the population is very homogeneous, with across the board, suggesting that language A and B co-exist in the population in an evenly distributed bilingual form. The second graph in this column shows the result from the game with the neighborhood size equal to 25 × 25, when the interactions among individual speakers are restricted. The equilibrium frequency of the population remains about the same, approximately equal to 0.75. However, the individual frequency becomes less constant. Some regions have higher individual frequencies than others, and local groups are formed with varying individual frequencies, as shown in the graph. The graph in the bottom of this column shows the result from the game with the neighborhood size further reduced to 5 × 5. While the population frequency is not significantly changed, the individual frequency shows even bigger variations, with even smaller local spots formed with higher or lower individual frequencies than average.

The graphs in the second column show similar results. Again, when the neighborhood size is set to 75 × 75, the equilibrium frequency of the population reaches approximately 0.4, which agrees with the direct prediction from the continuous model described in previous sections. The distribution of the individual frequency in the population is homogeneous, with across the board, as shown in the graph on the top. However, when the neighborhood size is reduced to 25 × 25, the individual frequency becomes less constant, although remains about the same. Some regions have higher individual frequencies than others, and local groups are formed with varying individual frequencies, as shown in the second graph. When the neighborhood size is further reduced to 5 × 5, the individual frequency varies in greater degrees in the population. Even smaller local spots are formed with higher or lower individual frequencies than average. The dynamic behaviors shown from these simulations agree with our intuition or experience in language development: Indeed, when communications are restricted to local groups, language variations will remain.

Discussion

The competition-only model as described in the Evolution with Competition section predicts that between two competing languages, if based only on competition, one would eventually die while the other takes over the whole population, and the co-existing state is unstable. This result is not so surprising, for it has already been discussed in previous studies, although from different perspectives. On the other hand, the intervention-only model as described in the Evolution with Intervention section shows that if controlled only by interventions, it is possible to prevent either language population from becoming too large or too small, and keep the population in a stable co-existing state. By combining the two, the model with both competition and intervention as described in the Evolution with Competition and Intervention section gives a more complete description on the evolution of multilingual populations when it is under the influences of both language competition and societal intervention. It predicts that languages may co-exist stably in multilingual forms if appropriate interventional strategies are employed. In addition, the interventional measures may not only be able to prevent language extinction but also direct populations to desired equilibrium states.

These predictions are based on the assumptions that the payoff functions P and P are increasing functions and and are decreasing functions, meaning that in some sense, competition favors large populations, while intervention subsidies small populations. For this reason, the functions in (16) and (17) are adopted for defining P, P, , and . Besides, P and P in (16) can be considered as a general set of functions that includes the one used in the Abrams-Strogatz model with α = 1.31 ± 0.25, which is well tested against real language data [13]. The functions for and in (17) are also related to some applications in genetic selection, where they are used to model genetic mutations with [6, 19]. However, the and values for genetic mutation are certainly in a much smaller scale than those for language conversion.

Note that this paper has not discussed in detail how the language is acquired, learned, and used, which may in fact relate directly to what language strategies really mean and how they can be changed. There are some obvious questions related to this issue. For example, what does it mean that an individual uses language A with 40% frequency? Is only reading language A without speaking counted as using language A? What happens if an individual can never learn a second language? How long does it take for an individual to change his/her language use from one frequency to another? These are legitimate concerns, but may require more linguistic characterizations on language learning and communication, and will hopefully be addressed in future work.

The work in this paper is to develop a general framework for modeling and simulating evolutionary dynamics of multilingual populations. Work remains to be done to analyze some real language groups, while several issues are yet to be resolved: First, some terms used in the proposed models need to be further discussed for exact implications such as “social or economic preferences”, which determine the parameters s and s, and “public policies”, “educational influences”, and “family influences”, which the parameters and depend on. Second, even if exact meanings of these terms are given, the corresponding parameters including α and are not easy to estimate. They may vary with varying populations and varying time periods. The language data is also difficult to collect across a meaningful period of time. Third, the factors that affect language evolution can be more than what a few parameters can cover. While the parameters used in the proposed models are key to determining general dynamic behaviors of multilingual populations, more parameters may be introduced, and additional mechanisms such as the reward and punishment schemes discussed in [23, 24] may be employed.

Computer simulation is a tool to extend the power of theoretical models. It can be used to not only validate the theory but also implement and test additional hypothesis such as the influences from the geographical or demographical structure of the population on multilingual evolution. The simulation done in this study is motivated by previous work on simulation of evolutionary dynamics such as Nowak and May 1992, Durrett and Levin 1997, and more recently, Chen, et al. 2008, Wang, et al. 2012, and Wang, et al. 2016 [25-29]. However, the simulation algorithms can be different for different models and even for the same model. Central to a simulation algorithm is the rule to update the current strategy. Different updating rules may result in very different outcomes. In this work, language competition is considered to be moderate. Therefore, if p > π, i.e., speaking A benepan class="Disease">fits more, then x will stay unchanged if it is already bigger than the average frequency y. It will increase to y only when it is less than y. The rule could be more aggressive such that x will increase by a certain amount even if it is bigger than y. The results would look different from what have been presented.

Research on language evolution has been pursued extensively in recent years including modeling and simulation. Work has focused on two different but related aspects of language evolution. One is on how a given language develops and evolves with respect to the “genetic changes” of its linguistic elements such as lexicon, syntax, semantics, pronunciations, etc [30-35]. The other is on how a mixed population of multilingual speakers changes dynamically in terms of the percentage of speakers of each language in the whole population [13–17, 36–40]. The work described in this paper belongs to the second category. Work on pan class="Disease">language death by Abrams and Strogatz 2003 has been followed and extended by several research groups [13–15, 17]. Other approaches have also been proposed, using Lofka-Volterra equations to characterize the evolutionary nature of language dynamics [36, 38], introducing reaction-diffusion models to include geographical influences on language competition [37, 39], and directly simulating the language dynamics with an agent-based or cellular automata models [16, 40]. The model proposed in this work shares some features of these previous approaches, but is focused more on the evolution of multilingualism and in other words, the co-existence of multiple languages in multilingual forms. It is intended to account for possible outside controls or interventions as well as inter-language competitions. Evolutionary games are defined with languages as competing strategies, and corresponding replicator equations are formed to characterize the dynamic behaviors of multilingual populations. The games are also recognized as a special class of potential games and therefore connected to a special class of potential maximization problems, which makes it possible to analyze and even control the equilibrium and stability conditions of given language populations.

To keep the description simple and easy to follow, the proposed models are introduced with and analyzed on bilingual populations only, but they are not limited to bilingual populations. They can actually be extended straightforwardly to general multilingual populations, as given in Appendix E. In order to make the paper more accessible, efforts have been made to keep the discussions from being too technical, with more formal definitions, theorems, and especially proofs all provided as appendices in the end of the paper. Nonetheless, these appendices are considered to be necessary contents as well for supporting the results and conclusions presented in the paper.

Appendix A: Evolutionary games [18, 19]

Let x, y ∈ S be the strategies of an individual and a given population, respectively, where S = {x ∈ R: Σ x = 1, x ≥ 0, i = 1, …, n}. Let π(x, y) = Σ x p(y) be the payoff function for the individual, where p is the payoff function for the ith pure strategist.

Definition 1 (Evolutionary Game). A Nash equilibrium for the evolutionary game with the payoff function π is a strategy x* such that Definition 2 (Evolutionary Stability). An equilibrium strategy x* is said to be evolutionarily stable if there is a positive number such that for all x ≠ x* and ,

Definition 3 (Replicator Equation). The following system of equations is called a system of replicator equations: where π(x, x) = Σ x p(x).

Theorem 1. The equilibrium strategy of an evolutionary game is a fixed point of the system of replicator equations. It is an asymptotically stable fixed point if it is an evolutionarily stable equilibrium strategy.

Theorem 2. An equilibrium strategy x* is evolutionarily stable if and only if there is a small pan class="Disease">neighborhood N of x* such that π(x*, x) > π(x, x) for all x ∈ N ∩ S, x ≠ x*.

Appendix B: Potential games [18, 20]

Definition 4 (Potential Game). A game defined by a payoff function π(x, y) = Σ x p(y) is called a potential game if there is a function f such that ∂f(y)/∂y = p(y) for all i = 1, …, n.

Definition 5 (Potential Maximization Problem). The following optimization problem is called a potential maximization problem: where S = {y ∈ R: Σ y = 1, y ≥ 0, i = 1, …, n} and ∂f(y)/∂y = p(y) for all i = 1, …, n.

Theorem 3. Let x* be a strategy for a given potential game. Then, x* is an equilibrium strategy for the game if and only if it is a KKT point of the corresponding potential maximization problem.

Theorem 4. Let x* be an equilibrium strategy for a given potential game. Then, x* is evolutionarily stable if and only if it is a strict local maximizer of the corresponding potential maximization problem.

Appendix C: Equilibrium strategies and stabilities

Cases for competition only

Consider the game in (2) with the payoff functions P and P being increasing functions as defined in (16). Let be an equilibrium strategy in one of the three cases given in (7).

Theorem 5. The equilibrium strategy is evolutionarily stable for case (a) and (b), but unstable for case (c).

Proof. For case (a), and . Let (y, y) be any strategy not equal to . Define . Then, . But , which is less than for ϵ sufficiently small. It follows that must be evolutionarily stable.

For case (b), and . Let (y, y) be any strategy not equal to . Define . Then, . But , which is less than for ϵ sufficiently small. It follows that must be evolutionarily stable.

For case (c), and . The only active constraint of the problem in (20) at this strategy is y + y = 1. The null space of the constraint at this strategy can be represented by a basis matrix z = (1, −1). Then, the projected Hessian on this space at is , which is always positive definite. Therefore, by the second order necessary conditions for constraint optimization, can never be a local maximizer of the problem in (20), and therefore, can never be evolutionarily stable.

Cases for intervention only

Consider the game in (5) with the payoff functions and being decreasing functions as defined in (17). Let be an equilibrium strategy in one of the three cases given in (8).

Theorem 6. The equilibrium strategy is evolutionarily unstable for case (a) and (b), but stable for case (c).

Proof. For case (a), and . Let (y, y) be any strategy not equal to . Define . Then, . But , which is greater than or equal to for ϵ sufficiently small. It follows that must be evolutionarily unstable.

For case (b), and . Let (y, y) be any strategy not equal to . Define . Then, . But , which is greater than or equal to for ϵ sufficiently small. It follows that must be evolutionarily unstable.

For case (c), and . The only active constraint of the problem in (23) at this strategy is y + y = 1. The null space of the constraint at this strategy can be represented by a basis matrix z = (1, −1). Then, the projected Hessian on this space at is , which is always negative definite. Therefore, by the second order sufficient conditions for constraint optimization, must be a strict local maximizer of the problem in (23), and must therefore be evolutionarily stable.

Appendix D: Dynamics of interventional strategies

Consider the game in (14) with the payoff functions and defined in (15), where P and P are given in (16) and and in (17). Let , , be an equilibrium strategy or a desired equilibrium strategy for the game. This equilibrium strategy may or may not be stable and unique depending on the choice of interventional strategies and .

Theorem 7 (Stability and Uniqueness for Strategy 1). Let and , 1 ≤ t ≤ min{1/s, 1/s}. Let , , be the corresponding equilibrium strategy. Assume that . Then, is evolutionarily stable and also unique.

Proof. Since , and . Then, the stability condition in (33) is satisfied at for all t, 1 ≤ t ≤ min{1/s, 1/s}. Therefore, is evolutionarily stable. Note that the same condition is also satisfied at all strategies (y, y).

Let ϕ be the function for the difference between and , . Then, . Note that . Then, if the condition in (33) is satisfied at all strategies (y, y), must be negative, and ϕ(y, y) is monotonically decreasing in y. It follows that must be unique.

Theorem 8 (Stability and Uniqueness for Strategy 2). Assume that . Let and , 1 ≤ t ≤ min{1/s, 1/s}. Let , , be the corresponding equilibrium strategy. Then, is unique and also evolutionarily stable.

Proof. Since is an equilibrium strategy, , i.e., This equation can be written equivalently to It can then be re-arranged to which implies that and hence , i.e., and . Reversely, it is also true.

Since there can only be one strategy such that and , is unique.

The uniqueness of can also be justified by noticing that for any strategy (y, y), can be written in the following form, which is equivalent to Note that , , and is negative if and positive if . Therefore, ϕ(y, y) is positive if and negative if , and is unique.

For stability, it suffices to show that the stability condition in (33) is satisfied at for and . Note that Since , It follows that . Similarly, Since , It follows that . The stability condition in (33) is then satisfied and is evolutionarily stable.

Theorem 9 (Stability for Strategy 3). Let , , be a given strategy. Let and . Assume that . Then, is an equilibrium strategy, and is also evolutionarily stable if it is selected such that or , where .

Proof. Given and , it is easy to verify that and , and , proving that is an equilibrium strategy.

Define a function ψ for the difference between the left and right hand sides of the stability condition in (33) without the constant terms and (α − 1): Let , , and substitute and into the function to obtain: The formula can be further simplified to Note that for . Then, since P(y) − P(y) is monotonically increasing in y, is positive if and negative if . Then, if is selected such that or , . Since , it follows that the stability condition in (33) is satisfied at and is evolutionarily stable.

Note that for Strategy 3, the equilibrium strategy , , may not always be unique, especially if is close to α − 1. For example, if , it is easy to verify that and , and , which means that and is also an equilibrium strategy of the game, in addition to the equilibrium strategy and . It can also be proved evolutionarily stable: For and , and if and , which can be simplified to Note that for . Then, since is monotonically decreasing in y, is negative if and positive if . Then, if is selected such that or , . Since , it follows that the stability condition in (33) is satisfied at and is also evolutionarily stable.

Appendix E: General models for multilingualism

Assume that there is a multilingual population of m languages. Let x = (x1, …, x) be the strategy vector for an individual speaker, where x is the frequency of the individual to speak language i, Σ x = 1. Let y = (y1, …, y) be the strategy vector for the population, with y being the average frequency of the population to speak language i, Σ y = 1. Thus, an x-speaker is an individual speaking language i with frequency x, and a y-population is a population speaking language i with average frequency y.

Assume that the evolution of the given population is influenced only by competition on population and social or economic preferences. Let P(y) be the payoff function for the individual who speaks only i when the average frequency of speaking i in the population is given by y. Assume that P, i = 1, …, m, are monotonically increasing functions. Then, the payoff function π for an x-speaker in y-population can be defined by the average payoff, π(x, y) = Σ x P(y), with which an evolutionary game can be formulated. A Nash equilibrium of this game is a strategy x* such that π(x*, x*) ≥ π(x, x*) for any strategy x. A corresponding system of replicator equations can be defined as

Similarly, assume that the evolution of the given population is influenced only by possible societal interventions. Let be the payoff function for the individual who speaks only i when the average frequency of speaking i in the population is given by y. Assume that , i = 1, …, m, are monotonically decreasing functions. Then, the payoff function for an x-speaker in y-population can be defined by the average payoff, , with which an evolutionary game can be formulated. A Nash equilibrium of this game is a strategy x* such that for any strategy x. A corresponding system of replicator equations can be defined as

Let , i = 1, …, m, and correspondingly, , where λ ∈ [0, 1]. Then, a general evolutionary game can be formulated for the population with the influences of both language competition and societal intervention. A Nash equilibrium of the game is a strategy x* such that for any strategy x. A corresponding system of replicator equations can be defined as The payoff functions P and can be defined more specifically such as , 1 < α ≤ 2 and , , where c and are scaling constants, s are indicators for social or economic preferences, and are language reversing rates due to interventions, .

17 Aug 2020

pan class="Chemical">PONE-D-20-05223

Why Multilingual, and How to Keep It -- An Evolutionary Dynamics Perspective

PLOS ONE

Dear Dr. Wu,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

Please submit your revised manuscript by Sep 19 2020 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

Please include the following items when submitting your revised manuscript:

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

We look forward to receiving your revised manuscript.

Kind regards,

Xiaojie Chen

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. Please consider referring to yourself at the first person of the singular ("I") rather than the plural ("we") in this single-authored paper.

3. We suggest you thoroughly copyedit your manuscript for language usage, spelling, and grammar. If you do not know anyone who can help you do this, you may wish to consider employing a professional scientific editing service.

Whilst you may use any professional scientific editing service of your choice, PLOS has partnered with both American Journal Experts (AJE) and Editage to provide discounted services to PLOS authors. Both organizations have experience helping authors meet PLOS guidelines and can provide language editing, translation, manuscript formatting, and figure formatting to ensure your manuscript meets our submission guidelines. To take advantage of our partnership with AJE, visit the AJE website (http://learn.aje.com/plos/) for a 15% discount off AJE services. To take advantage of our partnership with Editage, visit the Editage website (www.editage.com) and enter referral code PLOSEDIT for a 15% discount off Editage services.  If the PLOS editorial team finds any language issues in text that either AJE or Editage has edited, the service provider will re-edit the text for free.

Upon resubmission, please provide the following:

The name of the colleague or the details of the professional service that edited your manuscript

A copy of your manuscript showing your changes by either highlighting them or using track changes (uploaded as a *supporting information* file)

A clean copy of the edited manuscript (uploaded as the new *manuscript* file)

Additional Editor Comments (if provided):

The author should revise the manuscript according to the reviewers' comments.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Partly

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: N/A

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. pan class="Species">participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: In this MS, the author propose an evolutionary dynamic model to study the evolution of multilingualism. The model consists two different parts. Concretely, the first one relates to selection of languages based on competition and the second one relates to circumstances when selection of languages is altered. Through theoretical analysis and numerical simulation, the results show that the stable co-existence of languages is possible and extinction can be prevented by choosing appropriate interventional strategies. Furthermore, the author gives the stability conditions of equilibrium states for different interventional strategies. Overall, this is a nice and well written paper, tackling an interesting question in the evolution of cooperation. The mathematical analysis is done carefully (although I have to admit that I am unable to check all formulas) and the results are presented in an understandable way.

However, there are some remaining issues with the manuscript, requiring some answers.

1) The author respectively study the evolutionary dynamics in competition-only and intervention-only populations. Furthermore, the author also combines these two and explores the stability of the equilibrium points. He/she make a combination of the systems in (3) and (6) to obtain a new system. Can the author explain in detail why the system can be described by equation 9.

2) In the first part of the manuscript, the author studies the evolutionary dynamics of multilingualism in well-mixed populations, while in the latter part, he/she studies the evolutionary dynamics of multilingual on the 2D torus-shaped lattice. Here, I want to know the relationship between these two, because individuals in the lattice are not well-mixed.

3) The author study how multilingualism may evolve under the language competition and societal interventions. What does societal interventions usually mean here? What is the relationship between this kind of social intervention and some social incentive mechanisms? These social incentive mechanisms include punishment, reward, or exclusion such as Chen et al. PLoS Computational Biology, 14(7), e1006347 (2018), Liu et al. Nonlinear Dynamics, 97(1), 749-766, (2019).

4) The figures in the manuscript should be more exquisite. Concretely, some circles are incomplete. In addition, stable and unstable equilibrium points cannot be distinguished. What does the axis in Figure 5 represent?

5) There are some grammatical problems in the manuscript that need to be further checked. Besides, the format of some references needs to be adjusted, for example, [8], [31-32], [34], [36].

Reviewer #2: The paper investigates the evolution of bilingualism in the framework of evolutionary game theory. The author constructs an evolutionary game-theoretical model between two languages by considering not only the inner competition between them, as Abrams and Strogatz do (Nature 424: 900, 2003), but also the social interventions that are imposed by the external factors such as public policies, education, or family influences. By combining competition with intervention, they show both languages can be stably coexisted in the well-mixed populations. Besides, the author also performs computer simulations to explore the evolutionary dynamics on spatial lattices.

In general, I think this is a good paper. Particularly, it complements the results previously found by Abrams and Strogatz, that is, the evolutionary competition between two languages leads to the extinction of either language. The findings of this paper clearly show us that how bilingualism may evolve under the influences of both competition and intervention.

I didn't check the algebraic manipulations in detail, but the general methods employed seem to be sound and elegant.

I do have some suggestions as to how the style of presentation can be improved (see below). Once these suggestions are incorporated, I support publication.

(1) The title of the manuscript is "Why Multilingual, and How to Keep It -- An Evolutionary Dynamics Perspective". But I find that the author merely investigated the dynamic behaviors of bilingual population. If the author insists on usage of the present title, I think the author should extend the bilingual model to multilingual model, and mainly study the dynamic behaviors of multilingual population.

(2) It seems to me that the author assume that language competition and social intervention have equal impact in the evolution of bilingualism due to the same weights assigned to them in Eq. (9). This maybe a too strong assumption in my opinion. I think it is necessary to explore what happens if they have different impacts in the evolution of bilingualism? Is the present conclusion robust against such changes? If the answer is "No", what are the new results?

(3) With reference to the computer simulations on evolutionary game dynamics in spatial networks, I suggest citing the following papers: Phys. Rev. E 78, 051120 (2008) and Sci. Rep. 2, 740 (2012), wherein the authors used computer simulations for exploring the spatial dynamics of evolutionary games.

Reviewer #3: General evaluation:

My major concern with the paper is as follows: I think it is too theoretical to be published in a general journal such as PLOS ONE. In the current form, I would suggest to submit it to another journal where this kind of article type is common. Moreover, I think that there is a huge gap between empirical data and the model's predictions. I understand that a model with such a high degree of abstraction from reality is hard to be checked against empirical data. However, here even a preliminary attempt is missing (as the author mentioned himself in the Discussion section). Finally, I think there are also some at least questionable aspects of the model definition with respect to how languages are actually acquired, learned or competitively/cooperatively used in real societies. It is not clear to me that the model particularly defines language use. What are the language-specific aspects? It seems to me as if A and B could also be other cultural or social traits and customs, which can be adopted or abandoned, used to initiate competition or cooperation, and supported by societal interventions.

More detailed points:

line 24/25: “It assumes that the speakers are free to choose among available languages.” I think this is a unrealistic assumption. Speakers mostly learn a new language due to their (new) social environment, for example as a result of migration. They are mostly forced to learn the new local language to be able to participate in the social environment. Here it sounds as if speakers can change their language from one day to the other as a free choice.

line 81: “...that the larger the population percentage of a language, the more competitive the speaker of the language.” Why are they more competitive? It there empirical evidence for that? Shouldn't the payoff be dependent either i) on (mutual) intelligibility ii) or (probably more important here) on the social and economic benepan class="Disease">fits of using a particular language? I am not convinced that majority alone is sufficient for a language to be more competitive per se.

line 107/108: “...meaning that the smaller the population percentage of a language, the more incentive the speaker of the language receives.” What does the payoff actually stands for? Here it represents the incentive a speaker receives. Is this also the case for the competitive scenario? It is not clear at all what the payoff stands for in general. Maybe the author makes this point more clear to the reader.

Again, considering the high degree of abstraction of the model from the real world, it would be good to see how the model predicts real world phenomena. In the introduction, the author points to the many places in the world where multilingualism exists. But is there empirical evidence that interventions are the (unique) reason for its existence/perpetuation? Is it possible to get data about language competition and societal interventions, etc, that the model can be tested against?

Spelling:

Abstract: and playing -> and is playing

Abstract: accounts for selection -> accounts for the selection

In general: the article 'the' is dropped very often: (the) selection, (the) evolution, etc...

line 56: in home -> at home

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

Reviewer #3: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

18 Sep 2020

In the submitted files: Response to the Reviewers.

Submitted filename: PLOS Evolution of Multilingualism Response Letter.pdf

Click here for additional data file.

26 pan class="Gene">Oct 2020

Why multilingual, and how to keep it -- An evolutionary dynamics perspective

pan class="Chemical">PONE-D-20-05223R1

Dear Dr. Wu,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Xiaojie Chen

Academic Editor

PLOS ONE

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #2: All comments have been addressed

Reviewer #3: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: N/A

Reviewer #2: Yes

Reviewer #3: N/A

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. pan class="Species">participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: The author has carefully and thoughtfully revised this paper and it is now suitable for publication.

Reviewer #2: (No Response)

Reviewer #3: One major concern of my former review was the degree of abstraction of the current study and if it pan class="Disease">fits to a general journal such as PLOS ONE. The author pointed out covincingly that the topic in question - the modeling of multilingual populations - appears to be of high interest for readers in general journals such as PLOS ONE. Moreover, I think that the author treated all my further comments very carefully, which makes the paper much more comprehensible and gives it a much better flow of reading. Finally, the spelling is clearly improved.

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

Reviewer #3: No

28 pan class="Gene">Oct 2020

pan class="Chemical">PONE-D-20-05223R1

Why multilingual, and how to keep it – An evolutionary dynamics perspective

Dear Dr. Wu:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Professor Xiaojie Chen

Academic Editor

PLOS ONE