The stationary points of the potential energy function of the ${\phi }^4$ model on a two-dimensional square lattice with nearest-neighbor interactions are studied by means of two numerical methods: a numerical homotopy continuation method and a globally convergent Newton-Raphson method. We analyze the properties of the stationary points, in particular with respect to a number of quantities that have been conjectured to display signatures of the thermodynamic phase transition of the model. Although no such signatures are found for the nearest-neighbor ${\phi }^4$ model, our study illustrates the strengths and weaknesses of the numerical methods employed.

An equilibrium random surface multistep height model proposed in [Abraham and Newman, EPL 86, 16002 (2009)] is studied using a variant of the worm algorithm. In one limit, the model reduces to the two-dimensional Ising model in the height representation. When the Ising model constraint of single height steps is relaxed, the critical temperature and critical exponents are continuously varying functions of the parameter controlling height steps larger than one. Numerical estimates of the critical exponents can be mapped via a single parameter, the Coulomb gas coupling, to the exponents of the O$\left(n\right)$ loop model on the honeycomb lattice with $n⩽1$.

Irreversible thermodynamics provides interface conditions that yield temperature and chemical potential jumps at phase boundaries. The interfacial jumps allow unexpected transport phenomena, such as the inverted temperature profile [Pao, Phys. Fluids 14, 306 (1971)] and mass transfer from a cold to a warm liquid driven by a temperature difference across the vapor phase [Mills and Phillips, Chem. Phys. Lett. 372, 615 (2002)]. Careful evaluation of the thermodynamic laws has shown [Bedeaux et al., Physica A 169, 263 (1990)] that the inverted temperature profile is observed for processes with a high heat of vaporization. In this paper, we show that cold to warm mass transfer through the vapor from a cold to a warm liquid is only possible when the heat of evaporation is sufficiently small. A necessary criterium for the size of the mass transfer coefficient is given.

This paper studies the segregation behavior of binary granular particles with diameters at approximately 10:1 in a vertically vibrated container. An array of transitional separation patterns between reversed Brazilian nut (RBN) and Brazilian nut (BN) separations are observed, with their geometrical features carefully measured. The binary particle system develops into either a stable separation pattern when $f$ and Γ are relatively low or an oscillating pattern when $f$ and Γ are relatively high. We regard these patterns as different phases, in which the stable patterns can be divided into phases of RBN, RBN transitional (RBNT), BNT, and BN. A phase parameter λ between$-$1 and 1 is defined to describe the separation patterns based on the mass center height difference in large and small particles. By drawing $f$-Γ-λ phase diagrams, the system's tendency toward BN separation was found to increase with $f$ and decrease with Γ. Furthermore, the range of the tendency toward BN separation expands when the size of small particles rises. As the total mass of the small particles increases, the system's tendency toward RBN separation is enhanced. Abnormal points are also observed in the stable phase regions, and the oscillating phase shifts among the four stable phases with time. These stable phases can be explained via an analysis of the distribution of the dissipation energy, whereas the mechanism of the oscillating phase remains to be discovered.

Queuing models provide insight into the temporal inhomogeneity of human dynamics, characterized by the broad distribution of waiting times of individuals performing tasks. We theoretically study the queuing model of an agent trying to execute a task of interest, the priority of which may vary with time due to the agent's “state of mind.” However, its execution is disrupted by other tasks of random priorities. By considering the priority of the task of interest either decreasing or increasing algebraically in time, we analytically obtain and numerically confirm the bimodal and unimodal waiting time distributions with power-law decaying tails, respectively. These results are also compared to the updating time distribution of papers in arXiv.org and the processing time distribution of papers in Physical Review journals. Our analysis helps to understand human task execution in a more realistic scenario.

We study synchronization in a system of phase-only oscillators residing on the sites of a one-dimensional periodic lattice. The oscillators interact with a strength that decays as a power law of the separation along the lattice length and is normalized by a size-dependent constant. The exponent $\alpha$ of the power law is taken in the range $0⩽\alpha <1$. The oscillator frequency distribution is symmetric about its mean (taken to be zero) and is nonincreasing on $[0,\infty )$. In the continuum limit, the local density of oscillators evolves in time following the continuity equation that expresses the conservation of the number of oscillators of each frequency under the dynamics. This equation admits as a stationary solution the unsynchronized state uniform both in phase and over the space of the lattice. We perform a linear stability analysis of this state to show that when it is unstable, different spatial Fourier modes of fluctuations have different stability thresholds beyond which they grow exponentially in time with rates that depend on the Fourier modes. However, numerical simulations show that at long times all the nonzero Fourier modes decay in time, while only the zero Fourier mode (i.e., the “mean-field” mode) grows in time, thereby dominating the instability process and driving the system to a synchronized state. Our theoretical analysis is supported by extensive numerical simulations.

Modern theoretical methods for estimating the width of the chaotic layer in the presence of prominent marginal resonances are considered in the perturbed pendulum model of nonlinear resonance. The fields of applicability of these methods are explicitly and precisely formulated. The comparative accuracy is investigated in massive and long-run numerical experiments. It is shown that the methods are naturally subdivided in classes applicable for adiabatic and nonadiabatic cases of perturbation. It is explicitly shown that the pendulum approximation of marginal resonance works well in the nonadiabatic case. In this case, the role of marginal resonances in determining the total layer width is demonstrated to diminish with increasing main parameter $\lambda$ (equal to the ratio of the perturbation frequency to the frequency of small-amplitude phase oscillations on the resonance). Solely the “bending effect” is important in determining the total amplitude of the energy deviations of the near-separatrix motion at $\lambda \gtrsim 7$. In the adiabatic case, it is demonstrated that the geometrical form of the separatrix cell can be described analytically quite easily by means of using a specific representation of the separatrix map. It is shown that the nonadiabatic (and, to some extent, intermediary) case is most actual, in comparison with the adiabatic one, for the physical or technical applications that concern the energy jumps in the near-separatrix chaotic motion.

Time-dependent Poiseuille flow of a viscous compressible fluid confined between two planar walls is studied for a partial slip boundary condition with frequency-dependent slip length. After an initial uniform impulse parallel to the walls, the flow pattern quickly becomes nearly parabolic. For a narrow gap, a dynamic slip length can lead to damped oscillations of total fluid momentum.

We present an explicit analytic form for the two-breather solution of the nonlinear Schrödinger equation with imaginary eigenvalues. It describes various nonlinear combinations of Akhmediev breathers and Kuznetsov-Ma solitons. The degenerate case, when the two eigenvalues coincide, is quite involved. The standard inverse scattering technique does not generally provide an answer to this scenario. We show here that the solution can still be found as a special limit of the general second-order expression and appears as a mixture of polynomials with trigonometric and hyperbolic functions. A further restriction of this particular case, where the two eigenvalues are equal to $i$, produces the second-order rogue wave with two free parameters considered as differential shifts. The illustrations reveal a precarious dependence of wave profile on the degenerate eigenvalues and differential shifts. Thus we establish a hierarchy of second-order solutions, revealing the interrelated nature of the general case, the rogue wave, and the degenerate breathers.

The discrete system of equations for a chain consisting of a large number of spheres interacting via the Hertz force of index 3/2 in strain is examined in the very long wavelength limit, yielding an effective medium description. The resulting continuum second-order equation of motion possesses a subset of simple waves obeying a first-order equation of reduced index 5/4. These simple waves appear not to have examined before. For a given initial strain, the simple wave solution prescribes initial sphere centroid velocities. Together the initial strain and velocities are used in the second-order discrete system. Results for shock wave development compare very well between the second-order discrete system (minus physically valid oscillations) and the reduced first-order equation. A second-order simulation of colliding waves examines the ability of waves to pass through each other, with a phase advance accruing during the collision process. An arbitrary initial condition is shown to evolve toward a universal similarity solution proportional to ${(x/t)}^4$. A closed-form solution is given including the complete history of the waveform, shock location, and amplitude.

We consider the problem of homogenizing the Maxwell equations for periodic composites. The analysis is based on Bloch-Floquet theory. We calculate explicitly the reflection coefficient for a half space and derive and implement a computationally efficient continued-fraction expansion for the effective permittivity. Our results are illustrated by numerical computations for the case of two-dimensional systems. The homogenization theory of this paper is designed to predict various physically measurable quantities rather than to simply approximate certain coefficients in a partial differential equation.

Computations of the self-diffusion coefficient and viscosity in warm dense matter are presented with an emphasis on obtaining numerical convergence and a careful evaluation of the standard deviation. The transport coefficients are computed with the Green-Kubo relation and orbital-free molecular dynamics at the Thomas-Fermi-Dirac level. The numerical parameters are varied until the Green-Kubo integral is equal to a constant in the $t\to +\infty$ limit; the transport coefficients are deduced from this constant and not by extrapolation of the Green-Kubo integral. The latter method, which gives rise to an unknown error, is tested for the computation of viscosity; it appears that it should be used with caution. In the large domain of coupling constant considered, both the self-diffusion coefficient and viscosity turn out to be well approximated by simple analytical laws using a single effective atomic number calculated in the average-atom model.

We investigate the onset of “eigenstate thermalization” and the crossover to ergodicity in a system of one-dimensional fermions with increasing interaction. We analyze the fluctuations in the expectation values of most relevant few-body operators with respect to eigenstates. It turns out that these are intimately related to the inverse participation ratio of eigenstates displayed in the operator eigenbasis. Based on this observation, we find good evidence that eigenstate thermalization should set in even for vanishingly small perturbations in the thermodynamic limit.

The transport of $N$ globally coupled Brownian motors driven by a periodic force with phase disorder is investigated. An approximate theoretical analysis of the model is presented. The effects of the phase disorder and the driving strength of the periodic force on the transport of the coupled Brownian motors are discussed both theoretically and numerically. It is found that the increase of the periodical driving force decreases the average velocity, while the coupled particles may benefit from the phase disorder to enhance collective transport.

We study the energy minimization problem for an elastic interface in a random potential plus a quadratic well. As the position of the well is varied, the ground state undergoes jumps, called shocks or static avalanches. We introduce an efficient and systematic method to compute the statistics of avalanche sizes and manifold displacements. The tree-level calculation, i.e., mean-field limit, is obtained by solving a saddle-point equation. Graphically, it can be interpreted as the sum of all tree graphs. The 1-loop corrections are computed using results from the functional renormalization group. At the upper critical dimension the shock statistics is described by the Brownian force model (BFM), the static version of the so-called Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model in the nonequilibrium context of depinning. This model can itself be treated exactly in any dimension and its shock statistics is that of a Lévy process. Contact is made with classical results in probability theory on the Burgers equation with Brownian initial conditions. In particular we obtain a functional extension of an evolution equation introduced by Carraro and Duchon, which recursively constructs the tree diagrams in the field theory.

We analyze a recent experiment [Boudet, Cassagne, and Kellay, Phys. Rev. Lett. 103, 224501 (2009)] in which the shock created by the impact of a steel ball on a flowing monolayer of glass beads is studied quantitatively. We argue that radial momentum is conserved in the process and hence show that in two dimensions the shock radius increases in time $t$ as a power law $t^{1/3}$. This is confirmed in event driven simulations of an inelastic hard sphere system. The experimental data are compared with the theoretical prediction and are shown to compare well at intermediate times. At long times the experimental data exhibit a crossover to a different scaling behavior. We attribute this to the problem becoming effectively three dimensional due to the accumulation of particles at the shock front and propose a simple hard sphere model that incorporates this effect. Simulations of this model capture the crossover seen in the experimental data.

The diffusive growth rate of a polyhedral cell in dry three-dimensional foams depends on details of shape beyond cell topology, in contrast to the situation in two dimensions, where, by von Neumann's law, the growth rate depends only on the number of cell edges. We analyze the dependence of the instantaneous growth rate on the shape of single foam cells surrounded by uniform pressure; this is accomplished by supporting the cell with films connected to a wire frame and inducing cell distortions by deforming the wire frame. We consider three foam cells with a very simple topology; these are the Platonic foam cells, which satisfy Plateau's laws and are based on the trivalent Platonic solids (tetrahedron, cube, and dodecahedron). The Surface Evolver is used to model cell deformations induced through extension, compression, shear, and torsion of the wire frames. The growth rate depends on the deformation mode and frame size and can increase or decrease with increasing cell distortion. The cells have negative growth rates, in general, but dodecahedral cells subjected to torsion in small wire frames can have positive growth rates. The deformation of cubic cells is demonstrated experimentally.

We experimentally study the motion of optically driven colloidal particles on a circular path by varying their number $N$. Although an identical driving force is applied to each particle, their equally spaced configuration is hydrodynamically unstable, and a doublet configuration is spontaneously formed. In small-$N$ systems, the angular difference between neighboring particles exhibits oscillatory or nonoscillatory behavior. The number of oscillatory modes that appear depends on the maximum number of doublets that the system can contain. Frequent switching between different modes was observed with increasing $N$. The characteristic frequencies of the oscillatory modes are discussed theoretically by linear stability analysis of the equations that govern the motion of hydrodynamically coupled particles. The evaluated frequencies of the slowest modes exhibit reasonably good agreement with those of the mainly observed modes in experiments. The relationship between the characteristic frequencies and specific configurations is confirmed experimentally by setting a specific initial configuration for the particles. An increase in $N$ also enhances the mean angular velocity of the particles owing to the reduced effective viscosity in large-$N$ systems.

The wetting transition of the Blume-Capel model is studied by a finite-size scaling analysis of $L\times M$ lattices where competing boundary fields $\pm H_1$ act on the first or last row of the $L$ rows in the strip, respectively. We show that using the appropriate anisotropic version of finite-size scaling, critical wetting in $d=2$ is equivalent to a “bulk” critical phenomenon with exponents $\alpha =-1$, $\beta =0$, and $\gamma =3$. These concepts are also verified for the Ising model. For the Blume-Capel model, it is found that the field strength $H_{1c}\left(T\right)$ where critical wetting occurs goes to zero when the bulk second-order transition is approached, while $H_{1c}\left(T\right)$ stays nonzero in the region where in the bulk a first-order transition from the ordered phase, with nonzero spontaneous magnetization, to the disordered phase occurs. Interfaces between coexisting phases then show interfacial enrichment of a layer of the disordered phase which exhibits in the second-order case a finite thickness only. A tentative discussion of the scaling behavior of the wetting phase diagram near the tricritical point is also given.