We investigate synchronization and plateau splitting of coupled oscillators on a one-dimensional lattice with long-range interactions that decay over distance as a power law. We show that in the thermodynamic limit the dynamics of systems of coupled oscillators with power-law exponent $\alpha \le 1$ is identical to that of the all-to-all coupling case. For $\alpha >1$, oscillatory behavior of the phase coherence appears as a result of single plateau splitting into multiple plateaus. A coarse-graining method is used to investigate the onset of plateau splitting. We analyze a simple oscillatory state formed by two plateaus in detail and propose a systematic approach to predict the onset of plateau splitting. The prediction of breaking points of plateau splitting is in quantitatively good agreement with numerical simulations.

Defibrillation is the standard clinical treatment used to stop ventricular fibrillation. An electrical device delivers a controlled amount of electrical energy via a pair of electrodes in order to reestablish a normal heart rate. We propose a technique that is a combination of biphasic shocks applied with a four-electrode system rather than the standard two-electrode system. We use a numerical model of a one-dimensional ring of cardiac tissue in order to test and evaluate the benefit of this technique. We compare three different shock protocols, namely a monophasic and two types of biphasic shocks. The results obtained by using a four-electrode system are compared quantitatively with those obtained with the standard two-electrode system. We find that a huge reduction in defibrillation threshold is achieved with the four-electrode system. For the most efficient protocol (asymmetric biphasic), we obtain a reduction in excess of 80% in the energy required for a defibrillation success rate of 90%. The mechanisms of successful defibrillation are also analyzed. This reveals that the advantage of asymmetric biphasic shocks with four electrodes lies in the duration of the cathodal and anodal phase of the shock.

Self-propulsion of cellular microswimmers generates flow signatures, commonly classified as pusher and puller type, which characterize hydrodynamic interactions with other cells or boundaries. Using experimentally measured beat patterns, we compute that the flagellated green alga Chlamydomonas oscillates between pusher and puller, rendering it an approximately neutral swimmer, when averaging over its full beat cycle. Beyond a typical distance of $100\mu \mathrm{m}$ from the cell, inertia attenuates oscillatory microflows. We show that hydrodynamic interactions between cells oscillate in time and are of similar magnitude as stochastic swimming fluctuations. From our analysis, we also find that the rate of hydrodynamic dissipation varies in time, which implies that flagellar beat patterns are not optimized with respect to this measure.

An analysis has been performed for the osmotic pressure of ionic liquids in the electric double layer (EDL). By using the electromechanical approach, we first derive a differential equation that is valid for computing the osmotic pressure in the continuum limit of any incompressible fluid in EDL. Then a specific model for ionic liquids proposed by Bazant et al. [M. Z. Bazant, B. D. Storey, and A. A. Kornyshev, Phys. Rev. Lett. 106, 046102 (2011)] is adopted for more detailed computation of the osmotic pressure. Ionic liquids are characterized by the correlation and the steric effects of ions and their effects are analyzed. In the low voltage cases, the correlation effect is dominant and the problem becomes linear. For this low voltage limit, a closed form formula is derived for predicting the osmotic pressure in EDL with no overlapping. It is found that the osmotic pressure decreases as the correlation effect increases. The osmotic pressures at the nanoslit surface and nanoslit centerline are also obtained for the low voltage limit. For the cases of moderately high voltage with high correlation factor, approximate formulas are derived for estimating osmotic pressure values based on the concept of a condensed layer near the electrode. In order to corroborate the results predicted by analytical studies, the full nonlinear model has been solved numerically.

A mechanical analogy is used to analyze the interaction between the magnetic field, electric current, and deformation of interfaces in liquid metal batteries. In the framework of a low-mode, nondissipative, linear stability model, it is found that, during charging or discharging, a sufficiently large battery is prone to instabilities of two types. One is similar to the metal pad instability known to exist in the aluminum reduction cells. Another type is new. It is related to the destabilizing effect of the Lorentz force formed by the azimuthal magnetic field induced by the base current, and the current perturbations caused by the local variations of the thickness of the electrolyte layer.

In regimes far beyond the wave-breaking threshold of Raman amplification, we show that significant amplification can occur after the onset of wave breaking, before phase mixing destroys the coherent coupling between pump, probe, and plasma wave. Amplification in this regime is therefore a transient effect, with the higher-efficiency “coherent wave-breaking” (CWB) regime accessed by using a short, intense probe. Parameter scans illustrate the marked difference in behavior between below wave breaking, in which the energy-transfer efficiency is high but total energy transfer is low, wave breaking, in which efficiency is low, and CWB, in which moderate efficiencies allow the highest total energy transfer.

Hard-spin mean-field theory has recently been applied to Ising magnets, correctly yielding the absence and presence of an interface roughening transition respectively in $d=2$ and $d=3$ dimensions and producing the ordering-roughening phase diagram for isotropic and anisotropic systems. The approach has now been extended to the effects of quenched random pinning centers and missing bonds on the interface of isotropic and anisotropic Ising models in $d=3$. We find that these frozen impurities cause domain boundary roughening that exhibits consecutive thresholding transitions as a function of interaction anisotropy. For both missing-bond and pinning-center impurities, for moderately large values of the anisotropy, the systems saturate to the “solid-on-solid” limit, exhibiting a single universal curve for the domain boundary width as a function of impurity concentration.

We generalize a theory of diffusion of a massive particle by the way in which transport characteristics are described by analytical expressions that formally coincide with those for the overdamped massless case but contain a factor comprising the particle mass which can be calculated in terms of Risken's matrix continued fraction method (MCFM). Using this generalization, we aim to elucidate how large gradients of a periodic potential affect the current in a tilted periodic potential and the average current of adiabatically driven on-off flashing ratchets. For this reason, we perform calculations for a sawtooth potential of the period $L$ with an arbitrary sawtooth length ($l<L$) instead of the smooth potentials typically considered in MCFM-solvable problems. We find nonanalytic behavior of the transport characteristics calculated for the sharp extremely asymmetric sawtooth potential at $l\to 0$ which appears due to the inertial effect. Analysis of the temperature dependences of the quantities under study reveals the dominant role of inertia in the high-temperature region. In particular, we show, by the analytical strong-inertia approach developed for this region, that the temperature-dependent contribution to the mobility at zero force and to the related effective diffusion coefficient are proportional to $T^{-3/2}$ and $T^{-1/2}$, respectively, and have a logarithmic singularity at $l\to 0$.

We study phase transitions and the nature of order in a class of classical generalized $\text{O}\left(N\right)$ nonlinear $\sigma$ models (NLS) constructed by minimally coupling pure NLS with additional degrees of freedom in the form of (i) Ising ferromagnetic spins, (ii) an advective Stokesian velocity, and (iii) multiplicative noises. In examples (i) and (ii), and also (iii) with the associated multiplicative noise being not sufficiently long-ranged, we show that the models may display a class of unusual phase transitions between stiff and soft phases, where the effective spin stiffness respectively diverges and vanishes in the long wavelength limit at two dimensions (2D), unlike in pure NLS. In the stiff phase, in the thermodynamic limit the variance of the transverse spin (or, the Goldstone mode) fluctuations are found to scale with the system size $L$ in 2D as $lnlnL$ with a model-dependent amplitude, which is markedly weaker than the well-known $lnL$ dependence of the variance of the broken symmetry modes in models that display quasi-long-range order in 2D. Equivalently, for $N=2$ at 2D the equal-time spin-spin correlations decay in powers of inverse logarithm of the spatial separation with model-dependent exponents. These transitions are controlled by the model parameters those couple the $\text{O}\left(N\right)$ spins with the additional variables. In the presence of long-range noises in example (iii), true long-range order may set in 2D, depending upon the specific details of the underlying dynamics. Our results should be useful in understanding phase transitions in equilibrium and nonequilibrium low-dimensional systems with continuous symmetries in general.

We calculate the isentropic evolution of an instantaneously heated foil, assuming a van der Waals equation of state with the Maxwell construction. The analysis by Yuen and Barnard [Phys. Rev. E 92, 033019 (2015)] is extended for the particular case of three degrees of freedom. We assume heating to temperatures in the vicinity of the critical point. The self-similar profiles of the rarefaction waves describing the evolution of the foil display plateaus in density and temperature due to a phase transition from the single-phase to the two-phase regime. The hydrodynamic equations are expressed in a dimensionless form and the solutions form a set of universal curves, depending on a single parameter: the dimensionless initial entropy. We characterize the rarefaction waves by calculating how the plateau length, density, pressure, temperature, velocity, internal energy, and sound speed vary with dimensionless initial entropy.

Polymer-dispersed liquid-crystal (PDLC) composites long have been a focus of study for their unique electro-optical properties which have resulted in various applications such as switchable (transparent or translucent) windows. These composites are manufactured using desirable “bottom-up” techniques, such as phase separation of a liquid-crystal–polymer mixture, which enable production of PDLC films at very large scales. LC domains within PDLCs are typically spheroidal, as opposed to rectangular for an LCD panel, and thus exhibit substantially different behavior in the presence of an external field. The fundamental difference between spheroidal and rectangular nematic domains is that the former results in the presence of nanoscale orientational defects in LC order while the latter does not. Progress in the development and optimization of PDLC electro-optical properties has progressed at a relatively slow pace due to this increased complexity. In this work, continuum simulations are performed in order to capture the complex formation and electric field-driven switching dynamics of approximations of PDLC domains. Using a simplified elliptic cylinder (microcapillary) geometry as an approximation of spheroidal PDLC domains, the effects of geometry (aspect ratio), surface anchoring, and external field strength are studied through the use of the Landau–de Gennes model of the nematic LC phase.

We investigate a number of complex patterns driven by the electroconvection instability in a planarly aligned layer of a nematic liquid crystal. They are traced back to various secondary instabilities of the ideal roll patterns bifurcating at onset of convection, whereby the basic nematohydrodynamic equations are solved by common Galerkin expansion methods. Alternatively these equations are systematically approximated by a set of coupled amplitude equations. They describe slow modulations of the convection roll amplitudes, which are coupled to a flow field component with finite vorticity perpendicular to the layer and to a quasihomogeneous in-plane rotation of the director. It is demonstrated that the Galerkin stability diagram of the convection rolls is well reproduced by the corresponding one based on the amplitude equations. The main purpose of the paper is, however, to demonstrate that their direct numerical simulations match surprisingly well new experiments, which serves as a convincing test of our theoretical approach.

We perform a Monte Carlo study of $N$-step self-avoiding walks, attached to the corner of an impenetrable wedge in two dimensions ($d=2$), or the tip of an impenetrable cone in $d=3$, of sizes ranging up to $N={10}^6$ steps. We find that the critical exponent ${\gamma }_{\alpha }$, which determines the dependence of the number of available conformations on $N$ for a cone or wedge with opening angle $\alpha$, is in good agreement with the theory for $d=2$. We study the end-point distribution of the walks in the allowed space and find similarities to the known behavior of random walks (ideal polymers) in the same geometry. For example, the ratio between the mean square end-to-end distances of a polymer near the cone or wedge and a polymer in free space depends linearly on ${\gamma }_{\alpha }$, as is known for ideal polymers. We show that the end-point distribution of polymers attached to a wedge does not separate into a product of angular and radial functions, as it does for ideal polymers in the same geometry. The angular dependence of the end position of polymers near the wedge differs from theoretical predictions.

We study the evolution of high-amplitude stress pulses in periodic dissipative laminates taking into account the nonlinear constitutive equations of the components and their dissipative behavior. Aluminum-tungsten laminate was selected due to the large difference in acoustic impedances of components, the significant nonlinearity of the aluminum constitutive equation at the investigated range of stresses, and its possible practical applications. Laminates with different cell size, which controls the internal time scale, impacted by plates with different thicknesses that determine the incoming pulse duration, were investigated. It has been observed that the ratio of the duration of the incoming pulse to the internal characteristic time determines the nature of the high-amplitude dissipative propagating waves—a triangular oscillatory shock–like profile, a train of localized pulses, or a single localized pulse. These localized quasistationary waves resemble solitary waves even in the presence of dissipation: The similar pulses emerged from different initial conditions, indicating that they are inherent properties of the corresponding laminates; their characteristic length scale is determined by the scale of mesostructure, nonlinear properties of materials, and the stress amplitude; and a linear relationship exists between their speed and amplitude. They mostly recover their shapes after collision with phase shift. A theoretical description approximating the shape, length scale, and speed of these high-amplitude dissipative pulses was proposed based on the Korteweg–de Vries equation with a dispersive term determined by the mesostructure and a nonlinear term derived using Hugoniot curves of components.

Based on first- and second-order perturbation theory, we present a numerical study of the temporal buildup and decay of unsteady acoustic fields and acoustic streaming flows actuated by vibrating walls in the transverse cross-sectional plane of a long straight microchannel under adiabatic conditions and assuming temperature-independent material parameters. The unsteady streaming flow is obtained by averaging the time-dependent velocity field over one oscillation period, and as time increases, it is shown to converge towards the well-known steady time-averaged solution calculated in the frequency domain. Scaling analysis reveals that the acoustic resonance builds up much faster than the acoustic streaming, implying that the radiation force may dominate over the drag force from streaming even for small particles. However, our numerical time-dependent analysis indicates that pulsed actuation does not reduce streaming significantly due to its slow decay. Our analysis also shows that for an acoustic resonance with a quality factor $Q$, the amplitude of the oscillating second-order velocity component is $Q$ times larger than the usual second-order steady time-averaged velocity component. Consequently, the well-known criterion $v_1^{}\ll c_s^{}$ for the validity of the perturbation expansion is replaced by the more restrictive criterion $v_1^{}\ll c_s^{}/Q$. Our numerical model is available as supplemental material in the form of comsol model files and matlab scripts.

A general understanding of optimal control in nonequilibrium systems would illuminate the operational principles of biological and artificial nanoscale machines. Recent work has shown that a system driven out of equilibrium by a linear response protocol is endowed with a Riemannian metric related to generalized susceptibilities, and that geodesics on this manifold are the nonequilibrium control protocols with the lowest achievable dissipation. While this elegant mathematical framework has inspired numerous studies of exactly solvable systems, no description of the thermodynamic geometry yet exists when the metric cannot be derived analytically. Herein, we numerically construct the dynamic metric of the two-dimensional Ising model in order to study optimal protocols for reversing the net magnetization.

A recurring motif in the organization of biological tissues are networks of long, fibrillar protein strands effectively confined to cylindrical surfaces. Often, the fibers in such curved, quasi-two-dimensional (2D) geometries adopt a characteristic order: the fibers wrap around the central axis at an angle which varies with radius and, in several cases, is strongly bimodally distributed. In this Rapid Communication, we investigate the general problem of a 2D crosslinked network of semiflexible fibers confined to a cylindrical substrate, and demonstrate that in such systems the trade-off between bending and stretching energies, very generically, gives rise to crosshatched order. We discuss its general dependency on the radius of the confining cylinder, and present an intuitive model that illustrates the basic physical principle of curvature-induced order. Our findings shed new light on the potential origin of some curiously universal fiber orientational distributions in tissue biology, and suggests novel ways in which synthetic polymeric soft materials may be instructed or programmed to exhibit preselected macromolecular ordering.

Theoretical results for the extension of a polymer confined to a channel are usually derived in the limit of infinite contour length. But experimental studies and simulations of DNA molecules confined to nanochannels are not necessarily in this asymptotic limit. We calculate the statistics of the span and the end-to-end distance of a semiflexible polymer of finite length in the extended de Gennes regime, exploiting the fact that the problem can be mapped to a one-dimensional weakly self-avoiding random walk. The results thus obtained compare favorably with pruned-enriched Rosenbluth method (PERM) simulations of a three-dimensional discrete wormlike chain model of DNA confined in a nanochannel. We discuss the implications for experimental studies of linear $\lambda$-DNA confined to nanochannels at the high ionic strengths used in many experiments.

Research has suggested that the activity occurring in a variety of social, economic, and technological systems exhibits long-range fluctuations in time. Pronounced levels of rapidly occurring events are typically observed over short periods of time, followed by long periods of inactivity. Relatively few studies, however, have shed light on the degree to which inhomogeneous temporal processes can be detected at, and emerge from, different levels of analysis. Here we investigate patterns of human activity within an online forum in which communication can be assessed at three intertwined levels: the micro level of the individual users; the meso level of discussion groups and continuous sessions; and the macro level of the whole system. To uncover the relation between different levels, we conduct a number of numerical simulations of a zero-crossing model in which users' behavior is constrained by progressively richer and more realistic rules of social interaction. Results indicate that, when users are solipsistic, their bursty behavior is not sufficient for generating heavy-tailed interevent time distributions at a higher level. However, when users are socially interdependent, the power spectra and interevent time distributions of the simulated and real forums are remarkably similar at all levels of analysis. Social interaction is responsible for the aggregation of multiple bursty activities at the micro level into an emergent bursty activity pattern at a higher level. We discuss the implications of the findings for an emergentist account of burstiness in complex systems.