Dynamics of Public Opinion: Diverse Media and Audiences’ Choices

: Studies on the fundamental role of diverse media in the evolution of public opinion can protect us from the spreading of brainwashing, extremism, and terrorism. Many fear the information cocoon may result in polarization of the public opinion. Hence, in this work, we investigate how audiences’ choices among diverse media might influence public opinion. Specifically, we aim to figure out how peoples’ horizons (i.e., range of available media) and quantity, as well as the distribution of media, may shape the space of public opinion. We propose a novel model of opinion dynamics that considers different influences and horizons for every individual, and we carry out simulations using a real-world social network. Numerical simulations show that diversity in media can provide more choices to the people, although individuals only choose media within the bounds of their horizons, extreme opinions are more diluted, and no opinion polarizations emerge. Furthermore, we find that the distribution of media’s opinions can effectively influence the space for public opinion, but when the number of media grows to a certain level, its effect will reach a limitation. Finally, we show that the effect of campaigns for consciousness or education can be improved by constructing the opinion of media, which can provide a basis for the policy maker in the new media age.


Introduction
. In the modern age, media (especially new media) play significant roles in changing people's thoughts and ways of thinking. Human minds are complex and not rational in most situations, and their opinions are constantly influenced by the input of acquaintances and media such as TV, newspaper, YouTube, and Facebook. In the past, broadcast media featuring one-to-many information acted as gatekeepers of the information flow and primary vehicles for the messages of political elites or authorities to reach people (Quattrociocchi et al. ). Nowadays, with the continuous popularization of the internet, especially mobile internet, each individual has the ability to publish information and give opinions, which rapidly changes the means of spreading information and opinion. With the emergence of such plurality in media, the decentralization of media has become apparent, and new communication technologies continue to grow more attractive to the former audiences of traditional media (Jenkins ; Noll & Price ). Today, the media compete for people's attention, feeding (as much as possible) each individual a steady stream of content built for and seen only by them. The people of the real world may have become more deeply divided and fragmentary than ever before. .
Information cocoon is a well-known concept, it describes the phenomenon in which individuals only choose the information sources and contents they like, eventually a cocoon is formed to wrap themself (Sunstein ). Similar concepts are proposed by scholars continually, such as echo chamber (Jamieson & Cappella ), filter bubble (Pariser ), and information enclaves (Weeks et al. ). Those works are reflections of concerns about the opinion polarization brought by the diverse media. For instance, the High Level Expert Group on Media Diversity and Pluralism state the increasing filtering mechanisms may shake the foundations of democracy (Vīke-Freiberga et al. ).Hence, it is necessary for researchers to study public opinion formation when media have more diversity and people have more choices. .
Beyond these basic models of opinion dynamics, researchers have proposed several creative models. Based on the Social Judgment Theory, Jager & Amblard ( ) first introduced repulsions between agents to the traditional bounded confidence model. On this basis, Fan & Pedrycz ( , ) proposed social judgment-based opinion (SJBO) model. The continuous opinion and discrete actions (CODA) model (Martins ), meanwhile, is the name of a hybrid model in which agents have continuous opinions but their behaviors are discrete. Subsequently, by introducing the concept of trust to this model, the same author (Martins ) discussed the issue of trust in society. Although most people potentially enjoy the convenience of interaction that online technology provides, some people can only interact or get information o line, and there is also someone who chooses to do so o line. Hence, Dong et al. ( ) designed an online-o line model to unveil the interaction mechanism among the agents between the online and o line social networks. It reveals that online agents can smooth opinion changes in the opinion dynamics and decrease the number of opinion clusters. .
Additionally, as a useful tool for modeling and studying social opinions, opinion dynamics theories have been applied in various fields of scientific research. Topics include elections (Moya et al. ), the spread of extremist opinions (Fan & Pedrycz ), cyber violence (Liu et al. ), and even cognitive dissonance (Li et al. ).
. For several years now, researches have been focused on understanding the role of media (Hu & Zhu ; Martins et al. ; Pineda & Buendía ; Quattrociocchi et al. ) and the web (Ugander et al. ; Zhang et al. ) in the formation, evolution, and di usion of public opinions. Quattrociocchi et al. ( ) studied the formation of the agents' opinions under the influence of mass media and experts, based on an Italian political campaign in . Martins et al. ( ), who had introduced the repulsion behavior to the De uant-Weisbuch model, then studied the reaction of this model under external information, arriving at the conclusion that the repulsion behavior can promote people to reach consensus under the influence of external information. Hu & Zhu ( ) proposed an analyzable model, based on the multi-state voter model, to study the e ects of mass media on public opinions, with simulations showing that media influence is amplified by the interactions among agents within social networks. However, a majority of these works are focused on the impact of the mass media, and less attention has been dedicated to exploring how diverse media and audiences' choices might a ect the course of public opinions.

.
Considering previous researches by scholars from social and computer science, this paper aims to study the process by which opinions spread in a closed social network with various information sources (i.e., media) and allowing individuals choice on how to obtain information. The main objectives of this paper will be achieved by answering the following three questions.
• When individuals can choose their favorite media, will this leads to opinion polarization?
• How do the quantity and distribution patterns of the diverse media a ect the spread of the opinions?
• How can campaigns for consciousness or education be improved through specific media constructions?
. Therefore, we propose a novel model, social judgment and influence based opinion (SIBO), to describe how agents adjust their opinions a er interacting with other agents or with media. Multiple media holding various opinions are treated as external information sources bringing opinions to the closed social network, and individuals have the freedom to choose the media they like. The simulations are based on a real-life social network instead of on artificial data sets.
. The rest of this paper is organized as follows. In the following section, we present a new model of opinion dynamics and detailed simulation plans. Then, our simulation results and related discussions are provided, and the conclusion is summarized at the end.

Model
. We model the social network using agent-based modeling (ABM), proven for describing complicated sociotechnical systems such as Facebook and Twitter (Alvarez-Galvez ). In the model, each individual or medium is considered as an agent, and their connections are simplified as links or edges. The social network and information network are distinct two networks as modeled in our research, shown in Figure . As in the real world, the social network is not fully-connected, and agents are connected by their relationships whether online or o line. As Figure (b) shows, in today's media era, the form and content of the media have become more diverse. Considering that an individual chooses their information sources based on their own opinion (An et al. ), the individuals can and will choose specific media contents and forms. Consequently, in our study, the individuals' options are viewed as changing with e ects from two sources: acquaintances in the social network and media in the information network. .
Inspired by the SJBO model (Fan & Pedrycz , ), which can be seen as an extension of the widely used DW model, we present our model as the social judgment and influence based opinion (SIBO) model. A schematic diagram of our model is presented in Figure ; in the model, each agent's opinion is a continuous value in the range of [+1, −1]. In the opinion space [+1, −1], an opinion value of 0 means that the agent is entirely ignorant or neutral on one topic, while a positive value shows that the agent holds a positive or supportive opinion, and a negative value shows the opposite. A er agent i communicates with agent j from their social network at time t, the opinion of the agent i will shi according to: where x t i and x t j stand for the opinions of agent i and j at time t, respectively. Meanwhile, x t+1 i is the opinion of agent i at time t + 1 , µ ij , ε i and τ i are the influence coe icient, convergence threshold and repulsion threshold for agent i, respectively. From the Equation , we can see that the opinion of agent i is attracted to the persuader when their opinions' di erence is smaller than the convergence threshold. However, di erent from the bounded confidence model, when the opinion di erence is greater than the repulsion threshold, the opinion of agent i shi s away from the persuader's opinion. Except for the above two cases, the agent i sticks to its own opinion.

.
The influence coe icient µ ij is inspired by the concept of universal gravitation, which reflects the force between two agents that draws the two agents toward each other or push each other away. The value of µ ij is proportional to the influence of the persuader and inversely proportional to squared opinion di erence. Hence, the definition is given by: where m j stands for influence level of agent j, and its value is equal to the total number of edges for this agent (i.e., the total number of friends it has); the hyperbolic tangent is used as a saturation function to limit the value. In fact, the influence coe icient contains the information from the social network itself, which has gone ignored in previous models. By using this information, the agent with more influence can be more persuasive and has more impact on the opinion dissemination; such individuals are usually called opinion leaders.
. The initial definition of opinion leaders in social opinion dynamics was proposed in Katz and Lazarsfeld's book (Katz & Lazarsfeld ) as "the individuals who were likely to influence other persons in their immediate environment". This concept is a good indication of the asymmetry influence of the individuals in a social network. Many papers have appeared to analyze the function of opinion leaders in the context of opinion dynamics (Roch ; Watts & Dodds ; Chen et al. ) based on di erent opinion fusion rules. Meanwhile, the total number of an agent's edges, i.e., the degree of a node in graph theory, is the most intuitive information of an agent's influence within a network without considering the global structure of the graph (Bamakan et al. ). That is the reason why the value of µ ij and m j have a positive correlation. .
Similarly, when the agent i chooses a medium at time t, its opinion will be updated according to: where G k is the opinion of the medium k which as chosen by agent i. Although not all the media have the same degree of influence level, for simplicity, we set the coe icient to the maximum in Equation .
In addition, we introduce another parameter η i ∈ [0, 2] to the model, and it stands for the horizon of an agent. The agents will only choose media in terms of their horizon (i.e., their range of choices) in getting information according to our model: this phenomenon is becoming much more common in the new media era. In the model, the convergence threshold ε i always less than or equal to the repulsion threshold τ i , while the value of horizon η i has no such limitation, it can be less than ε i or greater than τ i . .
The repulsion threshold, as well as the convergence threshold, are inspired by the Social Judgment Theory (Sherif & Hovland ) which has been utilized in many articles. Those two confirmation bias related parameters determine how individuals change their opinion a er being confronted with another opinion. As Axelrod mentioned in his work (Axelrod ), in addition to beliefs, attitudes, and behavior, there are still more things over which interpersonal influence extends, such as language, art, etc. Therefore, the thresholds are representational parameters of internal emotion, and they are determined by factors such as education, growth environment, and so on. On the other hand, the parameter horizon we designed is to explore how the individuals' behavior of choice a ects the opinion dynamic. Hence, it represents a kind of active external behavior. Nowadays, people know a good deal about the media environments, and reflect upon how they use those environments, even provide a rational account of their actions (Webster ). For instance, one individual can try to use TikTok or Clubhouse just for curiosity, but this person may still feel hard to accept the information these new media provide due to its bias. Therefore, the repulsion threshold and horizon are conceptually di erent.

Simulation implementation
. Although the opinions in the model are continuous, the behavior of agents to obtain information and communicate with each other is discrete. Hence, in the simulation, the agents adjust their opinions in discrete timesteps, with the description of the implementation process being depicted in Figure . In each timestep T , every agent chooses a medium in terms of their horizon and then interacts with the medium first. A erward, considering the asynchronization in online and o line interactions (Ding et al. ), we consider that the agent would communicate with three random acquaintances (instead of one) in the social network. A er each interaction, the agents adjust their opinions according to Equations and . .
In this study, we use a real-life social network, rather than an artificially generated one, for simulation. The social network is a Facebook friendship network of Haverford College, Pennsylvania (Traud et al. ). As shown in Figure , this network has , agents and , edges, and the average degree is up to . (i.e., each individual has an average of . friends). Considering the inter-individual diversity, the agents' convergence thresholds ε i are uniformly distributed in the range of (0, 2), and the repulsion thresholds τ i are uniformly distributed in the range of (ε i , 2). Moreover, because agents change their opinions based only on di erences in opinion, it does not matter what or who the interaction object is. Therefore, the agents have the same threshold for convergence and repulsion toward other agents including media. When the public has a certain level of understanding of something, individuals generate all kinds of opinions. Hence, without loss of generality, the initial opinions of the agents are distributed uniformly over (−1, 1). In the simulation, the interaction cycle repeats times, i.e., T = 500. Considering the randomness of the whole process, the dynamic simulation will be carried out times with di erent possible initial configurations.

Simulation Results and Discussion
Simulation scenarios .
In order to answer the three questions given in the introduction, simulation scenarios based on di erent conditions are given in Table . Those scenarios are designed based on the one-factor-at-a-time methodology (ten Broeke et al. ), which reduces the influences on the result from other parameters by modifying as few parameters as possible.

Scenario Changing Parameters
Question to Answer Horizon of the agent How audiences' choices a ect opinion dynamics?
Opinion of one authoritarian medium How media opinion a ects opinion dynamics?
Quantity and distribution patterns of the media How the quantity and distribution patterns of the media a ect opinion dynamics?
Distribution patterns or changing opinion of the media How to guide the public opinion? It is clear enough in our times that a consensus of public opinion is rarely achieved. Thus, a metric is needed to reflect the degree of dispersion in public opinion. Although the number of opinion clusters is widely used, the algorithm that calculates this number o en requires iteration involving substantial computational costs. Therefore, we introduce a metric called the participation ratio (Nagel et al. ), borrowed from physics, and where x i is the final opinion of agent i and n is the total number of agents, while P is a scalar in the range of (0, 1]. P = 1 means that all the agents reach a consensus, and its value diminishes with increased diversity of discrete public opinions. In short, P −1 is approximately equal to the number of opinion clusters but requires less computation.

Scenario .
In this scenario, the influence of individual choice behavior is discussed. Hence, the number for media is only one: one authoritarian medium broadcasting an unchanged opinion of 0.6 (i.e., G = 0.6). The reason for choosing this opinion value will be given in our presentation of the next scenario. .
The horizon value η i is the same for all agents, set from 0 to 2 with intervals of 0.2. Besides, a condition has been added such that agents' horizon values are between each one's convergence and repulsion thresholds, which is more in line with practical situations. The simulation results are depicted in Figure , where H stands for the additional condition. To reveal more information from the random results, the box plot is utilized instead of the simple mean value plot. As demonstrated in Figure (a) and (b), the overall trend is that the average opinion increases and the participation ratio decreases with higher horizon values. The reason is not hard to understand: in a close-minded social network with a low horizon value, agents only get information from friends and media that have very similar opinions to theirs. In such a society, the agents with extreme opinions (−1 or +1) are more committed because they refuse to choose and obtain information from the media di erent from them. As a result, individuals have more potential to be attracted by extreme opinions, and they tend to bring out polarization as shown in Figure (c). To be clear, in order to reduce the randomness, the distribution graph is based on all the agents' opinions in the repeated simulations. . When η i < 0.6 , the final average opinion value contains greater uncertainty. Then, once η i reaches 0.6, the public opinion is influenced by the propaganda most of the time, but due to low horizon values, there remains a greater possibility that the propaganda repulses public opinion toward the negative side. This phenomenon can be seen from the crossed outlier in Figure (a). In addition, we can obtain a counterintuitive phenomenon whereby the public opinion decreases slightly a er η i > 0.8 . Analyzing opinion distribution, we find that this is because the value for the horizon is higher than most agents' repulsion threshold, which leaves some agents exposed to information they repulsive, pushing them to take attitudes of opposition. As for the additional condition H, with random horizon value between two thresholds, the average opinion of agents is relatively high. .
In summary, an open-minded social network with a high horizon improves the spreading of the media's opinions, but if the society is too open, the e ectiveness of authoritarian propaganda will be lower. Moreover, in this scenario, no matter whether society is close-minded or not, there are more than % of extremists exist.

Scenario
. Although many scholars have discussed the impact of mass media on public opinion, it is still necessary to study it in this paper because of the appearance of the audiences' choices. In this scenario, the opinion of the one medium varies from −1 to 1, and the agents' horizon values are between each one's convergence and repulsion thresholds. The results are shown in Figure . It is obvious that the medium has a strong impact on public opinion, as seen in Figure (a), but as the opinions have become extreme, the influence of the medium has declined slightly. When the opinion of the medium reaches 1 or −1, not only the average opinion can only reach 0.6 or −0.6, but opposed public opinions emerge in some cases. From Figure (b) it can be seen that the participation ratio decreases as the opinion's radicalness decreases, until the opinion is totally neutral. Even in the case where the medium is totally neutral, % agents hold extreme opinions (1 or −1), as shown in Figure  (d). .
The phenomenon of reduced media influence in (a) is related to agents' horizon values. In order to illustrate it, the simulation result of changing horizon with one medium hold an extreme opinion (G = 1) is given in (c). As observed in (c), the appearance of opposed public opinion diminishes as agents have wider horizons. In addition, compared with (a), extreme media opinions have a higher possibility of generating opposed public opinion; that is the reason for our choice of G = 0.6 in Scenario . .
In conclusion, the authoritarian media play an essential role in propaganda, which can change public opinion e ectively, and its influence can be mitigated with the development of open minds (i.e., high horizon values). Additionally, extreme media opinions lead to polarization, and those agents who were pushed to be opponents lowered the average opinion. Extreme media opinions also show greater potential to generate collective opposing public opinion in a close-minded society.

Scenario .
In this scenario, we explore how the quantity and distribution patterns of media a ect opinion dynamics. For this objective, the following simulations present multiple media with di erent distributions. We set the numbers for the media to , , , , and , accounting for %, %, %, %, and % of the total number of agents, respectively. To better reflect reality, we consider three di erent distributions of these media.

.
In the first case, the opinions of the media are evenly distributed in the space of public opinion with [−1, +1], which is highly idealized. Figure presents public opinion under di erent numbers of media with uniform distribution. It is important to note that a new metric called the Euclidean norm is introduced here. The Euclidean norm is related to the Euclidean distance used to describe the straight-line distance between two points. In our work, the Euclidean distance can indicate the distance between the public opinion and the average opinion of whereḠ is the average opinion of media, and in this caseḠ = 0. . Figure (a), although there is a great deal of randomness, the average opinion maintains values in the range of ±5 because the average opinion of the media is zero. In fact, the diversity of the media distracts the agents and leads to randomness and disorder, but the di erence in the quantity of the media makes little di erence. Furthermore, compared to the situation under one medium, both the average opinion and the Euclidean norm are more random under multiple media. Comparing the distribution in Figure (c) and (d) with Figure (d), it is evident that opinion is more evenly distributed on the opinion space because the agents have more choice among media by which to obtain information. In general, uniform distributed media bring strong randomness to the social network, hence breaking down the information monopoly of one authoritarian medium and also reducing the percentage of extreme opinions.

As observed in
. In case , the opinions of the media are normally distributed, with their mean set to 0 and the standard deviation at 0.3, i.e., G ∼ N 0, 0.3 2 . Although these parameters are carefully chosen, there is still a chance that the opinions will be greater than 1 or less than −1, so a saturation function should be used to maintain the opinions in [−1, +1]. Figure shows public opinion under di erent numbers of media with the normal distribution. It can be seen that, being di erent from case , normally distributed opinions of the media bring some level of orderliness to the outline of public opinion. First of all, average opinion is converged around ±0.2, more compact than in case , and the randomness is decreased with the increasing quantities of media.
. Figure (b), with more media available, the Euclidean norm is decreased significantly, but the trend slows down a er media have reached. In fact, due to the properties of the normal distribution, there is a high possibility that most opinions of the media are near neutral instead of extreme, leading to similar distributions in public opinion as shown in Figure (c) and (d) except with extreme opinions. Compared to case , this case not only reduces the extreme opinions, but also concentrates that of neutral opinion in a higher density, even more than in the case where only one medium is available (see Figure (d . In case , a special distribution is utilized in the simulation; we call it reverse-normal distribution (this is not a mathematical term). As the reverse of the normal distribution, most of the media hold extreme opinions instead of neutral opinions, and the opinions are obtained using the equation G = sgn(G) − G, G ∼ N 0, 0.3 2 . The e ect of this distribution on public opinion is depicted in Figure . It can be observed that although the average opinion is still around 0, the Euclidean norms mark a sharp increase with multiple media. Moreover, in Figure (c) and (d), the distribution of public opinion is similar to the opinion distribution characterizing the media's input to the system. Although extreme opinions predominate in social opinion, compared to the situation with one authoritarian medium that holds an extreme opinion (i.e., G = 1 or −1), various opinions are still existing and scattered on the opinion space in this case. . Combining all three cases above, our conclusions are obtained as follows. ) Compared with the situation under a single medium, a variety of media will bring choices for agents, thereby reducing the attraction of extreme opinions and leading to dispersed opinion distribution. It is a counterintuitive result, the diversity of the media is not only not conducive to opinion polarization, but also plays a role in discretizing public opinion. ) The distribution of the media's opinions plays an important role in influencing public opinion: with more than a single medium, opinions of diverse media with a gradual increase in density around a particular opinion can concentrate agents on the same opinion. The reason is that the media which hold opinions slightly away from the center act as agencies and bridges to guide greater numbers of agents to the center.

Moreover, as observed in
. Furthermore, it should be noted that the increased number of media is not always e ective, because similar media with high density attract agents with approximate opinions around them, isolating or repulsing other remaining agents. Similar phenomena have been observed in several works, such as Pineda & Buendía ( ), Pulick et al. ( ), and Gargiulo & Gandica ( ).In these papers, if media pressure is low, media are able to attract most of the agents but have a self-defeating e ect if the media pressure is too high. Meanwhile, compared with case , the e ectiveness of increasing media number is relatively poor in case , because the strong opinion and more media accentuate the extrusion e ect, which is in agreement with the result previously obtained in Gonzalez-Avella et al. ( ) and Gargiulo et al. ( ).

Scenario .
In real life, policy makers design their campaigns for consciousness or education to improve public, health, or environmental awareness in the general public. Take the public health crisis caused by the COVID-pandemic as an example. Governments worldwide try to propagandize the right behavior for protecting people from COVID-, from social distancing and hand washing to wearing medical masks. Combining the knowledge obtained from the above simulations and discussions, we present two methods to improve epidemic prevention knowledge in the general public. In this Scenario, opinion -stands for the agent distrusts the public health institution and panics about the epidemic, while opinion indicates that the agent trusts the public health institution and has obtained the correct knowledge of epidemic prevention from propaganda.Therefore, for the purpose of greater health awareness, the Euclidean norm is given by: The number of the media in this Scenario is set to .
. Method : The opinions of the media are generated by following arc-shaped function: where k = 1, 2, 3 · · · 217, p is the adjustable parameter which decides the average opinions of the media, and sat() is a saturation function as follows: Although the average opinion of the media decreases with larger p, average public opinion first increases then decreases as illustrated in Figure (a). It reflects the fact, as discussed in Scenario , that in guiding public opinion, the distribution pattern is more significant than the average opinion of the media. In fact, by constructing the space of media's opinions, it is possible to achieve much better results than through one authoritarian medium. This conclusion can be obtained by comparing Figure (a) and Figure (a): the average for public opinion can be as high as 0.7053 with Method , but it only reaches 0.5791 under one authoritarian medium in Scenario . However, Figure (b) shows that the Euclidean norm is smallest when p = 0.6, which means that public opinion is closest to the target. Comparing the opinion distribution of two situations (p = 0.3 and 0.6), it can be said that most agents have opinions equal to 1 when p = 0.3, but there are more agents close to the right side of the coordinate with fewer opponents when p = 0.6. In practice, policy makers can choose between those two situations or compromise according to their purpose.
. Method : A section of the media is utilized as a set of guide media to improve agents' health awareness, and those media change their opinions over time according to: As shown in Equation , the opinions of guide media, in timesteps, increase from −0.6 to 1. In this Method, the guide media are selected randomly while the opinions of other media are uniformly distributed in [−1, +1], and the percentage of guide media will increase from % up to %. The simulation results are given in Figure  . From Figure (a) and (b), it is evident that the average opinion increases and the Euclidean norm decrease with higher percentage of guide media, until it is over %. However, when all the media become guide media, the public opinion declines instead. The reason is similar to what we discussed in Scenario , because overpowering propaganda will increase antipathy and generate more opponents. In fact, compared with Method , this method not only has higher values (up to 0.8059) for average opinion, but most of the agents reach maximum when the time is over, with fewer opponents, as shown in Figure (  . Finally, by constructing opinions of media in space (Method ) or time (Method ), the e ect of campaigns for consciousness or education can be e icaciously improved. Nonetheless, too many media speaking with one voice may be counterproductive, thus generating more opponents. Of course, to achieve the best results (i.e., highest average public opinion), at least % of media in Method and % of media in Method need to broadcast the information with the same opinion, and this is unattainable in real life. Moreover, the research (Castaldo et al. ) shows that even if people's circadian rhythms change and they spend more time on social media, there is no evidence of information mass increasing during the lockdown. Therefore, too many media broadcasting similar information during lockdown time may lead to a backfired situation. Even so, a small proportion of media choosing opinions following our approach can still obtain certain achievements. Hence, the methods we proposed can o er insights that help policy makers to improve public, health, or environmental awareness in the general public.

Conclusions
. In this paper, we proposed a novel model of opinion dynamics, considering degrees of influence and free choice among social agents. Individuals in a social network interacted with their acquaintances and media they chose, and they attracted to each other's opinions when those opinions were similar enough, but pushing each other away when opinions were too di erent. Furthermore, we introduced diverse media and their representations of opinion, in di erent quantities and distributions, into the social network as external information sources. Finally, we proposed strategies to help policy makers improve the e ects of campaigns for consciousness or education.
. Our numerical simulations based on a real-world social network revealed the following conclusions. First, when given choices to make individually, people from open-minded societies are more susceptible to the media, while those from close-minded societies ignore the media and brew more extreme opinions. Furthermore, if every individual is open-minded, an overwhelming propaganda program could have the opposite of its intended e ect. Secondly, compared to the situation under one authoritarian medium, diverse media's presentation of more choices can result in scattered opinions and reduce the harm that extreme opinions could cause. Although the decentralized media will not lead to opinion polarization, we cannot a ord to ignore the threat the polarized media may bring. Finally, taking into account that the distribution of opinions can be a ected by the distribution of the media, we designed two methods of media programs with specific distribution and changing opinions. Simulations verified that the methods could e ectively enhance awareness among the public. .
In this work, we applied the model in a real-life social network to study the influence on public opinion brought by diverse media and audience choices. Nevertheless, our study has several limitations. In the future, we will introduce more realism to each aspect of the model and simulation, for instance, agents have prejudices on specific topics or feel loyal to certain media, and there are special individuals such as informed agents, inflexible agents, contrarians, and zealot agents. The simulations will not be limited to one network, while di erent networks from reality can be used, and it will be possible to study the e ect of the social networks' structure on public opinion. Moreover, contemporary media compete with each other and adjust their opinions to attract audiences, and this will also bring significant changes to the model. Nevertheless, the reality is full of complexity and randomness, and we cannot simulate all real-world elements in one enormous simulation; doing so would undoubtedly lead to undecipherable results. .
Considering that the study of public opinion is a combination of social science, computer science, and math, its future success will depend not only on modeling and simulation but also on creative and multidisciplinary methods. Therefore, new knowledge and technology, such as deep learning, can help us understand the real thoughts of the public, providing new perspectives in this field. Finally, such approaches will help us obtain more novel and practical explanations and conclusions.

Model Documentation
The agent-based model is written and simulated in Matlab. All related files containing the simulation model's code are accessible at https://www.comses.net/codebase-release/19cbafdc-15ca-4fe1-a818-05f3d422bec7https://w release/ cbafdc-ca-fe -a -f d bec . Please read the "readme.md" file for more details.

Appendix A: How the Opinion Leaders A ect Opinion Dynamic?
In this section, the role of individual influence in the formation of public opinion is discussed. Rewrite Equation as: the additional parameter s determines the distribution characteristics of individuals' influence. The increase of the parameter s represents the increase of the gap between the influence levels of di erent individuals. It also means that the influence of opinion leaders on public opinion will be more obvious. The value of s is set from increases to , and the agents' horizon values are between each one's convergence and repulsion thresholds. The other parameters' other settings are the same as in Scenario , and the results are shown in Figure . It can be seen from the Figure (a) that the average opinion of the public begins to decrease slightly with the increase of s, but the randomness will greatly increase a er s exceeds . Meanwhile, the trend of participation ratio is consistent with the increase of s, especially when the value of s exceeds , the participation ratio can be near . in some cases. It means that most individuals of the public have reached a certain degree of consensus under the leadership of opinion leaders. Figure (d) depicts the possible distribution of public opinions when s = 10000 (one simulation). Unlike the more scattered distribution in Figure (d), public opinions are clearly divided into two clusters, except for the % of the public who converge to the media opinion . , most of the rest are opponents with extreme opinions. Through observation Figure  (c), we can tell some key agents play vital roles in forming public opinion. For instance, the agent numbered has always insisted on its own opinion of . , and eventually attracted some followers.
In summary, as the value of s increases, the di erences in influence levels between individuals become larger. On the one hand, these opinion leaders with greater influence can help the media spread their opinions, but on the other hand, these leaders who hold obstinate opinions may also completely change the trend of public opinion and bring great randomness to the final public opinion. The main object of this paper is to explore how diverse media might influence public opinion. Therefore, the parameter s is set relatively small to reduce the randomness the opinion leader may bring.

Appendix B: Robustness Check with Di erent Repulsion Threshold Setting
In this Appendix section, we explore whether the distribution of τ i will a ect the general trends of simulation results. Di erent from the previous simulations, in this section, the repulsion thresholds τ i follow normal distribution N 1, 0.3 2 while satisfying τ i ≥ ε i . This setting limits the extreme value of τ i and makes the average value of τ i smaller. The simulations are taken based on similar settings of Scenario to except for the value of τ i , and the results are presented in Figures to .      Figure (c), the general trends are identical. However, through careful inspection, we can explore the following di erences compared with the previous simulation results. First, the average opinion values are smaller, especially when G = 1, only slightly bigger than G = 0.6 when η i ≥ 1.4. Secondly, the non-monotonic is still existing but not that obvious until the opinion of the medium is extreme enough. Both changes are not di icult to explain, the normally distributed τ i have limited extreme values and small mean, therefore, agents are more likely to repulse the opinion of the authoritarian medium. Meanwhile, the model is insensitive to the parameter η i due to more repulsive interactions, which leads to the inconspicuous of the non-monotonic.
In Figure , the simulation settings are the same as case and of Scenario except for the setting of τ i . The change patterns of the Euclidean norm are similar to the Figure (b) and (b). Nonetheless, the small τ i brings more randomness, the Euclidean norm even bigger with media than medium, as shown in Figure (a). Even then, the norm is continually decreasing as the number of media increases. The impact brought by the multiple media is still diminishing when the number increases, just as in Scenario .
The public opinions under two guide methods are depicted in Figures and . The curves' shapes are basically identical to Figures and , respectively; the di erences are that the average opinion is smaller and Euclidean norm is larger due to the intensive repulsive interaction. Moreover, in this situation, method has more possibilities of creating opposed public opinion, hence its e ect is worse than that of method . Because the smaller repulsive threshold setting makes it easier for the agents to revolt against the media broadcasting the same information. This result illuminates that the policy makers should consider both the appropriate method and the social environment simultaneously when designing the propaganda.
In this Appendix section, we have tested a di erent setting of the repulsive threshold, and obtained similar results as the previous simulations, proving the robustness of the model.