
Collective, especially groupbased, managerial decision making is crucial in
organizations. Using an evolutionary theoretic approach to collective decision
making, agentbased simulations were conducted to investigate how human
collective decision making would be affected by the agents' diversity in
problem understanding and/or behavior in discussion, as well as by their social
network structure. Simulation results indicated that groups with consistent
problem understanding tended to produce higher utility values of ideas and
displayed better decision convergence, but only if there was no grouplevel
bias in collective problem understanding. Simulation results also indicated the
importance of balance between selectionoriented (i.e., exploitative) and
variationoriented (i.e., explorative) behaviors in discussion to achieve
quality final decisions. Expanding the group size and introducing nontrivial
social network structure generally improved the quality of ideas at the cost of
decision convergence. Simulations with different social network topologies
revealed collective decision making on smallworld networks with high local
clustering tended to achieve highest decision quality more often than on random
or scalefree networks. Implications of this evolutionary theory and simulation
approach for future managerial research on collective, group, and multilevel
decision making are discussed.

We propose a new polynomialtime deterministic algorithm that produces an
approximated solution for the traveling salesperson problem. The proposed
algorithm ranks cities based on their priorities calculated using a power
function of means and standard deviations of their distances from other cities
and then connects the cities to their neighbors in the order of their
priorities. When connecting a city, a neighbor is selected based on their
neighbors' priorities calculated as another power function that additionally
includes their distance from the focal city to be connected. This repeats until
all the cities are connected into a single loop. The time complexity of the
proposed algorithm is $O(n^2)$, where $n$ is the number of cities. Numerical
evaluation shows that, despite its simplicity, the proposed algorithm produces
shorter tours with less time complexity than other conventional tour
construction heuristics. The proposed algorithm can be used by itself or as an
initial tour generator for other more complex heuristic optimization
algorithms.

Gene regulatory network (GRN)based morphogenetic models have recently gained
an increasing attention. However, the relationship between microscopic
properties of intracellular GRNs and macroscopic properties of morphogenetic
systems has not been fully understood yet. Here we propose a theoretical
morphogenetic model representing an aggregation of cells, and reveal the
relationship between criticality of GRNs and morphogenetic pattern formation.
In our model, the positions of the cells are determined by springmassdamper
kinetics. Each cell has an identical Kauffman's $NK$ random Boolean network
(RBN) as its GRN. We varied the properties of GRNs from ordered, through
critical, to chaotic by adjusting node indegree $K$. We randomly assigned four
cell fates to the attractors of RBNs for cellular behaviors. By comparing
diverse morphologies generated in our morphogenetic systems, we investigated
what the role of the criticality of GRNs is in forming morphologies. We found
that nontrivial spatial patterns were generated most frequently when GRNs were
at criticality. Our finding indicates that the criticality of GRNs facilitates
the formation of nontrivial morphologies in GRNbased morphogenetic systems.

We studied the longterm dynamics of evolutionary Swarm Chemistry by
extending the simulation length tenfold compared to earlier work and by
developing and using a new automated object harvesting method. Both macroscopic
dynamics and microscopic object features were characterized and tracked using
several measures. Results showed that the evolutionary dynamics tended to
settle down into a stable state after the initial transient period, and that
the extent of environmental perturbations also affected the evolutionary trends
substantially. In the meantime, the automated harvesting method successfully
produced a huge collection of spontaneously evolved objects, revealing the
system's autonomous creativity at an unprecedented scale.

Selforganization has been an important concept within a number of
disciplines, which Artificial Life (ALife) also has heavily utilized since its
inception. The term and its implications, however, are often confusing or
misinterpreted. In this work, we provide a minireview of selforganization and
its relationship with ALife, aiming at initiating discussions on this important
topic with the interested audience. We first articulate some fundamental
aspects of selforganization, outline its usage, and review its applications to
ALife within its soft, hard, and wet domains. We also provide perspectives for
further research.

Cooperation is ubiquitous at every level of living organisms. It is known
that spatial (network) structure is a viable mechanism for cooperation to
evolve. A recently proposed numerical metric, Average Gradient of Selection
(AGoS), a useful tool for interpreting and visualizing evolutionary dynamics on
networks, allows simulation results to be visualized on a onedimensional phase
space. However, stochastic mutation of strategies was not considered in the
analysis of AGoS. Here we extend AGoS so that it can analyze the evolution of
cooperation where mutation may alter strategies of individuals on networks. We
show that our extended AGoS correctly visualizes the final states of
cooperation with mutation in the individualbased simulations. Our analyses
revealed that mutation always has a negative effect on the evolution of
cooperation regardless of the payoff functions, fraction of cooperators, and
network structures. Moreover, we found that scalefree networks are the most
vulnerable to mutation and thus the dynamics of cooperation are altered from
bistability to coexistence on those networks, undergoing an imperfect pitchfork
bifurcation.

Due to the interdisciplinary nature of complex systems as a field, students
studying complex systems at University level have diverse disciplinary
backgrounds. This brings challenges (e.g. wide range of computer programming
skills) but also opportunities (e.g. facilitating interdisciplinary
interactions and projects) for the classroom. However, there is little
published regarding how these challenges and opportunities are handled in
teaching and learning Complex Systems as an explicit subject in higher
education, and how this differs in comparison to other subject areas. We seek
to explore these particular challenges and opportunities via an interviewbased
study of pioneering teachers and learners (conducted amongst the authors)
regarding their experiences. We compare and contrast those experiences, and
analyse them with respect to the educational literature. Our discussions
explored: approaches to curriculum design, how theories/models/frameworks of
teaching and learning informed decisions and experience, how diversity in
student backgrounds was addressed, and assessment task design. We found a
striking level of commonality in the issues expressed as well as the strategies
to handle them, for example a significant focus on problembased learning, and
the use of major studentled creative projects for both achieving and assessing
learning outcomes.

Whereas the relationship between criticality of gene regulatory networks
(GRNs) and dynamics of GRNs at a single cell level has been vigorously studied,
the relationship between the criticality of GRNs and system properties at a
higher level has remained unexplored. Here we aim at revealing a potential role
of criticality of GRNs at a multicellular level which are hard to uncover
through the singlecelllevel studies, especially from an evolutionary
viewpoint. Our model simulated the growth of a cell population from a single
seed cell. All the cells were assumed to have identical GRNs. We induced
genetic perturbations to the GRN of the seed cell by adding, deleting, or
switching a regulatory link between a pair of genes. From numerical
simulations, we found that the criticality of GRNs facilitated the formation of
nontrivial morphologies when the GRNs were critical in the presence of the
evolutionary perturbations. Moreover, the criticality of GRNs produced
topologically homogenous cell clusters by adjusting the spatial arrangements of
cells, which led to the formation of nontrivial morphogenetic patterns. Our
findings corresponded to an epigenetic viewpoint that heterogeneous and complex
features emerge from homogeneous and less complex components through the
interactions among them. Thus, our results imply that highly structured tissues
or organs in morphogenesis of multicellular organisms might stem from the
criticality of GRNs.

Many living and nonliving complex systems can be modeled and understood as
collective systems made of heterogeneous components that selforganize and
generate nontrivial morphological structures and behaviors. This chapter
presents a brief overview of our recent effort that investigated various
aspects of such morphogenetic collective systems. We first propose a
theoretical classification scheme that distinguishes four complexity levels of
morphogenetic collective systems based on the nature of their components and
interactions. We conducted a series of computational experiments using a
selfpropelled particle swarm model to investigate the effects of (1)
heterogeneity of components, (2) differentiation/redifferentiation of
components, and (3) local information sharing among components, on the
selforganization of a collective system. Results showed that (a) heterogeneity
of components had a strong impact on the system's structure and behavior, (b)
dynamic differentiation/redifferentiation of components and local information
sharing helped the system maintain spatially adjacent, coherent organization,
(c) dynamic differentiation/redifferentiation contributed to the development
of more diverse structures and behaviors, and (d) stochastic redifferentiation
of components naturally realized a selfrepair capability of selforganizing
morphologies. We also explored evolutionary methods to design novel
selforganizing patterns, using interactive evolutionary computation and
spontaneous evolution within an artificial ecosystem. These selforganizing
patterns were found to be remarkably robust against dimensional changes from 2D
to 3D, although evolution worked efficiently only in 2D settings.

Despite the recent advances in developing more effective thresholding methods
to convert weighted networks to unweighted counterparts, there are still
several limitations that need to be addressed. One such limitation is the
inability of the most existing thresholding methods to take into account the
topological properties of the original weighted networks during the
binarization process, which could ultimately result in unweighted networks that
have drastically different topological properties than the original weighted
networks. In this study, we propose a new thresholding method based on the
percolation theory to address this limitation. The performance of the proposed
method was validated and compared to the existing thresholding methods using
simulated and realworld functional connectivity networks in the brain.
Comparison of macroscopic and microscopic properties of the resulted unweighted
networks to the original weighted networks suggest that the proposed
thresholding method can successfully maintain the topological properties of the
original weighted networks.

As network science has matured as an established field of research, there are
already a number of courses on this topic developed and offered at various
higher education institutions, often at postgraduate levels. In those courses,
instructors adopted different approaches with different focus areas and
curricular designs. We collected information about 30 existing network science
courses from various online sources, and analyzed the contents of their syllabi
or course schedules. The topics and their curricular sequences were extracted
from the course syllabi/schedules and represented as a directed weighted graph,
which we call the topic network. Community detection in the topic network
revealed seven topic clusters, which matched reasonably with the concept list
previously generated by students and educators through the Network Literacy
initiative. The minimum spanning tree of the topic network revealed typical
flows of curricular contents, starting with examples of networks, moving onto
random networks and smallworld networks, then branching off to various
subtopics from there. These results illustrate the current state of consensus
formation (including variations and disagreements) among the network science
community on what should be taught about networks and how, which may also be
informative for K12 education and informal education.

We propose a slightly revised MillerHagberg (MH) algorithm that efficiently
generates a random network from a given expected degree sequence. The revision
was to replace the approximated edge probability between a pair of nodes with a
combinatorically calculated edge probability that better captures the
likelihood of edge presence especially where edges are dense. The computational
complexity of this combinatorial MH algorithm is still in the same order as the
original one. We evaluated the proposed algorithm through several numerical
experiments. The results demonstrated that the proposed algorithm was
particularly good at accurately representing highdegree nodes in dense,
heterogeneous networks. This algorithm may be a useful alternative of other
more established network randomization methods, given that the data are
increasingly becoming larger and denser in today's network science research.

We propose a novel computational method to extract information about
interactions among individuals with different behavioral states in a biological
collective from ordinary video recordings. Assuming that individuals are acting
as finite state machines, our method first detects discrete behavioral states
of those individuals and then constructs a model of their state transitions,
taking into account the positions and states of other individuals in the
vicinity. We have tested the proposed method through applications to two
realworld biological collectives: termites in an experimental setting and
human pedestrians in a university campus. For each application, a robust
tracking system was developed inhouse, utilizing interactive human
intervention (for termite tracking) or online agentbased simulation (for
pedestrian tracking). In both cases, significant interactions were detected
between nearby individuals with different states, demonstrating the
effectiveness of the proposed method.

This paper aims to establish theoretical foundations of graph product
multilayer networks (GPMNs), a family of multilayer networks that can be
obtained as a graph product of two or more factor networks. Cartesian, direct
(tensor), and strong product operators are considered, and then generalized. We
first describe mathematical relationships between GPMNs and their factor
networks regarding their degree/strength, adjacency, and Laplacian spectra, and
then show that those relationships can still hold for nonsimple and generalized
GPMNs. Applications of GPMNs are discussed in three areas: predicting epidemic
thresholds, modeling propagation in nontrivial space and time, and analyzing
higherorder properties of selfsimilar networks. Directions of future research
are also discussed.

This short article presents a summary of the NetSciEd (Network Science and
Education) initiative that aims to address the need for curricula, resources,
accessible materials, and tools for introducing K12 students and the general
public to the concept of networks, a crucial framework in understanding
complexity. NetSciEd activities include (1) the NetSci High educational
outreach program (since 2010), which connects high school students and their
teachers with regional university research labs and provides them with the
opportunity to work on network science research projects; (2) the NetSciEd
symposium series (since 2012), which brings network science researchers and
educators together to discuss how network science can help and be integrated
into formal and informal education; and (3) the Network Literacy: Essential
Concepts and Core Ideas booklet (since 2014), which was created collaboratively
and subsequently translated into 18 languages by an extensive group of network
science researchers and educators worldwide.

Cooperation is ubiquitous in every level of living organisms. It is known
that spatial (network) structure is a viable mechanism for cooperation to
evolve. Until recently, it has been difficult to predict whether cooperation
can evolve at a network (population) level. To address this problem, Pinheiro
et al. proposed a numerical metric, called Average Gradient of Selection (AGoS)
in 2012. AGoS can characterize and forecast the evolutionary fate of
cooperation at a population level. However, stochastic mutation of strategies
was not considered in the analysis of AGoS. Here we analyzed the evolution of
cooperation using AGoS where mutation may occur to strategies of individuals in
networks. Our analyses revealed that mutation always has a negative effect on
the evolution of cooperation regardless of the fraction of cooperators and
network structures. Moreover, we found that mutation affects the fitness of
cooperation differently on different social network structures.

A stock market is considered as one of the highly complex systems, which
consists of many components whose prices move up and down without having a
clear pattern. The complex nature of a stock market challenges us on making a
reliable prediction of its future movements. In this paper, we aim at building
a new method to forecast the future movements of Standard & Poor's 500 Index
(S&P 500) by constructing timeseries complex networks of S&P 500 underlying
companies by connecting them with links whose weights are given by the mutual
information of 60minute price movements of the pairs of the companies with the
consecutive 5,340 minutes price records. We showed that the changes in the
strength distributions of the networks provide an important information on the
network's future movements. We built several metrics using the strength
distributions and network measurements such as centrality, and we combined the
best two predictors by performing a linear combination. We found that the
combined predictor and the changes in S&P 500 show a quadratic relationship,
and it allows us to predict the amplitude of the one step future change in S&P
500. The result showed significant fluctuations in S&P 500 Index when the
combined predictor was high. In terms of making the actual index predictions,
we built ARIMA models. We found that adding the network measurements into the
ARIMA models improves the model accuracy. These findings are useful for
financial market policy makers as an indicator based on which they can
interfere with the markets before the markets make a drastic change, and for
quantitative investors to improve their forecasting models.

It is well known that cooperation cannot be an evolutionary stable strategy
for a noniterative game in a wellmixed population. In contrast, structured
populations favor cooperation since cooperators can benefit each other by
forming local clusters. Previous studies have shown that scalefree networks
strongly promote cooperation. However, little is known about the invasion
mechanism of cooperation in scalefree networks. To study microscopic and
macroscopic behaviors of cooperators' invasion, we conducted computational
experiments of the evolution of cooperation in scalefree networks where,
starting from all defectors, cooperators can spontaneously emerge by mutation.
Since the evolutionary dynamics are influenced by the definition of fitness, we
tested two commonly adopted fitness functions: accumulated payoff and average
payoff. Simulation results show that cooperation is strongly enhanced with the
accumulated payoff fitness compared to the average payoff fitness. However, the
difference between the two functions decreases as the average degree increases.
As the average degree increases, cooperation decreases with the accumulated
payoff fitness, while it increases with the average payoff fitness. Moreover,
with the average payoff fitness, lowdegree nodes play a more important role in
spreading cooperative strategies compared to the case of the accumulated payoff
fitness.

Calculating a product of multiple graphs has been studied in mathematics,
engineering, computer science, and more recently in network science,
particularly in the context of multilayer networks. One of the important
questions to be addressed in this area is how to characterize spectral
properties of a product graph using those of its factor graphs. While several
such characterizations have already been obtained analytically (mostly for
adjacency spectra), characterization of Laplacian spectra of direct product and
strong product graphs has remained an open problem. Here we develop practical
methods to estimate Laplacian spectra of direct and strong product graphs from
spectral properties of their factor graphs using a few heuristic assumptions.
Numerical experiments showed that the proposed methods produced reasonable
estimation with percentage errors confined within a +/10% range for most
eigenvalues.

Cooperation is ubiquitous in biological and social systems. Previous studies
revealed that a preference toward similar appearance promotes cooperation, a
phenomenon called tagmediated cooperation or communitarian cooperation. This
effect is enhanced when a spatial structure is incorporated, because space
allows agents sharing an identical tag to regroup to form locally cooperative
clusters. In spatially distributed settings, one can also consider migration of
organisms, which has a potential to further promote evolution of cooperation by
facilitating spatial clustering. However, it has not yet been considered in
spatial tagmediated cooperation models. Here we show, using computer
simulations of a spatial model of evolutionary games with organismal migration,
that tagbased segregation and homophilic cooperation arise for a wide range of
parameters. In the meantime, our results also show another evolutionarily
stable outcome, where a high level of heterophilic cooperation is maintained in
spatially wellmixed patterns. We found that these two different forms of
tagmediated cooperation appear alternately as the parameter for temptation to
defect is increased.

Networks have become increasingly relevant to everyday life as human society
has become increasingly connected. Attaining a basic understanding of networks
has thus become a necessary form of literacy for people (and for youths in
particular). At the NetSci 2014 conference, we initiated a yearlong process to
develop an educational resource that concisely summarizes essential concepts
about networks that can be used by anyone of school age or older. The process
involved several brainstorming sessions on one key question: "What should every
person living in the 21st century know about networks by the time he/she
finishes secondary education?" Different sessions reached diverse participants,
which included professional researchers in network science, educators, and
highschool students. The generated ideas were connected by the students to
construct a concept network. We examined community structure in the concept
network to group ideas into a set of important themes, which we refined through
discussion into seven essential concepts. The students played a major role in
this development process by providing insights and perspectives that were often
unrecognized by researchers and educators. The final result, "Network Literacy:
Essential Concepts and Core Ideas", is now available as a booklet in several
different languages from http://tinyurl.com/networkliteracy .

Recent research has established both a theoretical basis and strong empirical
evidence that effective social behavior plays a beneficial role in the
maintenance of physical and psychological wellbeing of people. To test whether
social behavior and wellbeing are also associated in online communities, we
studied the correlations between the recovery of patients with mental disorders
and their behaviors in online social media. As the source of the data related
to the social behavior and progress of mental recovery, we used PatientsLikeMe
(PLM), the world's first openparticipation research platform for the
development of patientcentered health outcome measures. We first constructed
an online social network structure based on patienttopatient ties among 200
patients obtained from PLM. We then characterized patients' online social
activities by measuring the numbers of "posts and views" and "helpful marks"
each patient obtained. The patients' recovery data were obtained from their
selfreported status information that was also available on PLM. We found that
some node properties (indegree, eigenvector centrality and PageRank) and the
two online social activity measures were significantly correlated with
patients' recovery. Furthermore, we recollected the patients' recovery data
two months after the first data collection. We found significant correlations
between the patients' social behaviors and the second recovery data, which were
collected two months apart. Our results indicated that social interactions in
online communities such as PLM were significantly associated with the current
and future recoveries of patients with mental disorders.

We study a mathematical model of social diffusion on a symmetric weighted
network where individual nodes' states gradually assimilate to local social
norms made by their neighbors' average states. Unlike physical diffusion, this
process is not state conservational and thus the global state of the network
(i.e., sum of node states) will drift. The asymptotic average node state will
be the average of initial node states weighted by their strengths. Here we show
that, while the global state is not conserved in this process, the inner
product of strength and state vectors is conserved instead, and perfect
positive correlation between node states and local averages of their
self/neighbor strength ratios always results in upward (or at least neutral)
global drift. We also show that the strength assortativity negatively affects
the speed of homogenization. Based on these findings, we propose an adaptive
link weight adjustment method to achieve the highest upward global drift by
increasing the strengthstate correlation. The effectiveness of the method was
confirmed through numerical simulations and implications for realworld social
applications are discussed.

The nature of concept learning is a core question in cognitive science.
Theories must account for the relative difficulty of acquiring different
concepts by supervised learners. For a canonical set of six category types, two
distinct orderings of classification difficulty have been found. One ordering,
which we call paradigmspecific, occurs when adult human learners classify
objects with easily distinguishable characteristics such as size, shape, and
shading. The general order occurs in all other known cases: when adult humans
classify objects with characteristics that are not readily distinguished (e.g.,
brightness, saturation, hue); for children and monkeys; and when categorization
difficulty is extrapolated from errors in identification learning. The
paradigmspecific order was found to be predictable mathematically by measuring
the logical complexity of tasks, i.e., how concisely the solution can be
represented by logical rules.
However, logical complexity explains only the paradigmspecific order but not
the general order. Here we propose a new difficulty measurement, information
complexity, that calculates the amount of uncertainty remaining when a subset
of the dimensions are specified. This measurement is based on Shannon entropy.
We show that, when the metric extracts minimal uncertainties, this new
measurement predicts the paradigmspecific order for the canonical six category
types, and when the metric extracts average uncertainties, this new measurement
predicts the general order. Moreover, for learning category types beyond the
canonical six, we find that the minimaluncertainty formulation correctly
predicts the paradigmspecific order as well or better than existing metrics
(Boolean complexity and GIST) in most cases.

We present NetSci High, our NSFfunded educational outreach program that
connects high school students who are underrepresented in STEM (Science
Technology Engineering and Mathematics), and their teachers, with regional
university research labs and provides them with the opportunity to work with
researchers and graduate students on teambased, yearlong network science
research projects, culminating in a formal presentation at a network science
conference. This short paper reports the content and materials that we have
developed to date, including lesson plans and tools for introducing high school
students and teachers to network science; empirical evaluation data on the
effect of participation on students' motivation and interest in pursuing STEM
careers; the application of professional development materials for teachers
that are intended to encourage them to use network science concepts in their
lesson plans and curriculum; promoting districtlevel interest and engagement;
best practices gained from our experiences; and the future goals for this
project and its subsequent outgrowth.