Chemical reactions subjected to time-varying external forces cannot generally be described through a fixed bottleneck near the transition-state barrier or dividing surface. A naive dividing surface attached to the instantaneous, but moving, barrier top also fails to be recrossing-free. We construct a moving dividing surface in phase space over a transition-state trajectory. This surface is recrossing-free for both Hamiltonian and dissipative dynamics. This is confirmed even for strongly anharmonic barriers using simulation. The power of transition-state theory is thereby applicable to chemical reactions and other activated processes even when the bottlenecks are time dependent and move across space.

In this paper, we investigate the evolution of cooperation in the spatial prisoner's dilemma game by incorporating partner choice into the framework of evolutionary game theory. Our research shows that the introduction of partner choice can notably promote the cooperative behavior in the prisoner's dilemma game. All the players are more likely to play the game with cooperators, which makes it easier for cooperators to form alliances. In particular, the system will be dominated completely by cooperators in a comfortable environment (i.e., with lower survival cost) because the cooperators can get more benefits than the defectors in this case due to their good reputation. In addition, we have found that the sustenance of cooperators improves notably as well in the snowdrift game and the stag-hunt game due to this introduction.

We have constructed a macroscopic driven system of chaotic Faraday waves whose statistical mechanics, we find, are surprisingly simple, mimicking those of a thermal gas. We use real-time tracking of a single floating probe, energy equipartition, and the Stokes-Einstein relation to define and measure a pseudotemperature and diffusion constant and then self-consistently determine a coefficient of viscous friction for a test particle in this pseudothermal gas. Because of its simplicity, this system can serve as a model for direct experimental investigation of nonequilibrium statistical mechanics, much as the ideal gas epitomizes equilibrium statistical mechanics.

We report a study of nonequilibrium relaxation in a two-dimensional random field Ising model at a nonzero temperature. We attempt to observe the coarsening from a different perspective with a particular focus on three dynamical quantities that characterize the kinetic coarsening. We provide a simple generalized scaling relation of coarsening supported by numerical results. The excellent data collapse of the dynamical quantities justifies our proposition. The scaling relation corroborates the recent observation that the average linear domain size satisfies different scaling behavior in different time regimes.

In olfactory search an immobile target emits chemical molecules at constant rate. The molecules are transported by the medium, which is assumed to be turbulent. Considering a searcher able to detect such chemical signals and whose motion follows the infotaxis strategy, we study the statistics of the first-passage time to the target when the searcher moves on a finite two-dimensional lattice of different geometries. Far from the target, where the concentration of chemicals is low, the direction of the searcher's first movement is determined by the geometry of the domain and the topology of the lattice, inducing strong fluctuations on the average search time with respect to the initial position of the searcher. The domain is partitioned in well-defined regions characterized by the direction of the first movement. If the search starts over the interface between two different regions, large fluctuations in the search time are observed.

When pulling a particle in a driven granular fluid with constant force $F_{\text{ex}}$, the probe particle approaches a steady-state average velocity $v$. This velocity and the corresponding friction coefficient of the probe $\zeta =F_{\text{ex}}/v$ are obtained within a schematic model of mode-coupling theory and compared to results from event-driven simulations. For small and moderate drag forces, the model describes the simulation results successfully for both the linear as well as the nonlinear region: The linear response regime (constant friction) for small drag forces is followed by shear thinning (decreasing friction) for moderate forces. For large forces, the model demonstrates a subsequent increasing friction in qualitative agreement with the data. The square-root increase of the friction with force found in [Fiege et al., Granul. Matter 14, 247 (2012)] is explained by a simple kinetic theory.

The elastic deformation of a soft solid induced by capillary forces crucially relies on the excess stress inside the solid-liquid interface. While for a liquid-liquid interface this “surface stress” is strictly identical to the “surface free energy,” the thermodynamic Shuttleworth equation implies that this is no longer the case when one of the phases is elastic. Here we develop a microscopic model that incorporates enthalpic interactions and entropic elasticity, based on which we explicitly compute as the surface stress and surface free energy. It is found that the compressibility of the interfacial region, through the Poisson ratio near the interface, determines the difference between surface stress and surface energy. We highlight the consequence of this finding by comparing with recent experiments and simulations on partially wetted soft substrates.

We determine the nematic-isotropic coexistence curve terminating at the critical point in a temperature-external field phase diagram for nematic liquid crystals with positive diamagnetic anisotropy, where the molecules are either perfectly uniaxial or biaxial using computer simulation of a lattice model. The coexistence curve is much steeper than that predicted by the standard Landau–de Gennes and Maier-Saupe mean-field theories. For the uniaxial system the critical magnetic field is estimated to be one order of magnitude lower than the mean-field estimate but of the same order of magnitude as the experimental measurement. Our study shows that molecular biaxiality could reduce the critical field strength significantly.

Using methods from random matrix theory researchers have recently calculated the full spectra of random networks with arbitrary degrees and with community structure. Both reveal interesting spectral features, including deviations from the Wigner semicircle distribution and phase transitions in the spectra of community structured networks. In this paper we generalize both calculations, giving a prescription for calculating the spectrum of a network with both community structure and an arbitrary degree distribution. In general the spectrum has two parts, a continuous spectral band, which can depart strongly from the classic semicircle form, and a set of outlying eigenvalues that indicate the presence of communities.

A numerical simulation of evolution of an identical interacting streamers array in semiconductors has been performed using the diffusion-drift approximation and taking into account the impact and tunnel ionization. It has been assumed that the external electric field $E_0$ is static and uniform, the background electrons and holes are absent, the initial avalanches start simultaneously from the nodes of the plane hexagonal lattice, which is perpendicular to the external field, but the avalanches and streamers are axially symmetric within a cylinder of radius $R$. It has been shown that under certain conditions, the interaction between the streamers leads finally either to the formation of two types of stationary ionization waves with corrugated front or to a stationary plane ionization wave. A diagram of different steady states of this type of waves in the plane of parameter $E_0,R$ has been presented, and a qualitative explanation of the plane partition into four different regions has been given. Characteristics of corrugated waves have been studied in detail and discussed in the region of $R$ and $E_0$ large values, in which the maximum field strength at the front is large enough for the tunnel ionization implementation. It has been shown that corrugated waves ionize semiconductors more efficiently than flat ones, especially in relatively weak external fields.

As a generalization of Newton's two body problem, we explore the dynamics of two massive line segments interacting gravitationally. The extension of each line segment or slash (/) provides extra degrees of freedom that enable the interplay between rotation and revolution in an especially simple example. This slash-slash (//) body problem can thereby elucidate the dynamics of nonspherical space structures, from asteroids to space stations. Fortunately, as we show, Newton's laws imply exact algebraic expressions for the force and torque between the slashes, and this greatly facilitates analysis. The diverse dynamics include a stable synchronous orbit, families of unstable periodic orbits, generic chaotic orbits, and spin-orbit coupling that can unbind the slashes. In particular, retrograde orbits where the slashes spin opposite to their orbits are stable, with regular dynamics and smooth parameter spaces, while prograde orbits are unstable, with chaotic dynamics and fractal parameter spaces.

We solve for the motion of charged particles in a helical time-periodic ABC (Arnold-Beltrami-Childress) magnetic field. The magnetic field lines of a stationary ABC field with coefficients $A=B=C=1$ are chaotic, and we show that the motion of a charged particle in such a field is also chaotic at late times with positive Lyapunov exponent. We further show that in time-periodic ABC fields, the kinetic energy of a charged particle can increase indefinitely with time. At late times the mean kinetic energy grows as a power law in time with an exponent that approaches unity. For an initial distribution of particles, whose kinetic energy is uniformly distributed within some interval, the probability density function of kinetic energy is, at late times, close to a Gaussian but with steeper tails.

To understand transition to turbulence via 3D disturbance growth, we report here results obtained from the solution of Navier-Stokes equation (NSE) to reproduce experimental results obtained by minimizing background disturbances and imposing deterministic excitation inside the shear layer. A similar approach was adopted in Sengupta and Bhaumik [Phys. Rev. Lett. 107, 154501 (2011)], where a route of transition from receptivity to fully developed turbulent stage was explained for 2D flow in terms of the spatio-temporal wave-front (STWF). The STWF was identified as the unit process of 2D turbulence creation for low amplitude wall excitation. Theoretical prediction of STWF for boundary layer was established earlier in Sengupta, Rao, and Venkatasubbaiah [Phys. Rev. Lett. 96, 224504 (2006)] from the Orr-Sommerfeld equation as due to spatiotemporal instability. Here, the same unit process of the STWF during transition is shown to be present for 3D disturbance field from the solution of governing NSE.

The influence of the evolution of numerous voids on macroscopic properties of materials is a multiscale problem that challenges computational research. A shock-wave compression model for brittle material, which can obtain both microscopic evolution and macroscopic shock properties, was developed using discrete element methods (lattice model). Using a model interaction-parameter-mapping procedure, qualitative features, as well as trends in the calculated shock-wave profiles, are shown to agree with experimental results. The shock wave splits into an elastic wave and a deformation wave in porous brittle materials, indicating significant shock plasticity. Void collapses in the deformation wave were the natural reason for volume shrinkage and deformation. However, media slippage and rotation deformations indicated by complex vortex patterns composed of relative velocity vectors were also confirmed as an important source of shock plasticity. With increasing pressure, the contribution from slippage deformation to the final plastic strain increased. Porosity was found to determine the amplitude of the elastic wave; porosity and shock stress together determine propagation speed of the deformation wave, as well as stress and strain on the final equilibrium state. Thus, shock behaviors of porous brittle material can be systematically designed for specific applications.

The mechanical response of granular beds under applied stresses is often characterized by repeated cycles of stick-slip. Using the discrete element method, we examine stick-slip from a concentrated force loading—imposed by a single grain that is drawn through a densely packed, periodic granular bed via a stiff virtual spring. Force chains continually form and collapse ahead of the intruder grain. A comprehensive characterization of the birth-death evolution of these load-bearing structures, along with their surrounding contact cycles, reveals a well-defined shear zone of around eight particle diameters from the intruder, encapsulating: (i) long force chains that form buttresses with the fixed bottom wall for support, (ii) a region where the collapse of the most stable, persistent three-cycles preferentially occur to the point where they are essentially depleted by the end of the first cycle of stick-slip, and (iii) an inner core where force chain buckling events concentrate. Dilatancy is greatest in this inner core, and in the region next to the free surface. During slip, secondary force chains briefly form behind the intruder: these transient force chains, most of which comprise only 3 particles, form in the direction that is roughly perpendicular to the intruder motion.

We investigate the local fluctuations in dense granular media subjected to athermal, quasistatic shearing, based on three-dimensional discrete element method simulations. By shearing granular assemblies of different polydispersities under constant-volume constraint, we quantify the characteristics of local structures (in terms of local volume and local anisotropy) and local deformation (using local shear strain and nonaffine displacement). The distribution of the local volume in a granular medium is found unchanged during the entire shearing process, which indicates a constant temperaturelike compactivity for the material. The compactivity is not, however, equilibrated among different particle groups in a polydisperse assembly. The local structures of a disordered granular assembly are inherently anisotropic. The fluctuations in local anisotropy can be well captured by a gamma or mixed-gamma distribution function, which is also unchanged during the shear. The local anisotropic orientation evolves towards the coaxial direction of the stress anisotropy with shear. The deformation characteristics of a jammed granular medium have their origins in the structural amorphousness. The local shear strain field depicts clear shear transformation zones which act as plasticity carriers. The spatial correlation of the local shear strains exhibits a fourfold pattern which is stronger in the stress deviatoric planes than in the stress isotropic plane. The fluctuations of nonaffine displacement suggest an isotropic granular temperature and an isotropic spatial correlation independent of the stress state. Both the local strain and the nonaffine displacement exhibit a power-law decayed distribution with a long-range correlation. We further modify the shear-transformation-zone theory to predict the pressure-dependent constitutive behavior of a sheared granular material and compare its prediction with our simulation data.

The morphology and the microscopic internal dynamics of a bidimensional gel formed by spontaneous aggregation of gold nanoparticles confined at the water surface are investigated by a suite of techniques, including grazing-incidence x-ray photon correlation spectroscopy (GI-XPCS). The range of concentrations studied spans across the percolation transition for the formation of the gel. The dynamical features observed by GI-XPCS are interpreted in view of the results of microscopic imaging; an intrinsic link between the mechanical modulus and internal dynamics is demonstrated for all the concentrations. Our work presents an example of a transition from a stretched to a compressed correlation function actively controlled by quasistatically varying the relevant thermodynamic variable. Moreover, by applying a model proposed some time ago by Duri and Cipelletti [Europhys. Lett. 76, 972 (2006)] we are able to build a master curve for the shape parameter, whose scaling factor allows us to quantify a “long-time displacement length.” This characteristic length is shown to converge, as the concentration is increased, to the “short-time localization length” determined by pseudo-Debye-Waller analysis of the initial contrast. Finally, the intrinsic dynamics of the system is then compared with that induced by means of a delicate mechanical perturbation applied to the interface.

We expand a previous study [Phys. Rev. E 86, 051611 (2012)] on the conditions for occurrence of strong anisotropy in the scaling properties of two-dimensional surfaces displaying generic scale invariance. In that study, a natural scaling ansatz was proposed for strongly anisotropic systems, which arises naturally when analyzing data from, e.g., thin-film production experiments. The ansatz was tested in Gaussian (linear) models of surface dynamics and in nonlinear models, like the Hwa-Kardar (HK) equation [Phys. Rev. Lett. 62, 1813 (1989)], which are susceptible of accurate approximations through the former. In contrast, here we analyze nonlinear equations for which such approximations fail. Working within generically scale-invariant situations, and as representative case studies, we formulate and study a generalization of the HK equation for conserved dynamics and reconsider well-known systems, such as the conserved and the nonconserved anisotropic Kardar-Parisi-Zhang equations. Through the combined use of dynamic renormalization group analysis and direct numerical simulations, we conclude that the occurrence of strong anisotropy in two-dimensional surfaces requires dynamics to be conserved. We find that, moreover, strong anisotropy is not generic in parameter space but requires, rather, specific forms of the terms appearing in the equation of motion, whose justification needs detailed information on the dynamical process that is being modeled in each particular case.