The adsorption of Ar on planar structureless substrates of alkali metals, alkaline-earth metal Mg, ${\text{CO}}_2$, and Au was analyzed by applying a density functional formalism which includes a recently proposed effective attractive pair potential conditioned to Ar. It is shown that this approach reproduces the experimental surface tension of the liquid-vapor interface over the entire bulk coexistence curve for temperatures $T$ spanning from the triple point $T_t$ up to the critical point $T_c$. The wetting properties were studied over the entire range temperatures $T_t\leftrightarrow T_c$. It was found that Ar wets all the investigated surfaces. The adsorption isotherms for alkali metals exhibit first-order phase transitions. Prewetting lines were resolved even for the less attractive surfaces. In the cases of Mg, ${\text{CO}}_2$, and Au a continuous growth for $T\ge T_t$ was obtained. A comparison with experimental data and other microscopic calculations is reported.

Autocatalytic reactions may propagate as solitary waves, namely, at a constant front velocity and with a stationary concentration profile, resulting from a balance between molecular diffusion and chemical reaction. When the reaction is exothermic, a thermal wave is linked to the chemical front. As the thermal diffusivity is nearly two orders of magnitude larger than the molecular one, the temperature profile spreads over length scales (mm) two orders of magnitude larger than the concentration one. Using an infrared camera, we measure the temperature profiles for a chlorite-tetrathionate autocatalytic reaction. The profiles are compared quantitatively to lattice Bhatnagar-Gross-Krook (BGK) numerical simulations. Our analysis also accounts for the lack of observation of the thermal wave for the iodate arsenous acid reaction.

We examine a naming game with two agents trying to establish a common vocabulary for $n$ objects. Such efforts lead to the emergence of language that allows for an efficient communication and exhibits some degree of homonymy and synonymy. Although homonymy reduces the communication efficiency, it seems to be a dynamical trap that persists for a long, and perhaps indefinite, time. On the other hand, synonymy does not reduce the efficiency of communication but appears to be only a transient feature of the language. Thus, in our model the role of synonymy decreases and in the long-time limit it becomes negligible. A similar rareness of synonymy is observed in present natural languages. The role of noise, that distorts the communicated words, is also examined. Although, in general, the noise reduces the communication efficiency, it also regroups the words so that they are more evenly distributed within the available “verbal” space.

The complex behavior that occurs when traffic lights are synchronized is studied for a row of interacting cars. The system is modeled through a cellular automaton. Two strategies are considered: all lights in phase and a “green wave” with a propagating green signal. It is found that the mean velocity near the resonant condition follows a critical scaling law. For the green wave, it is shown that the mean velocity scaling law holds even for random separation between traffic lights and is not dependent on the density. This independence on car density is broken when random perturbations are considered in the car velocity. Random velocity perturbations also have the effect of leading the system to an emergent state, where cars move in clusters, but with an average velocity which is independent of traffic light switching for large injection rates.

We study the evolution of cooperation in public goods games on different regular graphs as a function of the noise level underlying strategy adoptions. We focus on the effects that are brought about by different group sizes of public goods games in which individuals participate, revealing that larger groups of players may induce qualitatively different behavior when approaching the deterministic limit of strategy adoption. While by pairwise interactions an intermediate uncertainty by strategy adoptions may ensure optimal conditions for the survival of cooperators at a specific graph topology, larger groups warrant this only in the vicinity of the deterministic limit independently from the underlying graph. These discrepancies are attributed to the indirect linkage of otherwise not directly connected players, which is brought about by joint memberships within the larger groups. Thus, we show that increasing the group size may introduce an effective transition of the interaction topology, and that the latter shapes the noise dependence of the evolution of cooperation in case of pairwise interactions only.

Basic peculiarities of market price fluctuations are known to be well described by a recently developed random-walk model in a temporally deforming quadratic potential force whose center is given by a moving average of past price traces [M. Takayasu, T. Mizuno, and H. Takayasu, Physica A 370, 91 (2006)]. By analyzing high-frequency financial time series of exceptional events, such as bubbles and crashes, we confirm the appearance of higher-order potential force in the markets. We show statistical significance of its existence by applying the information criterion. This time series analysis is expected to be applied widely for detecting a nonstationary symptom in random phenomena.

We identify a type of pattern formation in spatially distributed active systems. We simulate one-dimensional two-component systems with predator-prey local interaction and pursuit-evasion taxis between the components. In a sufficiently large domain, spatially uniform oscillations in such systems are unstable with respect to small perturbations. This instability, through a transient regime appearing as spontaneous focal sources, leads to establishment of periodic traveling waves. The traveling wave regime is established even if boundary conditions do not favor such solutions. The stable wavelength is within a range bounded both from above and from below, and this range does not coincide with instability bands of the spatially uniform oscillations.

Laser cooling and crystallization of electron-ion plasma is studied using the Brownian dynamics simulation technique and taking into consideration the interaction of ions with the electron subsystem. It has been shown that the nonlinear dependence of laser friction force on the velocity of ions has to be taken into account in order to simulate in an adequate manner the cooling dynamics and obtain a correct estimate for minimum temperatures. It has been found that times required for formation of an ordered ionic structure can be much longer than the typical plasma cooling time.

The arrest of Langmuir wave collapse by quantum effects, first addressed by Haas and Shukla [Phys. Rev. E 79, 066402 (2009)] using a Rayleigh-Ritz trial function method is revisited, using rigorous estimates and systematic asymptotic expansions. The absence of blow up for the so-called quantum Zakharov equations is proved in two and three dimensions, whatever the strength of the quantum effects. The time-periodic behavior of the solution for initial conditions slightly in excess of the singularity threshold for the classical problem is established for various settings in two space dimensions. The difficulty of developing a consistent perturbative approach in three dimensions is also discussed and a semiphenomenological model is suggested for this case.

Inverse bremsstrahlung (IB) absorption and evolution of the electron distribution function (EDF) in a wide laser intensity range $\left({10}^{12}–{10}^{17}\text{W}/{\text{cm}}^2\right)$ have been studied systematically by a two velocity-dimension Fokker-Planck code. It is found that Langdon’s IB operator overestimates the absorption rate at high laser intensity, consequently with an overdistorted non-Maxwellian EDF. According to the small anisotropy of EDF in the oscillation frame, we introduce an IB operator which is similar to Langdon’s but without the low laser intensity limit. This operator is appropriate for self-consistently tackling the nonlinear effects of high laser intensity as well as non-Maxwellian EDF. Particularly, our operator is capable of treating IB absorption properly in the indirect and direct-drive inertial confinement fusion schemes with the National Ignition Facility and Laser MegaJoule laser parameters at focused laser intensity beyond ${10}^{15}\text{W}/{\text{cm}}^2$.

When thermal energies are weak, two-dimensional lamellar structures confined on a curved substrate display complex patterns arising from the competition between layer bending and compression in the presence of geometric constraints. We present broad design principles to engineer the geometry of the underlying substrate so that a desired lamellar pattern can be obtained by self-assembly. Two distinct physical effects are identified as key factors that contribute to the interaction between the shape of the underlying surface and the resulting lamellar morphology. The first is a local ordering field for the direction of each individual layer, which tends to minimize its curvature with respect to the three-dimensional embedding. The second is a nonlocal effect controlled by the intrinsic geometry of the surface that forces the normals to the (nearly incompressible) layers to lie on geodesics, leading to caustic formation as in optics. As a result, different surface morphologies with predominantly positive or negative Gaussian curvature can act as converging or diverging lenses, respectively. By combining these ingredients, as one would with different optical elements, complex lamellar morphologies can be obtained. This smectic optometry enables the manipulation of lamellar configurations for the design of materials.

It has been known for some time that some microorganisms can swim faster in high-viscosity gel-forming polymer solutions. These gel-like media come to mimic highly viscous heterogeneous environment that these microorganisms encounter in-vivo. The qualitative explanation of this phenomena first offered by Berg and Turner [Nature (London) 278, 349 (1979)], suggests that propulsion enhancement is a result of flagellum pushing on quasi-rigid loose polymer network formed in some polymer solutions. Inspired by these observations, inertia-less propulsion in a heterogeneous viscous medium composed of sparse array of stationary obstacles embedded into a incompressible Newtonian liquid is considered. It is demonstrated that for prescribed propulsion gaits, including propagating surface distortions and rotating helical filament, the propulsion speed is enhanced when compared to swimming in purely viscous solvent. It is also shown that the locomotion in heterogenous viscous media is characterized by improved hydrodynamic efficiency. The results of the rigorous numerical simulation of the rotating helical filament propelled through a random sparse array of stationary obstructions are in close agreement with predictions of the proposed resistive force theory based on effective media approximation.

We discuss the physics of embolic stroke using a minimal model of emboli moving through the cerebral arteries. Our model of the blood flow network consists of a bifurcating tree into which we introduce particles (emboli) that halt flow on reaching a node of similar size. Flow is weighted away from blocked arteries inducing an effective interaction between emboli. We justify the form of the flow weighting using a steady flow (Poiseuille) analysis and a more complicated nonlinear analysis. We discuss free flowing and heavily congested limits and examine the transition from free flow to congestion using numerics. The correlation time is found to increase significantly at a critical value and a finite-size scaling is carried out. An order parameter for nonequilibrium critical behavior is identified as the overlap of blockages’ flow shadows. Our work shows embolic stroke to be a feature of the cerebral blood flow network on the verge of a phase transition.

A standing difficulty in the problem of blood vessel tethering has been that only one of the two required boundary conditions can be fully specified, namely, that at the inner (endothelial) wall surface. The other, at the outer layer of the vessel wall, is not known except in the limiting case where the wall is fully tethered such that its outer layer is prevented from any displacement. In all other cases, where the wall is either free or partially tethered, a direct boundary condition is not available. We present a method of determining this missing boundary condition by considering the limiting case of a semi-infinite wall. The result makes it possible to define the degree of tethering imposed by surrounding tissue more accurately in terms of the displacement of the outer layer of the vessel wall, rather than in terms of equivalent added mass which has been done in the past. This new approach makes it possible for the first time to describe the effect of partial tethering in its full range, from zero to full tethering. The results indicate that high tethering leads to high stresses and low displacements within the vessel wall, while low tethering leads to low stresses and high displacements. Since both extremes would be damaging to wall tissue, particularly elastin, this suggest that moderate tethering would be optimum in the physiological setting.

We develop a model of the fluctuation dynamo in which the magnetic field is confined to thin flux ropes advected by a multiscale model of turbulence. Magnetic dissipation occurs only via reconnection of the flux ropes. This model can be viewed as an implementation of the asymptotic limit $R_{\text{m}}\to \infty$ for a continuous magnetic field, where magnetic dissipation is strongly localized to small regions of strong-field gradients. We investigate the kinetic-energy release into heat mediated by the dynamo action, both in our model and by solving the induction equation with the same flow. We find that a flux-rope dynamo is an order of magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic energy release in reconnections has a power-law form with the slope $-3$, consistent with the solar corona heating by nanoflares.

Embeddings are diffeomorphisms between some unseen physical attractor and a reconstructed image. Different embeddings may or may not be equivalent under isotopy. We regard embeddings as representations of the attractor, review the labels required to distinguish inequivalent representations for an important class of dynamical systems, and discuss the systematic ways inequivalent embeddings become equivalent as the embedding dimension increases until there is finally only one “universal” embedding in a suitable dimension.

We present microcanonical replica exchange molecular dynamics simulations as an alternative to canonical ones. Its advantage is the easily tunable high acceptance rate for replica exchange. We present the theory, comment on its actual implementation, and demonstrate its application for a common test case, the trp-cage protein.

In this paper we present the event-chain algorithms, which are fast Markov-chain Monte Carlo methods for hard spheres and related systems. In a single move of these rejection-free methods, an arbitrarily long chain of particles is displaced, and long-range coherent motion can be induced. Numerical simulations show that event-chain algorithms clearly outperform the conventional Metropolis method. Irreversible versions of the algorithms, which violate detailed balance, improve the speed of the method even further. We also compare our method with a recent implementations of the molecular-dynamics algorithm.

Using the parallel tempering algorithm and graphics processing unit accelerated techniques, we have performed large-scale Monte Carlo simulations of the Ising model on a square lattice with antiferromagnetic (repulsive) nearest-neighbor and next-nearest-neighbor interactions of the same strength and subject to a uniform magnetic field. Both transitions from the $\left(2\times 1\right)$ and row-shifted $\left(2\times 2\right)$ ordered phases to the paramagnetic phase are continuous. From our data analysis, re-entrance behavior of the $\left(2\times 1\right)$ critical line and a bicritical point which separates the two ordered phases at $\text{T}=0$ are confirmed. Based on the critical exponents we obtained along the phase boundary, Suzuki’s weak universality seems to hold.

We present a numerical technique, namely, triangular tessellation, to calculate the free energy associated with the adsorption of a colloidal particle at a flat interface. The theory and numerical scheme presented here are sufficiently general to handle nonconvex patchy colloids with arbitrary surface patterns characterized by a wetting angle, e.g., amphiphilicity. We ignore interfacial deformation due to capillary, electrostatic, or gravitational forces, but the method can be extended to take such effects into account. It is verified that the numerical method presented is accurate and sufficiently stable to be applied to more general situations than presented in this paper. The merits of the tessellation method prove to outweigh those of traditionally used semianalytic approaches, especially when it comes to generality and applicability.