Spindle-shaped nematic droplets (tactoids) form in solutions of rod-like molecules at the onset of the liquid crystalline phase. Their unique shape and internal structure result from the interplay of the elastic deformation of the nematic and anisotropic surface forces. The balance of these forces dictates that tactoids must display a continuous variation in aspect ratio and director-field configuration. Yet, such continuous transition has eluded observation for decades: tactoids have displayed either a bipolar configuration with particles aligned parallel to the droplet interface or a homogeneous configuration with particles aligned parallel to the long axis of the tactoid. Here, we report the first observation of the continuous transition in shape and director-field configuration of tactoids in true solutions of carbon nanotubes in chlorosulfonic acid. This observation is possible because the exceptional length of carbon nanotubes shifts the transition to a size range that can be visualized by optical microscopy. Polarization micrographs yield the interfacial and elastic properties of the system. Absorbance anisotropy measurements provide the highest nematic order parameter $\left(S=0.79\right)$ measured to date for a nematic phase of carbon nanotubes at coexistence with its isotropic phase.

The influence of networks topology on collective properties of dynamical systems defined upon it is studied in the thermodynamic limit. A network model construction scheme is proposed where the number of links and the average eccentricity are controlled. This is done by rewiring links of a regular one-dimensional chain according to a probability $p$ within a specific range $r$ that can depend on the number of vertices $N$. We compute the thermodynamical behavior of a system defined on the network, the $XY$-rotors model, and monitor how it is affected by the topological changes. We identify the network effective dimension $d$ as a crucial parameter: topologies with $d<2$ exhibit no phase transitions, while topologies with $d>2$ display a second-order phase transition. Topologies with $d=2$ exhibit states characterized by infinite susceptibility and macroscopic chaotic, turbulent dynamical behavior. These features are also captured by $d$ in the finite size context.

Social exclusion, keeping free riders from benefit sharing, plays an important role in sustaining cooperation in our world. Here we propose two different exclusion regimes, namely, peer exclusion and pool exclusion, to investigate the evolution of social exclusion in finite populations. In the peer exclusion regime, each excluder expels all the defectors independently, and thus bears the total cost on his own, while in the pool exclusion regime, excluders spontaneously form an institution to carry out rejection of the free riders, and each excluder shares the cost equally. In a public goods game containing only excluders and defectors, it is found that peer excluders outperform pool excluders if the exclusion costs are small, and the situation is converse once the exclusion costs exceed some critical points, which holds true for all the selection intensities and different update rules. Moreover, excluders can dominate the whole population under a suitable parameters range in the presence of second-order free riders (cooperators), showing that exclusion has prominent advantages over common costly punishment. More importantly, our finding indicates that the group exclusion mechanism helps the cooperative union to survive under unfavorable conditions. Our results may give some insights into better understanding the prevalence of such a strategy in the real world and its significance in sustaining cooperation.

The net steady state flow pattern of a distorting sphere is studied in the framework of the bilinear theory of swimming at low Reynolds number. It is argued that the starting point of a theory of interacting active particles should be based on such a calculation, since any arbitrarily chosen steady state flow pattern is not necessarily the result of a swimming motion. Furthermore, it is stressed that as a rule the phase of stroke is relevant in hydrodynamic interactions, so that the net flow pattern must be used with caution.

One of the intriguing questions in fluid dynamics is on the interrelation between dynamic singularities in the solutions of fluid dynamic equations – unboundedness of the velocity field in an appropriate norm – and the geometric ones – divergence of curvature at fluid interfaces. The present work focuses on two generic interfacial singularities – genuine cusps and cuspidal edges – found here in both two and three dimensions thus establishing a relation between real fluid interfaces and geometric singularity theory. The key finding is the necessary condition for the existence of geometric singularities, which is a variation of surface tension. It is also established here that the dynamic and geometric singularities entail each other only in the case of three-dimensional cusps. Explicit asymptotic solutions for the flow field and interface shape near steady-state singularities at fluid interfaces are developed as well. The practical motivation for the present study comes from the fundamental role interfacial singularities play in sustaining self-driven conversion of chemical into mechanical energy.

A theory for nonlinear excitations in quantum plasmas is presented for narrow-gap semiconductors by considering the combined effects of quantum and pseudorelativity. The system is governed by a coupled Klein-Gordon equation for the collective wave functions of the conduction electrons and Poisson's equation for the electrostatic potential. This gives a closed system, including the effects of charge separation, quantum tunneling, and pseudorelativity. By choosing the typical parameters of semiconductor InSb, the quasistationary soliton solution, which is a multipeaked dark soliton, is obtained numerically and shows depleted electron densities correlated with a localized potential. The dynamical simulation result shows that the dark soliton is stable and has a multipeaked profile, which is consistent with the quasistationary solution. The present model and results may be useful in understanding the nonlinear properties of semiconductor plasma on an ultrafast time scale.

To probe both the hydrodynamic nonequilibrium (HNE) and thermodynamic nonequilibrium (TNE) in the combustion process, a two-dimensional multiple-relaxation-time (MRT) version of lattice Boltzmann kinetic model (LBKM) for combustion phenomena is presented. The chemical energy released in the progress of combustion is dynamically coupled into the system by adding a chemical term to the LB kinetic equation. Aside from describing the evolutions of the conserved quantities, the density, momentum, and energy, which are what the Navier-Stokes model describes, the MRT-LBKM presents also a coarse-grained description on the evolutions of some nonconserved quantities. The current model works for both subsonic and supersonic flows with or without chemical reaction. In this model, both the specific-heat ratio and the Prandtl number are flexible, the TNE effects are naturally presented in each simulation step. The model is verified and validated via well-known benchmark tests. As an initial application, various nonequilibrium behaviors, including the complex interplays between various HNEs, between various TNEs, and between the HNE and TNE, around the detonation wave in the unsteady and steady one-dimensional detonation processes are preliminarily probed. It is found that the system viscosity (or heat conductivity) decreases the local TNE, but increases the global TNE around the detonation wave, that even locally, the system viscosity (or heat conductivity) results in two kinds of competing trends, to increase and to decrease the TNE effects. The physical reason is that the viscosity (or heat conductivity) takes part in both the thermodynamic and hydrodynamic responses.

The stresslet provides a simple model of the flow created by a small, freely swimming and neutrally buoyant aquatic organism and shows that the far field fluid disturbance created by such an organism in general decays as one over distance squared. Here we discuss a quieter swimming mode that eliminates the stresslet component of the flow and leads to a faster spatial decay of the fluid disturbance described by a force quadrupole that decays as one over distance cubed. Motivated by recent experimental results on fluid disturbances due to small aquatic organisms, we demonstrate that a three-Stokeslet model of a swimming organism which uses breast stroke type kinematics is an example of such a quiet swimmer. We show that the fluid disturbance in both the near field and the far field is significantly reduced by appropriately arranging the propulsion apparatus, and we find that the far field power laws are valid surprisingly close to the organism. Finally, we discuss point force models as a general framework for hypothesis generation and experimental exploration of fluid mediated predator-prey interactions in the planktonic world.

We consider a model of a stepping molecular motor consisting of two connected heads. Directional motion of the stepper takes place along a one-dimensional track. Each head is subject to a periodic potential without spatial reflection symmetry. When the potential for one head is switched on, it is switched off for the other head. Additionally, the system is subject to the influence of symmetric, white Lévy noise that mimics the action of external random forcing. The stepper exhibits motion with a preferred direction which is examined by analyzing the median of the displacement of a midpoint between the positions of the two heads. We study the modified dynamics of the stepper by numerical simulations. We find flux reversals as noise parameters are changed. Speed and direction appear to very sensitively depend on characteristics of the noise.

Evolutionary game theory is crucial to capturing the characteristic interaction patterns among selfish individuals. In a population of coordination games of two strategies, one of the central problems is to determine the fixation probability that the system reaches a state of networkwide of only one strategy, and the corresponding expectation times. The deterministic replicator equations predict the critical value of initial density of one strategy, which separates the two absorbing states of the system. However, numerical estimations of this separatrix show large deviations from the theory in finite populations. Here we provide a stochastic treatment of this dynamic process on complex networks of finite sizes as Markov processes, showing the evolutionary time explicitly. We describe analytically the effects of network structures on the intermediate fixations as observed in numerical simulations. Our theoretical predictions are validated by various simulations on both random and scale free networks. Therefore, our stochastic framework can be helpful in dealing with other networked game dynamics.

The eigenvalue spectrum of the matrix of directed weights defining a neural network model is informative of several stability and dynamical properties of network activity. Existing results for eigenspectra of sparse asymmetric random matrices neglect spatial or other constraints in determining entries in these matrices, and so are of partial applicability to cortical-like architectures. Here we examine a parameterized class of networks that are defined by sparse connectivity, with connection weighting modulated by physical proximity (i.e., asymmetric Euclidean random matrices), modular network partitioning, and functional specificity within the excitatory population. We present a set of analytical constraints that apply to the eigenvalue spectra of associated weight matrices, highlighting the relationship between connectivity rules and classes of network dynamics.

Extreme intensity pulses sharing statistical properties similar to rogue waves have been recently observed in a laser diode with phase-conjugate feedback [A. Karsaklian Dal Bosco, D. Wolfersberger, and M. Sciamanna, Opt. Lett. 38, 703 (2013)], but remain unexplained. We demonstrate here that a rate equation model of a laser diode that includes an instantaneous phase-conjugate feedback field reproduces qualitatively well the statistical features of these extreme events as identified in the experiment, i.e., the deviation of the intensity statistics to a Gaussian-shape statistics and the statistics of the time separating extreme events. The numerical simulations confirm the importance of the feedback strength in increasing the number of such extreme events and allow us to explain how extreme events emerge from a sequence of bifurcations on self-pulsating solutions, the so-called external cavity modes.

We report the inherently three-dimensional linear instabilities of a propulsive wake, produced by a flapping foil, mimicking the caudal fin of a fish or the wing of a flying animal. For the base flow, three sequential wake patterns appear as we increase the flapping amplitude: Bénard-von Kármán (BvK) vortex streets; reverse BvK vortex streets; and deflected wakes. Imposing a three-dimensional spanwise periodic perturbation, we find that the resulting Floquet multiplier $\left|\mu \right|$ indicates an unstable “short wavelength” mode at wave number $\beta =30$, or wavelength $\lambda =0.21$ (nondimensionalized by the chord length) at sufficiently high flow Reynolds number $\text{Re}=Uc/\nu \simeq 600$, where $U$ is the upstream flow velocity, $c$ is the chord length, and $\nu$ is the kinematic viscosity of the fluid. Another, “long wavelength” mode at $\beta =6$ ($\lambda =1.05$) becomes critical at somewhat higher Reynolds number, although we do not expect that this mode would be observed physically because its growth rate is always less than the short wavelength mode, at least for the parameters we have considered. The long wavelength mode has certain similarities with the so-called mode $A$ in the drag wake of a fixed bluff body, while the short wavelength mode appears to have a period of the order of twice that of the base flow, in that its structure seems to repeat approximately only every second cycle of the base flow. Whether it is appropriate to classify this mode as a truly subharmonic mode or as a quasiperiodic mode is still an open question however, worthy of a detailed parametric study with various flapping amplitudes and frequencies.

In Sec. II of the named paper, the authors study the planar model on a regular two-dimensional (square or triangular) lattice and find evidence of an ordering transition at finite temperatures in both cases. It is shown that their findings do not agree with, and appear to ignore, a number of mathematical results, which have been known in the literature for some decades to date, and entail orientational disorder at all finite temperatures, as well as the existence of the Berezinski–Kosterlitz-Thouless transition.

In Romano's Comment [Phys. Rev. E 91, 046501 (2015)] on Sec. II of our paper [Phys. Rev. E 90, 022502 (2014)], the author suggests that our findings concerning the nature of the ordering transition of our modified Lebwohl-Lasher model with two-dimensional planar rotators on a planar lattice are inconsistent with known mathematical results. We argue in this Reply that our findings are in fact in agreement with previous mathematical and simulation results and that the criticisms raised by Romano have no impact on the results presented in our paper.

Granger causality has long been a prominent method for inferring causal interactions between stochastic variables for a broad range of complex physical systems. However, it has been recognized that a moving average (MA) component in the data presents a serious confound to Granger causal analysis, as routinely performed via autoregressive (AR) modeling. We solve this problem by demonstrating that Granger causality may be calculated simply and efficiently from the parameters of a state-space (SS) model. Since SS models are equivalent to autoregressive moving average models, Granger causality estimated in this fashion is not degraded by the presence of a MA component. This is of particular significance when the data has been filtered, downsampled, observed with noise, or is a subprocess of a higher dimensional process, since all of these operations—commonplace in application domains as diverse as climate science, econometrics, and the neurosciences—induce a MA component. We show how Granger causality, conditional and unconditional, in both time and frequency domains, may be calculated directly from SS model parameters via solution of a discrete algebraic Riccati equation. Numerical simulations demonstrate that Granger causality estimators thus derived have greater statistical power and smaller bias than AR estimators. We also discuss how the SS approach facilitates relaxation of the assumptions of linearity, stationarity, and homoscedasticity underlying current AR methods, thus opening up potentially significant new areas of research in Granger causal analysis.

We generalize Morita's works [J. Phys. A 7, 289 (1974); J. Phys. A 7, 1613 (1974)] on ground states of Ising chains, for chains with a periodic structure and different spins, to any interaction order. The main assumption is translational invariance. The length of the irreducible blocks is a multiple of the period of the chain. If there is parity invariance, it restricts the length in general only in the diatomic case. There are degenerated states and under certain circumstances there could be nonregular ground states. We illustrate the results and give the ground state diagrams in several cases.

In this work we analyze the simultaneous emergence of diffusive energy transport and local thermalization in a nonequilibrium one-dimensional quantum system, as a result of integrability breaking. Specifically, we discuss the local properties of the steady state induced by thermal boundary driving in a $XXZ$ spin chain with staggered magnetic field. By means of efficient large-scale matrix product simulations of the equation of motion of the system, we calculate its steady state in the long-time limit. We start by discussing the energy transport supported by the system, finding it to be ballistic in the integrable limit and diffusive when the staggered field is finite. Subsequently, we examine the reduced density operators of neighboring sites and find that for large systems they are well approximated by local thermal states of the underlying Hamiltonian in the nonintegrable regime, even for weak staggered fields. In the integrable limit, on the other hand, this behavior is lost, and the identification of local temperatures is no longer possible. Our results agree with the intuitive connection between energy diffusion and thermalization.

In the usual Achlioptas processes the smallest clusters of a few randomly chosen ones are selected to merge together at each step. The resulting aggregation process leads to the delayed birth of a giant cluster and the so-called explosive percolation transition showing a set of anomalous features. We explore a process with the opposite selection rule, in which the biggest clusters of the randomly chosen ones merge together. We develop a theory of this kind of percolation based on the Smoluchowsky equation, find the percolation threshold, and describe the scaling properties of this continuous transition, namely, the critical exponents and amplitudes, and scaling functions. We show that, qualitatively, this transition is similar to the ordinary percolation one, though occurring in less connected systems.