Nonlinear dynamics and bifurcations of a planar undulating magnetic microswimmer.

Swimming microorganisms such as flagellated bacteria and sperm cells have fascinating locomotion capabilities. Inspired by their natural motion, there is an ongoing effort to develop artificial robotic nanoswimmers for potential in-body biomedical applications. A leading method for actuation of nanoswimmers is by applying a time-varying external magnetic field. Such systems have rich and nonlinear dynamics that call for simple fundamental models. A previous work studied forward motion of a simple two-link model with a passive elastic joint, assuming small-amplitude planar oscillations of the magnetic field about a constant direction. In this work, we found that there exists a faster, backward motion of the swimmer with very rich dynamics. By relaxing the small-amplitude assumption, we analyze the multiplicity of periodic solutions, as well as their bifurcations, symmetry breaking, and stability transitions. We have also found that the net displacement and/or mean swimming speed are maximized for optimal choices of various parameters. Asymptotic calculations are performed for the bifurcation condition and the swimmer's mean speed. The results may enable significantly improving the design aspects of magnetically actuated robotic microswimmers.

From Chaos to Ordering: New Studies in the Shannon Entropy of 2D Patterns.

Properties of the Voronoi tessellations arising from random 2D distribution points are reported. We applied an iterative procedure to the Voronoi diagrams generated by a set of points randomly placed on the plane. The procedure implied dividing the edges of Voronoi cells into equal or random parts. The dividing points were then used to construct the following Voronoi diagram. Repeating this procedure led to a surprising effect of the positional ordering of Voronoi cells, reminiscent of the formation of lamellae and spherulites in linear semi-crystalline polymers and metallic glasses. Thus, we can conclude that by applying even a simple set of rules to a random set of seeds, we can introduce order into an initially disordered system. At the same time, the Shannon (Voronoi) entropy showed a tendency to attain values that are typical for completely random patterns; thus, the Shannon (Voronoi) entropy does not distinguish the short-range ordering. The Shannon entropy and the continuous measure of symmetry of the patterns demonstrated the distinct asymptotic behavior, while approaching the close saturation values with the increase in the number of iteration steps. The Shannon entropy grew with the number of iterations, whereas the continuous measure of symmetry of the same patterns demonstrated the opposite asymptotic behavior. The Shannon (Voronoi) entropy is not an unambiguous measure of order in the 2D patterns. The more symmetrical patterns may demonstrate the higher values of the Shannon entropy.

Breather arrest, localization, and acoustic non-reciprocity in dissipative nonlinear lattices.

The effect of on-site damping on breather arrest, localization, and non-reciprocity in strongly nonlinear lattices is analytically and numerically studied. Breathers are localized oscillatory wavepackets formed by nonlinearity and dispersion. Breather arrest refers to breather disintegration over a finite "penetration depth" in a dissipative lattice. First, a simplified system of two nonlinearly coupled oscillators under impulsive excitation is considered. The exact relation between the number of beats (energy exchanges between oscillators), the excitation magnitude, and the on-site damping is derived. Then, these analytical results are correlated to those of the semi-infinite extension of the simplified system, where breather penetration depth is governed by a similar law to that of the finite beats in the simplified system. Finally, motivated by the experimental results of Bunyan, Moore, Mojahed, Fronk, Leamy, Tawfick, and Vakakis [Phys. Rev. E 97, 052211 (2018)], breather arrest, localization, and acoustic non-reciprocity in a non-symmetric, dissipative, strongly nonlinear lattice are studied. The lattice consists of repetitive cells of linearly grounded large-scale particles nonlinearly coupled to small-scale ones, and linear intra-cell coupling. Non-reciprocity in this lattice yields either energy localization or breather arrest depending on the position of excitation. The nonlinear acoustics governing non-reciprocity, and the surprising effects of existence of linear components in the coupling nonlinear stiffnesses, in the acoustics, are investigated.

Nonreciprocity in the dynamics of coupled oscillators with nonlinearity, asymmetry, and scale hierarchy.

In linear time-invariant dynamical and acoustical systems, reciprocity holds by the Onsager-Casimir principle of microscopic reversibility, and this can be broken only by odd external biases, nonlinearities, or time-dependent properties. A concept is proposed in this work for breaking dynamic reciprocity based on irreversible nonlinear energy transfers from large to small scales in a system with nonlinear hierarchical internal structure, asymmetry, and intentional strong stiffness nonlinearity. The resulting nonreciprocal large-to-small scale energy transfers mimic analogous nonlinear energy transfer cascades that occur in nature (e.g., in turbulent flows), and are caused by the strong frequency-energy dependence of the essentially nonlinear small-scale components of the system considered. The theoretical part of this work is mainly based on action-angle transformations, followed by direct numerical simulations of the resulting system of nonlinear coupled oscillators. The experimental part considers a system with two scales-a linear large-scale oscillator coupled to a small scale by a nonlinear spring-and validates the theoretical findings demonstrating nonreciprocal large-to-small scale energy transfer. The proposed study promotes a paradigm for designing nonreciprocal acoustic materials harnessing strong nonlinearity, which in a future application will be implemented in designing lattices incorporating nonlinear hierarchical internal structures, asymmetry, and scale mixing.

Heat conduction in a chain of colliding particles with a stiff repulsive potential.

One-dimensional billiards, i.e., a chain of colliding particles with equal masses, is a well-known example of a completely integrable system. Billiards with different particle masses is generically not integrable, but it still exhibits divergence of a heat conduction coefficient (HCC) in the thermodynamic limit. Traditional billiards models imply instantaneous (zero-time) collisions between the particles. We relax this condition of instantaneous impact and consider heat transport in a chain of stiff colliding particles with the power-law potential of the nearest-neighbor interaction. The instantaneous collisions correspond to the limit of infinite power in the interaction potential; for finite powers, the interactions take nonzero time. This modification of the model leads to a profound physical consequence-the probability of multiple (in particular triple) -particle collisions becomes nonzero. Contrary to the integrable billiards of equal particles, the modified model exhibits saturation of the heat conduction coefficient for a large system size. Moreover, the identification of scattering events with triple-particle collisions leads to a simple definition of the characteristic mean free path and a kinetic description of heat transport. This approach allows us to predict both the temperature and density dependencies for the HCC limit values. The latter dependence is quite counterintuitive-the HCC is inversely proportional to the particle density in the chain. Both predictions are confirmed by direct numerical simulations.

Localization in finite vibroimpact chains: Discrete breathers and multibreathers.

We explore the dynamics of strongly localized periodic solutions (discrete solitons or discrete breathers) in a finite one-dimensional chain of oscillators. Localization patterns with both single and multiple localization sites (breathers and multibreathers) are considered. The model involves parabolic on-site potential with rigid constraints (the displacement domain of each particle is finite) and a linear nearest-neighbor coupling. When the particle approaches the constraint, it undergoes an inelastic impact according to Newton's impact model. The rigid nonideal impact constraints are the only source of nonlinearity and damping in the system. We demonstrate that this vibro-impact model allows derivation of exact analytic solutions for the breathers and multibreathers with an arbitrary set of localization sites, both in conservative and in forced-damped settings. Periodic boundary conditions are considered; exact solutions for other types of boundary conditions are also available. Local character of the nonlinearity permits explicit derivation of a monodromy matrix for the breather solutions. Consequently, the stability of the derived breather and multibreather solutions can be efficiently studied in the framework of simple methods of linear algebra, and with rather moderate computational efforts. One reveals that that the finiteness of the chain fragment and possible proximity of the localization sites strongly affect both the existence and the stability patterns of these localized solutions.

Nonstationary heat conduction in one-dimensional chains with conserved momentum.

This Rapid Communication addresses the relationship between hyperbolic equations of heat conduction and microscopic models of dielectrics. Effects of the nonstationary heat conduction are investigated in two one-dimensional models with conserved momentum: Fermi-Pasta-Ulam (FPU) chain and chain of rotators (CR). These models belong to different universality classes with respect to stationary heat conduction. Direct numeric simulations reveal in both models a crossover from oscillatory decay of short-wave perturbations of the temperature field to smooth diffusive decay of the long-wave perturbations. Such behavior is inconsistent with parabolic Fourier equation of the heat conduction. The crossover wavelength decreases with increase in average temperature in both models. For the FPU model the lowest-order hyperbolic Cattaneo-Vernotte equation for the nonstationary heat conduction is not applicable, since no unique relaxation time can be determined.