We study a generalization of the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model of a particle in a Brownian force landscape, including retardation effects. We show that under monotonous driving the particle moves forward at all times, as it does in absence of retardation (Middleton's theorem). This remarkable property allows us to develop an analytical treatment. The model with an exponentially decaying memory kernel is realized in Barkhausen experiments with eddy-current relaxation and has previously been shown numerically to account for the experimentally observed asymmetry of Barkhausen pulse shapes. We elucidate another qualitatively new feature: the breakup of each avalanche of the standard ABBM model into a cluster of subavalanches, sharply delimited for slow relaxation under quasistatic driving. These conditions are typical for earthquake dynamics. With relaxation and aftershock clustering, the present model includes important ingredients for an effective description of earthquakes. We analyze quantitatively the limits of slow and fast relaxation for stationary driving with velocity $v>0$. The $v$-dependent power-law exponent for small velocities, and the critical driving velocity at which the particle velocity never vanishes, are modified. We also analyze nonstationary avalanches following a step in the driving magnetic field. Analytically, we obtain the mean avalanche shape at fixed size, the duration distribution of the first subavalanche, and the time dependence of the mean velocity. We propose to study these observables in experiments, allowing a direct measurement of the shape of the memory kernel and tracing eddy current relaxation in Barkhausen noise.

We study the dynamics of a tagged particle in an environment of point Brownian particles with hard-core interactions in an infinite one-dimensional channel (a single-file model). In particular, we examine the influence of initial conditions on the dynamics of the tagged particle. We compare two initial conditions: equal distances between particles and uniform density distribution. The effect is shown by the differences of mean-square-displacement and correlation function for the two ensembles of initial conditions. We discuss the violation of Einstein relation, and its dependence on the initial condition, and the difference between time and ensemble averaging. More specifically, using the Jepsen line, we will discuss how transport coefficients, like diffusivity, depend on the initial state. Our work shows that initial conditions determine the long time limit of the dynamics, and in this sense the system never forgets its initial state in complete contrast with thermal systems (i.e., a closed system that attains equilibrium independent of the initial state).

We examine a density-independent modification of the Vicsek model in which a particle interacts with neighbors defined by Delaunay triangulation. To feasibly simulate the model, an algorithm for repairing the triangulation over time was developed. This algorithm may also be applied to any time varying two-dimensional Delaunay triangulation. This model exhibits a continuous phase transition with noise, and a distinct set of critical exponents were measured which satisfy a hyperscaling relationship. The critical exponents are found to vary between a low and high velocity regime, but they are robust under the inclusion of a repulsive interaction. We present evidence that the correlation length approximately scales with the size of the system even in the ordered phase.

In our previous work [J. Chem. Phys. 136, 024502 (2012)], we reported a demixing phase transition of a quasi-two-dimensional, binary Heisenberg fluid mixture driven by the ferromagnetic interactions of the magnetic species. Here, we present a theoretical study for the time-dependent coarsening occurring within the two-phase region in the density-concentration plane, also known as spinodal decomposition. Our investigations are based on dynamical density functional theory (DDFT). The particles in the mixture are modeled as Gaussian soft spheres on a two-dimensional surface, where one component carries a classical spin of Heisenberg type. To investigate the two-phase region, we first present a linear stability analysis with respect to small, harmonic density perturbations. Second, to capture nonlinear effects, we calculate time-dependent structure factors by combining DDFT with Percus’ test particle method. For the growth of the average domain size $l$ during spinodal decomposition with time $t$, we observe a power-law behavior $l\propto t^{{\delta }_{\alpha }}$ with ${\delta }_m\simeq 0.333$ for the magnetic species and ${\delta }_n\simeq 0.323$ for the nonmagnetic species.

In this paper we extend the recent theory of shear localization in two-dimensional (2D) amorphous solids to three dimensions. In two dimensions the fundamental instability of shear localization is related to the appearance of a line of displacement quadrupoles that makes an angle of 45${}^{\circ }$ with the principal stress axis. In three dimensions the fundamental plastic instability is also explained by the formation of a lattice of anisotropic elastic inclusions. In the case of pure external shear stress, we demonstrate that this is a 2D triangular lattice of similar elementary events. It is shown that this lattice is arranged on a plane that, similarly to the 2D case, makes an angle of 45${}^{\circ }$ with the principal stress axis. This solution is energetically favorable only if the external strain exceeds a yield-strain value that is determined by the strain parameters of the elementary events and the Poisson ratio. The predictions of the theory are compared to numerical simulations and very good agreement is observed.

There is much interest in building deterministic continuum models from discrete agent-based models governed by local stochastic rules where an agent represents a biological cell. In developmental biology, cells are able to move and undergo cell division on and within growing tissues. A growing tissue is itself made up of cells which undergo cell division, thereby providing a significant transport mechanism for other cells within it. We develop a discrete agent-based model where domain agents represent tissue cells. Each agent has the ability to undergo a proliferation event whereby an additional domain agent is incorporated into the lattice. If a probability distribution describes the waiting times between proliferation events for an individual agent, then the total length of the domain is a random variable. The average behavior of these stochastically proliferating agents defining the growing lattice is determined in terms of a Fokker-Planck equation, with an advection and diffusion term. The diffusion term differs from the one obtained Landman and Binder [J. Theor. Biol. 259, 541 (2009)] when the rate of growth of the domain is specified, but the choice of agents is random. This discrepancy is reconciled by determining a discrete-time master equation for this process and an associated asymmetric nonexclusion random walk, together with consideration of synchronous and asynchronous updating schemes. All theoretical results are confirmed with numerical simulations. This study furthers our understanding of the relationship between agent-based rules, their implementation, and their associated partial differential equations. Since tissue growth is a significant cellular transport mechanism during embryonic growth, it is important to use the correct partial differential equation description when combining with other cellular functions.

Damage growth in composite materials is a complex process which is of interest in many fields of science and engineering. We consider this problem in a fiber bundle model where fibers undergo an aging process due to the accumulation of damage driven by the locally acting stress in a chemically active environment. By subjecting the bundle to a constant external load, fibers fail either when the load on them exceeds their individual intrinsic strength or when the accumulated internal damage exceeds a random threshold. We analyze the time evolution of the breaking process under low external loads where aging of fibers dominates. In the mean field limit, we show analytically that the aging system continuously accelerates in a way which can be characterized by an inverse power law of the event rate with a singularity that defines a failure time. The exponent is not universal; it depends on the details of the aging process. For localized load sharing, a more complex damage process emerges which is dominated by distinct spatial regions of the system with different degrees of stress concentration. Analytical calculations revealed that the final acceleration to global failure is preceded by a stationary accumulation of damage. When the disorder is strong, the accelerating phase has the same functional behavior as in the mean field limit. The analytical results are verified by computer simulations.

In the present study, we focus on evaluating the induced streaming electric field along the orthogonal directions in a narrow fluidic confinement in the presence of patterned surface wettability and modulated surface charges. We attempt to assess the implications of such modulations on the related important quantities and pinpoint the regimes of improved induced streaming potential field and the resulting anisotropy in the induced potential. Our results reveal that for certain combinations of the parameters characterizing the modulated slip, a significant amount of augmentation in the streaming electric field might be obtained, whereas in other cases the effects may lead to adverse consequences. We further demonstrate that the presence of anisotropic modulations on the channel walls give rise to considerable off-diagonal effects, which makes the streaming potential “disoriented” with the applied pressure gradient, when the same is not applied along one of the orthogonal directions. Our analysis also shows that one can remove such “mis-orientations” by finely tuning several relevant flow and geometric parameters, which may bear immense scientific and technological consequences towards an improved design of miniaturized energy conversion devices.

We investigated the cluster-size distribution of the residual nonwetting phase in a sintered glass-bead porous medium at two-phase flow conditions, by means of micro-computed-tomography ($\mu$CT) imaging with pore-scale resolution. Cluster-size distribution functions and cluster volumes were obtained by image analysis for a range of injected pore volumes under both imbibition and drainage conditions; the field of view was larger than the porosity-based representative elementary volume (REV). We did not attempt to make a definition for a two-phase REV but used the nonwetting-phase cluster-size distribution as an indicator. Most of the nonwetting-phase total volume was found to be contained in clusters that were one to two orders of magnitude larger than the porosity-based REV. The largest observed clusters in fact ranged in volume from 65% to 99% of the entire nonwetting phase in the field of view. As a consequence, the largest clusters observed were statistically not represented and were found to be smaller than the estimated maximum cluster length. The results indicate that the two-phase REV is larger than the field of view attainable by $\mu$CT scanning, at a resolution which allows for the accurate determination of cluster connectivity.

We study a number of filtering schemes for the reduction of the statistical error in nonadiabatic calculations by means of the quantum-classical Liouville equation. In particular, we focus on a scheme based on setting a threshold value on the sampling weights, so that when the threshold is overcome the value of the weight is reset, and on another approach which prunes the ensemble of the allowed nonadiabatic transitions according to a generalized sampling probability. Both methods have advantages and drawbacks, however, their combination drastically improves the performance of an algorithm known as the sequential short-time step propagation [MacKernan et al., J. Phys: Condens. Matter 14, 9069 (2002)], which is derived from a simple first order expansion of the quantum-classical propagator. Such an algorithm together with the combined filtering procedures produce results that compare very well with those obtained by means of numerically accurate path integral quantum calculations for the spin-boson model, even for intermediate and strong coupling regimes.

This paper studies the relaxation of the molecular spin angular velocity in the framework of generalized extended Navier-Stokes theory. Using molecular dynamics simulations, it is shown that for uncharged diatomic molecules the relaxation time decreases with increasing molecular moment of inertia per unit mass. In the regime of large moment of inertia the fast relaxation is wave-vector independent and dominated by the coupling between spin and the fluid streaming velocity, whereas for small inertia the relaxation is slow and spin diffusion plays a significant role. The fast wave-vector-independent relaxation is also observed for highly packed systems. The transverse and longitudinal spin modes have, to a good approximation, identical relaxation, indicating that the longitudinal and transverse spin viscosities have same value. The relaxation is also shown to be isomorphic invariant. Finally, the effect of the coupling in the zero frequency and wave-vector limit is quantified by a characteristic length scale; if the system dimension is comparable to this length the coupling must be included into the fluid dynamical description. It is found that the length scale is independent of moment of inertia but dependent on the state point.

Examples of self-propulsion in strongly fluctuating environments are abundant in nature, e.g., molecular motors and pumps operating in living cells. Starting from the Langevin equation of motion, we develop a stochastic thermodynamic description of noninteracting self-propelled particles using simple models of velocity-dependent forces. We derive fluctuation theorems for entropy production and a modified fluctuation-dissipation relation, characterizing the linear response in nonequilibrium steady states. We study these notions in a simple model of molecular motors, and in the Rayleigh-Helmholtz and energy-depot models of self-propelled particles.

A reciprocating quantum refrigerator is analyzed with the intention to study the limitations imposed by external noise. In particular we focus on the behavior of the refrigerator when it approaches the absolute zero. The cooling cycle is based on the Otto cycle with a working medium constituted by an ensemble of noninteracting harmonic oscillators. The compression and expansion segments are generated by changing an external parameter in the Hamiltonian. In this case the force constant of the harmonic oscillators $m{\omega }^2$ is modified from an initial to a final value. As a result, the kinetic and potential energy of the system do not commute causing frictional losses. By proper choice of scheduling function $\omega \left(t\right)$ frictionless solutions can be obtained in the noiseless case. We examine the performance of a refrigerator subject to noise. By expanding from the adiabatic limit we find that the external noise, Gaussian phase, and amplitude noises reduce the amount of heat that can be extracted but nevertheless the zero temperature can be approached.

In this article we address the problem of the nonlinear interaction of subdiffusive particles. We introduce the random walk model in which statistical characteristics of a random walker such as escape rate and jump distribution depend on the mean density of particles. We derive a set of nonlinear subdiffusive fractional master equations and consider their diffusion approximations. We show that these equations describe the transition from an intermediate subdiffusive regime to asymptotically normal advection-diffusion transport regime. This transition is governed by nonlinear tempering parameter that generalizes the standard linear tempering. We illustrate the general results through the use of the examples from cell and population biology. We find that a nonuniform anomalous exponent has a strong influence on the aggregation phenomenon.

Aiming at the mechanical properties of cross-linked biopolymers, we set up and analyze a model of two weakly bending wormlike chains subjected to a tensile force, with regularly spaced inter-chain bonds (cross-links) represented by harmonic springs. Within this model, we compute the force-extension curve and the differential stiffness exactly and discuss several limiting cases. Cross-links effectively stiffen the chain pair by reducing thermal fluctuations transverse to the force and alignment direction. The extra alignment due to cross-links increases both with growing number and with growing strength of the cross-links, and is most prominent for small force $f$. For large $f$, the additional, cross-link-induced extension is subdominant except for the case of linking the chains rigidly and continuously along their contour. In this combined limit, we recover asymptotically the elasticity of a weakly bending wormlike chain without constraints, stiffened by a factor of 4. The increase in differential stiffness can be as large as $100\%$ for small $f$ or large numbers of cross-links.

The sessile microbial communities known as biofilms exhibit varying architectures as environmental factors are varied, which for immersed biofilms includes the shear rate of the surrounding flow. Here we modify an established agent-based biofilm model to include affine flow and employ it to analyze the growth of surface roughness of single-species, three-dimensional biofilms. We find linear growth laws for surface geometry in both horizontal and vertical directions and measure the thickness of the active surface layer, which is shown to anticorrelate with roughness. Flow is shown to monotonically reduce surface roughness without affecting the thickness of the active layer. We argue that the rapid roughening is due to nonlocal surface interactions mediated by the nutrient field, which are curtailed when advection competes with diffusion. We further argue the need for simplified models to elucidate the underlying mechanisms coupling flow to growth.

Migratory birds can utilize the geomagnetic field for orientation and navigation through a widely accepted radical-pair mechanism. Although many theoretical works have been done, the available experimental results have not been fully considered, especially the temporary disorientation induced by the field which is increased by $30\%$ of the geomagnetic field and the disorientation of the very weak resonant field of 15 nT. In this paper, we consider the monotonicity of the singlet yield angular profile as the prerequisite of direction sensitivity, and find that for some optimal values of the hyperfine coupling parameters (that is, the order of ${10}^{-7}\sim {10}^{-6}$ meV) the experimental results available so far can be satisfied. We also investigate the effects of two decoherence environments and demonstrate that, in order to satisfy the available experimental results, the decoherence rate should be lower than the recombination rate. Finally, we investigate the effects of the fluctuating magnetic noises and find that the vertical noise destroys the monotonicity of the profile completely, but the parallel noise preserves the monotonicity perfectly and even can enhance the direction sensitivity.

Simple physical models based on fluid mechanics have long been used to understand the flow of vehicular traffic on freeways; analytically tractable models of flow on an urban grid, however, have not been as extensively explored. In an ideal world, traffic signals would be timed such that consecutive lights turned green just as vehicles arrived, eliminating the need to stop at each block. Unfortunately, this “green-wave” scenario is generally unworkable due to frustration imposed by competing demands of traffic moving in different directions. Until now this has typically been resolved by numerical simulation and optimization. Here, we develop a theory for the flow in an idealized system consisting of a long two-way road with periodic intersections. We show that optimal signal timing can be understood analytically and that there are counterintuitive asymmetric solutions to this signal coordination problem. We further explore how these theoretical solutions degrade as traffic conditions vary and automotive density increases.

We consider a system of two or four nonlinear sites coupled with binary chain waveguides. When a monochromatic wave is injected into the first (symmetric) propagation channel, the presence of cubic nonlinearity can lead to symmetry breaking, giving rise to emission of antisymmetric wave into the second (antisymmetric) propagation channel of the waveguides. We found that in the case of nonlinear plaquette, there is a domain in the parameter space where neither symmetry-preserving nor symmetry-breaking stable stationary solutions exit. As a result, injection of a monochromatic symmetric wave gives rise to emission of nonsymmetric satellite waves with energies differing from the energy of the incident wave. Thus, the response exhibits nonmonochromatic behavior.

We study chimeric states in a coupled phase oscillator system with piecewise linear nonlocal coupling. By modifying the details of the coupling, it is possible to obtain multiple chimeric states with a specified number of coherent regions and with specified phase relationships. The case of a two-component chimera is illustrated and the generalization to arbitrary chimeric configurations is discussed. The phase relations between the two clusters of phase oscillators is described in some detail.