We consider the multiple scattering of light by fluctuations of the director in a nematic liquid crystal. Using methods of numerical simulation the peak of the coherent backscattering and the coefficients of anisotropic diffusion are calculated. The calculations were carried out without simplifying assumptions on the properties of the liquid crystal. The process of multiple scattering was simulated as a random walk of photons in the medium. We investigated in detail the transition to the diffusion regime. The dependence of the diffusion coefficients on the applied magnetic field and the wavelength of light were studied. The results of simulation showed a nonmonotonic dependence of the diffusion coefficients on the external magnetic field. For calculation of the peak of the coherent backscattering we used the semianalytical approach as long as in nematic liquid crystals this peak is extremely narrow. The parameters of the backscattering peak and of diffusion coefficients which were found in numerical simulations were compared with the experimental data and the results of analytical calculation.

A number of biological processes, such as invasive plant species and cell migration, are composed of two key mechanisms: motility and reproduction. Due to the spatially exclusive interacting behavior of these processes a cellular automata (CA) model is specified to simulate a one-dimensional invasion process. Three (independence, Poisson, and 2D-Markov chain) approximations are considered that attempt to capture the average behavior of the CA. We show that our 2D-Markov chain approximation accurately predicts the state of the CA for a wide range of motility and reproduction rates.

Water flowing over a loose granular bed organizes into a braided river, a network of ephemeral and interacting channels. The temporal and spatial evolution of this network of braided channels is not yet quantitatively understood. In $\sim 1$ m-scale experiments, we found that individual channels exhibit a self-similar geometry and near-threshold transport conditions. Measurements of the rate of growth of topographic correlation length scales, the time scale of system-slope establishment, and the random spatial decorrelation of channel locations indicate together that the evolution of the braided river system may be diffusive in nature. This diffusion is due to the separation of scales between channel formation and network evolution, and the random motion of interacting channels when viewed at a coarse-grained scale.

We theoretically study the motion of surface-active Janus particles, driven by an effective slip velocity due to a nonuniform temperature or concentration field $\psi$. With parameters realized in thermal traps, we find that the torque exerted by the gradient $\nabla \psi$ inhibits rotational diffusion and favors alignment of the particle axes. In a swarm of active particles, this polarization adds a novel term to the drift velocity and modifies the collective behavior. Self-polarization in a nonuniform laser beam could be used for guiding hot particles along a given trajectory.

We introduce the violation fraction $\upsilon$ as the cumulative fraction of time that a mesoscopic system spends consuming entropy at a single trajectory in phase space. We show that the fluctuations of this quantity are described in terms of a symmetry relation reminiscent of fluctuation theorems, which involve a function $\Phi$, which can be interpreted as an entropy associated with the fluctuations of the violation fraction. The function $\Phi$, when evaluated for arbitrary stochastic realizations of the violation fraction, is odd upon the symmetry transformations that are relevant for the associated stochastic entropy production. This fact leads to a detailed fluctuation theorem for the probability density function of $\Phi$. We study the steady-state limit of this symmetry in the paradigmatic case of a colloidal particle dragged by optical tweezers through an aqueous solution. Finally, we briefly discuss possible applications of our results for the estimation of free-energy differences from single-molecule experiments.

Recent experiments have demonstrated that identical material samples can charge one another after being brought into symmetric contact. The mechanism for this charging is not known. In this article, we use a simplified one-dimensional lattice model to analyze charging in the context of agitated particles. We find that the electric field from a single weakly polarized grain can feed back on itself by polarizing its neighbors, leading to an exponential growth in polarization. We show that, by incorporating partial neutralization between neighboring polarized particles, either uniform alignment of dipoles or complex charge and polarization waves can be produced. We reproduce a polarized state experimentally using identical colliding particles and raise several issues for future study.

The hydrodynamic equation governing the homogeneous time evolution of the temperature in a model of confined granular gas is studied by means of the Enskog equation. The existence of a normal solution of the kinetic equation is assumed as a condition for hydrodynamics. Dimensional analysis implies a scaling of the distribution function that is used to determine it in the first Sonine approximation, with a coefficient that evolves in time through its dependence on the temperature. The theoretical predictions are compared to numerical results obtained by the direct simulation Monte Carlo method and a good agreement is found. The relevance of the normal homogeneous distribution function to derive inhomogeneous hydrodynamic equations, for instance using the Champan-Enskog algorithm, is indicated.

The role of interfacial rheology on the bulk linear viscoelastic moduli of low concentration bovine albumin solutions is probed. Previously reported soft gel properties of these systems were attributed to either protein aggregation or organization within the bulk. Instead, these behaviors are shown to be attributable to the measurement error caused by interfacial rheology due to adsorption of bovine serum albumin to the air and water interface. Even at low bulk concentrations, fast interfacial adsorption results in erroneous measurements. When these effects are removed, the solutions are viscous dominated with a dynamic viscosity slightly larger than water.

Self-assembly of particles, suspended in a drying droplet, were studied by the Monte Carlo method. The Brownian diffusion of particles was simulated accounting for the effect of hydrodynamic retardation and interparticle interactions. The model allowed for explaining formation of the “coffee ring” patterns even without accounting for the radial flows towards the three-phase contact line. Morphologies of the drying patterns and their dependence on interparticle interactions and concentration of particles are discussed.

We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We focus on the common case where the bidirectional coupling of intracellular calcium concentration and membrane voltage dynamics produces calcium and voltage alternans that are temporally in phase. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolarization and calcium transient (CaT) alternans, similarly to the well-studied case of voltage-driven alternans. In contrast, further away from the bifurcation, the amplitude of CaT alternans jumps discontinuously at the nodes separating out-of-phase regions, while the amplitude of repolarization alternans remains smooth. We identify universal dynamical features of SDA pattern formation and evolution in the presence of those jumps. We show that node motion of discontinuous SDA patterns is strongly hysteretic even in homogeneous tissue due to the novel phenomenon of “unidirectional pinning”: node movement can only be induced towards, but not away from, the pacing site in response to a change of pacing rate or physiological parameter. In addition, we show that the wavelength of discontinuous SDA patterns scales linearly with the conduction velocity restitution length scale, in contrast to the wavelength of smooth patterns that scales sublinearly with this length scale. Those results are also shown to be robust against cell-to-cell fluctuations due to the property that unidirectional node motion collapses multiple jumps accumulating in nodal regions into a single jump. Amplitude equation predictions are in good overall agreement with ionic model simulations. Finally, we briefly discuss physiological implications of our findings. In particular, we suggest that due to the tendency of conduction blocks to form near nodes, the presence of unidirectional pinning makes calcium-driven alternans potentially more arrhythmogenic than voltage-driven alternans.

Fluctuations in gene expression give identical cells access to a spectrum of phenotypes that can serve as a transient, nongenetic basis for natural selection by temporarily increasing drug resistance. In this study, we demonstrate using mathematical modeling and simulation that certain gene regulatory network motifs, specifically coherent feedforward loop motifs, can facilitate the development of nongenetic resistance by increasing cell-to-cell variability and the time scale at which beneficial phenotypic states can be maintained. Our results highlight how regulatory network motifs enabling transient, nongenetic inheritance play an important role in defining reproductive fitness in adverse environments and provide a selective advantage subject to evolutionary pressure.

We discuss the exit probability of the one-dimensional $q$-voter model and present tools to obtain estimates about this probability, both through simulations in large networks (around ${10}^7$ sites) and analytically in the limit where the network is infinitely large. We argue that the result $E\left(\rho \right)=\frac{{\rho }^q}{{\rho }^q+{(1-\rho )}^q}$, that was found in three previous works [F. Slanina, K. Sznajd-Weron, and P. Przybyła, Europhys. Lett. 82, 18006 (2008); R. Lambiotte and S. Redner, Europhys. Lett. 82, 18007 (2008), for the case $q=2;$ and P. Przybyła, K. Sznajd-Weron, and M. Tabiszewski, Phys. Rev. E 84, 031117 (2011), for $q>2$] using small networks (around ${10}^3$ sites), is a good approximation, but there are noticeable deviations that appear even for small systems and that do not disappear when the system size is increased (with the notable exception of the case $q=2$). We also show that, under some simple and intuitive hypotheses, the exit probability must obey the inequality $\frac{{\rho }^q}{{\rho }^q+(1-\rho )}\le E\left(\rho \right)\le \frac{\rho }{\rho +{(1-\rho )}^q}$ in the infinite size limit. We believe this settles in the negative the suggestion made [S. Galam and A. C. R. Martins, Europhys. Lett. 95, 48005 (2001)] that this result would be a finite size effect, with the exit probability actually being a step function. We also show how the result that the exit probability cannot be a step function can be reconciled with the Galam unified frame, which was also a source of controversy.

A one-component bistable reaction-diffusion system with asymmetric nonlocal coupling is derived as a limiting case of a two-component activator-inhibitor reaction-diffusion model with differential advection. The effects of asymmetric nonlocal couplings in such a bistable reaction-diffusion system are then compared to the previously studied case of a system with symmetric nonlocal coupling. We carry out a linear stability analysis of the spatially homogeneous steady states of the model and numerical simulations of the model to show how the asymmetric nonlocal coupling controls and alters the steady states and the front dynamics in the system. In a second step, a third fast reaction-diffusion equation is included which induces the formation of more complex patterns. A linear stability analysis predicts traveling waves for asymmetric nonlocal coupling, in contrast to a stationary Turing patterns for a system with symmetric nonlocal coupling. These findings are verified by direct numerical integration of the full equations with nonlocal coupling.

Spatially localized structures are key components of turbulence and other spatiotemporally chaotic systems. From a dynamical systems viewpoint, it is desirable to obtain corresponding exact solutions, though their existence is not guaranteed. A damping filter method is introduced to obtain variously localized solutions and adapted in two typical cases. This method introduces a spatially selective damping effect to make a good guess at the exact solution, and we can obtain an exact solution through a continuation with the damping amplitude. The first target is a steady solution to the Swift-Hohenberg equation, which is a representative of bistable systems in which localized solutions coexist and a model for spanwise-localized cases. Not only solutions belonging to the well-known snaking branches but also those belonging to isolated branches known as “isolas” are found with continuation paths between them in phase space extended with the damping amplitude. This indicates that this spatially selective excitation mechanism has an advantage in searching spatially localized solutions. The second target is a spatially localized traveling-wave solution to the Kuramoto-Sivashinsky equation, which is a model for streamwise-localized cases. Since the spatially selective damping effect breaks Galilean and translational invariances, the propagation velocity cannot be determined uniquely while the damping is active, and a singularity arises when these invariances are recovered. We demonstrate that this singularity can be avoided by imposing a simple condition, and a localized traveling-wave solution is obtained with a specific propagation speed.

The influence of gas-surface interactions of a dilute gas confined between two parallel walls on the heat flux predictions is investigated using a combined Monte Carlo (MC) and molecular dynamics (MD) approach. The accommodation coefficients are computed from the temperature of incident and reflected molecules in molecular dynamics and used as effective coefficients in Maxwell-like boundary conditions in Monte Carlo simulations. Hydrophobic and hydrophilic wall interactions are studied, and the effect of the gas-surface interaction potential on the heat flux and other characteristic parameters like density and temperature is shown. The heat flux dependence on the accommodation coefficient is shown for different fluid-wall mass ratios. We find that the accommodation coefficient is increasing considerably when the mass ratio is decreased. An effective map of the heat flux depending on the accommodation coefficient is given and we show that MC heat flux predictions using Maxwell boundary conditions based on the accommodation coefficient give good results when compared to pure molecular dynamics heat predictions. The accommodation coefficients computed for a dilute gas for different gas-wall interaction parameters and mass ratios are transferred to compute the heat flux predictions for a dense gas. Comparison of the heat fluxes derived using explicit MD, MC with Maxwell-like boundary conditions based on the accommodation coefficients, and pure Maxwell boundary conditions are discussed. A map of the heat flux dependence on the accommodation coefficients for a dense gas, and the effective accommodation coefficients for different gas-wall interactions are given. In the end, this approach is applied to study the gas-surface interactions of argon and xenon molecules on a platinum surface. The derived accommodation coefficients are compared with values of experimental results.

We show that for a large class of stationary Markov processes the total variation distance between the final equilibrium distribution and that at a given time is a strongly monotonic vanishing function. We illustrate this for basic paradigmatic processes and discuss how, for systems susceptible to a canonical description, this can be interpreted as a statistical arrow of time that exists besides the standard decrease of free energy.

Evaluating the entropy production (EP) along a stochastic trajectory requires the knowledge of the system probability distribution, an ensemble quantity notoriously difficult to measure. In this paper we show that the EP of nonautonomous systems in contact with multiple reservoirs can be expressed solely in terms of physical quantities measurable at the single-trajectory level with a suitable preparation of the initial condition. As a result, we identify universal energy and particle fluctuation relations valid for any measurement time. We apply our findings to an electronic junction model, which may be used to verify our prediction experimentally.

We study the center-of-mass motion in systems of trapped interacting particles with space- and velocity-dependent friction and anharmonic traps. Our approach, based on a dynamical ansatz assuming a fixed density profile, allows us to obtain information at once for a wide range of binary interactions and interaction strengths, at linear and nonlinear levels. Our findings are first tested on different simple models by comparison with direct numerical simulations. Then, we apply the method to characterize the motion of the center of mass of a magneto-optical trap and its dependence on the number of trapped atoms. Our predictions are compared with experiments performed on a large ${\text{Rb}}^{85}$ magneto-optical trap.

Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segregation criterion based on Navier-Stokes granular hydrodynamics is written in terms of the tracer diffusion transport coefficients, whose dependence on the parameters of the system (masses, sizes, and coefficients of restitution) is explicitly determined from a solution of the inelastic Boltzmann equation. The theoretical predictions are validated by means of Monte Carlo and molecular dynamics simulations, showing that Navier-Stokes hydrodynamics produces accurate segregation criteria even under strong shearing and/or inelasticity.