We analyze the effect of the Basset history force on the sedimentation or rising of inertial particles in a two-dimensional convection flow. When memory effects are neglected, the system exhibits rich dynamics, including periodic, quasiperiodic, and chaotic attractors. Here we show that when the full advection dynamics is considered, including the history force, both the nature and the number of attractors change, and a fractalization of their basins of attraction appears. In particular, we show that the history force significantly weakens the horizontal diffusion and changes the speed of sedimentation or rising. The influence of the history force is dependent on the size of the advected particles, being stronger for larger particles.

Standing waves are studied as solutions of a complex Ginzburg-Landau equation subjected to local and global time-delay feedback terms. The onset is described as an instability of the uniform oscillations with respect to spatially periodic perturbations. The solution of the standing wave pattern is given analytically and studied through simulations.

We show that exploding dissipative solitons can arise in a reaction-diffusion system for a range of parameters. As a function of a vorticity parameter, we observe a sequence of transitions from oscillatory localized states via meandering dissipative solitons to exploding dissipative solitons propagating in one direction for long times followed by the reverse cascade back to oscillatory localized states. While exploding dissipative solitons are known from the cubic-quintic complex Ginzburg-Landau (CGL) equation, propagating exploding dissipative solitons appear to require for their existence a system of lower symmetry such as the reaction-diffusion model studied here.

Textures (i.e., smooth space nonuniform distributions of the order parameter) in biaxial nematics turned out to be much more complex and interesting than expected. Scanning the literature we find only a very few publications on this topic. Thus, the immediate motivation of the present paper is to develop a systematic procedure to study, classify, and visualize possible textures in biaxial nematics. Based on the elastic energy of a biaxial nematic (written in the most simple form that involves the least number of phenomenological parameters) we derive and solve numerically the Lagrange equations of the first kind. It allows one to visualize the solutions and offers a deep insight into their geometrical and topological features. Performing Fourier analysis we find some particular textures possessing two or more characteristic space periods (we term such solutions quasiperiodic ones because the periods are not necessarily commensurate). The problem is not only of intellectual interest but also of relevance to optical characteristics of the liquid-crystalline textures.

We study the driven flow of an unentangled glassy polymer film with a free upper surface and supported below by a substrate using nonequilibrium molecular dynamics simulations based on a bead-spring model. Above the glass transition temperature $T_g$, simple Poiseuille laminar flow is observed with the film mobility defined as the flow current density per unit pressure gradient scaling as $h^3$ with the film thickness $h$. Below $T_g$, the film mobility becomes independent of $h$, signifying surface transport. This is in full agreement with recent experiments on the time evolution of capillary waves in polystyrene films supported by silica. A mobile layer is found responsible for the surface transport, as previously conjectured. Our result also shows that it has a velocity profile decaying exponentially into the bulk.

Exact equations of motion for the microscopically defined collective density $\widehat{\rho }(\mathbf{x},t)$ and the momentum density $\widehat{\mathbf{g}}(\mathbf{x},t)$ of a fluid have been obtained in the past starting from the corresponding Langevin equations representing the dynamics of the fluid particles. In the present work we average these exact equations of microscopic dynamics over the local equilibrium distribution to obtain stochastic partial differential equations for the coarse-grained densities with smooth spatial and temporal dependence. In particular, we consider Dean's exact balance equation for the microscopic density of a system of interacting Brownian particles to obtain the basic equation of the dynamic density functional theory with noise. Our analysis demonstrates that on thermal averaging the dependence of the exact equations on the bare interaction potential is converted to dependence on the corresponding thermodynamic direct correlation functions in the coarse-grained equations.

The normal modes for three to seven particle two-dimensional (2D) dust clusters in a complex plasma are investigated using an $N$-body simulation. The ion wakefield downstream of each particle is shown to induce coupling between horizontal and vertical modes. The rules of mode coupling are investigated by classifying the mode eigenvectors employing the Bessel and trigonometric functions indexed by order integers ($m$, $n$). It is shown that coupling only occurs between two modes with the same $m$ and that horizontal modes having a higher shear contribution exhibit weaker coupling. Three types of resonances are shown to occur when two coupled modes have the same frequency. Discrete instabilities caused by both the first and third type of resonances are verified and instabilities caused by the third type of resonance are found to induce melting. The melting procedure is observed to go through a two-step process with the solid-liquid transition closely obeying the Lindemann criterion.

The dynamics of particles in solution or suspension is influenced by thermal fluctuations and hydrodynamic interactions. Several mesoscale methods exist to account for these solvent-induced effects such as Brownian dynamics with hydrodynamic interactions and hybrid molecular dynamics-stochastic rotation dynamics methods. Here we compare two ways of coupling solutes to the solvent with stochastic rotation dynamics (SRD) to Brownian dynamics with and without explicit hydrodynamic interactions. In the first SRD scheme [SRD with collisional coupling (CC)] the solutes participate in the collisional step with the solvent and in the second scheme [SRD with central force coupling (CFC)] the solutes interact through direct forces with the solvent, generating slip boundary conditions. We compare the transport coefficients of neutral and charged solutes in a model system obtained by these simulation schemes. Brownian dynamics without hydrodynamic interactions is used as a reference to quantify the influence of hydrodynamics on the transport coefficients as modeled by the different methods. We show that, in the dilute range, the SRD CFC method provides results similar to those of Brownian dynamics with hydrodynamic interactions for the diffusion coefficients and for the electrical conductivity. The SRD CC scheme predicts diffusion coefficients close to those obtained by Brownian dynamic simulations without hydrodynamic interactions, but accounts for part of the influence of hydrodynamics on the electrical conductivity.

The role of dynamical topologies in the evolution of cooperation has received considerable attention, as some studies have demonstrated that dynamical networks are much better than static networks in terms of boosting cooperation. Here we study a dynamical model of evolution of cooperation on stochastic dynamical networks in which there are no permanent partners to each agent. Whenever a new link is created, its duration is randomly assigned without any bias or preference. We allow the agent to adaptively adjust the duration of each link during the evolution in accordance with the feedback from game interactions. By Monte Carlo simulations, we find that cooperation can be remarkably promoted by this adaptive dynamical linking mechanism both for the game of pairwise interactions, such as the Prisoner's Dilemma game (PDG), and for the game of group interactions, illustrated by the public goods game (PGG). And the faster the adjusting rate, the more successful the evolution of cooperation. We also show that in this context weak selection favors cooperation much more than strong selection does. What is particularly meaningful is that the prosperity of cooperation in this study indicates that the rationality and selfishness of a single agent in adjusting social ties can lead to the progress of altruism of the whole population.

A boundary condition for the Boltzmann equation (kinetic boundary condition, KBC) at the vapor-liquid interface of argon is constructed with the help of molecular dynamics (MD) simulations. The KBC is examined at a constant liquid temperature of 85 K in a wide range of nonequilibrium states of vapor. The present investigation is an extension of a previous one by Ishiyama, Yano, and Fujikawa [Phys. Rev. Lett. 95, 084504 (2005)] and provides a more complete form of the KBC. The present KBC includes a thermal accommodation coefficient in addition to evaporation and condensation coefficients, and these coefficients are determined in MD simulations uniquely. The thermal accommodation coefficient shows an anisotropic behavior at the interface for molecular velocities normal versus tangential to the interface. It is also found that the evaporation and condensation coefficients are almost constant in a fairly wide range of nonequilibrium states. The thermal accommodation coefficient of the normal velocity component is almost unity, while that of the tangential component shows a decreasing function of the density of vapor incident on the interface, indicating that the tangential velocity distribution of molecules leaving the interface into the vapor phase may deviate from the tangential parts of the Maxwell velocity distribution at the liquid temperature. A mechanism for the deviation of the KBC from the isotropic Maxwell KBC at the liquid temperature is discussed in terms of anisotropic energy relaxation at the interface. The liquid-temperature dependence of the present KBC is also discussed.

We perform phase-field simulations of unsteady crystal growth in a three-dimensional capillary. Motivated by the appearance of chirality-symmetry breaking periodic states in our preceding study, we here focus on more general dynamic states. Most of these are obtained in the limit of isotropic surface tension, but we test genericity by looking at a few cases with weak anisotropy. Whereas steady states are similar for all channel shapes studied so far, including channels with circular, hexagonal, quadratic, and triangular cross sections, there is a stronger dependence on the cross section for time-dependent states. Various oscillatory modes are identified and discussed, including rotating and swinging patterns as well as pulsating modes containing one, two, and four fingers, respectively.

This work analyzes methods for the identification and the stability under perturbation of a territorial community structure with specific reference to transportation networks. We considered networks of commuters for a city and an insular region. In both cases, we have studied the distribution of commuters’ trips (i.e., home-to-work trips and vice versa). The identification and stability of the communities’ cores are linked to the land-use distribution within the zone system, and therefore their proper definition may be useful to transport planners.

We study the influence of hydrophobic effects on streaming potential mediated flow through a narrow confinement. In a clear departure from the approach used in prior works, we use a phase-field model to capture the hydrophobicity-induced depletion in the near wall region, and express the variation of viscosity and permittivity across the interfacial layer in terms of the phase-field variable. We then use these in the determination of the flow velocity, and highlight the sensitive interplay between the intrinsic length scale of the electrical double layer and that of the depletion in terms of the variations of an effective normalized viscosity that captures the electroviscous effect. We expect that this work will be an important step forward in the realistic continuum modeling of interfacial physics in the particular context of streaming potential mediated flows.

Nonlinear dynamics of an intense circularly polarized laser beam interacting with a hot magnetized plasma is investigated. Using a relativistic fluid model, a modified nonlinear Schrödinger equation is derived based on a quasineutral approximation, which is valid for hot plasma. Using a three-dimensional model, spatial-temporal development of the laser pulse is investigated. The occurrence of some nonlinear phenomena such as self-focusing, self-modulation, light trapping, and filamentation of the laser pulse is discussed. Also the effect of polarization and external magnetic field on the nonlinear evolution of these phenomena is studied.

Phase separation in binary mixtures in the presence of Janus particles has been studied in terms of a Cahn-Hilliard model coupled to the Langevin equations describing the particle dynamics. We demonstrate that the phase separation process is arrested leading to unexpected regular stripe patterns in the concentration field. The underlying pattern forming mechanism has been elucidated: The twofold absorption properties on the surface of Janus particles with respect to the two components of a binary mixture trigger in their neighborhood spatial concentration variations. They result in an effective interaction between the particles mediated by the binary mixture. Our findings open a route to design composite materials with nanoscale lamellar morphologies where the pattern wavelength can be tuned by changing the wetting properties of the Janus particles.

The dynamics of edge dislocations with parallel Burgers vectors, moving in the same slip plane, is mapped onto Dyson's model of a two-dimensional Coulomb gas confined in one dimension. We show that the tail distribution of the velocity of dislocations is power law in form, as a consequence of the pair interaction of nearest neighbors in one dimension. In two dimensions, we show the presence of a pairing phase transition in a system of interacting dislocations with parallel Burgers vectors. The scaling exponent of the velocity distribution at effective temperatures well below this pairing transition temperature can be derived from the nearest-neighbor interaction, while near the transition temperature, the distribution deviates from the form predicted by the nearest-neighbor interaction, suggesting the presence of collective effects.

The central limit theorem ensures that a sum of random variables tends to a Gaussian distribution as their total number tends to infinity. However, for a class of positive random variables, we find that the sum tends faster to a log-normal distribution. Although the sum tends eventually to a Gaussian distribution, the distribution of the sum is always close to a log-normal distribution rather than to any Gaussian distribution if the summands are numerous enough. This is in contrast to the current consensus that any log-normal distribution is due to a product of random variables, i.e., a multiplicative process, or equivalently to nonlinearity of the system. In fact, the log-normal distribution is also observable for a sum, i.e., an additive process that is typical of linear systems. We show conditions for such a sum, an analytical example, and an application to random scalar fields such as those of turbulence.

We consider the domino tilings of an Aztec diamond with a cut-off corner of macroscopic square shape and given size and address the bulk properties of tilings as the size is varied. We observe that the free energy exhibits a third-order phase transition when the cut-off square, increasing in size, reaches the arctic ellipse—the phase separation curve of the original (unmodified) Aztec diamond. We obtain this result by studying the thermodynamic limit of a certain nonlocal correlation function of the underlying six-vertex model with domain wall boundary conditions, the so-called emptiness formation probability (EFP). We consider EFP in two different representations: as a $\tau$ function for Toda chains and as a random matrix model integral. The latter has a discrete measure and a linear potential with hard walls; the observed phase transition shares properties with both Gross-Witten-Wadia and Douglas-Kazakov phase transitions.

It is often stated that heat baths with finite degrees of freedom i.e., finite baths, are sources of Tsallis distributions for classical Hamiltonian systems. By using well-known fundamental statistical mechanics expressions, we rigorously show that Tsallis distributions with fat tails are possible only for finite baths with constant negative heat capacity, while constant positive heat capacity finite baths yield decays with sharp cutoff with no fat tails. However, the correspondence between Tsallis distributions and finite baths holds at the expense of violating the equipartition theorem for finite classical systems at equilibrium. We comment on the implications of the finite bath for the recent attempts towards a $q$-generalized central limit theorem.

When particle flux is regulated by multiple factors such as particle supply and varying transport rate, it is important to identify the respective dominant regimes. We extend the well-studied totally asymmetric simple exclusion model to investigate the interplay between a controlled entrance and a local defect site. The model mimics cellular transport phenomena where there is typically a finite particle pool and nonuniform moving rates due to biochemical kinetics. Our simulations reveal regions where, despite an increasing particle supply, the current remains constant while particles redistribute in the system. Exploiting a domain wall approach with mean-field approximation, we provide a theoretical ground for our findings. The results in steady-state current and density profiles provide quantitative insights into the regulation of the transcription and translation process in bacterial protein synthesis.