Local changes in the topology of electricity grids can cause overloads far away from the disturbance [D. Witthaut and M. Timme, Eur. Phys. J. B 86, 377 (2013)], making the prediction of the robustness against changes in the topology—for example, caused by power outages or grid extensions—a challenging task. The impact of single-line additions on the long-range response of dc electricity grids has recently been studied [D. Labavić, R. Suciu, H. Meyer-Ortmanns, and S. Kettemann, Eur. Phys. J.: Spec. Top. 223, 2517 (2014)]. By solving the real part of the static ac load flow equations, we conduct a similar investigation for ac grids. In a regular two-dimensional grid graph with cyclic boundary conditions, we find a power law decay for the change of power flow as a function of distance to the disturbance over a wide range of distances. The power exponent increases and saturates for large system sizes. By applying the same analysis to the German transmission grid topology, we show that also in real-world topologies a long-ranged response can be found.

We propose a physical model for the capsids of tailed archaeal viruses as viscoelastic membranes under tension. The fluidity is generated by thermal motion of scarlike structures that are an intrinsic feature of the ground state of large particle arrays covering surfaces with nonzero Gauss curvature. The tension is generated by a combination of the osmotic pressure of the enclosed genome and an extension force generated by filamentous structure formation that drives the formation of the tails. In continuum theory, the capsid has the shape of a surface of constant mean curvature: an unduloid. Particle arrays covering unduloids are shown to exhibit pronounced subdiffusive and diffusive single-particle transport at temperatures that are well below the melting temperature of defect-free particle arrays on a surface with zero Gauss curvature.

During cell division, sister chromatids are segregated by the mitotic spindle, a bipolar assembly of interdigitating antiparallel polar filaments called microtubules. The spindle contains the midzone, a stable region of overlapping antiparallel microtubules, that is essential for maintaining bipolarity. Although a lot is known about the molecular players involved, the mechanism underlying midzone formation and maintenance is still poorly understood. We study the interaction of polar filaments that are cross-linked by molecular motors moving directionally and by passive cross-linkers diffusing along microtubules. Using a particle-based stochastic model, we find that the interplay of motors and passive cross-linkers can generate a stable finite overlap between a pair of antiparallel polar filaments. We develop a mean-field theory to study this mechanism in detail and investigate the influence of steric interactions between motors and passive cross-linkers on the overlap dynamics. In the presence of interspecies steric interactions, passive cross-linkers mimic the behavior of molecular motors and stable finite overlaps are generated even for non-cross-linking motors. Finally, we develop a mean-field theory for a bundle of aligned polar filaments and show that they can self-organize into a spindlelike pattern. Our work suggests possible ways as to how cells can generate spindle midzones and control their extensions.

Acoustic signals generated by filamentation of ultrashort terawatt laser pulses in water are characterized experimentally. Measurements reveal a strong influence of input pulse duration on the shape and intensity of the acoustic wave. Numerical simulations of the laser pulse nonlinear propagation and the subsequent water hydrodynamics and acoustic wave generation show that the strong acoustic emission is related to the mechanism of superfilamention in water. The elongated shape of the plasma volume where energy is deposited drives the far-field profile of the acoustic signal, which takes the form of a radially directed pressure wave with a single oscillation and a very broad spectrum.

Cell motility and tissue morphogenesis depend crucially on the dynamic remodeling of actomyosin networks. An actomyosin network consists of an actin polymer network connected by cross-linker proteins and motor protein myosins that generate internal stresses on the network. A recent discovery shows that for a range of experimental parameters, actomyosin networks contract to clusters with a power-law size distribution [J. Alvarado, Nat. Phys. 9, 591 (2013)]. Here, we argue that actomyosin networks can exhibit a robust critical signature without fine-tuning because the dynamics of the system can be mapped onto a modified version of percolation with trapping (PT), which is known to show critical behavior belonging to the static percolation universality class without the need for fine-tuning of a control parameter. We further employ our PT model to generate experimentally testable predictions.

Nonlinear plastic modes (NPMs) are collective displacements that are indicative of imminent plastic instabilities in elastic solids. In this work we formulate the atomistic theory that describes the reversible evolution of NPMs and their associated stiffnesses under external deformations. The deformation dynamics of NPMs is compared to those of the analogous observables derived from atomistic linear elastic theory, namely, destabilizing eigenmodes of the dynamical matrix and their associated eigenvalues. The key result we present and explain is that the dynamics of NPMs and of destabilizing eigenmodes under external deformations follow different scaling laws with respect to the proximity to imminent instabilities. In particular, destabilizing modes vary with a singular rate, whereas NPMs exhibit no such singularity. As a result, NPMs converge much earlier than destabilizing eigenmodes to their common final form at plastic instabilities. This dynamical difference between NPMs and linear destabilizing eigenmodes underlines the usefulness of NPMs for predicting the locus and geometry of plastic instabilities, compared to their linear-elastic counterparts.

We revisit the classical problem of flow of electrolyte solutions through charged capillary nanopores or nanotubes as described by the capillary pore model (also called “space charge” theory). This theory assumes very long and thin pores and uses a one-dimensional flux-force formalism which relates fluxes (electrical current, salt flux, and fluid velocity) and driving forces (difference in electric potential, salt concentration, and pressure). We analyze the general case with overlapping electric double layers in the pore and a nonzero axial salt concentration gradient. The $3\times 3$ matrix relating these quantities exhibits Onsager symmetry and we report a significant new simplification for the diagonal element relating axial salt flux to the gradient in chemical potential. We prove that Onsager symmetry is preserved under changes of variables, which we illustrate by transformation to a different flux-force matrix given by Gross and Osterle [J. Chem. Phys. 49, 228 (1968)]. The capillary pore model is well suited to describe the nonlinear response of charged membranes or nanofluidic devices for electrokinetic energy conversion and water desalination, as long as the transverse ion profiles remain in local quasiequilibrium. As an example, we evaluate electrical power production from a salt concentration difference by reverse electrodialysis, using an efficiency versus power diagram. We show that since the capillary pore model allows for axial gradients in salt concentration, partial loops in current, salt flux, or fluid flow can develop in the pore. Predictions for macroscopic transport properties using a reduced model, where the potential and concentration are assumed to be invariant with radial coordinate (“uniform potential” or “fine capillary pore” model), are close to results of the full model.

Depletion interaction potentials between micron-sized colloidal particles are induced by nanometer-scale surfactant micelles composed of hexaethylene glycol monododecyl ether $\left({\mathrm{C}}_{12}{\mathrm{E}}_6\right)$, and they are measured by video microscopy. The strength and range of the depletion interaction is revealed to arise from variations in shape anisotropy of the surfactant micelles. This shape anisotropy increases with increasing sample temperature. By fitting the colloidal interaction potentials to theoretical models, we extract micelle length and shape anisotropy as a function of temperature. This work introduces shape anisotropy tuning as a means to control interparticle interactions in colloidal suspensions, and it shows how the interparticle depletion potentials of micron-scale objects can be employed to probe the shape and size of surrounding macromolecules at the nanoscale.

Experiments are conducted in a two-dimensional monolayer vibrofluidized bed of glass beads, with a goal to understand the transition scenario and the underlying microstructure and dynamics in different patterned states. At small shaking accelerations ($\mathrm{\Gamma }=A{\omega }^2/g<1$, where $A$ and $\omega =2\pi f$ are the amplitude and angular frequency of shaking and $g$ is the gravitational acceleration), the particles remain attached to the base of the vibrating container; this is known as the solid bed (SB). With increasing $\mathrm{\Gamma }$ (at large enough shaking amplitude $A/d$) and/or with increasing $A/d$ (at large enough $\mathrm{\Gamma }$), the sequence of transitions/bifurcations unfolds as follows: SB (“solid bed”) to BB (“bouncing bed”) to LS (“Leidenfrost state”) to “2-roll convection” to “1-roll convection” and finally to a gas-like state. For a given length of the container, the coarsening of multiple convection rolls leading to the genesis of a “single-roll” structure (dubbed the multiroll transition) and its subsequent transition to a granular gas are two findings of this work. We show that the critical shaking intensity (${\mathrm{\Gamma }}_{\mathrm{BB}}^{\mathrm{LS}}$) for the $\text{BB}\to \text{LS}$ transition has a power-law dependence on the particle loading ($F=h_0/d$, where $h_0$ is the number of particle layers at rest and $d$ is the particle diameter) and the shaking amplitude ($A/d$). The characteristics of BB and LS states are studied by calculating (i) the coarse-grained density and temperature profiles and (ii) the pair correlation function. It is shown that while the contact network of particles in the BB state represents a hexagonal-packed structure, the contact network within the “floating cluster” of the LS resembles a liquid-like state. An unsteadiness of the Leidenfrost state has been uncovered wherein the interface (between the floating cluster and the dilute collisional layer underneath) and the top of the bed are found to oscillate sinusoidally, with the oscillation frequency closely matching the frequency of external shaking. Therefore, the granular Leidenfrost state is a period-1 wave as is the case for the BB state.

We generalize the telegrapher's equation to allow for anomalous transport. We derive the space-time fractional telegrapher's equation using the formalism of the persistent random walk in continuous time. We also obtain the characteristic function of the space-time fractional process and study some particular cases and asymptotic approximations. Similarly to the ordinary telegrapher's equation, the time-fractional equation also presents distinct behaviors for different time scales. Specifically, transitions between different subdiffusive regimes or from superdiffusion to subdiffusion are shown by the fractional equation as time progresses.

We present high-resolution measurements of the pair interactions between dielectric spheres dispersed in a fluid medium with a low dielectric constant. Despite the absence of charge control agents or added organic salts, these measurements reveal strong and long-ranged repulsions consistent with substantial charges on the particles whose interactions are screened by trace concentrations of mobile ions in solution. The dependence of the estimated charge on the particles' radii is consistent with charge renormalization theory and, thus, offers insights into the charging mechanism in this interesting class of model systems. The measurement technique, based on optical-tweezer manipulation and artifact-free particle tracking, makes use of optimal statistical methods to reduce measurement errors to the femtonewton frontier while covering an extremely wide range of interaction energies.

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for noncylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. We provide error estimates of our numerical solution schemes and systematically assess their convergence properties. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.

Self-assembly is ubiquitous in nature, particularly in biology, where it underlies the formation of protein quaternary structure and protein aggregation. Quaternary structure assembles deterministically and performs a wide range of important functions in the cell, whereas protein aggregation is the hallmark of a number of diseases and represents a nondeterministic self-assembly process. Here we build on previous work on a lattice model of deterministic self-assembly to investigate nondeterministic self-assembly of single lattice tiles and mixtures of two tiles at varying relative concentrations. Despite limiting the simplicity of the model to two interface types, which results in 13 topologically distinct single tiles and 106 topologically distinct sets of two tiles, we observe a wide variety of concentration-dependent behaviors. Several two-tile sets display critical behaviors in the form of a sharp transition from bound to unbound structures as the relative concentration of one tile to another increases. Other sets exhibit gradual monotonic changes in structural density, or nonmonotonic changes, while again others show no concentration dependence at all. We catalog this extensive range of behaviors and present a model that provides a reasonably good estimate of the critical concentrations for a subset of the critical transitions. In addition, we show that the structures resulting from these tile sets are fractal, with one of two different fractal dimensions.

We explore the spectra and localization properties of the $N$-site banded one-dimensional non-Hermitian random matrices that arise naturally in sparse neural networks. Approximately equal numbers of random excitatory and inhibitory connections lead to spatially localized eigenfunctions and an intricate eigenvalue spectrum in the complex plane that controls the spontaneous activity and induced response. A finite fraction of the eigenvalues condense onto the real or imaginary axes. For large $N$, the spectrum has remarkable symmetries not only with respect to reflections across the real and imaginary axes but also with respect to ${90}^{\circ }$ rotations, with an unusual anisotropic divergence in the localization length near the origin. When chains with periodic boundary conditions become directed, with a systematic directional bias superimposed on the randomness, a hole centered on the origin opens up in the density-of-states in the complex plane. All states are extended on the rim of this hole, while the localized eigenvalues outside the hole are unchanged. The bias-dependent shape of this hole tracks the bias-independent contours of constant localization length. We treat the large-$N$ limit by a combination of direct numerical diagonalization and using transfer matrices, an approach that allows us to exploit an electrostatic analogy connecting the “charges” embodied in the eigenvalue distribution with the contours of constant localization length. We show that similar results are obtained for more realistic neural networks that obey “Dale's law” (each site is purely excitatory or inhibitory) and conclude with perturbation theory results that describe the limit of large directional bias, when all states are extended. Related problems arise in random ecological networks and in chains of artificial cells with randomly coupled gene expression patterns.

We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.

Using a highly viscous magnetic fluid, the dynamics in the aftermath of the Rosensweig instability can be slowed down by more than 2000 times. In this way we expand the regime where the growth rate is predicted to scale linearly with the bifurcation parameter by six orders of magnitude, while this regime is tiny for standard ferrofluids and cannot be resolved experimentally there. We measure the growth of the pattern by means of a two-dimensional imaging technique, and find that the slopes of the growth and decay rates are not the same—a qualitative discrepancy with respect to the theoretical predictions. We solve this discrepancy by taking into account a viscosity which is assumed to be different for the growth and decay. This may be a consequence of the measured shear thinning of the ferrofluid.

Glasses display a wide array of nonlinear acoustic phenomena at temperatures $T\lesssim 1$ K. This behavior has traditionally been explained by an ensemble of weakly coupled, two-level tunneling states, a theory that is also used to describe the thermodynamic properties of glasses at low temperatures. One of the most striking acoustic signatures in this regime is the existence of phonon echoes, a feature that has been associated with two-level systems with the same formalism as spin echoes in NMR. Here we report the existence of a distinctly different type of acoustic echo in classical models of glassy materials. Our simulations consist of finite-ranged, repulsive spheres and also particles with attractive forces using Lennard-Jones interactions. We show that these echoes are due to anharmonic, weakly coupled vibrational modes and perhaps provide an alternative explanation for the phonon echoes observed in glasses at low temperatures.

Neurites, the long cellular protrusions that form the routes of the neuronal network, are capable of actively extending during early morphogenesis or regenerating after trauma. To perform this task, they rely on their cytoskeleton for mechanical support. In this paper, we present a three-component active gel model that describes neurites in the three robust mechanical states observed experimentally: collapsed, static, and motile. These states arise from an interplay between the physical forces driven by the growth of the microtubule-rich inner core of the neurite and the acto-myosin contractility of its surrounding cortical membrane. In particular, static states appear as a mechanical balance between traction and compression of these two parallel structures. The model predicts how the response of a neurite to a towing force depends on the force magnitude and recovers the response of neurites to several drug treatments that modulate the cytoskeleton active and passive properties.

Microfibers with their elongated shape and translation symmetry can act as important components in various soft materials, notably for their mechanics on the microscopic level. Here we demonstrate the mechanical response of a micro-object to imposed chirality, in this case, the tilt of disclination rings in an achiral nematic medium caused by the chiral surface anchoring on an immersed microfiber. This coupling between chirality and mechanical response, used to demonstrate sensing of chirality of electrospun cellulose microfibers, is revealed in the optical micrographs due to anisotropy in the elastic response of the host medium. We provide an analytical explanation of the chirogyral effect supported with numerical simulations and perform an experiment to test the effect of the cell confinement and fiber size. We controllably twist the microfibers and demonstrate the response of the nematic medium. More generally the demonstrated study provides means for experimental discrimination of surface properties and allows mechanical control over the shape of disclination rings.

We report experiments concerning the melting of ice disks (85 mm in diameter and 14 mm in height) at the surface of a thermalized water bath. During the melting, the ice disks undergo translational and rotational motions. In particular, the disks rotate. The rotation speed has been found to increase with the bath temperature. We investigated the flow under the bottom face of the ice disks by a particle image velocimetry technique. We find that the flow goes downwards and also rotates horizontally, so that a vertical vortex is generated under the ice disk. The proposed mechanism is the following. In the vicinity of the bottom face of the disk, the water eventually reaches the temperature of $4{}^{\circ }\mathrm{C}$ for which the water density is maximum. The $4{}^{\circ }\mathrm{C}$ water sinks and generates a downwards plume. The observed vertical vorticity results from the flow in the plume. Finally, by viscous entrainment, the horizontal rotation of the flow induces the solid rotation of the ice block. This mechanism seems generic: any vertical flow that generates a vortex will induce the rotation of a floating object.