By an appropriate choice for the phase of the complex nucleonic field and going over to the traveling coordinate, we reduce the coupled Higgs equations to the Hamiltonian form and treat the resulting equation using the dynamical system theory. We present a phase-space analysis of its stable points. The results of our study demonstrate that the equation can support both traveling- and standing-wave solutions. The traveling-wave solution appears in the form of a soliton and resides in the midst of doubly periodic standing-wave solutions.

How high should be the rate of immigration into a stochastic population in order to significantly reduce the probability of observing the population become extinct? Is there any relation between the population size distributions with and without immigration? Under what conditions can one justify the simple patch occupancy models, which ignore the population distribution and its dynamics in a patch, and treat a patch simply as either occupied or empty? We answer these questions by exactly solving a simple stochastic model obtained by adding a steady immigration to a variant of the Verhulst model: a prototypical model of an isolated stochastic population.

It is well known that the $q$-state clock model can exhibit a Kosterlitz-Thouless (KT) transition if $q$ is equal to or greater than a certain threshold, which has been believed to be five. However, recent numerical studies indicate that helicity modulus does not vanish in the high-temperature phase of the five-state clock model as predicted by the KT scenario. By performing Monte Carlo calculations under the fluctuating twist boundary condition, we show that it is because the five-state clock model does not have the fully continuous $U\left(1\right)$ symmetry even in the high-temperature phase while the six-state clock model does. We suggest that the upper transition of the five-state clock model is actually a weaker cousin of the KT transition so that it is $q\ge 6$ that exhibits the genuine KT behavior.

Recent numerical studies on an infinite number of identical spherical beads in Hertzian contact showed the presence of frequency bands [Jayaprakash, Starosvetsky, Vakakis, Peeters, and Kerschen, Nonlinear Dyn. 63, 359 (2011)]. These bands, denoted here as propagation and attenuation bands (PBs and ABs), are typically present in linear or weakly nonlinear periodic media; however, their counterparts are not intuitive in essentially nonlinear periodic media where there is a complete lack of classical linear acoustics, i.e., in “sonic vacua.” Here, we study the effects of PBs and ABs on the forced dynamics of ordered, uncompressed granular systems. Through numerical and experimental techniques, we find that the dynamics of these systems depends critically on the frequency and amplitude of the applied harmonic excitation. For fixed forcing amplitude, at lower frequencies, the oscillations are large in amplitude and governed by strongly nonlinear and nonsmooth dynamics, indicating PB behavior. At higher frequencies the dynamics is weakly nonlinear and smooth, in the form of compressed low-amplitude oscillations, indicating AB behavior. At the boundary between the PB and the AB large-amplitude oscillations due to resonance occur, giving rise to collisions between beads and chaotic dynamics; this renders the forced dynamics sensitive to initial and forcing conditions, and hence unpredictable. Finally, we study asymptotically the near field standing wave dynamics occurring for high frequencies, well inside the AB.

Rotating spiral waves organize spatial patterns in chemical, physical, and biological excitable systems. Factors affecting their dynamics, such as spatiotemporal drift, are of great interest for particular applications. Here, we propose a quantitative description for spiral wave dynamics on curved surfaces which shows that for a wide class of systems, including the Belousov-Zhabotinsky reaction and anisotropic cardiac tissue, the Ricci curvature scalar of the surface is the main determinant of spiral wave drift. The theory provides explicit equations for spiral wave drift direction, drift velocity, and the period of rotation. Depending on the parameters, the drift can be directed to the regions of either maximal or minimal Ricci scalar curvature, which was verified by direct numerical simulations.

The rogue wave solutions (rational multibreathers) of the nonlinear Schrödinger equation (NLS) are tested in numerical simulations of weakly nonlinear and fully nonlinear hydrodynamic equations. Only the lowest order solutions from 1 to 5 are considered. A higher accuracy of wave propagation in space is reached using the modified NLS equation, also known as the Dysthe equation. This numerical modeling allowed us to directly compare simulations with recent results of laboratory measurements in Chabchoub et al. [Phys. Rev. E 86, 056601 (2012)]. In order to achieve even higher physical accuracy, we employed fully nonlinear simulations of potential Euler equations. These simulations provided us with basic characteristics of long time evolution of rational solutions of the NLS equation in the case of near-breaking conditions. The analytic NLS solutions are found to describe the actual wave dynamics of steep waves reasonably well.

We report experiments on a millimeter-sized pendant drop of ferrofluid in a horizontal magnetic field. The initial drop size is chosen just below the breakup threshold under gravity. As the magnetic field is increased, the drop tilts in order to align with the direction of the total volume force that is exerted on it: weight plus magnetic force. The breakup is controlled by a generalized Bond number based on this total force and on the radius of the neck of the drop. The evolution of drop shape turns out to be a complex process governed by many parameters such as the angle between the total force and the needle, the drop size relative to the needle radius, and the wettability of the liquid on the needle material. This suggests a certain universality, that a single value of the critical Bond number is found regardless of magnetic fluid properties and whether the force is inclined or not.

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. With the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g., in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic and also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a lattice Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account in the time-dependent diffusion coefficient in a more complex porous medium.

Extreme events taking place on networks are not uncommon. We show that it is possible to manipulate the extreme event occurrence probabilities and distribution over the nodes of a scale-free network by tuning the nodal capacity. This can be used to reduce the number of extreme event occurrences. However, monotonic nodal capacity enhancements, beyond a point, do not lead to any substantial reduction in the number of extreme events. We point out the practical implication of this result for network design in the context of reducing extreme event occurrences.

A bundle of fibers has been considered here as a model for composite materials, where breaking of the fibers occur due to a combined influence of applied load (stress) and external noise. Through numerical simulation and a mean-field calculation we show that there exists a robust phase boundary between continuous (no waiting time) and intermittent fracturing regimes. In the intermittent regime, throughout the entire rupture process avalanches of different sizes are produced and there is a waiting time between two consecutive avalanches. The statistics of waiting times follows a $\Gamma$ distribution and the avalanche distribution shows power-law scaling, similar to what has been observed in the case of earthquake events and bursts in fracture experiments. We propose a prediction scheme that can tell when the system is expected to reach the continuous fracturing point from the intermittent phase.

Uniaxial compaction and breakage of rodlike particle packing has been studied using a discrete element method simulation. A scaling relationship between the applied stress, the number of breakage events, and the number-mean particle length has been derived and compared with computational experiments. Based on results for a wide range of intrinsic particle strengths and initial particle lengths, it seems that a single universal relation can be used to describe the incidence of breakage events during compaction of rodlike particle layers.

The photoinduced motion of an oil droplet on an aqueous solution under local irradiation by a green laser is reported. The results showed that a repulsive force is generated on pure water, while an attractive force is observed with an aqueous solution containing a surfactant. The driving force is discussed in terms of a thermal Marangoni effect. The switching on the photothermal effect is interpreted by taking into account the advection caused by the spatial gradient of the surface tension under local heating by a laser. A numerical model revealed that the geometrical profile of the surface tension around the droplet determines the mode of advection around the droplet and causes switching in the direction of migrations.

The zeros of the canonical partition functions for flexible square-well polymer chains have been approximately computed for chains up to length 256 for a range of square-well diameters. We have previously shown that such chain molecules can undergo a coil-globule and globule-crystal transition as well as a direct coil-crystal transition. Here we show that each of these transitions has a well-defined signature in the complex-plane map of the partition function zeros. The freezing transitions are characterized by nearly circular rings of uniformly spaced roots, indicative of a discontinuous transition. The collapse transition is signaled by the appearance of an elliptical horseshoe segment of roots that pinches down towards the positive real axis and defines a boundary to a root-free region of the complex plane. With increasing chain length, the root density on the circular ring and in the space adjacent to the elliptical boundary increases and the leading roots move towards the positive real axis. For finite-length chains, transition temperatures can be obtained by locating the intersection of the ellipse and/or circle of roots with the positive real axis. A finite-size scaling analysis is used to obtain transition temperatures in the long-chain (thermodynamic) limit. The collapse transition is characterized by crossover and specific-heat exponents of $\varphi \approx 0.76\left(2\right)$ and $\alpha \approx 0.66\left(2\right)$, respectively, consistent with a second-order phase transition.

Asymmetric patchy particle models have recently been shown to describe the crystallization of small globular proteins with near-quantitative accuracy. Here, we investigate how asymmetry in patch geometry and bond energy generally impacts the phase diagram and nucleation dynamics of this family of soft matter models. We find the role of the geometry asymmetry to be weak, but the energy asymmetry to markedly interfere with the crystallization thermodynamics and kinetics. These results provide a rationale for the success and occasional failure of the proposal of George and Wilson for protein crystallization conditions as well as physical guidance for developing more effective protein crystallization strategies.

Using an improved version of an evolutionary algorithm originally proposed by François and Hakim [Proc. Natl. Acad. Sci. USA 101, 580 (2004)], we generated small gene regulatory networks in which the concentration of a target protein oscillates in time. These networks may serve as candidates for oscillatory modules to be found in larger regulatory networks and protein interaction networks. The algorithm was run for ${10}^5$ times to produce a large set of oscillating modules, which were systematically classified and analyzed. The robustness of the oscillations against variations of the kinetic rates was also determined, to filter out the least robust cases. Furthermore, we show that the set of evolved networks can serve as a database of models whose behavior can be compared to experimentally observed oscillations. The algorithm found three smallest (core) oscillators in which nonlinearities and number of components are minimal. Two of those are two-gene modules: the mixed feedback loop, already discussed in the literature, and an autorepressed gene coupled with a heterodimer. The third one is a single gene module which is competitively regulated by a monomer and a dimer. The evolutionary algorithm also generated larger oscillating networks, which are in part extensions of the three core modules and in part genuinely new modules. The latter includes oscillators which do not rely on feedback induced by transcription factors, but are purely of post-transcriptional type. Analysis of post-transcriptional mechanisms of oscillation may provide useful information for circadian clock research, as recent experiments showed that circadian rhythms are maintained even in the absence of transcription.

This paper studies international migration from a complex-network perspective. We define the international migration network (IMN) as the weighted-directed graph where nodes are world countries and links account for the stock of migrants originated in a given country and living in another country at a given point in time. We characterize the binary and weighted architecture of the network and its evolution over time in the period 1960–2000. We find that the IMN is organized around a modular structure with a small-world binary pattern displaying disassortativity and high clustering, with power-law distributed weighted-network statistics. We also show that a parsimonious gravity model of migration can account for most of observed IMN topological structure. Overall, our results suggest that socioeconomic, geographical, and political factors are more important than local-network properties in shaping the structure of the IMN.

To motivate more people to participate in vaccination campaigns, various subsidy policies are often supplied by government and the health sectors. However, these external incentives may also alter the vaccination decisions of the broader public, and hence the choice of incentive needs to be carefully considered. Since human behavior and the networking-constrained interactions among individuals significantly impact the evolution of an epidemic, here we consider the voluntary vaccination on human contact networks. To this end, two categories of typical subsidy policies are considered: (1) under the free subsidy policy, the total amount of subsidy is distributed to a certain fraction of individual and who are vaccinated without personal cost, and (2) under the partial-offset subsidy policy, each vaccinated person is offset by a certain amount of subsidy. A vaccination decision model based on evolutionary game theory is established to study the effects of these different subsidy policies on disease control. Simulations suggest that, because the partial-offset subsidy policy encourages more people to take vaccination, its performance is significantly better than that of the free subsidy policy. However, an interesting phenomenon emerges in the partial-offset scenario: with limited amount of total subsidy, a moderate subsidy rate for each vaccinated individual can guarantee the group-optimal vaccination, leading to the maximal social benefits, while such an optimal phenomenon is not evident for the free subsidy scenario.

We study sign changes and scaling laws in the Cartesian components of the velocity and vorticity of rotating turbulence, in the helicity, and in the components of vertically averaged fields. Data for the analysis are provided by high-resolution direct numerical simulations of rotating turbulence with different forcing functions, with up to ${1536}^3$ grid points, with Reynolds numbers between $\approx 1100$ and $\approx 5100$, and with moderate Rossby numbers between $\approx 0.06$ and $\approx 8$. When rotation is negligible, all Cartesian components of the velocity show similar scaling, in agreement with the expected isotropy of the flow. However, in the presence of rotation, only the vertical components of the fields show clear scaling laws, with evidence of possible sign singularity in the limit of an infinite Reynolds number. Horizontal components of the velocity are smooth and do not display rapid fluctuations for arbitrarily small scales. The vertical velocity and vorticity, as well as the vertically averaged vertical velocity and vorticity, show the same scaling within error bars, in agreement with theories that predict that these quantities have the same dynamical equation for very strong rotation.

We investigate the transport of periodic trains of droplets through microfluidic networks having one inlet, one outlet, and nodes consisting of $\mathsf{T}$ junctions. Variations of the dilution of the trains, i.e., the distance between drops, reveal the existence of various hydrodynamic regimes characterized by the number of preferential paths taken by the drops. As the dilution increases, this number continuously decreases until only one path remains explored. Building on a continuous approach used to treat droplet traffic through a single asymmetric loop, we determine selection rules for the paths taken by the drops and we predict the variations of the fraction of droplets taking these paths with the parameters at play including the dilution. Our results show that as dilution decreases, the paths are selected according to the ascending order of their hydrodynamic resistance in the absence of droplets. The dynamics of these systems controlled by time-delayed feedback is complex: We observe a succession of periodic regimes separated by a wealth of bifurcations as the dilution is varied. In contrast to droplet traffic in single asymmetric loops, the dynamical behavior in networks of loops is sensitive to initial conditions because of extra degrees of freedom.