Cilia and flagella are hairlike organelles that propel cells through fluid. The active motion of the axoneme, the motile structure inside cilia and flagella, is powered by molecular motors of the axonemal dynein family. These motors generate forces and torques that slide and bend the microtubule doublets within the axoneme. To create regular waveforms, the activities of the dyneins must be coordinated. It is thought that coordination is mediated by stresses due to radial, transverse, or sliding deformations, and which build up within the moving axoneme and feed back on dynein activity. However, which particular components of the stress regulate the motors to produce the observed waveforms of the many different types of flagella remains an open question. To address this question, we describe the axoneme as a three-dimensional bundle of filaments and characterize its mechanics. We show that regulation of the motors by radial and transverse stresses can lead to a coordinated flagellar motion only in the presence of twist. We show that twist, which could arise from torque produced by the dyneins, couples curvature to transverse and radial stresses. We calculate emergent beating patterns in twisted axonemes resulting from regulation by transverse stresses. The resulting waveforms are similar to those observed in flagella of Chlamydomonas and sperm. Due to the twist, the waveform has nonplanar components, which result in swimming trajectories such as twisted ribbons and helices, which agree with observations.

We study elastic interactions between colloidal particles and deformation landscapes of undulations in a cholesteric liquid crystal under an electric field applied normal to cholesteric layers. The onset of undulation instability is influenced by the presence of colloidal inclusions and, in turn, layers’ undulations mediate the spatial patterning of particle locations. We find that the bending of cholesteric layers around a colloidal particle surface prompts the local nucleation of an undulations lattice at electric fields below the well-defined threshold known for liquid crystals without inclusions, and that the onset of the resulting lattice is locally influenced, both dimensionally and orientationally, by the initial arrangements of colloids defined using laser tweezers. Spherical particles tend to spatially localize in the regions of strong distortions of the cholesteric layers, while colloidal nanowires exhibit an additional preference for multistable alignment offset along various vectors of the undulations lattice. Magnetic rotation of superparamagnetic colloidal particles couples with the locally distorted helical axis and undulating cholesteric layers in a manner that allows for a controlled three-dimensional translation of these particles. These interaction modes lend insight into the physics of liquid crystal structure-colloid elastic interactions, as well as point the way towards guided self-assembly of reconfigurable colloidal composites with potential applications in diffraction optics and photonics.

Brownian particles interacting via repulsive soft-core potentials can spontaneously aggregate, despite repelling each other, and form periodic crystals of particle clusters. We study this phenomenon in low-dimensional situations (one and two dimensions) at two levels of description: by performing numerical simulations of the discrete particle dynamics and by linear and nonlinear analysis of the corresponding Dean-Kawasaki equation for the macroscopic particle density. Restricting to low dimensions and neglecting fluctuation effects, we gain analytical insight into the mechanisms of the instability leading to clustering which turn out to be the interplay among diffusion, the intracluster forces, and the forces between neighboring clusters. We show that the deterministic part of the Dean-Kawasaki equation provides a good description of the particle dynamics, including width and shape of the clusters and over a wide range of parameters, and analyze with weakly nonlinear techniques the nature of the pattern-forming bifurcation in one and two dimensions. Finally, we briefly discuss the case of attractive forces.

We propose a new topological characterization of ribonucleic acid (RNA) secondary structures with pseudoknots based on two topological invariants. Starting from the classic arc representation of RNA secondary structures, we consider a model that couples both (i) the topological genus of the graph and (ii) the number of crossing arcs of the corresponding primitive graph. We add a term proportional to these topological invariants to the standard free energy of the RNA molecule, thus obtaining a novel free-energy parametrization that takes into account the abundance of topologies of RNA pseudoknots observed in RNA databases.

We study the dynamics of a probe particle driven by a constant force through a colloidal glass of hard spheres. This nonequilibrium and anisotropic problem is investigated using a new implementation of the mode-coupling approximation with multiple relaxation channels and Langevin dynamics simulations. A force threshold is found, below which the probe remains localized, while above it the probe acquires a finite velocity. We focus on the localized regime, comparing theory and simulations concerning the dynamics in the length scale of the cage and the properties of the transition to the delocalized regime, such as the critical power-law decay of the probe correlation function. Probe van Hove functions predicted by the theory show exponential tails reminiscent of an intermittent dynamics of the probe. This scenario is microscopically supported by simulations.

Recent development of single-molecule manipulation technologies has made it possible to exert constant force and torque on individual DNA biopolymers to probe their elastic characteristics and structural stability. It has been previously shown that depending on the nature of applied mechanical constraints, DNA can exist in several forms including B-, L-, and P-DNA. However, there is still a lack of understanding of how structural heterogeneity of DNA, which may naturally arise due to sequence-dependent DNA properties, protein binding, or DNA damage, influences local stability of the above DNA states. To provide a more complete and detailed description of the DNA mechanics, we developed a theoretical framework based on transfer-matrix calculations and demonstrated how it can be used to predict the DNA behavior upon application of a wide range of force and torque constraints. The resulting phase diagram shows DNA structural transitions that are in good agreement with previous experimental and theoretical studies. We further discuss how the constructed formalism can be extended to include local inhomogeneities in the DNA physical properties, thus making it possible to investigate the effect of DNA sequence as well as protein binding on DNA structural stability.

Based on the explicit knowledge of a Hamiltonian of mean force, the classical statistical mechanics and equilibrium thermodynamics of open systems in contact with a thermal environment at arbitrary interaction strength can be formulated. Yet, even though the Hamiltonian of mean force uniquely determines the equilibrium phase space probability density of a strongly coupled open system, the knowledge of this probability density alone is insufficient to determine the Hamiltonian of mean force, needed in constructing the underlying statistical mechanics and thermodynamics. We demonstrate that under the assumption that the Hamiltonian of mean force is known, an extension of thermodynamic structures from the level of averaged quantities to fluctuating objects (i.e., a stochastic thermodynamics) is possible. However, such a construction undesirably also involves a vast ambiguity. This situation is rooted in the eminent lack of a physical guiding principle allowing us to distinguish a physically meaningful theory out of a multitude of other equally conceivable ones.

The dynamics of two weakly coupled parametric oscillators are studied in the neighborhood of the primary subharmonic instability. The nature of both primary and secondary instabilities depends in a critical way on the permutation symmetries, if any, that remain after coupling is considered, and this depends on the relative phases of the parametric forcing terms. Detailed bifurcation sets, revealing a complex series of transitions organized in part by Bogdanov-Takens points, are calculated for representative sets of parameters. In the particular case of out-of-phase forcing the predictions of the coupled oscillator model are compared with direct numerical simulations and with recent experiments on modulated cross waves. Both the initial Hopf bifurcation and the subsequent saddle-node heteroclinic bifurcation are confirmed.

We consider a homogeneous population of stochastic neurons that are driven by weak common noise (stimulus). To capture and analyze the joint firing events within the population, we introduce the partial synchronous output of the population. This is a time series defined by the events that at least a fixed fraction $\gamma \in [0,1]$ of the population fires simultaneously within a small time interval. For this partial synchronous output we develop two analytical approaches to the correlation statistics. In the Gaussian approach we represent the synchronous output as a nonlinear transformation of the summed population activity and approximate the latter by a Gaussian process. In the combinatorial approach the synchronous output is represented by products of box-filtered spike trains of the single neurons. In both approaches we use linear-response theory to derive approximations for statistical measures that hold true for weak common noise. In particular, we calculate the mean value and power spectrum of the synchronous output and the cross-spectrum between synchronous output and common noise. We apply our results to the leaky integrate-and-fire neuron model and compare them to numerical simulations. The combinatorial approach is shown to provide a more accurate description of the statistics for small populations, whereas the Gaussian approximation yields compact formulas that work well for a sufficiently large population size. In particular, in the Gaussian approximation all statistical measures reveal a symmetry in the synchrony threshold $\gamma$ around the mean value of the population activity. Our results may contribute to a better understanding of the role of coincidence detection in neural signal processing.

In many important systems exhibiting crackling noise—an intermittent avalanchelike relaxation response with power-law and, thus, self-similar distributed event sizes—the “laws” for the rate of activity after large events are not consistent with the overall self-similar behavior expected on theoretical grounds. This is particularly true for the case of seismicity, and a satisfying solution to this paradox has remained outstanding. Here, we propose a generalized description of the aftershock rates which is both self-similar and consistent with all other known self-similar features. Comparing our theoretical predictions with high-resolution earthquake data from Southern California we find excellent agreement, providing particularly clear evidence for a unified description of aftershocks and foreshocks. This may offer an improved framework for time-dependent seismic hazard assessment and earthquake forecasting.

Glow discharge polymer hydrocarbon (GDP-CH) is used as the ablator material in inertial confinement fusion (ICF) capsules for the Laser Mégajoule and National Ignition Facility. Due to its fabrication process, GDP-CH chemical composition and structure differ from commercially available plastics and detailed knowledge of its properties in the warm dense matter regime is needed to achieve accurate design of ICF capsules. First-principles ab initio simulations of the GDP-CH principal Hugoniot up to 8 Mbar were performed using the quantum molecular dynamics (QMD) code abinit and showed that atomic bond dissociation has an effect on the compressibility. Results from these simulations are used to parametrize a quantum semiempirical model in order to generate a tabulated equation of state that includes dissociation. Hugoniot measurements obtained from an experiment conducted at the LULI2000 laser facility confirm QMD simulations as well as EOS modeling. We conclude by showing the EOS model influence on shock timing in a hydrodynamic simulation.

We present a simulation and data analysis technique to investigate first-order phase transitions and the associated transition barriers. The simulation technique is based on the real microcanonical ensemble where the sum of kinetic and potential energy is kept constant. The method is tested for the droplet condensation-evaporation transition in a Lennard-Jones system with up to 2048 particles at fixed density, using simple Metropolis-like sampling combined with a replica-exchange scheme. Our investigation of the microcanonical ensemble properties reveals that the associated transition barrier is significantly lower than in the canonical counterpart. Along the line of investigating the microcanonical ensemble behavior, we develop a framework for general ensemble evaluations. This framework is based on a clear separation between system-related and ensemble-related properties, which can be exploited to specifically tailor artificial ensembles suitable for first-order phase transitions.

Many single-stranded (ss) ribonucleic acid (RNA) viruses self-assemble from capsid protein subunits and the nucleic acid to form an infectious virion. It is believed that the electrostatic interactions between the negatively charged RNA and the positively charged viral capsid proteins drive the encapsidation, although there is growing evidence that the sequence of the viral RNA also plays a role in packaging. In particular, the sequence will determine the possible secondary structures that the ssRNA will take in solution. In this work, we use a mean-field theory to investigate how the secondary structure of the RNA combined with electrostatic interactions affects the efficiency of assembly and stability of the assembled virions. We show that the secondary structure of RNA may result in negative osmotic pressures while a linear polymer causes positive osmotic pressures for the same conditions. This may suggest that the branched structure makes the RNA more effectively packaged and the virion more stable.

We derive the hydrodynamic equations for nematic liquid crystals lying on curved substrates. We invoke the Lagrange-Rayleigh variational principle to adapt the Ericksen-Leslie theory to two-dimensional nematics in which a degenerate anchoring of the molecules on the substrate is enforced. The only constitutive assumptions in this scheme concern the free-energy density, given by the two-dimensional Frank potential, and the density of dissipation which is required to satisfy appropriate invariance requirements. The resulting equations of motion couple the velocity field, the director alignment, and the curvature of the shell. To illustrate our findings, we consider the effect of a simple shear flow on the alignment of a nematic lying on a cylindrical shell.

A recent experiment by Kwak et al. [Sci. Rep. 5, 9010 (2015)] demonstrated the relevance of resonance-assisted tunneling for optical microcavities where resonance chains emerge in phase space due to boundary deformations. In this paper we adapt the perturbative description of resonance-assisted tunneling to calculate optical modes and the imaginary part of their complex wavenumber which determines the lifetime of the mode. We demonstrate our method at three example cavity shapes and compare our results to numerical data and perturbation theory for weakly deformed microdisk cavities.

The ratchet effect of a sine-Gordon kink is investigated in the absence of any external force while the symmetry of the field potential at every time instant is maintained. The directed motion appears by a time shift of the sine-Gordon potential through a time-dependent additional phase. A symmetry analysis provides the necessary conditions for the existence of net motion. It is also shown analytically, by using a collective coordinate theory, that the novel physical mechanism responsible for the appearance of the ratchet effect is the coupled dynamics of the kink width with the background field. Biharmonic and dichotomic periodic variations of the additional phase of the sine-Gordon potential are considered. The predictions established by the symmetry analysis and the collective coordinate theory are verified by means of numerical simulations. Inversion and maximization of the resulting current as a function of the system parameters are investigated.

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.