Search

From length scale distributions characterizing frozen amorphous domains, we relate the s ensemble method with standard cooling and aging protocols for formin0 glass. We show that in a class of models where space-time scaling is in harmony with that of experiment, the spatial distributions of excitations obtained with the .v ensemble are identical to those obtained through cooling or aging, but the computational effort for applying the s ensemble is generally many orders of magnitude smaller than that of straightforward numerical simulation of cooling or aging. We find that in contrast to the equilibrium ergodic state, a nonequilibrium length scale characterizes the anticorrelation between excitations and encodes the preparation history of glass states. [ABSTRACT FROM AUTHOR]

We present a thermodynamic theory for a generic population of M individuals distributed into N groups (clusters). We construct the ensemble of all distributions with fixed M and N, introduce a selection functional that embodies the physics that governs the population, and obtain the distribution that emerges in the scaling limit as the most probable among all distributions consistent with the given physics. We develop the thermodynamics of the ensemble and establish a rigorous mapping to regular thermodynamics. We treat the emergence of a so-called giant component as a formal phase transition and show that the criteria for its emergence are entirely analogous to the equilibrium conditions in molecular systems. We demonstrate the theory by an analytic model and confirm the predictions by Monte Carlo simulation. [ABSTRACT FROM AUTHOR]

We perform numerical measurements of the moments of the position of a tracer particle in a two-dimensional periodic billiard model (Lorentz gas) with infinite corridors. This model is known to exhibit a weak form of superdiffusion, in the sense that there is a logarithmic correction to the linear growth in time of the meansquared displacement. We show numerically that this expected asymptotic behavior is easily overwhelmed by the subleading linear growth throughout the time range accessible to numerical simulations. We compare our simulations to analytical results for the variance of the anomalously rescaled limiting normal distributions. [ABSTRACT FROM AUTHOR]

Diffusion in nonhomogeneous media is described by a dynamical process driven by a general Lévy noise and subordinated to a random time; the subordinator depends on the position. This problem is approximated by a multiplicative process subordinated to a random time: it separately takes into account effects related to the medium structure and the memory. Density distributions and moments are derived from the solutions of the corresponding Langevin equation and compared with the numerical calculations for the exact problem. Both subdiffusion and enhanced diffusion are predicted. Distribution of the process satisfies the fractional Fokker-Planck equation. [ABSTRACT FROM AUTHOR]

Probing deeper into the existing issues regarding the exit probability (EP) in one-dimensional dynamical models, we consider several models where the states are represented by Ising spins and the information flows inwards. At zero temperature, these systems evolve to either of two absorbing states. The EP, E(x), which is the probability that the system ends up with all up spins starting with x fraction of up spins, is found to have the general form E(x) = xα/[xα + (1 -- x)α]. The EP exponent a strongly depends on r, the range of interaction, the symmetry of the model, and the induced fluctuation. Even in a nearest-neighbor model, a nonlinear form of the EP can be obtained by controlling the fluctuations, and for the same range, different models give different results for a. Nonuniversal behavior of the EP is thus clearly established and the results are compared to those of existing studies in models with outflow dynamics to distinguish the two dynamical scenarios. [ABSTRACT FROM AUTHOR]

The ergodicity breaking parameter is a measure for the heterogeneity among different trajectories of one ensemble. In this report, this parameter is calculated for fractional Brownian motion with a random change of time scale, often called "subordination." We show that this quantity is the same as the known continuous time random walks case. [ABSTRACT FROM AUTHOR]

The adaptive voter model is a paradigmatic model in the study of opinion formation. Here we propose an extension for this model, in which conflicts are resolved by obtaining another opinion, and analytically study the time required for consensus to emerge. Our results shed light on the rich phenomenology of both the original and extended adaptive voter models, including a dynamical phase transition in the scaling behavior of the mean time to consensus. [ABSTRACT FROM AUTHOR]

Trapping processes constitute a primary problem of random walks, which characterize various other dynamical processes taking place on networks. Most previous works focused on the case of binary networks, while there is much less related research about weighted networks. In this paper, we propose a general framework for the trapping problem on a weighted network with a perfect trap fixed at an arbitrary node. By utilizing the spectral graph theory, we provide an exact formula for mean first-passage time (MFPT) from one node to another, based on which we deduce an explicit expression for average trapping time (ATT) in terms of the eigenvalues and eigenvectors of the Laplacian matrix associated with the weighted graph, where ATT is the average of MFPTs to the trap over all source nodes. We then further derive a sharp lower bound for the ATT in terms of only the local information of the trap node, which can be obtained in some graphs. Moreover, we deduce the ATT when the trap is distributed uniformly in the whole network. Our results show that network weights play a significant role in the trapping process. To apply our framework, we use the obtained formulas to study random walks on two specific networks: trapping in weighted uncorrelated networks with a deep trap, the weights of which are characterized by a parameter, and Levy random walks in a connected binary network with a trap distributed uniformly, which can be looked on as random walks on a weighted network. For weighted uncorrelated networks we show that the ATT to any target node depends on the weight parameter, that is, the ATT to any node can change drastically by modifying the parameter, a phenomenon that is in contrast to that for trapping in binary networks. For Levy random walks in any connected network, by using their equivalence to random walks on a weighted complete network, we obtain the optimal exponent characterizing Levy random walks, which have the minimal average of ATTs taken over all target nodes. [ABSTRACT FROM AUTHOR]

SLE(k,p) is a variant of Schramm-Loewner Evolution (SLE) which describes the curves which are not conformal invariant, but are self-similar due to the presence of some other preferred points on the boundary. In this paper we study the left passage probability (LPP) of SLE(fc,p) through field theoretical framework and find the differential equation governing this probability. This equation is numerically solved for the special case K = 2 and hp = 0 in which hp is the conformal weight of the boundary changing (bcc) operator. It may be referred to loop erased random walk (LERW) and Abelian sandpile model (ASM) with a sink on its boundary. For the curve which starts from fo and conditioned by a change of boundary conditions at x0, we find that this probability depends significantly on the factor xo -- £o- We also present the perturbative general solution for large x0. As a prototype, we apply this formalism to SLE0k,k -- 6) which governs the curves that start from and end on the real axis. [ABSTRACT FROM AUTHOR]

We consider multiple diseases spreading in a static configuration model network. We make standard assumptions that infection transmits from neighbor to neighbor at a disease-specific rate and infected individuals recover at a disease-specific rate. Infection by one disease confers immediate and permanent immunity to infection by any disease. Under these assumptions, we find a simple, low-dimensional ordinary differential equations model which captures the global dynamics of the infection. The dynamics depend strongly on initial conditions. Although we motivate this Rapid Communication with infectious disease, the model may be adapted to the spread of other infectious agents such as competing political beliefs, or adoption of new technologies if these are influenced by contacts. As an example, we demonstrate how to model an infectious disease which can be prevented by a behavior change. [ABSTRACT FROM AUTHOR]

Percolation describes the sudden emergence of large-scale connectivity as edges are added to a lattice or random network. In the Bohman-Frieze-Wormald model (BFW) of percolation, edges sampled from a random graph are considered individually and either added to the graph or rejected provided that the fraction of accepted edges is never smaller than a decreasing function with asymptotic value of α, a constant. The BFW process has been studied as a model system for investigating the underlying mechanisms leading to discontinuous phase transitions in percolation. Here we focus on the regime α ∈ [0.6,0.95] where it is known that only one giant component, denoted C1, initially appears at the discontinuous phase transition. We show that at some point in the supercritical regime C1 stops growing and eventually a second giant component, denoted emerges in a continuous percolation transition. The delay between the emergence of Q and C2 and their asymptotic sizes both depend on the value of a and we establish by several techniques that there exists a bifurcation point αc = 0.763 ± 0.002. For α ∈ [0.6,αc), Q stops growing the instant it emerges and the delay between the emergence of Q and C2 decreases with increasing a. For a ∈ (αc,0.95], in contrast, Q continues growing into the supercritical regime and the delay between the emergence of C1 and C2 increases with increasing α. As we show, αc marks the minimal delay possible between the emergence of C1 and C2 (i.e., the smallest edge density for which C2 can exist). We also establish many features of the continuous percolation of C2 including scaling exponents and relations. [ABSTRACT FROM AUTHOR]

The chemical potentials of multicomponent fluids are derived in terms of the pair correlation functions for arbitrary number of components, interaction potentials, and dimensionality. The formally exact result is particularized to hard-sphere mixtures with zero or positive nonadditivity. As a simple application, the chemical potentials of three-dimensional additive hard-sphere mixtures are derived from the Perçus-Yevick theory and the associated equation of state is obtained. This Perçus-Yevick chemical-route equation of state is shown to be more accurate than the virial equation of state. An interpolation between the chemical-potential and compressibility routes exhibits a better performance than the well-known Boublik-Mansoori-Carnahan-Starling-Leland equation of state. [ABSTRACT FROM AUTHOR]

We aim at quantitatively determining transport parameters like conductivity, mean free path, etc., for simple models of spatially completely disordered quantum systems, comparable to 'the systems which are sometimes referred to as Lifshitz models. While some low-energy eigenstates in such models always show Anderson localization, we focus on models for which most states of the full spectrum are delocalized, i.e., on the metallic regime. For the latter we determine transport parameters in the limit of high temperatures and low fillings using linear response theory. The Einstein relation (proportionality of conductivity and diffusion coefficient) is addressed numerically and analytically and found to hold. Furthermore, we find the transport behavior for some models to be in accord with a Boltzmann equation, i.e., rather long mean free paths, exponentially decaying currents, while this does not apply to other models even though they are also almost completely delocalized. [ABSTRACT FROM AUTHOR]

The paper authored by Lima et al. [Phys. Rev. E 61, 3260 (2000)] has discussed the dispersion relation and Landau damping of a Langmuir wave in the context of the nonextensive statistics proposed by Tsallis. However, the results obtained in this paper are not appropriate. In this comment on the paper we shall derive the correct analytic formulas for both the dispersion relation and Landau damping in the Tsallis formalism. We hope that this comment will be useful in providing the correct results. [ABSTRACT FROM AUTHOR]

The principle of mixed convolved action provides a new rigorous weak variational formalism for a broad range of initial value problems in mathematical physics and mechanics. Here, the focus is initially on classical single-degree-of-freedom oscillators incorporating either Kelvin-Voigt or Maxwell dissipative elements and then, subsequently, on systems that utilize fractional-derivative constitutive models. In each case, an appropriate mixed convolved action is formulated, and a corresponding weak form is discretized in time using temporal shape functions to produce an algorithm suitable for numerical solution. Several examples are considered to validate the mixed convolved action principles and to investigate the performance of the numerical algorithms. For undamped systems, the algorithm is found to be symplectic and unconditionally stable with respect to the time step. In the case of dissipative systems, the approach is shown to be robust and to be accurate with good convergence characteristics for both classical and fractional-derivative based models. As part of the derivations, some interesting results in the calculus of Caputo fractional derivatives also are presented. [ABSTRACT FROM AUTHOR]

We experimentally investigate the impact of a liquid jet on a soap film. We observe that the jet never breaks the film and that two qualitatively different steady regimes may occur. The first one is a refractionlike behavior obtained at small incidence angles when the jet crosses the film and is deflected by the film-jet interaction. For larger incidence angles, the jet is absorbed by the film, giving rise to a new class of flows in which the jet undulates along the film with a characteristic wavelength. Besides its fundamental interest, this paper presents a different way to guide a micrometric flow of liquid in the inertial regime and to probe foam stability submitted to violent perturbations at the soap film scale. [ABSTRACT FROM AUTHOR]

The analytical two- and three-dimensional vortex solitons with arbitrary values of vorticity are constructed in the cubic defocusing media with spatially inhomogeneous nonlinearity. The values of the nonlinearity coefficients are zero near the center and increase rapidly toward the periphery. In addition to the analytical ones, a number of vortex solitons are found numerically. It is shown that analytical vortex solitons are stable. Also, the stability region of the numerically constructed vortex solitons are given [ABSTRACT FROM AUTHOR]

The influences of spatial frequency distributions on complete amplitude death are explored by studying an array of diffusively coupled oscillators. We found that with all possible sets of spatial frequency distributions, the two critical coupling strengths ∈c1 (lower-bounded value) and ∈c2 (upper-bounded value) needed to get complete amplitude death exhibit a universal power law and a log-normal distribution respectively, which has long tails in both cases. This is significant for dynamics control, since large variations of ∈c1 and ∈c2 are possible for some spatial arrangements. Moreover, we explore optimal spatial distributions with the smallest (largest) ∈c1 or ∈c1. [ABSTRACT FROM AUTHOR]

Many natural and artificial networks evolve in time. Nodes and connections appear and disappear at various time scales, and their dynamics has profound consequences for any processes in which they are involved. The first empirical analysis of the temporal patterns characterizing dynamic networks are still recent, so that many questions remain open. Here, we study how random walks, as a paradigm of dynamical processes, unfold on temporally evolving networks. To this aim, we use empirical dynamical networks of contacts between individuals, and characterize the fundamental quantities that impact any general process taking place upon them. Furthermore, we introduce different randomizing strategies that allow us to single out the role of the different properties of the empirical networks. We show that the random walk exploration is slower on temporal networks than it is on the aggregate projected network, even when the time is properly rescaled. In particular, we point out that a fundamental role is played by the temporal correlations between consecutive contacts present in the data. Finally, we address the consequences of the intrinsically limited duration of many real world dynamical networks. Considering the fundamental prototypical role of the random walk process, we believe that these results could help to shed light on the behavior of more complex dynamics on temporally evolving networks. [ABSTRACT FROM AUTHOR]

We demonstrate numerically bidirectional sorting of flocking particles interacting with an array of V-shaped asymmetric barriers. Each particle aligns with the average swimming direction of its neighbors according to the Vicsek model and experiences additional steric interactions as well as repulsion from the fixed barriers. We show that particles preferentially localize to one side of the barrier array over time and that the direction of this rectification can be reversed by adjusting the particle-particle exclusion radius or the noise term in the equations of motion. These results provide a conceptual basis for isolation and sorting of single-cell and multicellular organisms that move collectively according to flocking-type interaction rules. [ABSTRACT FROM AUTHOR]

Elements of networks interact in many ways, so modeling them with graphs requires multiple types of edges (or network layers). Here we show that such multiplex networks are generically more vulnerable to global cascades than simplex networks. We generalize the threshold cascade model [Watts, Proc. Natl. Acad. Sei. USA 99, 5766 (2002)] to multiplex networks, in which a node activates if a sufficiently large fraction of neighbors in any layer are active. We show that both combining layers (i.e., realizing other interactions play a role) and splitting a network into layers (i.e., recognizing distinct kinds of interactions) facilitate cascades. Notably, layers unsusceptible to global cascades can cooperatively achieve them if coupled. On one hand, this suggests fundamental limitations on predicting cascades without full knowledge of a system's multiplexity; on the other hand, it offers feasible means to control cascades by introducing or removing sparse layers in an existing network. [ABSTRACT FROM AUTHOR]

The propagation of elastic waves in a box under direct shear, filled with glass beads and being sheared at constant rates, is studied experimentally and theoretically. The respective velocities are shown to be essentially unchanged from that in a static granular system under the same pressure and shear stress but without a shear band. Influence of shear band on sound behaviors are also briefly discussed. [ABSTRACT FROM AUTHOR]

Diffusive properties of a monodisperse system of interacting particles confined to a quasi-one-dimensional channel are studied using molecular dynamics simulations. We calculate numerically the mean-squared displacement (MSD) and investigate the influence of the width of the channel (or the strength of the confinement potential) on diffusion in finite-size channels of different shapes (i.e., straight and circular). The transition from single-file diffusion to the two-dimensional diffusion regime is investigated. This transition [regarding the calculation of the scaling exponent (&agr;) of the MSD (&Dgr;x2(t)) ∝ t&agr;] as a function of the width of the channel is shown to change depending on the channel's confinement profile. In particular, the transition can be either smooth (i.e., for a parabolic confinement potential) or rather sharp (i.e., for a hard-wall potential), as distinct from infinite channels where this transition is abrupt. This result can be explained by qualitatively different distributions of the particle density for the different confinement potentials. [ABSTRACT FROM AUTHOR]

In this paper, we study many geometrical properties of contour loops to characterize the morphology of synthetic multifractal rough surfaces, which are generated by multiplicative hierarchical cascading processes. To this end, two different classes of multifractal rough surfaces are numerically simulated. As the first group, singular measure multifractal rough surfaces are generated by using the p model. The smoothened multifractal rough surface then is simulated by convolving the first group with a so-called Hurst exponent, H*. The generalized multifractal dimension of isoheight lines (contours), D(q), correlation exponent of contours, xl, cumulative distributions of areas, &xgr;, and perimeters, n, are calculated for both synthetic multifractal rough surfaces. Our results show that for both mentioned classes, hyperscaling relations for contour loops are the same as that of monofractal systems. In contrast to singular measure multifractal rough surfaces, H* plays a leading role in smoothened multifractal rough surfaces. All computed geometrical exponents for the first class depend not only on its Hurst exponent but also on the set of p values. But in spite of multifractal nature of smoothened surfaces (second class), the corresponding geometrical exponents are controlled by H*, the same as what happens for monofractal rough surfaces. [ABSTRACT FROM AUTHOR]

Diffusion coefficients are intrinsically random in subdiffusion attributable to power-law trapping. Using deterministic biased and unbiased diffusion models, we investigate the Einstein relation for single trajectories in subdiffusion. The difference in the generalized Lyapunov exponent between biased and unbiased deterministic diffusions is related to the velocity under a bias. By Hopfs ergodic theorem, the ratios between the velocities and the Lyapunov exponents for single trajectories converge to a universal constant, which is proportional to the strength of the bias. Based on a certain transport coefficient obtained from a single trajectory, we provide a relation for the transport coefficients divided by the Lyapunov exponent and generalize the Einstein relation for single trajectories. [ABSTRACT FROM AUTHOR]

We study theoretically the effect of the variation of the elastic constants on the stability of cholesteric blue phases (BPs). We demonstrate that the stability of BPs is greatly enhanced when the bend elastic constant Ä33 is smaller, in agreement with recent experimental findings. Larger splay (K11) and twist (K22) elastic constants also stabilize BPs. Our argument sheds light on the subtle effect of the variation of the elastic constants on the stability of BPs, and could guide the synthesis of liquid crystals with a wider temperature range of stable BPs. [ABSTRACT FROM AUTHOR]

This paper investigates the dynamics of the flapping regime of a filament placed in a two-dimensional soap-film flow for different filament lengths and flow speeds. It was found that the onset of flapping is quasiperiodic, with the main flapping amplitude and frequency modulated by low-amplitude, low-frequency oscillation. At higher flow velocities, the oscillation becomes chaotic. The transition to chaos occurs via the quasiperiodic route to chaos. A new bistability phenomenon was discovered in which the system alternates between the stretched-straight and oscillatory states, which is here referred to as "switching oscillation." Unlike some previously reported forms of bistability, in this case the system alternates between the two states continuously, without any external perturbation. [ABSTRACT FROM AUTHOR]

A coupling phase is deemed to be crucial in stabilizing behavior in nonlinear systems. In this paper, we study how the coupling phase influences the delay-induced oscillation death (OD) in coupled oscillators. The OD boundaries are identified analytically even in the presence of the coupling phase. We find that OD only occurs for a coupling phase belonging to a certain interval. The optimal coupling phase, under which the largest OD island forms, is characterized well by a power law scaling with respect to the frequency. The coupling phase turns out to be a key parameter that determines a delay-induced OD. Furthermore, the controlling role of the coupling phase generally is proved to hold fairly for networked delay-coupled oscillators. [ABSTRACT FROM AUTHOR]

We investigate the existence and stability of traveling-wave solutions in a continuum field of nonlocally coupled identical phase oscillators with distance-dependent propagation delays. A comprehensive stability diagram in the parametric space of the system is presented that shows a rich structure of multistable regions and illuminates the relative influences of time delay, the nonlocality parameter and the intrinsic oscillator frequency on the dynamics of these states. A decrease in the intrinsic oscillator frequency leads to a break-up of the stability domains of the traveling waves into disconnected regions in the parametric space. These regions exhibit a tongue structure for high connectivity, whereas they submerge into the stable region of the synchronous state for low connectivity. One finding is the existence of forbidden regions in the parametric space where no phase-locked solutions are possible. We also discover a new class of nonstationary breather states for this model system that are characterized by periodic oscillations of the complex order parameter. [ABSTRACT FROM AUTHOR]

Robustness of two coupled networks systems has been studied separately only for dependency coupling [Buldyrev et al., Nature (London) 464, 1025 (2010)] and only for connectivity coupling [Leicht and D'Souza, e-print arXiv:0907.0894], Here we study, using a percolation approach, a more realistic coupled networks system where both interdependent and interconnected links exist. We find rich and unusual phase-transition phenomena including hybrid transition of mixed first and second order, i.e., discontinuities like in a first-order transition of the giant component followed by a continuous decrease to zero like in a second-order transition. Moreover, we find unusual discontinuous changes from second-order to first-order transition as a function of the dependency coupling between the two networks. [ABSTRACT FROM AUTHOR]

We investigate the aggregation kinetics of a simulated telechelic polymer gel. In the hybrid molecular dynamics (MD)/Monte Carlo (MC) algorithm, aggregates of associating end groups form and break according to MC rules, while the position of the polymers in space is dictated by MD. As a result, the aggregate sizes change over time. In order to describe this aggregation process, we employ master equations. They define changes in the number of aggregates of a certain size in terms of reaction rates. These reaction rates indicate the likelihood that two aggregates combine to form a large one, or that a large aggregate splits into two smaller parts. The reaction rates are obtained from the simulations for a range of temperatures. Our results indicate that the rates are not only temperature dependent, but also a function of the sizes of the aggregates involved in the reaction. Using the measured rates, solutions to the master equations are shown to be stable and in agreement with the aggregate size distribution, as obtained directly from simulation data. Furthermore, we show how temperature-induced variations in these rates give rise to the observed changes in the aggregate distribution that characterizes the sol-gel transition. [ABSTRACT FROM AUTHOR]

The nonlinear dynamics of a nonisothermal thin liquid film on a horizontal cylindrical surface subjected to axial harmonic vibration is investigated in this paper. It is found that the capillary instability of an axisymmetric film on a still cylinder can be saturated via the Marangoni effect by heating the film at the gas side, similar to the saturation of the Rayleigh-Taylor instability in the planar case. If the capillary instability is not saturated by the Marangoni effect, the combined capillary-Marangoni instability can be saturated by the harmonic axial forcing provided that the forcing amplitude exceeds a certain critical value depending on the rest of parameters. The critical amplitude is demonstrated to increase with the Marangoni number when other parameters are fixed. The critical forcing amplitude at a given forcing frequency is shown to be well approximated by a piecewise linear function of the Marangoni number. [ABSTRACT FROM AUTHOR]

When a voltage is applied across a thin layer of cholesteric liquid crystal, fingers of cholesteric alignment can form and propagate in the layer. In computer simulation, based on experimental laboratory results, we demonstrate that these cholesteric fingers can solve selected problems of computational geometry, logic, and arithmetics. We show that branching fingers approximate a planar Voronoi diagram, and nonbranching fingers produce a convex subdivision of concave polygons. We also provide a detailed blueprint and simulation of a one-bit half-adder functioning on the principles of collision-based computing, where the implementation is via collision of liquid crystal fingers with obstacles and other fingers. [ABSTRACT FROM AUTHOR]

A three-dimensional symmetric binary fluid is studied, as a function of temperature, in the two-phase (liquid-liquid) coexistence region via Monte Carlo simulations. Particular focus has been in the understanding of curvature-dependent interfacial tension, which is observed to vary as &sgr;(R) = &sgr;(∞)/[1 +2(£/R)2|, implying that a Tolman length is zero in the limit R → &∞;. The length ℓ is found to have a critical divergence the same as the correlation length, but its amplitude is significantly larger (ℓ ⋍ 4&xgr;). Our findings hence imply that the barrier against homogeneous nucleation is significantly reduced (in comparison with the classical nucleation theory) in the critical region. We also report results for the critical behavior of the flat interfacial tension &sgr;(∞) and the concentration susceptibility, as well as the amplitude ratios involving these thermodynamic quantities. Noting that the interatomic potential in our model is described by the Lennard-Jones form that decays faster that 1/r3, all of our results for critical phenomena are expectedly consistent with the Ising universality class of three spatial dimensions. [ABSTRACT FROM AUTHOR]

We consider a nanodimer in solution with asymmetric thermal properties that shows self-propelled motion. One monomer of the nanodimer can be heated to a fixed temperature producing a radially symmetric temperature gradient. The thermophoretic properties of the second monomer produce then a propulsion against or toward the heated particle, such that the nanodimer becomes a puller or pusher nanoswimmer. We combine our simulation measurements with a theoretical analysis that satisfactorily characterizes the self-propelled velocity with the temperature gradient, and the thermophoretic properties of the bead. [ABSTRACT FROM AUTHOR]

We reexamine the convolution approximation commonly used in the mode-coupling theory (MCT) of nonergodic states of classical fluids. This approximation concerns the static correlation functions used as input in the MCT treatment of the dynamics. Besides the hard-sphere model, we consider interaction potentials that present a short-range tail, either attractive or repulsive, beyond the hard core. By using accurate static correlation functions obtained from the fundamental measures functional for hard spheres, we show that the role of three-body direct correlations can be more significant than what is inferred from previous simple ansatzs for pure hard spheres. This may in particular impact the location of the glass transition line and the nonergodicity parameter. [ABSTRACT FROM AUTHOR]

The elastocapillary instability of a flexible plate plunged in a liquid bath is analyzed theoretically. We show that the plate can bend due to two separate destabilizing mechanisms, when the liquid is partially wetting the solid. For contact angles &thgr;e > π/2, the capillary forces acting tangential to the surface are compressing the plate and can induce a classical buckling instability. However, a second mechanism appears due to capillary forces normal to surface. These induce a destabilizing torque that tends to bend the plate for any value of the contact angle &thgr;e > 0. We denote these mechanisms as "buckling" and "bending," respectively and identify the two corresponding dimensionless parameters that govern the elastocapillary stability. The onset of instability is determined analytically and the different bifurcation scenarios are worked out for experimentally relevant conditions. [ABSTRACT FROM AUTHOR]

We report a detailed study, using state-of-the-art simulation and theoretical methods, of the effective (depletion) potential between a pair of big hard spheres immersed in a reservoir of much smaller hard spheres, the size disparity being measured by the ratio of diameters q ≡ &sgr;s/&sgr;b. Small particles are treated grand canonically. their influence being parameterized in terms of their packing fraction in the reservoir &eegr;rs. Two Monte Carlo simulation schemes-the geometrical cluster algorithm, and staged particle insertion-are deployed to obtain accurate depletion potentials for a number of combinations of q ≼ 0.1 and. After applying corrections for simulation finite-size effects, the depletion potentials are compared with the prediction of new density functional theory (DFT) calculations based on the insertion trick using the Rosenfeld functional and several subsequent modifications. While agreement between the DFT and simulation is generally good, significant discrepancies are evident at the largest reservoir packing fraction accessible to our simulation methods, namely, &eegr;rs = 0.35. These discrepancies are, however, small compared to those between simulation and the much poorer predictions of the Derjaguin approximation at this &eegr;rs. The recently proposed morphometric approximation performs better than Derjaguin but is somewhat poorer than DFT for the size ratios and small-sphere packing fractions that we consider. The effective potentials from simulation, DFT, and the morphometric approximation were used to compute the second virial coefficient B2 as a function of &eegr;rs. Comparison of the results enables an assessment of the extent to which DFT can be expected to correctly predict the propensity toward fluid-fluid phase separation in additive binary hard-sphere mixtures with q ≼ 0.1. In all, the new simulation results provide a fully quantitative benchmark for assessing the relative accuracy of theoretical approaches for calculating depletion potentials in highly size-asymmetric mixtures. [ABSTRACT FROM AUTHOR]

Long-lived quasistationary states, associated with stationary stable solutions of the Vlasov equation, are found in systems with long-range interactions. Studies of the relaxation time in a model of N globally coupled particles moving on a ring, the Hamiltonian mean-field model (HMF), have shown that it diverges as NY for large N, with &ggr; -≃1-7 for some initial conditions with homogeneously distributed particles. We propose a method for identifying exact inhomogeneous steady states in the thermodynamic limit, based on analyzing models of uncoupled particles moving in an external field. For the HMF model, we show numerically that the relaxation time of these states diverges with N&ggr; with the exponent &ggr; ≃ 1. The method, applicable to other models with globally coupled particles, also allows an exact evaluation of the stability I imit of homogeneous steady states. In some cases, it provides a good approximation for the correspondence between the initial condition and the final steady state. [ABSTRACT FROM AUTHOR]

Salt transport in bulk electrolytes is limited by diffusion and advection, but in microstructures with charged surfaces (e.g., microfluidic devices, porous media, soils, or biological tissues) surface conduction and electro-osmotic flow also contribute to ionic fluxes. For small applied voltages, these effects lead to well known linear electrokinetic phenomena. In this paper, we predict some surprising nonlinear dynamics that can result from the competition between bulk and interfacial transport at higher voltages. When counterions are selectively removed by a membrane or electrode, a "deionization shock" can propagate through the microstructure, leaving in its wake an ultrapure solution, nearly devoid of coions and colloidal impurities. We elucidate the basic physics of deionization shocks and develop a mathematical theory of their existence, structure, and stability, allowing for slow variations in surface charge or channel geometry. Via asymptotic approximations and similarity solutions, we show that deionization shocks accelerate and sharpen in narrowing channels, while they decelerate and weaken, and sometimes disappear, in widening channels. These phenomena may find applications in separations (deionization, decontamination, biological assays) and energy storage (batteries, supercapacitors) involving electrolytes in microstructures. [ABSTRACT FROM AUTHOR]

We report on the dynamics of a model frictional system submitted to minute external perturbations. The system consists of a chain of sliders connected through elastic springs that rest on an incline. By introducing cyclic expansions and contractions of the springs we observe a reptation of the chain. We account for the average reptation velocity theoretically. The velocity of small systems exhibits a series of plateaus as a function of the incline angle. Due to elastic effects, there exists a critical amplitude below which the reptation is expected to cease. However, rather than a full stop of the creep, we observe in numerical simulations a transition between a continuous-creep and an irregular-creep regime when the critical amplitude is approached. The latter transition is reminiscent of the transition between the continuous and the irregular compaction of granular matter submitted to periodic temperature changes. [ABSTRACT FROM AUTHOR]

This paper investigates the lift force exerted on an elliptical obstacle immersed in a granular flow through analytical calculations and computer simulations. The results are shown as a function of the obstacle si/e. orientation with respect to the flow direction (tilt angle), the restitution coefficient and ellipse eccentricity. The theoretical argument, based on the force exerted on the obstacle due to inelastic, frictionless collisions of a very dilute flow, captures the qualitative features of the lift, but fails to reproduce the data quantitatively. The reason behind this disagreement is that the dilute flow assumption on which this argument is built breaks down as a granular shock wave forms in front of the obstacle. More specifically, the shock wave changes the grains impact velocity at the obstacle, decreasing the overall net lift obtained from a very dilute flow. [ABSTRACT FROM AUTHOR]

We determine the complexity of several constraint satisfaction problems using the quantum adiabatic algorithm in its simplest implementation. We do so by studying the size dependence of the gap to the first excited state of "typical" instances. We find that, at large sizes N, the complexity increases exponentially for all models that we study. We also compare our results against the complexity of the analogous classical algorithm WalkSAT and show that the harder the problem is for the classical algorithm, the harder it is also for the quantum adiabatic algorithm. [ABSTRACT FROM AUTHOR]

Game-theoretical models where the rules of the game and the interaction structure both coevolve with the game dynamics-multiadaptive games-capture very flexible situations where cooperation among selfish agents can emerge. In this work, we will discuss a multiadaptive model presented in a recent Letter [Phys. Rev. Lett. 106, 028702 (2011)] as well as generalizations of it. The model captures a nonequilibrium situation where social unrest increases the incentive to cooperate and, simultaneously, agents are partly free to influence with whom they interact. First, we investigate the details of how feedback from the behavior of agents determines the emergence of cooperation and hierarchical contact structures. We also study the stability of the system to different types of noise, and find that different regions of parameter space show very different response. Some types of noise can destroy an all-cooperator (C) state. If, on the other hand, hubs are stable, then so is the all-C state. Finally, we investigate the dependence of the ratio between the time scales of strategy updates and the evolution of the interaction structure. We find that a comparatively fast strategy dynamics is a prerequisite for the emergence of cooperation. [ABSTRACT FROM AUTHOR]

For a broad class of random walks with anisotropic scattering kernels and absorption, we derive explicit formulas that allow expressing the moments of the collision number nv performed in a volume V as a function of the particle equilibrium distribution. Our results apply to arbitrary domains V and boundary conditions, and allow assessing the hitting statistics for systems where the typical displacements are comparable to the domain size, so that the diffusion limit is possibly not attained. An example is discussed for one-dimensional random flights with exponential displacements, where analytical calculations can be carried out. [ABSTRACT FROM AUTHOR]

We consider continuous-time random walks (CTRWs) in which the walkers have a finite probability to alter the waiting-time and/or step-length transport properties of their environment, resulting in possibly transient anomalous diffusion. We refer to these CTRWs as transmogrifying continuous-time random walks (TCTRWs) to emphasize that they change the form of the transport properties of their environment, and in a possibly strange way. The particular case in which the CTRW waiting-time density has a finite probability to be permanently altered at a given site, following a visitation by a walker, is considered in detail. Master equations for the probability density function of transmogrifying random walkers are derived, and results are compared with Monte Carlo simulations. An interesting finding is that TCTRWs can generate transient subdiffusion or transient superdiffusion without invoking truncated or tempered power law densities for either the waiting times or the step lengths. The transient subdiffusion or transient superdiffusion arises in TCTRWs with Gaussian step-length densities and exponential waiting-time densities when the altered average waiting time is greater than or less than, respectively, the original average waiting time. [ABSTRACT FROM AUTHOR]

Existence of a nonequilibrium tricritical point has been revealed in the Hamiltonian mean-field model by a nonequilibrium statistical mechanics. This statistical mechanics gives a distribution function containing unknown parameters, and the parameters are determined by solving simultaneous equations depending on a given initial state. Due to difficulty in solving these equations, pointwise numerical detection of the tricritical point has been unavoidable on a plane characterizing a family of initial states. In order to look into the tricritical point, we expand the simultaneous equations with respect to the order parameter and reduce them to one algebraic equation. The tricritical point is precisely identified by analyzing coefficients of the reduced equation. Reentrance to an ordered phase in a high-energy region is revisited around the obtained tricritical point. [ABSTRACT FROM AUTHOR]

The Boltzmann equation subject to a general boundary condition is expanded in a power series with respect to a thermodynamic force disturbing a gaseous system. Recurrence relations between the terms of the expansion are obtained using the main properties of the collision integral and of the gas-surface interaction kernel. The reciprocal relation for nonlinear irreversible phenomena, i.e. a relation between the terms of different orders, is obtained. The relations can be used to estimate the range of applicability of linearized solutions and to predict nonlinear phenomena in gaseous systems. [ABSTRACT FROM AUTHOR]

In an asymmetric Szilard engine model of Maxwell's demon, we show the equivalence between information theoretical and thermodynamic entropies when the demon erases information optimally. The work gain by the engine can be exactly canceled out by the work necessary to reset the demon's memory after optimal data compression in the manner of Shannon before the erasure. [ABSTRACT FROM AUTHOR]

Motivated by the phenomenology of transport through the Golgi apparatus of cells, we study a multispecies model with boundary injection of one species of particle, interconversion between the different species of particle, and driven diffusive movement of particles through the system by chipping of a single particle from a site. The model is analyzed in one dimension using equations for particle currents. It is found that, depending on the rates of various processes and the asymmetry in the hopping, the system may exist either in a steady phase, in which the average mass at each site attains a time-independent value, or in a "growing" phase, in which the total mass grows indefinitely in time, even in a finite system. The growing phases have interesting spatial structure. In particular, we find phases in which some spatial regions of the system have a constant average mass, while other regions show unbounded growth. [ABSTRACT FROM AUTHOR]