We study a switching synchronization phenomenon taking place in one-dimensional memristive networks when the memristors switch from the high- to low-resistance state. It is assumed that the distributions of threshold voltages and switching rates of memristors are arbitrary. Using the Laplace transform, a set of nonlinear equations describing the memristors dynamics is solved exactly, without any approximations. The time dependencies of memristances are found, and it is shown that the voltage falls across memristors are proportional to their threshold voltages. A compact expression for the network switching time is derived.

Filamentous semiflexible networks define the mechanical and physical properties of many materials such as cytoskeleton. In the absence of a distinct unit cell, the Mikado fiber network model is commonly used algorithm for representing the microstructure of these networks in numerical models. Nevertheless, certain types of filamentous structures such as collagenous tissues, at early stages of their development, are assembled by growth of individual fibers from random nucleation sites. In this work, we develop a computational model to investigate the mechanical response of such networks by characterizing their nonaffine behavior. We show that the deformation of these networks is nonaffine at all length scales. Furthermore, similar to Mikado networks, the degree of nonaffinity in these structures decreases with increasing the probing length scale, the network fiber density, and/or the bending stiffness of constituting filaments. Nevertheless, despite the lower coordination number of these networks, their deformation field is more affine than that of the Mikado networks with the same fiber density and fiber mechanical properties.

We formulate a numerical method for predicting the tensorial linear response of a rigid, asymmetrically charged body to an applied electric field. This prediction requires calculating the response of the fluid to the Stokes drag forces on the moving body and on the countercharges near its surface. To determine the fluid's motion, we represent both the body and the countercharges using many point sources of drag known as Stokeslets. Finding the correct flow field amounts to finding the set of drag forces on the Stokeslets that is consistent with the relative velocities experienced by each Stokeslet. The method rigorously satisfies the condition that the object moves with no transfer of momentum to the fluid. We demonstrate that a sphere represented by 1999 well-separated Stokeslets on its surface produces flow and drag force like a solid sphere to 1% accuracy. We show that a uniformly charged sphere with 3998 body and countercharge Stokeslets obeys the Smoluchowski prediction [F. Morrison, J. Colloid Interface Sci. 34, 210 (1970)] for electrophoretic mobility when the countercharges lie close to the sphere. Spheres with dipolar and quadrupolar charge distributions rotate and translate as predicted analytically to 4% accuracy or better. We describe how the method can treat general asymmetric shapes and charge distributions. This method offers promise as a way to characterize and manipulate asymmetrically charged colloid-scale objects from biology (e.g., viruses) and technology (e.g., self-assembled clusters).

We investigate the dependence of the displacements of a molecular motor embedded inside a glassy material on its folding characteristic time ${\tau }_f$. We observe two different time regimes. For slow foldings (regime I) the diffusion evolves very slowly with ${\tau }_f$, while for rapid foldings (regime II) the diffusion increases strongly with ${\tau }_f$($D\approx {\tau }_f^{-2}$), suggesting two different physical mechanisms. We find that in regime I the motor's displacement during the folding process is counteracted by a reverse displacement during the unfolding, while in regime II this counteraction is much weaker. We notice that regime I behavior is reminiscent of the scallop theorem that holds for larger motors in a continuous medium. We find that the difference in the efficiency of the motor's motion explains most of the observed difference between the two regimes. For fast foldings the motor trajectories differ significantly from the opposite trajectories induced by the following unfolding process, resulting in a more efficient global motion than for slow foldings. This result agrees with the fluctuation theorems expectation for time reversal mechanisms. In agreement with the fluctuation theorems we find that the motors are unexpectedly more efficient when they are generating more entropy, a result that can be used to increase dramatically the motor's motion.

Kinetic theory (KT) has been successfully used to model rapid granular flows in which particle interactions are frictionless and near elastic. However, it fails when particle interactions become frictional and inelastic. For example, the KT is not able to accurately predict the free cooling process of a vibrated granular medium that consists of inelastic frictional particles under microgravity. The main reason that the classical KT fails to model these flows is due to its inability to account for the particle surface friction and its inelastic behavior, which are the two most important factors that need be considered in modeling collisional granular flows. In this study, we have modified the KT model that is able to incorporate these two factors. The inelasticity of a particle is considered by establishing a velocity-dependent expression for the restitution coefficient based on many experimental studies found in the literature, and the particle friction effect is included by using a tangential restitution coefficient that is related to the particle friction coefficient. Theoretical predictions of the free cooling process by the classical KT and the improved KT are compared with the experimental results from a study conducted on an airplane undergoing parabolic flights without the influence of gravity [Y. Grasselli, G. Bossis, and G. Goutallier, Europhys. Lett. 86, 60007 (2009)]. Our results show that both the velocity-dependent restitution coefficient and the particle surface friction are important in predicting the free cooling process of granular flows; the modified KT model that integrates these two factors is able to improve the simulation results and leads to better agreement with the experimental results.

We characterize the water repartition within the partially saturated (two-phase) zone (PSZ) during evaporation from mixed wettable porous media by controlling the wettability of glass beads, their sizes, and as well the surrounding relative humidity. Here, capillary numbers are low and under these conditions, the percolating front is stabilized by gravity. Using experimental and numerical analyses, we find that the PSZ saturation decreases with the Bond number, where packing of smaller particles have higher saturation values than packing made of larger particles. Results also reveal that the extent (height) of the PSZ, as well as water saturation in the PSZ, both increase with wettability. We also numerically calculate the saturation exclusively contained in connected liquid films and results show that values are less than the expected PSZ saturation. These results strongly reflect that the two-phase zone is not solely made up of connected capillary networks but also made of disconnected water clusters or pockets. Moreover, we also find that global saturation (PSZ + full wet zone) decreases with wettability, confirming that greater quantity of water is lost via evaporation with increasing hydrophilicity. These results show that connected liquid films are favored in more-hydrophilic systems while disconnected water pockets are favored in less-hydrophilic systems.

We propose a Monte Carlo algorithm designed to simulate quantum as well as classical systems at equilibrium, bridging the algorithmic gap between quantum and classical thermal simulation algorithms. The method is based on a decomposition of the quantum partition function that can be viewed as a series expansion about its classical part. We argue that the algorithm not only provides a theoretical advancement in the field of quantum Monte Carlo simulations, but is optimally suited to tackle quantum many-body systems that exhibit a range of behaviors from “fully quantum” to “fully classical,” in contrast to many existing methods. We demonstrate the advantages, sometimes by orders of magnitude, of the technique by comparing it against existing state-of-the-art schemes such as path integral quantum Monte Carlo and stochastic series expansion. We also illustrate how our method allows for the unification of quantum and classical thermal parallel tempering techniques into a single algorithm and discuss its practical significance.

A polymer repelled by unfavorable interactions with a uniform flat surface may still be pinned to attractive edges and corners. This is demonstrated by considering adsorption of a two-dimensional ideal polymer to an attractive corner of a repulsive wedge. The well-known mapping between the statistical mechanics of an ideal polymer and the quantum problem of a particle in a potential is then used to analyze the singular behavior of the unbinding transition of the polymer. The divergence of the localization length is found to be governed by an exponent that varies continuously with the angle (when reflex). Numerical treatment of the discrete (lattice) version of such an adsorption problem confirms this behavior.

We report the phenomenon of self-organized emergence of hierarchical multilayered structures and chimera states in dynamical networks with adaptive couplings. This process is characterized by a sequential formation of subnetworks (layers) of densely coupled elements, the size of which is ordered in a hierarchical way, and which are weakly coupled between each other. We show that the hierarchical structure causes the decoupling of the subnetworks. Each layer can exhibit either a two-cluster state, a periodic traveling wave, or an incoherent state, and these states can coexist on different scales of subnetwork sizes.

A metapopulation structure in landscape ecology comprises a group of interacting spatially separated subpopulations or patches of the same species that may experience several local extinctions. This makes the investigation of survivability (in the form of global oscillation) of a metapopulation on top of diverse dispersal topologies extremely crucial. However, among various dispersal topologies in ecological networks, which one can provide higher metapopulation survivability under local extinction is still not well explored. In this article, we scrutinize the robustness of an ecological network consisting of prey-predator patches having Holling type I functional response, against progressively extinct population patches. We present a comprehensive study on this while considering global, small-world, and scale-free dispersal of the subpopulations. Furthermore, we extend our work in enhancing survivability in the form of sustained global oscillation by introducing asymmetries in the dispersal rates of the considered species. Our findings affirm that the asynchrony among the patches plays an important role in the survivability of a metapopulation. In order to demonstrate the model independence of the observed phenomenon, we perform a similar analysis for patches exhibiting Holling type II functional response. On the grounds of the obtained results, our work is expected to provide a better perception of the influence of dispersal arrangements on the global survivability of ecological networks.

Flexible magnetic filaments have garnered considerable attention as prospective materials for the creation of different microdevices. We describe a theoretical model of a ferromagnetic filament and derive its equations of motion by variational techniques. The numerical algorithm used to solve the filament dynamics in magnetic fields of different configurations is described. It is found that in a rotating field the filament transitions between synchronous and asynchronous regimes with respect to the rotating field, similarly to a rigid magnetic dipole. The mean angular velocity of the filament is well described by a relation valid for a rigid magnetic dipole with quantitative differences attributable to the flexibility of the filament.

The work is focused on obtaining boundary conditions for a one-sided numerical model of thermoconvective instabilities in evaporating pinned sessile droplets of ethanol on heated substrates. In the one-sided model, appropriate boundary conditions for heat and mass transfer equations are required at the droplet surface. Such boundary conditions are obtained in the present work based on a derived semiempirical theoretical formula for the total droplet's evaporation rate, and on a two-parametric nonisothermal approximation of the local evaporation flux. The main purpose of these boundary conditions is to be applied in future three-dimensional (3D) one-sided numerical models in order to save a lot of computational time and resources by solving equations only in the droplet domain. Two parameters, needed for the nonisothermal approximation of the local evaporation flux, are obtained by fitting computational results of a 2D two-sided numerical model. Such model is validated here against parabolic flight experiments and the theoretical value of the total evaporation rate. This study combines theoretical, experimental, and computational approaches in convective evaporation of sessile droplets. The influence of the gravity level on evaporation rate and contributions of different mechanisms of vapor transport (diffusion, Stefan flow, natural convection) are shown. The qualitative difference (in terms of developing thermoconvective instabilities) between steady-state and unsteady numerical approaches is demonstrated.

A substantial viscosity increase by the addition of a low dose of nanoparticles to the base fluids can well influence the dynamics of viscous fingering. There is a lack of detailed theoretical studies that address the effect of the presence of nanoparticles on unstable miscible displacements. In this study, the impact of nonreactive nanoparticle presence on the stability and subsequent mixing of an originally unstable binary system is examined using linear stability analysis (LSA) and pseudospectral-based direct numerical simulations (DNS). We have parametrized the role of both nondepositing and depositing nanoparticles on the stability of miscible displacements using the developed static and dynamic parametric analyses. Our results show that nanoparticles have the potential to weaken the instabilities of an originally unstable system. Our LSA and DNS results also reveal that nondepositing nanoparticles can be used to fully stabilize an originally unstable front while depositing particles may act as temporary stabilizers whose influence diminishes in the course of time. In addition, we explain the existing inconsistencies concerning the effect of the nanoparticle diffusion coefficient on the dynamics of the system. This study provides a basis for further research on the application of nanoparticles for control of viscosity-driven instabilities.

We propose a quantum heat engine composed of two superconducting transmission line resonators interacting with each other via an optomechanical-like coupling. One resonator is periodically excited by a thermal pump. The incoherently driven resonator induces coherent oscillations in the other one due to the coupling. A limit cycle, indicating finite power output, emerges in the thermodynamical phase space. The system implements an all-electrical analog of a photonic piston. Instead of mechanical motion, the power output is obtained as a coherent electrical charging in our case. We explore the differences between the quantum and classical descriptions of our system by solving the quantum master equation and classical Langevin equations. Specifically, we calculate the mean number of excitations, second-order coherence, as well as the entropy, temperature, power, and mean energy to reveal the signatures of quantum behavior in the statistical and thermodynamic properties of the system. We find evidence of a quantum enhancement in the power output of the engine at low temperatures.

We study long-range power-law correlated disorder on square and cubic lattices. In particular, we present high-precision results for the percolation thresholds and the fractal dimension of the largest clusters as a function of the correlation strength. The correlations are generated using a discrete version of the Fourier filtering method. We consider two different metrics to set the length scales over which the correlations decay, showing that the percolation thresholds are highly sensitive to such system details. By contrast, we verify that the fractal dimension $d_{\mathrm{f}}$ is a universal quantity and unaffected by the choice of metric. We also show that for weak correlations, its value coincides with that for the uncorrelated system. In two dimensions we observe a clear increase of the fractal dimension with increasing correlation strength, approaching $d_{\mathrm{f}}\to 2$. The onset of this change does not seem to be determined by the extended Harris criterion.

We study the behavior of a classical two-component ionic plasma made up of nonadditive hard disks with additional logarithmic Coulomb interactions between them. Due to the Coulomb repulsion, long-wavelength total density fluctuations are suppressed and the system is globally hyperuniform. Short-range volume effects lead to phase separation or to heterocoordination for positive or negative nonadditivities, respectively. These effects compete with the hidden long-range order imposed by hyperuniformity. As a result, the critical behavior of the mixture is modified, with long-wavelength concentration fluctuations partially damped when the system is charged. It is also shown that the decrease of configurational entropy due to hyperuniformity originates from contributions beyond the two-particle level. Finally, despite global hyperuniformity, we show that in our system the spatial configuration associated with each component separately is not hyperuniform, i.e., the system is not “multihyperuniform.”

The Ising model in two dimensions with special toroidal boundary conditions is analyzed. These boundary conditions, which we call duality-twisted boundary conditions, may be interpreted as inserting a specific defect line (“seam”) in the system, along noncontractible circles of the cylinder, before closing it into a torus. We derive exact expressions for the eigenvalues of a transfer matrix for the critical ferromagnetic Ising model on the $M\times N$ square lattice wrapped on the torus with a specific defect line. As a result we have obtained analytically the partition function for the Ising model with such boundary conditions. In the case of infinitely long cylinders of circumference L with duality-twisted boundary conditions we obtain the asymptotic expansion of the free energy and the inverse correlation lengths. We find that the ratio of subdominant finite-size correction terms in the asymptotic expansion of the free energy and the inverse correlation lengths should be universal. We verify such universal behavior in the framework of a perturbating conformal approach by calculating the universal structure constant $C_{n1n}$ for descendent states generated by the operator product expansion of the primary fields. For such states the calculations of an universal structure constants is a difficult task, since it involves knowledge of the four-point correlation function, which in general is not fixed by conformal invariance except for some particular cases, including the Ising model.

We consider the advection-diffusion transport of tracers in a one-parameter family of plane periodic flows where the patterns of streamlines feature regions of confined circulation in the shape of “cat's eyes,” separated by meandering jets with ballistic motion inside them. By varying the parameter, we proceed from the regular two-dimensional lattice of eddies without jets to the sinusoidally modulated shear flow without eddies. When a weak thermal noise is added, i.e., at large Péclet numbers, several intermediate time scales arise, with qualitatively and quantitatively different transport properties: depending on the parameter of the flow, the initial position of a tracer, and the aging time, motion of the tracers ranges from subdiffusive to superballistic. We report on results of extensive numerical simulations of the mean-squared displacement for different initial conditions in ordinary and aged situations. These results are compared with a theory based on a Lévy walk that describes the intermediate-time ballistic regime and gives a reasonable description of the behavior for a certain class of initial conditions. The interplay of the walk process with internal circulation dynamics in the trapped state results at intermediate time scales in nonmonotonic characteristics of aging not captured by the Lévy walk model.

The problem of zero crossings is of great historical prevalence and promises extensive application. The challenge is to establish precisely how the autocorrelation function or power spectrum of a one-dimensional continuous random process determines the density function of the intervals between the zero crossings of that process. This paper investigates the case where periodicities are incorporated into the autocorrelation function of a smooth process. Numerical simulations, and statistics about the number of crossings in a fixed interval, reveal that in this case the zero crossings segue between a random and deterministic point process depending on the relative time scales of the periodic and nonperiodic components of the autocorrelation function. By considering the Laplace transform of the density function, we show that incorporating correlation between successive intervals is essential to obtaining accurate results for the interval variance. The same method enables prediction of the density function tail in some regions, and we suggest approaches for extending this to cover all regions. In an ever-more complex world, the potential applications for this scale of regularity in a random process are far reaching and powerful.

Brownian systems often surmount energy barriers by absorbing and emitting heat to and from their local environment. Usually, the temperature gradients created by this heat exchange are assumed to dissipate instantaneously. Here we relax this assumption to consider the case where Brownian dynamics on a time-independent potential can lead to self-induced temperature gradients. In the same way that externally imposed temperature gradients can cause directed motion, these self-induced gradients affect the dynamics of the Brownian system. The result is a coupling between the local environment and the Brownian subsystem. We explore the resulting dynamics and thermodynamics of these coupled systems and develop a robust method for numerical simulation. In particular, by focusing on one-dimensional situations, we show that self-induced temperature gradients reduce barrier-crossing rates. We also consider a heat engine and a heat pump based on temperature gradients induced by a Brownian system in a nonequilibrium potential.