Invasion and interaction determine population composition in an open evolving ecological system.

It is well-known that interactions between species determine the population composition in an ecosystem. Conventional studies have focused on fixed population structures to reveal how interactions shape population compositions. However, interaction structures are not fixed but change over time due to invasions. Thus, invasion and interaction play an important role in shaping communities. Despite its importance, however, the interplay between invasion and interaction has not been well explored. Here, we investigate how invasion affects the population composition with interactions in open evolving ecological systems considering generalized Lotka-Volterra-type dynamics. Our results show that the system has two distinct regimes. One is characterized by low diversity with abrupt changes of dominant species in time, appearing when the interaction between species is strong and invasion slowly occurs. On the other hand, frequent invasions can induce higher diversity with slow changes in abundances despite strong interactions. It is because invasion happens before the system reaches its equilibrium, which drags the system from its equilibrium all the time. All species have similar abundances in this regime, which implies that fast invasion induces regime shift. Therefore, whether invasion or interaction dominates determines the population composition.

Micropyramidal Flexible Ion Gel Sensor for Multianalyte Discrimination and Strain Compensation.

A highly sensitive and flexible gas sensor that can detect a wide range of chemicals is crucial for wearable applications. However, conventional single resistance-based flexible sensors face challenges in maintaining chemical sensitivity under mechanical stress and can be affected by interfering gases. This study presents a versatile approach for fabricating a micropyramidal flexible ion gel sensor, which accomplishes sub-ppm sensitivity (<80 ppb) at room temperature and discrimination capability between various analytes, including toluene, isobutylene, ammonia, ethanol, and humidity. The discrimination accuracy of our flexible sensor is as high as 95.86%, enhanced by using machine learning-based algorithms. Moreover, its sensing capability remains stable with only a 2.09% change from the flat state to a 6.5 mm bending radius, further amplifying its universal usage for wearable chemical sensing. Therefore, we envision that a micropyramidal flexible ion gel sensor platform assisted by machine learning-based algorithms will provide a new strategy toward next-generation wearable sensing technology.

Degree distributions under general node removal: Power-law or Poisson?

Perturbations made to networked systems may result in partial structural loss, such as a blackout in a power-grid system. Investigating the resulting disturbance in network properties is quintessential to understand real networks in action. The removal of nodes is a representative disturbance, but previous studies are seemingly contrasting about its effect on arguably the most fundamental network statistic, the degree distribution. The key question is about the functional form of the degree distributions that can be altered during node removal or sampling. The functional form is decisive in the remaining subnetwork's static and dynamical properties. In this work, we clarify the situation by utilizing the relative entropies with respect to the reference distributions in the Poisson and power-law form, to quantify the distance between the subnetwork's degree distribution and either of the reference distributions. Introducing general sequential node removal processes with continuously different levels of hub protection to encompass a series of scenarios including uniform random removal and preferred or protective (i.e., biased random) removal of the hub, we classify the altered degree distributions starting from various power-law forms by comparing two relative entropy values. From the extensive investigation in various scenarios based on direct node-removal simulations and by solving the rate equation of degree distributions, we discover in the parameter space two distinct regimes, one where the degree distribution is closer to the power-law reference distribution and the other closer to the Poisson distribution.

Power-grid stability predictions using transferable machine learning.

Complex network analyses have provided clues to improve power-grid stability with the help of numerical models. The high computational cost of numerical simulations, however, has inhibited the approach, especially when it deals with the dynamic properties of power grids such as frequency synchronization. In this study, we investigate machine learning techniques to estimate the stability of power-grid synchronization. We test three different machine learning algorithms-random forest, support vector machine, and artificial neural network-training them with two different types of synthetic power grids consisting of homogeneous and heterogeneous input-power distribution, respectively. We find that the three machine learning models better predict the synchronization stability of power-grid nodes when they are trained with the heterogeneous input-power distribution rather than the homogeneous one. With the real-world power grids of Great Britain, Spain, France, and Germany, we also demonstrate that the machine learning algorithms trained on synthetic power grids are transferable to the stability prediction of the real-world power grids, which implies the prospective applicability of machine learning techniques on power-grid studies.

Synchronization of active rotators interacting with environment.

Multiple organs in a living system respond to environmental changes, and the signals from the organs regulate the physiological environment. Inspired by this biological feedback, we propose a simple autonomous system of active rotators to explain how multiple units are synchronized under a fluctuating environment. We find that the feedback via an environment can entrain rotators to have synchronous phases for specific conditions. This mechanism is markedly different from the simple entrainment by a common oscillatory external stimulus that is not interacting with systems. We theoretically examine how the phase synchronization depends on the interaction strength between rotators and environment. Furthermore, we successfully demonstrate the proposed model by realizing an analog electric circuit with microelectronic devices. This bioinspired platform can be used as a sensor for monitoring varying environments and as a controller for amplifying signals by their feedback-induced synchronization.

On structural and dynamical factors determining the integrated basin instability of power-grid nodes.

In electric power systems delivering alternating current, it is essential to maintain its synchrony of the phase with the rated frequency. The synchronization stability that quantifies how well the power-grid system recovers its synchrony against perturbation depends on various factors. As an intrinsic factor that we can design and control, the transmission capacity of the power grid affects the synchronization stability. Therefore, the transition pattern of the synchronization stability with the different levels of transmission capacity against external perturbation provides the stereoscopic perspective to understand the synchronization behavior of power grids. In this study, we extensively investigate the factors affecting the synchronization stability transition by using the concept of basin stability as a function of the transmission capacity. For a systematic approach, we introduce the integrated basin instability, which literally adds up the instability values as the transmission capacity increases. We first take simple 5-node motifs as a case study of building blocks of power grids, and a more realistic IEEE 24-bus model to highlight the complexity of decisive factors. We find that both structural properties such as gate keepers in network topology and dynamical properties such as large power input/output at nodes cause synchronization instability. The results suggest that evenly distributed power generation and avoidance of bottlenecks can improve the overall synchronization stability of power-grid systems.

In-depth data on the network structure and hourly activity of the Central Chilean power grid.

Network science enables us to improve the performance of complex systems such as traffic, communication, and power grids. To do so, it is necessary to use a well-constructed flawless network dataset associated with the system of interest. In this study, we present the dataset of the Chilean power grid. We harmonized data from three diverse sources to generate a unified dataset. Through an intensive review on the raw data, we filter out inconsistent errors and unrealistic faults, making the data more trustworthy. In contrast to other network dataset for power grids, we especially focus on preserving the physical structure of nodes' connection incorporating the 'tap' structure. As a result, we provide three different versions of the dataset: 'with-tap', 'without-tap', and 'reduced versions'. Along with structure, we incorporate various attributes of the nodes and edges such as the geo-coordinates, voltage of transmission lines, and the time series data of generation or consumption. These data are useful for network scientists to analyze the performance and dynamic stability of power grids.

Laser Writing Block Copolymer Self-Assembly on Graphene Light-Absorbing Layer.

Recent advance of high-power laser processing allows for rapid, continuous, area-selective material fabrication, typically represented by laser crystallization of silicon or oxides for display applications. Two-dimensional materials such as graphene exhibit remarkable physical properties and are under intensive development for the manufacture of flexible devices. Here we demonstrate an area-selective ultrafast nanofabrication method using low intensity infrared or visible laser irradiation to direct the self-assembly of block copolymer films into highly ordered manufacturing-relevant architectures at the scale below 12 nm. The fundamental principles underlying this light-induced nanofabrication mechanism include the self-assembly of block copolymers to proceed across the disorder-order transition under large thermal gradients, and the use of chemically modified graphene films as a flexible and conformal light-absorbing layers for transparent, nonplanar, and mechanically flexible surfaces.

Nash equilibrium and evolutionary dynamics in semifinalists' dilemma.

We consider a tournament among four equally strong semifinalists. The players have to decide how much stamina to use in the semifinals, provided that the rest is available in the final and the third-place playoff. We investigate optimal strategies for allocating stamina to the successive matches when players' prizes (payoffs) are given according to the tournament results. From the basic assumption that the probability to win a match follows a nondecreasing function of stamina difference, we present symmetric Nash equilibria for general payoff structures. We find three different phases of the Nash equilibria in the payoff space. First, when the champion wins a much bigger payoff than the others, any pure strategy can constitute a Nash equilibrium as long as all four players adopt it in common. Second, when the first two places are much more valuable than the other two, the only Nash equilibrium is such that everyone uses a pure strategy investing all stamina in the semifinal. Third, when the payoff for last place is much smaller than the others, a Nash equilibrium is formed when every player adopts a mixed strategy of using all or none of its stamina in the semifinals. In a limiting case that only last place pays the penalty, this mixed-strategy profile can be proved to be a unique symmetric Nash equilibrium, at least when the winning probability follows a Heaviside step function. Moreover, by using this Heaviside step function, we study the tournament by using evolutionary replicator dynamics to obtain analytic solutions, which reproduces the corresponding Nash equilibria on the population level and gives information on dynamic aspects.

Anomalous rapid defect annihilation in self-assembled nanopatterns by defect melting.

Molecular self-assembly commonly suffers from dense structural defect formation. Spontaneous defect annihilation in block copolymer (BCP) self-assembly is particularly retarded due to significant energy barrier for polymer chain diffusion and structural reorganization. Here we present localized defect melting induced by blending short neutral random copolymer chain as an unusual method to promote the defect annihilation in BCP self-assembled nanopatterns. Chemically neutral short random copolymer chains blended with BCPs are specifically localized and induce local disordered states at structural defect sites in the self-assembled nanopatterns. Such localized "defect melting" relieves the energy penalty for polymer diffusion and morphology reorganization such that spontaneous defect annihilation by mutual coupling is anomalously accelerated upon thermal annealing. Interestingly, neutral random copolymer chain blending also causes morphology-healing self-assembly behavior that can generate large-area highly ordered 10 nm scale nanopattern even upon poorly defined defective prepatterns. Underlying mechanisms of the unusual experimental findings are thoroughly investigated by three-dimensional self-consistent field theory calculation.

Large-scale quantitative analysis of painting arts.

Scientists have made efforts to understand the beauty of painting art in their own languages. As digital image acquisition of painting arts has made rapid progress, researchers have come to a point where it is possible to perform statistical analysis of a large-scale database of artistic paints to make a bridge between art and science. Using digital image processing techniques, we investigate three quantitative measures of images - the usage of individual colors, the variety of colors, and the roughness of the brightness. We found a difference in color usage between classical paintings and photographs, and a significantly low color variety of the medieval period. Interestingly, moreover, the increment of roughness exponent as painting techniques such as chiaroscuro and sfumato have advanced is consistent with historical circumstances.

Wrinkle-directed self-assembly of block copolymers for aligning of nanowire arrays.

Highly aligned metal nanowire arrays with feature sizes approaching 10 nm are fabricated. This is made possible by the self-assembly of block copolymers (BCPs) on graphene-wrinkle arrays. Thickness-modulated BCP films confined on the wrinkled reduced graphene oxide (rGO) surface promote the strict alignment of the self-assembled BCP lamellae in the direction of the film thickness gradient.

Costly bilingualism model in a population with one zealot.

We consider a costly bilingualism model in which one can take two strategies in parallel. We investigate how a single zealot triggers the cascading behavior and how the compatibility of the two strategies affects when interacting patterns change. First, the role of the interaction range in the cascading is studied by increasing the range from local to global. We find that people sometimes do not favor taking the superior strategy even though its payoff is higher than that of the inferior one. This is found to be caused by the local interactions rather than the global ones. Applying this model to social networks, we find that the location of the zealot is also important for larger cascading in heterogeneous networks.

Phase-shift inversion in oscillator systems with periodically switching couplings.

A system's response to external periodic changes can provide crucial information about its dynamical properties. We investigate the synchronization transition, an archetypical example of a dynamic phase transition, in the framework of such a temporal response. The Kuramoto model under periodically switching interactions has the same type of phase transition as the original mean-field model. Furthermore, we see that the signature of the synchronization transition appears in the relative delay of the order parameter with respect to the phase of oscillating interactions as well. Specifically, the phase shift becomes significantly larger as the system gets closer to the phase transition, so that the order parameter at the minimum interaction density can even be larger than that at the maximum interaction density, counterintuitively. We argue that this phase-shift inversion is caused by the diverging relaxation time, in a similar way to the resonance near the critical point in the kinetic Ising model. Our result, based on exhaustive simulations on globally coupled systems as well as scale-free networks, shows that an oscillator system's phase transition can be manifested in the temporal response to the topological dynamics of the underlying connection structure.

Cells transformed by PLC-gamma 1 overexpression are highly sensitive to clostridium difficile toxin A-induced apoptosis and mitotic inhibition.

Phospholipase C-γl (PLC-γl) expression is associated with cellular transformation. Notably, PLC-gamma is up-regulated in colorectal cancer tissue and breast carcinoma. Because exotoxins released by Clostridium botulinum have been shown to induce apoptosis and promote growth arrest in various cancer cell lines, we examined here the potential of Clostridium difficile toxin A to selectively induce apoptosis in cells transformed by PLC-γl overexpression. We found that PLC-γl-transformed cells, but not vectortransformed (control) cells, were highly sensitive to C. difficile toxin A-induced apoptosis and mitotic inhibition. Moreover, expression of the proapoptotic Bcl2 family member, Bim, and activation of caspase-3 were significantly up-regulated by toxin A in PLC-γl-transformed cells. Toxin A-induced cell rounding and paxillin dephosphorylation were also significantly higher in PLC-γl-transformed cells than in control cells. These findings suggest that C. difficile toxin A may have potential as an anticancer agent against colorectal cancers and breast carcinomas in which PLC-γl is highly up-regulated.

Exact solutions for mass-dependent irreversible aggregations.

We consider the mass-dependent aggregation process (k+1)X→X, given a fixed number of unit mass particles in the initial state. One cluster is chosen proportional to its mass and is merged into one, either with k neighbors in one dimension, or--in the well-mixed case--with k other clusters picked randomly. We find the same combinatorial exact solutions for the probability to find any given configuration of particles on a ring or line, and in the well-mixed case. The mass distribution of a single cluster exhibits scaling laws and the finite-size scaling form is given. The relation to the classical sum kernel of irreversible aggregation is discussed.

Percolation transitions are not always sharpened by making networks interdependent.

We study a model for coupled networks introduced recently by Buldyrev et al., [Nature (London) 464, 1025 (2010)], where each node has to be connected to others via two types of links to be viable. Removing a critical fraction of nodes leads to a percolation transition that has been claimed to be more abrupt than that for uncoupled networks. Indeed, it was found to be discontinuous in all cases studied. Using an efficient new algorithm we verify that the transition is discontinuous for coupled Erdös-Rényi networks, but find it to be continuous for fully interdependent diluted lattices. In 2 and 3 dimensions, the order parameter exponent ß is larger than in ordinary percolation, showing that the transition is less sharp, i.e., further from discontinuity, than for isolated networks. Possible consequences for spatially embedded networks are discussed.

Explosive percolation is continuous, but with unusual finite size behavior.

We study four Achlioptas-type processes with "explosive" percolation transitions. All transitions are clearly continuous, but their finite size scaling functions are not entirely holomorphic. The distributions of the order parameter, i.e., the relative size s(max)/N of the largest cluster, are double humped. But-in contrast to first-order phase transitions-the distance between the two peaks decreases with system size N as N(-η) with η>0. We find different positive values of ß (defined via (s(max)/N)∼(p-p(c))ß for infinite systems) for each model, showing that they are all in different universality classes. In contrast, the exponent Θ (defined such that observables are homogeneous functions of (p-p(c))N(Θ)) is close to-or even equal to-1/2 for all models.

Enhancing synchronization by directionality in complex networks.

We propose a method called the residual edge-betweenness gradient (REBG) to enhance the synchronizability of networks by assigning the link direction while keeping the topology and link weights unchanged. Direction assignment has been shown to improve the synchronizability of undirected networks in general, but we find that in some cases incommunicable components emerge and networks fail to synchronize. We show that the REBG method improves the residual degree gradient (RDG) method by effectively avoiding the synchronization failure. Further experiments show that the REBG method enhances the synchronizability in networks with a community structure compared with the RDG method.