In this work, the ratchet dynamics of Brownian particles driven by an external sinusoidal (harmonic) force is investigated. The gating ratchet effect is observed when another harmonic is used to modulate the spatially symmetric potential in which the particles move. For small amplitudes of the harmonics, it is shown that the current (average velocity) of particles exhibits a sinusoidal shape as a function of a precise combination of the phases of both harmonics. By increasing the amplitudes of the harmonics beyond the small-limit regime, departures from the sinusoidal behavior are observed and current reversals can also be induced. These current reversals persist even for the overdamped dynamics of the particles.

In a manner similar to the molecular chaos that underlies the stable thermodynamics of gases, a neuronal system may exhibit microscopic instability in individual neuronal dynamics while a macroscopic order of the entire population possibly remains stable. In this study, we analyze the microscopic stability of a network of neurons whose macroscopic activity obeys stable dynamics, expressing either monostable, bistable, or periodic state. We reveal that the network exhibits a variety of dynamical states for microscopic instability residing in a given stable macroscopic dynamics. The presence of a variety of dynamical states in such a simple random network implies more abundant microscopic fluctuations in real neural networks which consist of more complex and hierarchically structured interactions.

We investigate the effects of inherent fluctuations and system size in the dynamics of domain between uniform symmetric states. In the case of monotonous kinks, this dynamics is characterized by exhibiting nonsymmetric random walks, being attracted to the system borders. For nonmonotonous interface, the dynamics is replaced by a hopping dynamic. Based on bistable universal models, we characterize the origin of these unexpected dynamics through use of the stochastic kinematic laws for the interface position and the survival probability. Numerical simulations show a quite good agreement with the theoretical predictions.

In this paper a conjugate interface method is developed by performing extrapolations along the normal direction. Compared to other existing conjugate models, our method has several technical advantages, including the simple and straightforward algorithm, accurate representation of the interface geometry, applicability to any interface-lattice relative orientation, and availability of the normal gradient. The model is validated by simulating the steady and unsteady convection-diffusion system with a flat interface and the steady diffusion system with a circular interface, and good agreement is observed when comparing the lattice Boltzmann results with respective analytical solutions. A more general system with unsteady convection-diffusion process and a curved interface, i.e., the cooling process of a hot cylinder in a cold flow, is also simulated as an example to illustrate the practical usefulness of our model, and the effects of the cylinder heat capacity and thermal diffusivity on the cooling process are examined. Results show that the cylinder with a larger heat capacity can release more heat energy into the fluid and the cylinder temperature cools down slower, while the enhanced heat conduction inside the cylinder can facilitate the cooling process of the system. Although these findings appear obvious from physical principles, the confirming results demonstrates the application potential of our method in more complex systems. In addition, the basic idea and algorithm of the counter-extrapolation procedure presented here can be readily extended to other lattice Boltzmann models and even other computational technologies for heat and mass transfer systems.

Hierarchical organization is a cornerstone of complexity and multifractality constitutes its central quantifying concept. For model uniform cascades the corresponding singularity spectra are symmetric while those extracted from empirical data are often asymmetric. Using selected time series representing such diverse phenomena as price changes and intertransaction times in financial markets, sentence length variability in narrative texts, Missouri River discharge, and sunspot number variability as examples, we show that the resulting singularity spectra appear strongly asymmetric, more often left sided but in some cases also right sided. We present a unified view on the origin of such effects and indicate that they may be crucially informative for identifying the composition of the time series. One particularly intriguing case of this latter kind of asymmetry is detected in the daily reported sunspot number variability. This signals that either the commonly used famous Wolf formula distorts the real dynamics in expressing the largest sunspot numbers or, if not, that their dynamics is governed by a somewhat different mechanism.

In subdivided populations, migration acts together with selection and genetic drift and determines their evolution. Building upon a recently proposed method, which hinges on the emergence of a time scale separation between local and global dynamics, we study the fixation properties of subdivided populations in the presence of balancing selection. The approximation implied by the method is accurate when the effective selection strength is small and the number of subpopulations is large. In particular, it predicts a phase transition between species coexistence and biodiversity loss in the infinite-size limit and, in finite populations, a nonmonotonic dependence of the mean fixation time on the migration rate. In order to investigate the fixation properties of the subdivided population for stronger selection, we introduce an effective coarser description of the dynamics in terms of a voter model with intermediate states, which highlights the basic mechanisms driving the evolutionary process.

Burgers-Kardar-Parisi-Zhang (KPZ) scaling has recently (re-) surfaced in a variety of physical contexts, ranging from anharmonic chains to quantum systems such as open superfluids, in which a variety of random forces may be encountered and/or engineered. Motivated by these developments, we here provide a generalization of the KPZ universality class to situations with long-ranged temporal correlations in the noise, which purposefully break the Galilean invariance that is central to the conventional KPZ solution. We compute the phase diagram and critical exponents of the KPZ equation with $1/f$ noise $\left({\mathrm{KPZ}}_{1/f}\right)$ in spatial dimensions $1\le d<4$ using the dynamic renormalization group with a frequency cutoff technique in a one-loop truncation. Distinct features of ${\mathrm{KPZ}}_{1/f}$ are (i) a generically scale-invariant, rough phase at high noise levels that violates fluctuation-dissipation relations and exhibits hyperthermal statistics even in $d=1$, (ii) a fine-tuned roughening transition at which the flow fulfills an emergent thermal-like fluctuation-dissipation relation, that separates the rough phase from (iii) a massive phase in $1<d<4$ (in $d=1$ the interface is always rough). We point out potential connections to nonlinear hydrodynamics with a reduced set of conservation laws and noisy quantum liquids.

We consider an extension of classical stochastic reaction-diffusion (RD) dynamics to open quantum systems. We study a class of models of hard-core particles on a one-dimensional lattice whose dynamics is generated by a quantum master operator. Particle hopping is coherent while reactions, such as pair annihilation or pair coalescence, are dissipative. These are quantum open generalizations of the $A+A\to ⌀$ and $A+A\to A$ classical RD models. We characterize the relaxation of the state towards the stationary regime via a decomposition of the system Hilbert space into transient and recurrent subspaces. We provide a complete classification of the structure of the recurrent subspace (and the nonequilibrium steady states) in terms of the dark states associated to the quantum master operator and its general spectral properties. We also show that, in one dimension, relaxation towards these absorbing dark states is slower than that predicted by a mean-field analysis due to fluctuation effects, in analogy with what occurs in classical RD systems. Numerical simulations of small systems suggest that the decay of the density in one dimension, in both the open quantum $A+A\to ⌀$ and $A+A\to A$ systems, behaves asymptotically as $t^{-b}$ with $1/2<b<1$.

We develop a model of transport in multiscale porous media which accounts for adsorption in the different porosity scales. This model employs statistical mechanics to upscale molecular simulation and describe adsorption and transport at larger time and length scales. Using atom-scale simulations, which capture the changes in adsorption and transport with temperature, pressure, pore size, etc., this approach does not assume any adsorption or flow type. Moreover, by relating the local chemical potential $\mu \left(\mathbf{r}\right)$ and density $\rho \left(\mathbf{r}\right)$, the present model accounts for adsorption effects and possible changes in the confined fluid state upon transport. This model constitutes a bottom-up framework of adsorption and transport in multiscale materials as it (1) describes the adsorption-transport interplay, (2) accounts for the hydrodynamics breakdown at the nm scale, and (3) is multiscale.

A set of experiments was carried out on the motion of granular materials in a vertically vibrated annular system with a sawtooth-shaped base. We observed the coexistence of granular flow and surface patterns such as periodic subharmonic waves and kink pairs. Different patterns can transit each other as control parameters vary. The flow varies with space and time, and oppositely directed flows can occur at different levels and different moments. The magnitude and direction of the flow depend on the parameters defining the system in a complex manner. The motion of the patterns relates to the granular flow in a different way from that in which the wave in an ordinary fluid relates to the moving fluid. A preliminary explanation is given to our experimental findings.

Nuclear magnetic resonance (NMR) relaxation experimentation is an effective technique for probing the dynamics of proton spins in porous media, but interpretation requires the application of appropriate spin-diffusion models. Molecular dynamics (MD) simulations of porous silicate-based systems containing a quasi-two-dimensional water-filled pore are presented. The MD simulations suggest that the residency time of the water on the pore surface is in the range 0.03–12 ns, typically 2–5 orders of magnitude less than values determined from fits to experimental NMR measurements using the established surface-layer (SL) diffusion models of Korb and co-workers [Phys. Rev. E 56, 1934 (1997)]. Instead, MD identifies four distinct water layers in a tobermorite-based pore containing surface ${\mathrm{Ca}}^{2+}$ ions. Three highly structured water layers exist within 1 nm of the surface and the central region of the pore contains a homogeneous region of bulklike water. These regions are referred to as layer 1 and 2 (L1, L2), transition layer (TL), and bulk (B), respectively. Guided by the MD simulations, a two-layer (2L) spin-diffusion NMR relaxation model is proposed comprising two two-dimensional layers of slow- and fast-moving water associated with L2 and layers TL+B, respectively. The 2L model provides an improved fit to NMR relaxation times obtained from cementitious material compared to the SL model, yields diffusion correlation times in the range 18–75 ns and 28–40 ps in good agreement with MD, and resolves the surface residency time discrepancy. The 2L model, coupled with NMR relaxation experimentation, provides a simple yet powerful method of characterizing the dynamical properties of proton-bearing porous silicate-based systems such as porous glasses, cementitious materials, and oil-bearing rocks.

Using molecular dynamics simulations, the influence of the surface charge density of a nanotube on the static dielectric permittivity $\epsilon$ of confined water was reported. Whereas the dielectric anisotropy between the radial and axial directions of water confined in hydrophilic and hydrophobic membranes and the increase in axial dielectric permittivity with respect to the bulk value have previously been described, we found that an increase in the surface charge density leads to a drastic decrease in $\epsilon$ into the axial direction. The decrease in $\epsilon$ is accompanied by a strong slowdown in the rotational dynamics of water molecules. We show that this effect is due to the strong orientation of water molecules induced by the surface charge. Thus, by controlling the surface charge in nanotubes and nanocavities, it is possible to tune the dielectric permittivity of confined fluids at the nanoscale.

The effect of a local shear transformation on plastic deformation of a three-dimensional amorphous solid is studied using molecular dynamics simulations. We consider a spherical inclusion, which is gradually transformed into an ellipsoid of the same volume and then converted back into the sphere. It is shown that at sufficiently large strain amplitudes, the deformation of the material involves localized plastic events that are identified based on the relative displacement of atoms before and after the shear transformation. We find that the density profiles of cage jumps decay away from the inclusion, which correlates well with the radial dependence of the local deformation of the material. At the same strain amplitude, the plastic deformation becomes more pronounced in the cases of weakly damped dynamics or large time scales of the shear transformation. We show that the density profiles can be characterized by the universal function of the radial distance multiplied by a dimensionless factor that depends on the friction coefficient and the time scale of the shear event.

We investigate properties of evolving linguistic networks defined by the word-adjacency relation. Such networks belong to the category of networks with accelerated growth but their shortest-path length appears to reveal the network size dependence of different functional form than the ones known so far. We thus compare the networks created from literary texts with their artificial substitutes based on different variants of the Dorogovtsev-Mendes model and observe that none of them is able to properly simulate the novel asymptotics of the shortest-path length. Then, we identify the local chainlike linear growth induced by grammar and style as a missing element in this model and extend it by incorporating such effects. It is in this way that a satisfactory agreement with the empirical result is obtained.

We introduce two-dimensional (2D) linear and nonlinear Talbot effects. They are produced by propagating periodic 2D diffraction patterns and can be visualized as 3D stacks of Talbot carpets. The nonlinear Talbot effect originates from 2D rogue waves and forms in a bulk 3D nonlinear medium. The recurrences of an input rogue wave are observed at the Talbot length and at the half-Talbot length, with a $\pi$ phase shift; no other recurrences are observed. Differing from the nonlinear Talbot effect, the linear effect displays the usual fractional Talbot images as well. We also find that the smaller the period of incident rogue waves, the shorter the Talbot length. Increasing the beam intensity increases the Talbot length, but above a threshold this leads to a catastrophic self-focusing phenomenon which destroys the effect. We also find that the Talbot recurrence can be viewed as a self-Fourier transform of the initial periodic beam that is automatically performed during propagation. In particular, linear Talbot effect can be viewed as a fractional self-Fourier transform, whereas the nonlinear Talbot effect can be viewed as the regular self-Fourier transform. Numerical simulations demonstrate that the rogue-wave initial condition is sufficient but not necessary for the observation of the effect. It may also be observed from other periodic inputs, provided they are set on a finite background. The 2D effect may find utility in the production of 3D photonic crystals.

When repelling particles are confined by a transverse potential in quasi-one-dimensional geometry, the straight line equilibrium configuration becomes unstable at small confinement, in favor of a staggered row that may be inhomogeneous or homogeneous. This conformational phase transition is a pitchfork bifurcation called the zigzag transition. We study the zigzag transition in infinite and periodic finite systems with short-range interactions. We provide numerical evidence that in this case the bifurcation is subcritical since it exhibits phase coexistence and hysteretic behavior. The physical mechanism responsible for the change in the bifurcation character is the nonlinear coupling between the transverse soft mode at the transition and the longitudinal Goldstone mode linked to the translational or rotational invariance of the zigzag pattern. An asymptotic analysis, near the bifurcation threshold and assuming an infinite system, gives an explicit expression for the normal form of the bifurcation. We establish the subcriticality, and we describe with excellent precision the inhomogeneous zigzag patterns observed in the simulations. A direct test of the physical mechanism responsible for the bifurcation character evidences a quantitative agreement.

Chimera states consisting of domains of coherently and incoherently oscillating nonlocally coupled phase oscillators in systems with spatial inhomogeneity are studied. The inhomogeneity is introduced through the dependence of the oscillator frequency on its location. Two types of spatial inhomogeneity, localized and spatially periodic, are considered and their effects on the existence and properties of multicluster and traveling chimera states are explored. The inhomogeneity is found to break up splay states, to pin the chimera states to specific locations, and to trap traveling chimeras. Many of these states can be studied by constructing an evolution equation for a complex order parameter. Solutions of this equation are in good agreement with the results of numerical simulations.

The motility of swimming bacteria near solid surfaces has implications in a wide range of scenarios, including water treatment facilities, microfluidics, and biomedical implants. Using the boundary element method to numerically solve the equations of low Reynolds number fluid flow, we investigate the dynamics of a model swimmer propelled by rotating a single helical flagellum. Building on previous simulation results for swimmers near a single plane boundary, we introduce a second, parallel boundary and show that the bacterial trajectories change as the two plates are brought closer together. Analysis of this dynamical system shows that the configuration in the center of the channel and parallel to the walls is an unstable equilibrium state for large plate separations, but it becomes the only stable position for swimmers when the plate separation is reduced to three to four times the cell width. Our model also predicts that transient trajectories, i.e., those not at steady states, can exhibit curvature in the opposite sense to that expected from the well-known explanation for circular bacterial paths near a single wall.

The conventional fluid description of multicomponent plasma, supplemented by an appropriate equation of state for the macroparticle component, is used to evaluate the longitudinal sound velocity of Yukawa fluids. The obtained results are in very good agreement with those obtained earlier employing the quasilocalized charge approximation and molecular dynamics simulations in a rather broad parameter regime. Thus, a simple yet accurate tool to estimate the sound velocity across coupling regimes is proposed, which can be particularly helpful in estimating the dust-acoustic velocity in strongly coupled dusty (complex) plasmas. It is shown that, within the present approach, the sound velocity is completely determined by particle-particle correlations and the neutralizing medium (plasma), apart from providing screening of the Coulomb interaction, has no other effect on the sound propagation. The ratio of the actual sound velocity to its “ideal gas” (weak coupling) scale only weakly depends on the coupling strength in the fluid regime but exhibits a pronounced decrease with the increase of the screening strength. The limitations of the present approach in applications to real complex plasmas are briefly discussed.