Granular surfaces subjected to repeated passage often spontaneously develop a corrugated pattern. In this study, we scrutinized the growth dynamics of surface corrugation in a self-rotating sand bed that is traced by the edge of a vertically oscillating arm. We found that both the rotation velocity and the frequency of the oscillator strongly affect the occurrence rate of corrugation as well as the time evolution in the corrugation pattern, due to the intermittent collision between the oscillator and sand bed. We also discovered that the growth dynamics involves two distinct collective modes that describe the translational motion of the corrugation pattern on the sand bed.

In this paper, we investigate experimentally the dynamics of particles in dense granular suspensions when both shear and external vibrations are applied. We study in detail how vibrations affect particle reorganization at the local scale and modify the apparent rheology. The nonlocal nature of the rheology when no vibrations are applied is evidenced, in agreement with previous numerical studies from the literature. It is also shown that vibrations induce structural reorganizations, which tend to homogenize the system and cancel the nonlocal properties.

By means of the lubrication approximation we obtain the full family of static bidimensional profiles of a liquid resting on a substrate under partial-wetting conditions imposed by a disjoining-conjoining pressure. We show that for a set of quite general disjoining-conjoining pressure potentials, the free surface can adopt only five nontrivial static patterns; in particular, we find solutions when the height goes to zero which describe satisfactorily the complete free surface for a finite amount of fluid deposited on a substrate. To test the extension of the applicability of our solutions, we compare them with those obtained when the lubrication approximations are not employed and under conditions where the lubrication hypothesis are not strictly valid, and also with axisymmetric solutions. For a given disjoining-conjoining potential, we report a new analytical solution that accounts for all the five possible solutions.

We construct a boundary integral representation for the low-Reynolds-number flow in a channel in the presence of freely suspended particles (or droplets) of arbitrary size and shape. We demonstrate that lubrication theory holds away from the particles at horizontal distances exceeding the channel height and derive a multipole expansion of the flow which is dipolar to the leading approximation. We show that the dipole moment of an arbitrary particle is a weighted integral of the stress and the flow at the particle surface, which can be determined numerically. We introduce the equation of motion that describes hydrodynamic interactions between arbitrary, possibly different, distant particles, with interactions determined by the product of the mobility matrix and the dipole moment. Further, the problem of three identical interacting spheres initially aligned in the streamwise direction is considered and the experimentally observed “pair exchange” phenomenon is derived analytically and confirmed numerically. For nonaligned particles, we demonstrate the formation of a configuration with one particle separating from a stable pair. Our results suggest that in a dilute initially homogenous particulate suspension flowing in a channel the particles will eventually separate into singlets and pairs.

The numerical solution of the general coupled nonlinear Schrödinger equations on unbounded domains is considered by applying the artificial boundary method in this paper. In order to design the local absorbing boundary conditions for the coupled nonlinear Schrödinger equations, we generalize the unified approach previously proposed [J. Zhang et al., Phys. Rev. E 78, 026709 (2008)]. Based on the methodology underlying the unified approach, the original problem is split into two parts, linear and nonlinear terms, and we then achieve a one-way operator to approximate the linear term to make the wave out-going, and finally we combine the one-way operator with the nonlinear term to derive the local absorbing boundary conditions. Then we reduce the original problem into an initial boundary value problem on the bounded domain, which can be solved by the finite difference method. The stability of the reduced problem is also analyzed by introducing some auxiliary variables. Ample numerical examples are presented to verify the accuracy and effectiveness of our proposed method.

The present work shows a method for stable simulations via the lattice Boltzmann (LB) model for electromagnetic waves (EM) transiting homogeneous media. LB models for such media were already presented in the literature, but they suffer from numerical instability when the media transitions are sharp. We use one of these models in the limit of pure vacuum derived from Liu and Yan [Appl. Math. Model. 38, 1710 (2014)] and apply an extension that treats the effects of polarization and magnetization separately. We show simulations of simple examples in which EM waves travel into media to quantify error scaling, stability, accuracy, and time scaling. For conductive media, we use the Strang splitting and check the simulations accuracy at the example of the skin effect. Like pure EM propagation, the error for the static limits, which are constructed with a current density added in a first-order scheme, can be less than $1\%$. The presented method is an easily implemented alternative for the stabilization of simulation for EM waves propagating in spatially complex structured media properties and arbitrary transitions.

Permutation entropy (PE) is a statistic used widely for the detection of structure within a time series. Embedding delay times at which the PE is reduced are characteristic timescales for which such structure exists. Here, a generalized scheme is investigated where embedding delays are represented by vectors rather than scalars, permitting PE to be calculated over a $(D-1)$-dimensional space, where $D$ is the embedding dimension. This scheme is applied to numerically generated noise, sine wave and logistic map series, and experimental data sets taken from a vertical-cavity surface emitting laser exhibiting temporally localized pulse structures within the round-trip time of the laser cavity. Results are visualized as PE maps as a function of embedding delay, with low PE values indicating combinations of embedding delays where correlation structure is present. It is demonstrated that vector embedding delays enable identification of structure that is ambiguous or masked, when the embedding delay is constrained to scalar form.

We propose a computationally effective approach that builds on Landau mean-field theory in combination with modern nonequilibrium statistical mechanics to model and interpret protein dynamics and structure formation in small- to wide-angle x-ray scattering (S/WAXS) experiments. We develop the methodology by analyzing experimental data in the case of Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data qualitatively with a good precision and over an extended temperature range. We explain experimental observations in terms of protein phase structure, and we make predictions for future experiments and for how to analyze data at different ambient temperature values. We conclude that the approach we propose has the potential to become a highly accurate, computationally effective, and predictive tool for analyzing S/WAXS data. For this, we compare our results with those obtained previously in an all-atom molecular dynamics simulation.

The evaporation of droplets can give rise to a wide range of interesting phenomena in which the dynamics of the evaporation are crucial. In this work, we find simple scaling laws for the evaporation dynamics of axisymmetric droplets pinned on millimeter-sized pillars. Different laws are found depending on whether evaporation is limited by the diffusion of vapor molecules or by the transfer rate across the liquid-vapor interface. For the diffusion-limited regime, we find that a mass-loss rate equal to 3/7 of that of a free-standing evaporating droplet brings a good balance between simplicity and physical correctness. We also find a scaling law for the evaporation of multicomponent solutions. The scaling laws found are validated against experiments of the evaporation of droplets of (1) water, (2) blood plasma, and (3) a mixture of water and polyethylene glycol, pinned on acrylic pillars of different diameters. These results shed light on the macroscopic dynamics of evaporation on pillars as a first step towards the understanding of other complex phenomena that may be taking place during the evaporation process, such as particle transport and chemical reactions.

We develop a gray-scale lattice Boltzmann (LB) model to study fluid flow combined with heat transfer for flow through porous media where voxels may be partially solid (or void). Heat transfer in rocks may lead to deformation, which in turn can modulate the fluid flow and so has significant contribution to rock permeability. The LB temperature field is compared to a finite difference solution of the continuum partial differential equations for fluid flow in a channel. Excellent quantitative agreement is found for both Poiseuille channel flow and Brinkman flow. The LB model is then applied to sample porous media such as packed beds and also more realistic sandstone rock sample, and both the convective and diffusive regimes are recovered when varying the thermal diffusivity. It is found that while the rock permeability can be comparatively small (order milli-Darcy), the temperature field can show significant variation depending on the thermal convection of the fluid. This LB method has significant advantages over other numerical methods such as finite and boundary element methods in dealing with coupled fluid flow and heat transfer in rocks which have irregular and nonsmooth pore spaces.

In the context of a discrete analog of the classical Heisenberg model, we investigate the critical behavior of the icosahedron model, where the interaction energy is defined as the inner product of neighboring vector spins of unit length pointing to the vertices of the icosahedron. The effective correlation length and magnetization of the model are calculated by means of the corner-transfer-matrix renormalization group (CTMRG) method. A scaling analysis with respect to the cutoff dimension $m$ in CTMRG reveals a second-order phase transition characterized by the exponents $\nu =1.62\pm 0.02$ and $\beta =0.12\pm 0.01$. We also extract the central charge from the classical analog of entanglement entropy as $c=1.90\pm 0.02$, which cannot be explained by the minimal series of conformal field theory.

This work deals with thermostated fluctuating systems subjected to driven transformations of the internal energetics. The main focus is on generally multidimensional systems with continuous configurational degrees of freedom over which overdamped Markovian fluctuations take place (diffusive regime of the motion). Mutual bounds are established between the average energy dissipation, the deviation between nonequilibrium probability density and underlying equilibrium distribution due to the system's lag, and the statistical properties of the components of the directed flow induced by the transformation itself. The directed flow is here expressed in terms of time-dependent “drift velocity” associated with the probability current in a advection-like formulation of the nonstationary Fokker-Planck equation. Consideration of the drift makes that the bounds achieved here extend the inequality derived by Vaikuntanathan and Jarzynski [Europhys. Lett. 87, 60005 (2009)] involving only dissipation and lag. The key relations are then specified for the so-called stochastic pumps, i.e., systems that reach a periodic steady state in response of cyclic transformations and that are prototypes of nonautonomous dissipative converters of input energy into directed motion; a one-dimensional case model is adopted to illustrate the main features. Complementary results concerning bounds between the evolution rates of dissipation and lag, valid for both overdamped and underdamped dynamics, are also presented.

We analyze a family of exact finite energy solutions to Maxwell's equations. These solutions are a subset of the modified-power-spectrum solutions found by Ziolkowski [Phys. Rev. A 39, 2005 (1989)]. There are three characteristic parameters in the solutions: $q_1,q_2$, and $k_0$. $q_1$ and $q_2$ are related to the frequency bandwidth of the solution. In the parameter space of $k_0{q}_1\gg 1$ and $k_0{q}_2\gg 1$, they represent quasimonochromatic continuous wave fields with the main angular frequency $k_{0}c$ and energy localized in the transverse directions. Under the restriction of $q_1\ll q_2$, the beam propagates mainly in the $+z$ direction with velocity $c$ and limited diffraction.

We provide a classification of the possible flows of two-component Bose-Einstein condensates evolving from initially discontinuous profiles. We consider the situation where the dynamics can be reduced to the consideration of a single polarization mode (also denoted as “magnetic excitation”) obeying a system of equations equivalent to the Landau-Lifshitz equation for an easy-plane ferromagnet. We present the full set of one-phase periodic solutions. The corresponding Whitham modulation equations are obtained together with formulas connecting their solutions with the Riemann invariants of the modulation equations. The problem is not genuinely nonlinear, and this results in a non-single-valued mapping of the solutions of the Whitham equations with physical wave patterns as well as the appearance of interesting elements—contact dispersive shock waves—that are absent in more standard, genuinely nonlinear situations. Our analytic results are confirmed by numerical simulations.

We investigate the dynamics of coupled identical chaotic Lorenz oscillators just above the subcritical Hopf bifurcation. In the absence of coupling, the motion is on a strange chaotic attractor and the fixed points of the system are all unstable. With the coupling, the unstable fixed points are converted into chaotic attractors, and the system can exhibit a multiplicity of coexisting attractors. Depending on the strength of the coupling, the motion of the individual oscillators can be synchronized (both in and out of phase) or desynchronized and in addition there can be mixed phases. We find that the basins have a complex structure: the state that is asymptotically reached shows extreme sensitivity to initial conditions. The basins of attraction of these different states are characterized using a variety of measures and depending on the strength of the coupling, they are intermingled or quasiriddled.

We study the effect of common noise on coupled active rotators. While such a noise always facilitates synchrony, coupling may be attractive (synchronizing) or repulsive (desynchronizing). We develop an analytical approach based on a transformation to approximate angle-action variables and averaging over fast rotations. For identical rotators, we describe a transition from full to partial synchrony at a critical value of repulsive coupling. For nonidentical rotators, the most nontrivial effect occurs at moderate repulsive coupling, where a juxtaposition of phase locking with frequency repulsion (anti-entrainment) is observed. We show that the frequency repulsion obeys a nontrivial power law.

Previously we reported [A. Wynveen et al., Phys. Rev. E 89, 022725 (2014)] that requiring that the systems regarded as lifelike be out of chemical equilibrium in a model of abstracted polymers undergoing ligation and scission first introduced by Kauffman [S. A. Kauffman, The Origins of Order (Oxford University Press, New York, 1993), Chap. 7] implied that lifelike systems were most probable when the reaction network was sparse. The model was entirely statistical and took no account of the bond energies or other energetic constraints. Here we report results of an extension of the model to include effects of a finite bonding energy in the model. We studied two conditions: (1) A food set is continuously replenished and the total polymer population is constrained but the system is otherwise isolated and (2) in addition to the constraints in (1) the system is in contact with a finite-temperature heat bath. In each case, detailed balance in the dynamics is guaranteed during the computations by continuous recomputation of a temperature [in case (1)] and of the chemical potential (in both cases) toward which the system is driven by the dynamics. In the isolated case, the probability of reaching a metastable nonequilibrium state in this model depends significantly on the composition of the food set, and the nonequilibrium states satisfying lifelike condition turn out to be at energies and particle numbers consistent with an equilibrium state at high negative temperature. As a function of the sparseness of the reaction network, the lifelike probability is nonmonotonic, as in our previous model, but the maximum probability occurs when the network is less sparse. In the case of contact with a thermal bath at a positive ambient temperature, we identify two types of metastable nonequilibrium states, termed locally and thermally alive, and locally dead and thermally alive, and evaluate their likelihood of appearance, finding maxima at an optimal temperature and an optimal degree of sparseness in the network. We use a Euclidean metric in the space of polymer populations to distinguish these states from one another and from fully equilibrated states. The metric can be used to characterize the degree and type of chemical equilibrium in observed systems, as we illustrate for the proteome of the ribosome.

Nanoscale dynamics with cycles of receptor diffusion and immobilization by cell-external-or-internal factors is a key process in living cell adhesion phenomena at the origin of a plethora of signal transduction pathways. Motivated by modern correlation microscopy approaches, the receptor correlation functions in physical models based on diffusion-influenced reaction is studied. Using analytical and stochastic modeling, this paper focuses on the hybrid regime where diffusion and reaction are not truly separable. The time receptor autocorrelation functions are shown to be indexed by different time scales and their asymptotic expansions are given. Stochastic simulations show that this analysis can be extended to situations with a small number of molecules. It is also demonstrated that this analysis applies when receptor immobilization is coupled to environmental noise.

Protein and lipid nanodomains are prevalent on the surface of mammalian cells. In particular, it has been recently recognized that ion channels assemble into surface nanoclusters in the soma of cultured neurons. However, the interactions of these molecules with surface nanodomains display a considerable degree of heterogeneity. Here, we investigate this heterogeneity and develop statistical tools based on the recurrence of individual trajectories to identify subpopulations within ion channels in the neuronal surface. We specifically study the dynamics of the ${\mathrm{K}}^{+}$ channel Kv1.4 and the ${\mathrm{Na}}^{+}$ channel Nav1.6 on the surface of cultured hippocampal neurons at the single-molecule level. We find that both these molecules are expressed in two different forms with distinct kinetics with regards to surface interactions, emphasizing the complex proteomic landscape of the neuronal surface. Further, the tools presented in this work provide new methods for the analysis of membrane nanodomains, transient confinement, and identification of populations within single-particle trajectories.

Micro-organisms expend energy moving through complex media. While propulsion speed is an important property of locomotion, efficiency is another factor that may determine the swimming gait adopted by a micro-organism in order to locomote in an energetically favorable manner. The efficiency of swimming in a Newtonian fluid is well characterized for different biological and artificial swimmers. However, these swimmers often encounter biological fluids displaying shear-thinning viscosities. Little is known about how this nonlinear rheology influences the efficiency of locomotion. Does the shear-thinning rheology render swimming more efficient or less? How does the swimming efficiency depend on the propulsion mechanism of a swimmer and rheological properties of the surrounding shear-thinning fluid? In this work, we address these fundamental questions on the efficiency of locomotion in a shear-thinning fluid by considering the squirmer model as a general locomotion model to represent different types of swimmers. Our analysis reveals how the choice of surface velocity distribution on a squirmer may reduce or enhance the swimming efficiency. We determine optimal shear rates at which the swimming efficiency can be substantially enhanced compared with the Newtonian case. The nontrivial variations of swimming efficiency prompt questions on how micro-organisms may tune their swimming gaits to exploit the shear-thinning rheology. The findings also provide insights into how artificial swimmers should be designed to move through complex media efficiently.