The concept of coherency loss is proposed to understand wrinkle branching as a pathway toward hierarchical wrinkling pattern formation in a compressed film-substrate system with gradient stiffness of the film or substrate. A simple model indicates that the wrinkle branching arises when the characteristic length of the stiffness inhomogeneity zone is larger than the coherency persistent length, which depends on the amplitude of the stiffness inhomogeneity. Numerical simulations of nonlinear wrinkles based on the model of the Föppl–von Kármán plate on compliant substrates show how regulating the size and amplitude of the stiffness inhomogeneities results in branched wrinkles in striking agreement with the existing observations. The paper reveals the origin of such kinds of branched wrinkles and may provide a guideline for controllable hierarchical wrinkles by patterning the stiffness gradient.

We study the propulsion of a one-dimensional (1D) polymer chain under sinusoidal external forces in the overdamped (low Reynolds number) regime. We show that, when hydrodynamical interactions are included, the polymer presents directional motion which depends on the phase differences of the external force applied along the chain. Moreover, the velocity shows a maximum as a function of the frequency. We discuss the relevance of all these results in light of recent nanotechnology experiments.

We present an analytical study of the heat conduction in a quantum harmonic chain of oscillators with alternate particle masses. We find that a thermal-insulating effect already observed in some classical models consisting of chains with alternate masses also holds in the quantum version of the one-dimensional harmonic model with thermal reservoirs at the boundaries. Namely, we show that if the heat conduction is to be avoided, then alternating small and large masses in a wire is more effective than having all of the masses equally large. The presence of a similar effect in different models indicates that the phenomenon may be rather general, with possible applications in the heat-flow control.

With a simple phenomenological metapopulation model, which characterizes the invasion process of an influenza pandemic from a source to a subpopulation at risk, we compare the efficiency of inter- and intrapopulation interventions in delaying the arrival of an influenza pandemic. We take travel restriction and patient isolation as examples, since in reality they are typical control measures implemented at the inter- and intrapopulation levels, respectively. We find that the intrapopulation interventions, e.g., patient isolation, perform better than the interpopulation strategies such as travel restriction if the response time is small. However, intrapopulation strategies are sensitive to the increase of the response time, which might be inevitable due to socioeconomic reasons in practice and will largely discount the efficiency.

We present a phenomenological dynamical model able to describe the stretching features of the curve of DNA length vs applied force. As concerns the chain, the model is based on the discrete wormlike chain model with elastic modifications, which properly describes the elongation features at low and intermediate forces. The dynamics is developed under a double-well potential with a linear term, which, at high forces, accounts for the narrow transition present in the DNA elongation (overstretching). A quite good agreement between simulation and experiment is obtained.

We propose an efficient strategy to suppress epidemic explosion in heterogeneous metapopulation networks, wherein each node represents a subpopulation with any number of individuals and is assigned a curing rate that is proportional to $k^{\alpha }$ with the node degree $k$ and an adjustable parameter $\alpha$. We perform stochastic simulations of the dynamical reaction-diffusion processes associated with the susceptible-infected-susceptible model in scale-free networks. We find that the epidemic threshold reaches a maximum when $\alpha$ is tuned at ${\alpha }_{\mathrm{opt}}\simeq 1.3$. This nontrivial phenomenon is robust to the change of the network size and the average degree. In addition, we carry out a mean field analysis to further validate our scheme, which also demonstrates that epidemic explosion follows different routes for $\alpha$ larger or less than ${\alpha }_{\mathrm{opt}}$. Our work suggests that in order to efficiently suppress epidemic spreading on heterogeneous complex networks, subpopulations with higher degrees should be allocated more resources than just being linearly dependent on $k$.

We study the dynamics of the classical and quantum mechanical scattering of a wave packet from an oscillating barrier. Our main focus is on the dependence of the transmission coefficient on the initial energy of the wave packet for a wide range of oscillation frequencies. The behavior of the quantum transmission coefficient is affected by tunneling phenomena, resonances, and kinematic effects emanating from the time dependence of the potential. We show that when kinematic effects dominate (mainly in intermediate frequencies), classical mechanics provides very good approximation of quantum results. In that frequency region, the classical and quantum transmission coefficients are in optimal agreement. Moreover, the transmission threshold (i.e., the energy above which the transmission coefficient becomes larger than a specific small threshold value) is found to exhibit a minimum. We also consider the form of the transmitted wave packet and we find that for low values of the frequency the incoming classical and quantum wave packet can be split into a train of well-separated coherent pulses, a phenomenon that admits purely classical kinematic interpretation.

It has been recently demonstrated that self-defocusing (SDF) media with cubic nonlinearity, whose local coefficient grows from the center to the periphery fast enough, support stable bright solitons without the use of any linear potential. Our objective is to test the genericity of this mechanism for other nonlinearities, by applying it to one- and two-dimensional (1D and 2D) quintic SDF media. The models may be implemented in optics (in particular, in colloidal suspensions of nanoparticles), and the 1D model may be applied to the description of the Tonks-Girardeau gas of ultracold bosons. In 1D, the nonlinearity-modulation function is taken as $g_0+{\sinh }^2\left(\beta x\right)$. This model admits a subfamily of exact solutions for fundamental solitons. Generic soliton solutions are constructed in a numerical form and also by means of the Thomas-Fermi and variational approximations (TFA and VA). In particular, a new ansatz for the VA is proposed, in the form of “raised $\mathrm{sech},$” which provides for an essentially better accuracy than the usual Gaussian ansatz. The stability of all the fundamental (nodeless) 1D solitons is established through the computation of the corresponding eigenvalues for small perturbations and also verified by direct simulations. Higher-order 1D solitons with two nodes have a limited stability region, all the modes with more than two nodes being unstable. It is concluded that the recently proposed inverted Vakhitov-Kolokolov stability criterion for fundamental bright solitons in systems with SDF nonlinearities holds here too. Particular exact solutions for 2D solitons are produced as well.

Ab initio potentials are implemented into the direct simulation Monte Carlo (DSMC) method. Such an implementation allows us to model transport phenomena in rarefied gases without any fitting parameter of intermolecular collisions usually extracted from experimental data. Applying the method proposed by Sharipov and Strapasson [Phys. Fluids 24, 011703 (2012)], the use of ab initio potentials in the DSMC requires the same computational efforts as the widely used potentials such as hard spheres, variable hard sphere, variable soft spheres, etc. At the same time, the ab initio potentials provide more reliable results than any other one. As an example, the transport coefficients of a binary mixture He-Ar, viz., viscosity, thermal conductivity, and thermal diffusion factor, have been calculated for several values of the mole fraction.

The study of the effects of spatially uniform fields on the steady-state properties of Axelrod's model has yielded plenty of counterintuitive results. Here, we reexamine the impact of this type of field for a selection of parameters such that the field-free steady state of the model is heterogeneous or multicultural. Analyses of both one- and two-dimensional versions of Axelrod's model indicate that the steady state remains heterogeneous regardless of the value of the field strength. Turning on the field leads to a discontinuous decrease on the number of cultural domains, which we argue is due to the instability of zero-field heterogeneous absorbing configurations. We find, however, that spatially nonuniform fields that implement a consensus rule among the neighborhood of the agents enforce homogenization. Although the overall effects of the fields are essentially the same irrespective of the dimensionality of the model, we argue that the dimensionality has a significant impact on the stability of the field-free homogeneous steady state.

The present study describes wave propagation characteristics in a weakly disordered two-dimensional granular media composed of a square array of spheres accommodating interstitial cylindrical intruders. Previous investigations, performed experimentally as well as numerically, emphasized that wave-front shapes in similar systems are tunable via choice of material combinations. Here, we investigate the effects of statistical variation in the particle diameters and compare the effects of the resulting disorder in experiments and numerical simulations, finding good agreement.

We determine the speed of a crystallization (or, more generally, a solidification) front as it advances into the uniform liquid phase after the system has been quenched into the crystalline region of the phase diagram. We calculate the front speed by assuming a dynamical density functional theory (DDFT) model for the system and applying a marginal stability criterion. Our results also apply to phase field crystal (PFC) models of solidification. As the solidification front advances into the unstable liquid phase, the density profile behind the advancing front develops density modulations and the wavelength of these modulations is a dynamically chosen quantity. For shallow quenches, the selected wavelength is precisely that of the crystalline phase and so well-ordered crystalline states are formed. However, when the system is deeply quenched, we find that this wavelength can be quite different from that of the crystal, so the solidification front naturally generates disorder in the system. Significant rearrangement and aging must subsequently occur for the system to form the regular well-ordered crystal that corresponds to the free energy minimum. Additional disorder is introduced whenever a front develops from random initial conditions. We illustrate these findings with simulation results obtained using the PFC model.

We investigate experimentally the spatial distributions of heavy and neutrally buoyant particles of finite size in a fully turbulent flow. Because their Stokes number (i.e., the ratio of the particle viscous relaxation time to a typical flow time scale) is close to unity, one may expect both classes of particles to aggregate in specific flow regions. This is not observed. Using a Voronoï analysis we show that neutrally buoyant particles sample turbulence homogeneously, whereas heavy particles do cluster. These results show that several dimensionless numbers are needed in the modeling (and understanding) of the behavior of particles entrained by turbulent motions.

Parameter space of a driven damped oscillator in a double well potential presents either a chaotic trajectory with sign oscillating amplitude or a nonchaotic trajectory with a fixed sign amplitude. A network of such delay coupled damped oscillators is shown to present chaotic dynamics while the sign amplitude of each damped oscillator is randomly frozen. This phenomenon of random broken global symmetry of the network simultaneous with random freezing of each degree of freedom is accompanied by the existence of exponentially many randomly frozen chaotic attractors with the size of the network. Results are exemplified by a network of modified Duffing oscillators with infinite range pseudoinverse delayed interactions.

Combining experiments and theory, we investigate the break-up dynamics of deformable objects, such as drops and bubbles, against a linear micro-obstacle. Our experiments bring the role of the viscosity contrast $\Delta \eta$ between dispersed and continuous phases to light: the evolution of the critical capillary number to break a drop as a function of its size is either nonmonotonic ($\Delta \eta >0$) or monotonic ($\Delta \eta \le 0$). In the case of positive viscosity contrasts, experiments and modeling reveal the existence of an unexpected critical object size for which the critical capillary number for breakup is minimum. Using simple physical arguments, we derive a model that well describes observations, provides diagrams mapping the four hydrodynamic regimes identified experimentally, and demonstrates that the critical size originating from confinement solely depends on geometrical parameters of the obstacle.

Different nonkinetic approaches are adopted in this paper towards theoretically predicting the experimentally observed phenomenon of enhanced mass flow rates accompanying pressure-driven rarefied gas flows through microchannels. Our analysis utilizes a full set of mechanically consistent volume-diffusion hydrodynamic equations, allowing complete, closed-form, analytical solutions to this class of problems. As an integral part of the analysis, existing experimental data pertaining to the subatmospheric pressure dependence of viscosity were analyzed. The several nonkinetic approaches investigated were (1) pressure-dependent viscosity exponent model, (2) slip-velocity models, and (3) volume diffusion model. We explored the ability to predict the gas's mass flow rate over the full range of Knudsen numbers, including furnishing a physically sound interpretation of the well-known Knudsen minimum observed in the mass flow rate. Matching of a pressure-dependent viscosity model, one that follows the standard temperature-viscosity power law and its supporting single momentum diffusion mechanism, did not allow an accurate interpretation of the data. Rather, matching of this model with the flow rate was found to mismatch the experimental pressure dependence of the viscosity. An additional transport mechanism model, one based on volume diffusion, offered a comprehensive understanding of the Knudsen minimum, while also resulting in excellent agreement with experimental data well into the transition regime (up to a Knudsen number of 5).

Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that are representative of quasi-two-dimensional behavior of the flow. It is known that when energy is injected into the fluid at an intermediate scale $L_f$, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the inverse cascade range at a small but fixed Rossby number, ${\mathrm{Ro}}_f\approx 0.05$. Several numerical simulations with helical and nonhelical forcing functions are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade range with reasonably large Reynolds number, the analysis is based on large eddy simulations which include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of energy in two-dimensional (2D) turbulence with a $\sim k_{\perp }^{-5/3}$ scaling, and the other that corresponds to a steeper $\sim k_{\perp }^{-3}$ spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to the 2D modes. The spectrum that emerges depends on the anisotropy of the forcing function, the former solution prevailing for forcings in which more energy is injected into the 2D modes while the latter prevails for isotropic forcing. In the case of anisotropic forcing, whence the energy goes from the 2D to the 3D modes at low wave numbers, large-scale shear is created, resulting in a time scale ${\tau }_{sh}$, associated with shear, thereby producing a $\sim k^{-1}$ spectrum for the total energy with the horizontal energy of the 2D modes still following a $\sim k_{\perp }^{-5/3}$ scaling.

It is shown that the dynamics of amplitude-modulated compressional dispersive Alfvénic (CDA) waves in a collisional megnetoplasma is governed by a complex Ginzburg-Landau (CGL) equation. The nonlinear dispersion relation for the modulational instability of the CDA waves is derived and investigated numerically. It is found that the growth rate of the modulational instability decreases (increases) with the increase of the normalized electron-ion collision frequency $\alpha$ (the plasma $\beta$). The modulational instability criterion for the CGL equation is defined precisely and investigated numerically. The region of the modulational instability becomes narrower with the increase of $\alpha$ and $\beta$, indicating that the system dissipates the wave energy by collisions, and a stable CDA wave envelope packet in the form of a hole will be a dominant localized pulse. For a collisionless plasma, i.e., $\alpha =0$, the CGL equation reduces to the standard nonlinear Schrödinger (NLS) equation. The latter is used to investigate the modulational (in)stability region for the CDA waves in a collisionless magnetoplasma. It is shown that, within unstable regions, a random set of nonlinearly interacting CDA perturbations leads to the formation of CDA rogue waves. In order to demonstrate that the characteristics of the CDA rogue waves are influenced by the plasma $\beta$, the relevant numerical analysis of the appropriate nonlinear solution of the NLS equation is presented. The application of our investigation to space and laboratory magnetoplasmas is discussed.

According to Harada and Sasa [Phys. Rev. Lett. 95, 130602 (2005)], heat production generated in a nonequilibrium steady state can be inferred from measuring response and correlation functions. In many colloidal systems, however, it is a nontrivial task to determine response functions, whereas details about spatial steady state trajectories are easily accessible. Using a simple conditional averaging procedure, we show how this fact can be exploited to reliably evaluate average heat production. We test this method using Brownian dynamics simulations, and apply it to experimental data of an interacting driven colloidal system.

We study the effects of scattering lengths on Lévy walks in quenched, one-dimensional random and fractal quasilattices, with scatterers spaced according to a long-tailed distribution. By analyzing the scaling properties of the random-walk probability distribution, we show that the effect of the varying scattering length can be reabsorbed in the multiplicative coefficient of the scaling length. This leads to a superscaling behavior, where the dy-namical exponents and also the scaling functions do not depend on the value of the scattering length. Within the scaling framework, we obtain an exact expression for the multiplicative coefficient as a function of the scattering length both in the annealed and in the quenched random and fractal cases. Our analytic results are compared with numerical simulations, with excellent agreement, and are supposed to hold also in higher dimensions.