We consider subdiffusion of a particle in a one-dimensional system with a thin partially permeable and partially absorbing wall. The system with the wall can be used to filter diffusing particles. Passing through the wall, the particle can be absorbed with a certain probability. Knowing the Green's functions we derive boundary conditions at the wall. The boundary conditions take a specific form in which fractional time derivatives are involved. The temporal evolution of the probability that a diffusing particle has not been absorbed is also considered.

In this paper we have investigated through the numerical solution of the basic equation as well as through the dynamic model the influence of higher-order correction terms to the nonlinear amplification (absorption) and to the nonlinear refractive index on the self-frequency shift of Raman dissipative solitons. We have found a nonlinear dependence of the self-frequency shift of Raman dissipative solitons on the parameter describing intrapulse Raman scattering in the presence of the saturation of the nonlinear gain. With the increase of the absolute value of the saturation of the nonlinear gain, the maximum absolute value of the frequency shift decreases and its position moves to larger values of the parameter describing intrapulse Raman scattering. The increase in the value of the nonlinear gain leads to an increase in the maximum absolute value of the frequency shift, without changing its position. We have also observed the nonlinear dependence of the absolute value of the frequency shift on the parameter describing intrapulse Raman scattering in the presence of higher-order correction term to the nonlinear refractive index. The discovered nonlinear dependence of the self-frequency shift on the value of the saturation of the nonlinear gain as well as on the higher-order correction term to the nonlinear refractive index can be used for the better understanding and control of the spectral characteristics of Raman dissipative solitons. The dynamic model correctly describes all the features of the observed phenomena.

Although many previous studies have focused on the Brazil nut effect, segregation in a self-gravitating circular aggregate remains relatively unexplored. In this paper, size segregation in a two-dimensional assembly of grains in a circular geometry is studied through discrete element method (DEM) numerical simulations. We show that radial segregation within an asteroid submitted to periodic perturbations is not limited to the surface but also occurs in its core. The characteristic time and the overall efficiency of the segregation mechanism are studied as the intensity of the perturbation, the frictional properties, and rotational freedom of individual grains are varied.

Quantitative, three-dimensional (3D) spatially resolved magnetic resonance flow imaging (flow MRI) methods are presented to characterize structure-flow correlations in a 4-mm-diameter plug of Ketton limestone rock using undersampled k- and q-space data acquisition methods combined with compressed sensing (CS) data reconstruction techniques. The acquired MRI data are coregistered with an X-ray microcomputed tomography (μCT) image of the same rock sample, allowing direct correlation of the structural features of the rock with local fluid transport characteristics. First, 3D velocity maps acquired at 35 μm isotropic spatial resolution showed that the flow was highly heterogeneous, with ∼10% of the pores carrying more than 50% of the flow. Structure-flow correlations were found between the local flow velocities through pores and the size and topology (coordination number) associated with these pores. These data show consistent trends with analogous data acquired for flow through a packing of 4-mm-diameter spheres, which may be due to the microstructure of Ketton rock being a consolidation of approximately spherical grains. Using two-dimensional and 3D visualization of coregistered μCT images and velocity maps, complex pore-scale flow patterns were identified. Second, 3D spatially resolved propagators were acquired at 94 μm isotropic spatial resolution. Flow dispersion within the rock was examined by analyzing each of the 331 776 local propagators as a function of observation time. Again, the heterogeneity of flow within the rock was shown. Quantification of the mean and standard deviation of each of the local propagators showed enhanced mixing occurring within the pore space at longer observation times. These spatially resolved measurements also enable investigation of the length scale of a representative elementary volume. It is shown that for a 4-mm-diameter plug this length scale is not reached.

The transport properties of the weakly nonlinear (WNL) two-dimensional (2D) quasilongitudinal dust lattice mode is studied in an experimentally realized highly viscous, strongly coupled, weakly ionized plasma [V. E. Fortov et al., Phys. Rev. Lett. 109, 055002 (2012)]. The WNL dynamics is found to be described by a 2D dissipative-dispersive nonlinear partial differential equation. The analytical and computational (for gas discharge plasma parameters) results predict strong viscosity induced Shilnikov homoclinic chaos, which, in turn, can cause a phase transition.

The pore size distribution of a porous medium is often estimated from the retention curve or the invading fluid flow rate curve using simple relationships more or less explicitly based on the consideration that the porous medium is made of a bundle of cylindrical parallel tubes. This type of determination is tested using pore network simulations. Starting from two- or three-dimensional networks, the characteristics of which are known a priori, the estimation of the throat size distribution (TSD) is performed using the standard methods in the case of drainage. Results show a significant discrepancy with the input data. The disagreement is more pronounced when the fluid flow rate curve is employed together with the parallel tubes assumption. The physical origins of these shortcomings are identified. A method, based on pore network simulations combined with a genetic algorithm and the hill climbing algorithm, is then designed, which makes simultaneous use of the nonwetting fluid flow rate curve and the retention curve of the medium. Very significant improvement is achieved in the estimation of the TSD using this procedure.

Diffusion-limited aggregation (DLA) has served for 40 years as a paradigmatic example for the creation of fractal growth patterns. In spite of thousands of references, no exact result for the fractal dimension $D$ of DLA is known. In this Letter we announce an exact result for off-lattice DLA grown on a line embedded in the plane $D=3/2$. The result relies on representing DLA with iterated conformal maps, allowing one to prove self-affinity, a proper scaling limit, and a well-defined fractal dimension. Mathematical proofs of the main results are available in N. Berger, E. B. Procaccia, and A. Turner, Growth of stationary Hastings-Levitov, arXiv:2008.05792.

The environment of an open quantum system is usually modelled as a large many-body quantum system. However, when an isolated quantum system itself is a many-body quantum system, the question of how large and complex it must be to generate internal equilibration is an open key-point in the literature. In this work, by monitoring the degree of equilibration of a single spin through its purity degradation, we are able to sense the chaotic behavior of the generic spin chain to which it is coupled. Quite remarkably, this holds even in the case of extremely short spin chains composed of three spins, where we can also reproduce the whole integrable to chaos transition. Finally, we discuss implications on quantum control experiments and show that quantum chaos reigns over the best degree of control achieved, even in small chains.

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modeling approach for the fluid phases, we derive the governing equations for the dynamics of the ternary fluid flow and membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method for the ternary fluid (Eulerian description) and a finite difference representation of the membrane (Lagrangian description). The ternary fluid and membrane solvers are coupled through the immersed boundary method. For validation purposes, we consider the relaxation dynamics of a two-dimensional elastic capsule placed at a fluid-fluid interface. The capsule shapes, resulting from the balance of surface tension and elastic forces, are compared with equilibrium numerical solutions obtained by surface evolver. Furthermore, the Galilean invariance of the proposed model is proven. The proposed approach is versatile, allowing for the simulation of a wide range of geometries. To demonstrate this, we address the problem of a capillary bridge formed between two deformable capsules.

Understanding the global dynamical behavior of a network of coupled oscillators has been a topic of immense research in many fields of science and engineering. Various factors govern the resulting dynamical behavior of such networks, including the number of oscillators and their coupling schemes. Although these factors are seldom significant in large populations, a small change in them can drastically affect the global behavior in small populations. In this paper, we perform an experimental investigation on the effect of these factors on the coupled behavior of a minimal network of candle-flame oscillators. We observe that strongly coupled oscillators exhibit the global behavior of in-phase synchrony and amplitude death, irrespective of the number and the topology of oscillators. However, when they are weakly coupled, their global behavior exhibits the intermittent occurrence of multiple stable states in time. We report the experimental discovery of partial amplitude death in a network of candle-flame oscillators, in addition to the observation of other dynamical states including clustering, chimera, and weak chimera. We also show that closed-loop networks tend to hold global synchronization for longer duration as compared to open-loop networks.

When a network is inferred from data, two types of errors can occur: false positive and false negative conclusions about the presence of links. We focus on the influence of local network characteristics on the probability $\alpha$ of false positive conclusions, and on the probability $\beta$ of false negative conclusions, in the case of networks of coupled oscillators. We demonstrate that false conclusion probabilities are influenced by local connectivity measures such as the shortest path length and the detour degree, which can also be estimated from the inferred network when the true underlying network is not known $a$ priori. These measures can then be used for quantification of the confidence level of link conclusions, and for improving the network reconstruction via advanced concepts of link weights thresholding.

Optimizing traffic flow is essential for easing congestion. However, even when globally optimal, coordinated, and individualized routes are provided, users may choose alternative routes which offer lower individual costs. By analyzing the impact of selfish route choices on performance using the cavity method, we find that a small ratio of selfish route choices improves the global performance of uncoordinated transportation networks but degrades the efficiency of optimized systems. Remarkably, compliant users always gain in the former and selfish users may gain in the latter, under some parameter conditions. The theoretical results are in good agreement with large-scale simulations. Iterative route switching by a small fraction of selfish users leads to Nash equilibria close to the globally optimal routing solution. Our theoretical framework also generalizes the use of the cavity method, originally developed for the study of equilibrium states, to analyze iterative game-theoretical problems. These results shed light on the feasibility of easing congestion by route coordination when not all vehicles follow the coordinated routes.

The dynamics of networks on Heider's balance theory moves toward reducing the tension by constantly reevaluating the interactions to achieve a state of balance. Conflict of interest, however, is inherent in most complex systems; frequently, there are multiple ideals or states of balance, and moving towards one could work against another. In this paper, by introducing the competitive balance theory, we study the evolution of balance in the presence of conflicts of interest. In our model, the assumption is that different states of balance compete in the evolution process to dominate the system. We ask whether, through these interactions, different states of balance compete to prevail their own ideals or a set of coexisting ideals in a balanced condition is a possible outcome. The results show that although there is a symmetry in the type of balance the system either evolves towards a symmetry breaking, where one of the states of balance dominates the system, or, less frequently, the competing states of balance coexist in a jammed state.

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions—referred to as “glass stickiness”—can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, negative stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history dependence and what it is possibly affected by are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles—referred to here as glass “stickiness.” We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses; the latter are therefore deemed “thermomechanically inannealable.” The implications of our findings and future research directions are discussed.

Organic compounds with bent-core (BC) molecules usually form the layered smectic liquid crystals with tilted molecules and polarization (P) which lies in the plane of the layer. A few such compounds have been found in which P itself tilts out of the plane of the layer, and the medium with general tilt ($\mathrm{Sm}{C}_g$) of the molecules has the low chiral triclinic symmetry. We discuss the geometric constraints of molecular packing in this structure to show that projecting groups attached to one of the arms of the BC molecules favors the formation of the $\mathrm{Sm}{C}_g$ phase. We also extend our model for the modulated phases exhibited by BC molecules to show that the stripe structure made of bilayers shown by a few BC compounds, which is a signature of $\mathrm{Sm}{C}_g$ layers, is preferentially formed by rotation of the BC molecules about their long axes, rather than about the layer normal. The theoretical results based on the unified model which describes all the modulated phases exhibited by the BC molecules broadly reflect experimental trends.

The phase behavior and adsorption kinetics of hard-core particles on a honeycomb lattice are studied by means of random sequential adsorption with surface diffusion. We concentrate on reversible adsorption by introducing a desorption process into our previous model and varying the equilibrium rate constant as a control parameter. We find that an exact prediction of the temporal evolution of fractional surface coverage and the surface pressure dynamics of reversible adsorption can be achieved by use of the blocking function of a system with irreversible adsorption of highly mobile particles. For systems out of equilibrium we observe several features of glassy dynamics, such as slow relaxation dynamics, the memory effect, and aging. In particular, the analysis of our system in the limit of small desorption probability shows simple aging behavior with a power-law decay. A detailed discussion of Gibbs adsorption isotherm for nonequilibrium adsorption is given, which exhibits a hysteresis between this system and its equilibrium counterpart.

Euglena gracilis is a unicellular organism that swims by beating a single anterior flagellum. We study the nonplanar waveforms spanned by the flagellum during a swimming stroke and the three-dimensional flows that they generate in the surrounding fluid. Starting from a small set of time-indexed images obtained by optical microscopy on a swimming Euglena cell, we construct a numerical interpolation of the stroke. We define an optimal interpolation (which we call synthetic stroke) by minimizing the discrepancy between experimentally measured velocities (of the swimmer) and those computed by solving numerically the equations of motion of the swimmer driven by the trial interpolated stroke. The good match we obtain between experimentally measured and numerically computed trajectories provides a first validation of our synthetic stroke. We further validate the procedure by studying the flow velocities induced in the surrounding fluid. We compare the experimentally measured flow fields with the corresponding quantities computed by solving numerically the Stokes equations for the fluid flow, in which the forcing is provided by the synthetic stroke, and find good matching. Finally, we use the synthetic stroke to derive a coarse-grained model of the flow field resolved in terms of a few dominant singularities. The far field is well approximated by a time-varying Stresslet, and we show that the average behavior of Euglena during one stroke is that of an off-axis puller. The reconstruction of the flow field closer to the swimmer body requires a more complex system of singularities. A system of two Stokeslets and one Rotlet, that can be loosely associated with the force exerted by the flagellum, the drag of the body, and a torque to guarantee rotational equilibrium, provides a good approximation.

This paper reports on the mechanism of the hysteresis in the transition between regular and Mach shock wave reflections. We disclose that, for a given inflow Mach number, a stable reflection configuration should maintain the minimal dissipation. As the wedge angle varies, the set of the minimal dissipation points forms the valley lines in the dissipation landscape, and these valley lines compose the hysteresis loop. The saddle-nodes, intersections of the ridge line, and the valley lines are actually the transition points. Additionally, the predicted reflection configurations agree well with the experimental and numerical results, validating this theory.

It is shown theoretically and by simulation that a Gaussian laser beam of relativistic intensity interacting with a uniform-thickness plasma slab of azimuthally varying density can acquire orbital angular momentum (OAM). During the interaction, the laser ponderomotive force and the charge-separation force impose a torque on the plasma particles. The affected laser light and plasma ions gain oppositely directed axial OAM, but the plasma electrons remain almost OAM free. High OAM conversion efficiency is achieved due to the strong azimuthal electromagnetic energy flow during the laser phase modulation. The present scheme should provide useful reference for applications requiring relativistic-intense vortex light.