Asymmetry of magnetic objects in a fluid under an oscillating magnetic field leads to a wealth of nonequilibrium dynamics phenomena including a novel ratchet effect without an asymmetric substrate. These nonlinear dynamics are explained in the framework of the Stokes’ model by a drag coefficient, which depends on the direction of motion. This approach is general and is independent of the physical mechanism responsible for this directional dependence of the drag coefficient as well as the size of the object. The theoretical model is experimentally verified for two systems, a nonrigid magnetic microcoil and a chiral magnetic macroobject immersed in a bounded fluid.

We investigate a simplified model of two fully connected magnetic systems maintained at different temperatures by virtue of being connected to two independent thermal baths while simultaneously being interconnected with each other. Using generating functional analysis, commonly used in statistical mechanics, we find exactly soluble expressions for their individual magnetization that define a two-dimensional nonlinear map, the equations of which have the same form as those obtained for densely connected equilibrium systems. Steady states correspond to the fixed points of this map, separating the parameter space into a rich set of nonequilibrium phases that we analyze in asymptotically high and low (nonequilibrium) temperature limits. The theoretical formalism is shown to revert to the classical nonequilibrium steady state problem for two interacting systems with a nonzero heat transfer between them that catalyzes a phase transition between ambient nonequilibrium states.

Transitions to absorbing states are of fundamental importance in nonequilibrium physics as well as ecology. In ecology, absorbing states correspond to the extinction of species. We here study the spatial population dynamics of three cyclically interacting species. The interaction scheme comprises both direct competition between species as in the cyclic Lotka-Volterra model, and separated selection and reproduction processes as in the May-Leonard model. We show that the dynamic processes leading to the transient maintenance of biodiversity are closely linked to attractors of the nonlinear dynamics for the overall species’ concentrations. The characteristics of these global attractors change qualitatively at certain threshold values of the mobility and depend on the relative strength of the different types of competition between species. They give information about the scaling of extinction times with the system size and thereby the stability of biodiversity. We define an effective free energy as the negative logarithm of the probability to find the system in a specific global state before reaching one of the absorbing states. The global attractors then correspond to minima of this effective energy landscape and determine the most probable values for the species’ global concentrations. As in equilibrium thermodynamics, qualitative changes in the effective free energy landscape indicate and characterize the underlying nonequilibrium phase transitions. We provide the complete phase diagrams for the population dynamics and give a comprehensive analysis of the spatio-temporal dynamics and routes to extinction in the respective phases.

We study effect of weak noise on the dynamics of a hair bundle model near the excitability threshold and near a subcritical Hopf bifurcation. We analyze numerically noise-induced structural changes in the probability density and the power spectral density of the model. In particular, we show that weak noise can induce oscillations with two distinct frequencies in both excitable and limit-cycle regimes. We then applied a recently developed technique of stochastic sensitivity functions which allows us to estimate threshold values of noise intensity corresponding to these transitions.

The invasion of air into porous systems in drying processes is often localized in soft materials, such as colloidal suspensions and granular pastes, and it typically develops in the form of cracks before ordinary drying begins. To investigate such processes, we construct an invasion percolation model on a deformable lattice for cohesionless elastic systems, and with this model we derive the condition under which cracking occurs. A Griffith-like condition characterized by a dimensionless parameter is proposed, and its validity is checked numerically. This condition indicates that the ease with which cracking occurs increases as the particles composing the material become smaller, as the rigidity of the system increases and as the degree of heterogeneity characterizing the drying processes decreases.

We study, both analytically and numerically, the dynamics of mechanical oscillators kept in motion by a feedback force, which is generated electronically from a signal produced by the oscillators themselves. This kind of self-sustained systems may become standard in the design of frequency-control devices at microscopic scales. Our analysis is thus focused on their synchronization properties under the action of external forces and on the joint dynamics of two to many coupled oscillators. Existence and stability of synchronized motion are assessed in terms of the mechanical properties of individual oscillators, namely, their natural frequencies and damping coefficients, and synchronization frequencies are determined. Similarities and differences with synchronization phenomena in other coupled oscillating systems are emphasized.

Complex dynamics associated with multistability have been studied extensively in the past but mostly for low-dimensional nonlinear dynamical systems. A question of fundamental interest is whether multistability can arise in high-dimensional physical systems. Motivated by the ever increasing widespread use of nanoscale systems, we investigate a prototypical class of nanoelectromechanical systems: electrostatically driven Si nanowires, mathematically described by a set of driven, nonlinear partial differential equations. We develop a computationally efficient algorithm to solve the equations. Our finding is that multistability and complicated structures of basins of attraction are common types of dynamics, and the latter can be attributed to extensive transient chaos. Implications of these phenomena to device operations are discussed.

Since Kolmogorov proposed his phenomenological theory of hydrodynamic turbulence in 1941, the description of the mechanism leading to the energy cascade and anomalous scaling remains an open problem in fluid mechanics. Soon after, in 1949, Onsager noticed that the scaling properties in the inertial range imply nondifferentiability of the velocity field in the limit of vanishing viscosity. This observation suggests that the turbulence mechanism may be related to a finite-time singularity (blowup) of incompressible Euler equations. However, the existence of such blowup is still an open problem too. In this paper, we show that the blowup indeed represents the driving mechanism of the inertial range for a simplified (shell) model of turbulence. Here, blowups generate coherent structures (instantons), which travel through the inertial range in finite time and are described by universal self-similar statistics. The anomaly (deviation of scaling exponents of velocity moments from the Kolmogorov theory) is related analytically to the process of instanton creation using the large deviation principle. The results are confirmed by numerical simulations.

Horseshoe distributions are shell-like particle distributions that can arise in space and laboratory plasmas when particle beams propagate into increasing magnetic fields. The present paper studies the stability and the dynamics of wave packets interacting resonantly with electrons presenting low-energy horseshoe or shell-type velocity distributions in a magnetized plasma. The linear instability growth rates are determined as a function of the ratio of the plasma to the cyclotron frequencies, of the velocity and the opening angle of the horseshoe, and of the relative thickness of the shell. The nonlinear stage of the instability is investigated numerically using a symplectic code based on a three-dimensional Hamiltonian model. Simulation results show that the dynamics of the system is mainly governed by wave-particle interactions at Landau and normal cyclotron resonances and that the high-order normal cyclotron resonances play an essential role. Specific features of the dynamics of particles interacting simultaneously with two or more waves at resonances of different natures and orders are discussed, showing that such complex processes determine the main characteristics of the wave spectrum's evolution. Simulations with wave packets presenting quasicontinuous spectra provide a full picture of the relaxation of the horseshoe distribution, revealing two main phases of the evolution: an initial stage of wave energy growth, characterized by a fast filling of the shell, and a second phase of slow damping of the wave energy, accompanied by final adjustments of the electron distribution. The influence of the density inhomogeneity along the horseshoe on the wave-particle dynamics is also discussed.

We present extensive Monte Carlo simulations on a two-dimensional $XY$ model with a modified form of interaction potential. Thermodynamic quantities other than energy, specific heat, etc. (such as magnetization, susceptibility, and fourth-order cumulant of magnetization) are obtained using multiple histogram reweighting of the data obtained from the simulations. We employ an approach which eliminates the need to construct two-dimensional histograms. This approach makes judicious use of computer memory as well as CPU time. Lee-Kosterlitz's method of finite size scaling for a first-order transition and analysis using Binder's cumulant method allow us to make an accurate determination of the transition temperature.

We investigate non-Brownian particles suspended in a periodic shear-flow using simulations. Following Metzger and Butler [Phys. Rev. E 82, 051406 (2010)], we show that the chaotic dynamics arising from lubrication interactions are too weak to generate an observable particle dispersion. The irreversibility observed in periodic flow is dominated by contact interactions. Nonetheless, we show that lubrication interactions must be included in the calculation to obtain results that agree with experiments.

We study both theoretically and experimentally the electro-optical properties of vertically aligned deformed helix ferroelectric liquid crystals (VADHFLC) with subwavelength pitch that are governed by the electrically induced optical biaxiality of the smectic helical structure. The key theoretical result is that the principal refractive indices of homogenized VADHFLC cells exhibit the quadratic nonlinearity and such behavior might be interpreted as an orientational Kerr effect caused by the electric-field-induced orientational distortions of the FLC helix. In our experiments, it has been observed that, for sufficiently weak electric fields, the magnitude of biaxiality is proportional to the square of electric field in good agreement with our theoretical results for the effective dielectric tensor of VADHFLCs. Under certain conditions, the 2$\pi$ phase modulation of light, which is caused by one of the induced refractive indices, is observed without changes in ellipticity of incident light.

Long-distance bidirectional transport of organelles depends on the coordinated motion of various motor proteins on the cytoskeleton. Recent quantitative live cell imaging in the elongated hyphal cells of Ustilago maydis has demonstrated that long-range motility of motors and their endosomal cargo occurs on unipolar microtubules (MTs) near the extremities of the cell. These MTs are bundled into antipolar bundles within the central part of the cell. Dynein and kinesin-3 motors coordinate their activity to move early endosomes (EEs) in a bidirectional fashion where dynein drives motility towards MT minus ends and kinesin towards MT plus ends. Although this means that one can easily assign the drivers of bidirectional motion in the unipolar section, the bipolar orientations in the bundle mean that it is possible for either motor to drive motion in either direction. In this paper we use a multilane asymmetric simple exclusion process modeling approach to simulate and investigate phases of bidirectional motility in a minimal model of an antipolar MT bundle. In our model, EE cargos (particles) change direction on each MT with a turning rate $\Omega$ and there is switching between MTs in the bundle at the minus ends. At these ends, particles can hop between MTs with rate $q_1$ on passing from a unipolar to a bipolar section (the obstacle-induced switching rate) or $q_2$ on passing in the other direction (the end-induced switching rate). By a combination of numerical simulations and mean-field approximations, we investigate the distribution of particles along the MTs for different values of these parameters and of $\Theta$, the overall density of particles within this closed system. We find that even if $\Theta$ is low, the system can exhibit a variety of phases with shocks in the density profiles near plus and minus ends caused by queuing of particles. We discuss how the parameters influence the type of particle that dominates active transport in the bundle.

Airplane boarding process is an example where disorder properties of the system are relevant to the emergence of universality classes. Based on a simple model, we present a systematic analysis of finite-size effects in boarding time, and propose a comprehensive view of the role of sequential disorder in the scaling behavior of boarding time against the plane size. Using numerical simulations and mathematical arguments, we find how the scaling behavior depends on the number of seat columns and the range of sequential disorder. Our results show that new scaling exponents can arise as disorder is localized to varying extents.

The robustness of a network of networks (NON) under random attack has been studied recently [Gao et al., Phys. Rev. Lett. 107, 195701 (2011)]. Understanding how robust a NON is to targeted attacks is a major challenge when designing resilient infrastructures. We address here the question how the robustness of a NON is affected by targeted attack on high- or low-degree nodes. We introduce a targeted attack probability function that is dependent upon node degree and study the robustness of two types of NON under targeted attack: (i) a tree of $n$ fully interdependent Erdős-Rényi or scale-free networks and (ii) a starlike network of $n$ partially interdependent Erdős-Rényi networks. For any tree of $n$ fully interdependent Erdős-Rényi networks and scale-free networks under targeted attack, we find that the network becomes significantly more vulnerable when nodes of higher degree have higher probability to fail. When the probability that a node will fail is proportional to its degree, for a NON composed of Erdős-Rényi networks we find analytical solutions for the mutual giant component $P_{\infty }$ as a function of $p$, where $1-p$ is the initial fraction of failed nodes in each network. We also find analytical solutions for the critical fraction $p_c$, which causes the fragmentation of the $n$ interdependent networks, and for the minimum average degree ${\overline{k}}_{\min }$ below which the NON will collapse even if only a single node fails. For a starlike NON of $n$ partially interdependent Erdős-Rényi networks under targeted attack, we find the critical coupling strength $q_c$ for different $n$. When $q>q_c$, the attacked system undergoes an abrupt first order type transition. When $q\le q_c$, the system displays a smooth second order percolation transition. We also evaluate how the central network becomes more vulnerable as the number of networks with the same coupling strength $q$ increases. The limit of $q=0$ represents no dependency, and the results are consistent with the classical percolation theory of a single network under targeted attack.

We study the quantum dynamics of the central-depleted-well (CDW) regime in a three-mode Bose Hubbard model subject to a confining parabolic potential. By introducing a suitable set of momentum-like modes we identify the microscopic variables involved in the quantization process and the dynamical algebra of the model. We describe the diagonalization procedure showing that the model reduces to a double oscillator. Interestingly, we find that the parameter-space domain where this scheme entails a discrete spectrum well reproduces the two regions where the classical trimer excludes unstable oscillations. Spectral properties are examined in different limiting cases together with various delocalization effects. These are shown to characterize quantum states of the CDW regime in the proximity of the borderline with classically unstable domains.

The growth behavior of stationary axisymmetric vortices and of oscillatory, nonaxisymmetric spiral vortices in Taylor-Couette flow of a ferrofluid in between differentially rotating cylinders is analyzed using a numerical linear stability analysis. The investigation is done as a function of the inner and outer cylinder's rotation rates, the axial wave number of the vortex flows, and the magnitude of an applied homogeneous axial magnetic field. In particular, the consequences of incorporating elongational flow effects in the magnetization balance equation on the marginal control parameters that separate growth from decay behavior are determined. That is done for several values of the transport coefficient that measures the strength of these effects.

Self-assembly of structures from vertically aligned, charged dust particle bundles within a glass box placed on the lower, powered electrode of a Gaseous Electronics Conference rf reference cell were produced and examined experimentally. Self-organized formation of one-dimensional vertical chains, two-dimensional zigzag structures, and three-dimensional helical structures of triangular, quadrangular, pentagonal, hexagonal, and heptagonal symmetries are shown to occur. System evolution is shown to progress from a one-dimensional chain structure, through a zigzag transition to a two-dimensional, spindlelike structure, and then to various three-dimensional, helical structures exhibiting multiple symmetries. Stable configurations are found to be dependent upon the system confinement, ${\gamma }^2={\left({\omega }_{0h}/{\omega }_{0v}\right)}^2$ (where ${\omega }_{0h,v}$ are the horizontal and vertical dust resonance frequencies), the total number of particles within a bundle, and the rf power. For clusters having fixed numbers of particles, the rf power at which structural phase transitions occur is repeatable and exhibits no observable hysteresis. The critical conditions for these structural phase transitions as well as the basic symmetry exhibited by the one-, two-, and three-dimensional structures that subsequently develop are in good agreement with the theoretically predicted configurations of minimum energy determined employing molecular dynamics simulations for charged dust particles confined in a prolate, spheroidal potential as presented theoretically by Kamimura and Ishihara [Kamimura and Ishihara, Phys. Rev. E 85, 016406 (2012)].

This work compares computational methods for laser pulse propagation in hollow waveguides filled with rare gases at high pressures, with applications in extreme nonlinear optics in the midinfrared wavelength region. As the wavelength of light $\lambda =2\pi /k$ increases with respect to the transverse size $R$ of a leaky waveguide, the loss of light out of the waveguide upon propagation, in general, increases. The now standard numerical approach for studying such structures is based on expansion of the propagating field into approximate leaky waveguide modes. We compare this approach to an improved method that resolves the electric field in real space and correctly captures the energy loss through the waveguide wall. The comparison reveals that the expansion-based approach overestimates losses that occur in nonlinearly reshaped pulsed waveforms. For a modest increase in computational effort, the alternate method offers a physically more accurate model to describe phenomena (e.g., extreme pulse-selfcompression) in waveguides with smaller values of $kR$.