We study numerically and analytically the instabilities associated with phase separation in a solid layer on which an external material flux is imposed. The first instability is localized within a boundary layer at the exposed free surface by a process akin to spinodal decomposition. In the limiting static case, when there is no material flux, the coherent spinodal decomposition is recovered. In the present problem, stability analysis of the time-dependent and nonuniform base states as well as numerical simulations of the full governing equations are used to establish the dependence of the wavelength and onset of the instability on parameter settings and its transient nature as the patterns eventually coarsen into a flat moving front. The second instability is related to the Mullins-Sekerka instability in the presence of elasticity and arises at the moving front between the two phases when the flux is reversed. Stability analyses of the full model and the corresponding sharp-interface model are carried out and compared. Our results demonstrate how interface and bulk instabilities can be analyzed within the same framework which allows us to identify and distinguish each of them clearly. The relevance for a detailed understanding of both instabilities and their interconnections in a realistic setting is demonstrated for a system of equations modeling the lithiation and delithiation processes within the context of lithium ion batteries.

In view of the results achieved in a previously related work [A. Altieri, S. Franz, and G. Parisi, J. Stat. Mech. (2016) 093301], regarding a Plefka-like expansion of the free energy up to the second order in the perceptron model, we improve the computation here focusing on the role of third-order corrections. The perceptron model is a simple example of constraint satisfaction problem, falling in the same universality class as hard spheres near jamming and hence allowing us to get exact results in high dimensions for more complex settings. Our method enables to define an effective potential (or Thouless-Anderson-Palmer free energy), namely a coarse-grained functional, which depends on the generalized forces and the effective gaps between particles. The analysis of the third-order corrections to the effective potential reveals that, albeit irrelevant in a mean-field framework in the thermodynamic limit, they might instead play a fundamental role in considering finite-size effects. We also study the typical behavior of generalized forces and we show that two kinds of corrections can occur. The first contribution arises since the system is analyzed at a finite distance from jamming, while the second one is due to finite-size corrections. We nevertheless show that third-order corrections in the perturbative expansion vanish in the jamming limit both for the potential and the generalized forces, in agreement with the isostaticity argument proposed by Wyart and coworkers. Finally, we analyze the relevant scaling solutions emerging close to the jamming line, which define a crossover regime connecting the control parameters of the model to an effective temperature.

The existence and exact form of the continuum expression of the discrete nonlogarithmic $q$-entropy is an important open problem in generalized thermostatistics, since its possible lack implies that nonlogarithmic $q$-entropy is irrelevant for the continuous classical systems. In this work, we show how the discrete nonlogarithmic $q$-entropy in fact converges in the continuous limit and the negative of the $q$-entropy with continuous variables is demonstrated to lead to the (Csiszár type) $q$-relative entropy just as the relation between the continuous Boltzmann-Gibbs expression and the Kullback-Leibler relative entropy. As a result, we conclude that there is no obstacle for the applicability of the $q$-entropy to the continuous classical physical systems.

The transport properties of a bosonic chain have been calculated by placing the ends of the chain in contact with thermal and particle reservoirs at different temperatures and chemical potentials. The contact with the reservoirs is described by the use of a quantum Fokker-Planck-Kramers equation, which is a canonical quantization of the classical Fokker-Planck-Kramers equation. From the quantum equation we obtain equations for the covariances of the creation and annihilation boson operators and solve them in the stationary state for small interactions. From the covariances we determine the Onsager coefficients and in particular the conductance, which was found to be finite for any chain size leading to an infinite conductivity and the absence of Fourier's law.

We investigate dynamic magneto-optic effects in a ferromagnetic nematic liquid crystal experimentally and theoretically. Experimentally we measure the magnetization and the phase difference of the transmitted light when an external magnetic field is applied. As a model we study the coupled dynamics of the magnetization, $\mathbf{M}$, and the director field, $\mathbf{n}$, associated with the liquid crystalline orientational order. We demonstrate that the experimentally studied macroscopic dynamic behavior reveals the importance of a dynamic cross-coupling between $\mathbf{M}$ and $\mathbf{n}$. The experimental data are used to extract the value of the dissipative cross-coupling coefficient. We also make concrete predictions about how reversible cross-coupling terms between the magnetization and the director could be detected experimentally by measurements of the transmitted light intensity as well as by analyzing the azimuthal angle of the magnetization and the director out of the plane spanned by the anchoring axis and the external magnetic field. We derive the eigenmodes of the coupled system and study their relaxation rates. We show that in the usual experimental setup used for measuring the relaxation rates of the splay-bend or twist-bend eigenmodes of a nematic liquid crystal one expects for a ferromagnetic nematic liquid crystal a mixture of at least two eigenmodes.

The response of a thin layer of granular material to an external pure shear imposed at its base is investigated. The experiments show that, even for noncohesive materials, the resulting deformation of the material is inhomogeneous. Indeed, a novel smooth pattern, consisting of a periodic modulation of the shear deformation of the free surface, is revealed by an image-correlation technique. These observations are in contrast with the previous observation of the fracture pattern in cohesive granular materials subjected to stretching. For cohesive materials, the instability is due to the weakening of the material which results from the rupture of capillary bridges that bond the grains to one another. For noncohesive materials, the rupture of the capillary bridges cannot be invoked anymore. We show that the instability results from the decrease of friction on shearing. PACS: 89.75.Kd: Pattern formation in complex systems; 83.60.Uv: Rheology: fracture; 45.70.Qj: Pattern formation in granular matter

We report an electrohydrodynamic instability in a premixed stagnation flame under manipulations of a dc electric field. This instability occurs when the electric field strength is at a certain value below the breakdown threshold, which is 0.75 kV/cm in the experimental setup. Above this value the flame front suddenly transits from a substrate-stabilized near-flat shape to a nozzle-stabilized conical surface, accompanied by a jump in the electric current through the flame field. At the transition moment, the flame spontaneously propagates upstream to the nozzle while the flow velocity at the upstream of the flame front decreases to zero, as revealed by high-speed photographs and PIV measurements. These phenomena indicate a transient balance between the fluid inertia and the electric body force around the instability threshold. A quantitative model suggests that the flame instability can be explained by a positive feedback loop, where the electric field applies a nonuniform electric body force, pulling the flame front upstream, and the pulled flame front in turn enhances the local electric body force. The electrohydrodynamic instability occurs when the electric pulling is strong enough and both the growth rates and the magnitudes of the electric body force on flame exceed those of the fluid dynamic pressure.

According to the acoustic fluidization hypothesis, elastic waves at a characteristic frequency form inside seismic faults even in the absence of an external perturbation. These waves are able to generate a normal stress which contrasts the confining pressure and promotes failure. Here, we study the mechanisms responsible for this wave activation via numerical simulations of a granular fault model. We observe the particles belonging to the percolating backbone, which sustains the stress, to perform synchronized oscillations over ellipticlike trajectories in the fault plane. These oscillations occur at the characteristic frequency of acoustic fluidization. As the applied shear stress increases, these oscillations become perpendicular to the fault plane just before the system fails, opposing the confining pressure, consistently with the acoustic fluidization scenario. The same change of orientation can be induced by external perturbations at the acoustic fluidization frequency.

In this paper, subordinated stochastic processes are considered, where the renewal process acting as the operational time. It is assumed that the observation of the process begins at a certain time after the start of the renewal process. A recurrence formula was derived for calculating the multipoint probability density functions of the aged renewal process. Two-point correlation functions for certain subordinated stochastic processes, particularly for the generalized Ornstein-Uhlenbeck process, were calculated. A model of relaxation in a disordered medium with traps and obstacles is proposed.

We consider the steady-state behavior of pairs of active particles having different persistence times and diffusivities. To this purpose we employ the active Ornstein-Uhlenbeck model, where the particles are driven by colored noises with exponential correlation functions whose intensities and correlation times vary from species to species. By extending Fox's theory to many components, we derive by functional calculus an approximate Fokker-Planck equation for the configurational distribution function of the system. After illustrating the predicted distribution in the solvable case of two particles interacting via a harmonic potential, we consider systems of particles repelling through inverse power-law potentials. We compare the analytic predictions to computer simulations for such soft-repulsive interactions in one dimension and show that at linear order in the persistence times the theory is satisfactory. This work provides the toolbox to qualitatively describe many-body phenomena, such as demixing and depletion, by means of effective pair potentials.

Jerky flow in solids results from collective dynamics of dislocations which gives rise to serrated deformation curves and a complex evolution of the strain heterogeneity. A rich example of this phenomenon is the Portevin–Le Chatelier effect in alloys. The corresponding spatiotemporal patterns showed some universal features which provided a basis for a well-known phenomenological classification. Recent studies revealed peculiar features in both the stress serration sequences and the kinematics of deformation bands in Al-based alloys containing fine microstructure elements, such as nanosize precipitates and/or submicron grains. In the present work, jerky flow of an AlMgScZr alloy is studied using statistical analysis of stress serrations and the accompanying acoustic emission. As in the case of coarse-grained binary AlMg alloys, the amplitude distributions of acoustic events obey a power-law scaling which is usually considered as evidence of avalanchelike dynamics. However, the scaling exponents display specific dependences on the strain and strain rate for the investigated materials. The observed effects bear evidence to a competition between the phenomena of synchronization and randomization of dislocation avalanches, which may shed light on the mechanisms leading to a high variety of jerky flow patterns observed in applied alloys.

In this article, multiple eigen-systems including linear growth rates and eigen-functions have been discovered for the Rayleigh-Taylor instability (RTI) by numerically solving the Sturm-Liouville eigen-value problem in the case of two-dimensional plane geometry. The system called the first mode has the maximal linear growth rate and is just extensively studied in literature. Higher modes have smaller eigen-values, but possess multi-peak eigen-functions which bring on multiple pairs of vortices in the vorticity field. A general fitting expression for the first four eigen-modes is presented. Direct numerical simulations show that high modes lead to appearances of multi-layered spike-bubble pairs, and lots of secondary spikes and bubbles are also generated due to the interactions between internal spikes and bubbles. The present work has potential applications in many research and engineering areas, e.g., in reducing the RTI growth during capsule implosions in inertial confinement fusion.

Motivated by the recent experiment by Bordia et al. [Nat. Phys. 13, 460 (2017)], we study the single particle delocalization phenomena of the Aubry-André (AA) model subjected to periodic drives. In two distinct cases we construct an equivalent classical description to illustrate that the drive-induced delocalization phenomena stems from an instability and the onset of chaos in the underlying dynamics. In the first case we analyze the delocalization and the thermalization in a time modulated AA potential with respect to driving frequency and demonstrate that there exists a threshold value of the amplitude of the drive. In the next example, we show that the periodic modulation of the phase of the hopping amplitude induced by a gauge field leads to an unusual effect on delocalization with a nonmonotonic dependence on the driving frequency. Within a window of such a driving frequency a delocalized Floquet band with a mobility edge appears, exhibiting multifractality in the spectrum as well as in the Floquet eigenfunctions. Finally, we explore the effect of interaction and discuss how the results of the present analysis can be tested experimentally.

Optimal sequential inference, or filtering, for the state of a deterministic dynamical system requires simulation of the Frobenius-Perron operator, that can be formulated as the solution of a continuity equation. For low-dimensional, smooth systems, the finite-volume numerical method provides a solution that conserves probability and gives estimates that converge to the optimal continuous-time values, while a Courant-Friedrichs-Lewy-type condition assures that intermediate discretized solutions remain positive density functions. This method is demonstrated in an example of nonlinear filtering for the state of a simple pendulum, with comparison to results using the unscented Kalman filter, and for a case where rank-deficient observations lead to multimodal probability distributions.

We prove a generalized fluctuation-dissipation theorem for a certain class of out-of-time-ordered correlators (OTOCs) with a modified statistical average, which we call bipartite OTOCs, for general quantum systems in thermal equilibrium. The difference between the bipartite and physical OTOCs defined by the usual statistical average is quantified by a measure of quantum fluctuations known as the Wigner-Yanase skew information. Within this difference, the theorem describes a universal relation between chaotic behavior in quantum systems and a nonlinear-response function that involves a time-reversed process. We show that the theorem can be generalized to higher-order $n$-partite OTOCs as well as in the form of generalized covariance.

Scale-invariant power-law distributions for acoustic emission signals are ubiquitous in several plastically deforming materials. However, power-law distributions for acoustic emission energies are reported in distinctly different plastically deforming situations such as hcp and fcc single and polycrystalline samples exhibiting smooth stress-strain curves and in dilute metallic alloys exhibiting discontinuous flow. This is surprising since the underlying dislocation mechanisms in these two types of deformations are very different. So far, there have been no models that predict the power-law statistics for discontinuous flow. Furthermore, the statistics of the acoustic emission signals in jerky flow is even more complex, requiring multifractal measures for a proper characterization. There has been no model that explains the complex statistics either. Here we address the problem of statistical characterization of the acoustic emission signals associated with the three types of the Portevin–Le Chatelier bands. Following our recently proposed general framework for calculating acoustic emission, we set up a wave equation for the elastic degrees of freedom with a plastic strain rate as a source term. The energy dissipated during acoustic emission is represented by the Rayleigh-dissipation function. Using the plastic strain rate obtained from the Ananthakrishna model for the Portevin–Le Chatelier effect, we compute the acoustic emission signals associated with the three Portevin–Le Chatelier bands and the Lüders-like band. The so-calculated acoustic emission signals are used for further statistical characterization. Our results show that the model predicts power-law statistics for all the acoustic emission signals associated with the three types of Portevin–Le Chatelier bands with the exponent values increasing with increasing strain rate. The calculated multifractal spectra corresponding to the acoustic emission signals associated with the three band types have a maximum spread for the type C bands and decreasing with types B and A. We further show that the acoustic emission signals associated with Lüders-like band also exhibit a power-law distribution and multifractality.

Amplitude-modulated bursting (AMB), characterized by oscillations appearing in the envelope of the active phase of bursting, is a novel class of bursting rhythms reported recently. The present paper aims to report a simple and effective method, i.e., the multiple-frequency slow parametric modulation (MFSPM) method, for obtaining such a bursting pattern. We show that the MFSPM can be well controlled so that it may exhibit multiple continuous ups and downs in the active area. Then, the amplitude of the traced active state alternates between increases and decreases accordingly, which leads to oscillations in the envelope of the active phase, and AMB is thus created. Based on this, the route to AMB by the MFSPM is presented. The validity of the approach is demonstrated by several examples. The proposed approach does not depend on specific systems or bifurcations and thus is a general method.

Contact processes form a large and highly interesting class of dynamic processes on networks, including epidemic and information-spreading networks. While devising stochastic models of such processes is relatively easy, analyzing them is very challenging from a computational point of view, particularly for large networks appearing in real applications. One strategy to reduce the complexity of their analysis is to rely on approximations, often in terms of a set of differential equations capturing the evolution of a random node, distinguishing nodes with different topological contexts (i.e., different degrees of different neighborhoods), such as degree-based mean-field (DBMF), approximate-master-equation (AME), or pair-approximation (PA) approaches. The number of differential equations so obtained is typically proportional to the maximum degree $k_{\max }$ of the network, which is much smaller than the size of the master equation of the underlying stochastic model, yet numerically solving these equations can still be problematic for large $k_{\max }$. In this paper, we consider AME and PA, extended to cope with multiple local states, and we provide an aggregation procedure that clusters together nodes having similar degrees, treating those in the same cluster as indistinguishable, thus reducing the number of equations while preserving an accurate description of global observables of interest. We also provide an automatic way to build such equations and to identify a small number of degree clusters that give accurate results. The method is tested on several case studies, where it shows a high level of compression and a reduction of computational time of several orders of magnitude for large networks, with minimal loss in accuracy.

We study the relaxation dynamics of fully clustered networks (maximal number of triangles) to an unclustered state under two different edge dynamics—the double-edge swap, corresponding to degree-preserving randomization of the configuration model, and single edge replacement, corresponding to full randomization of the Erdős-Rényi random graph. We derive expressions for the time evolution of the degree distribution, edge multiplicity distribution and clustering coefficient. We show that under both dynamics networks undergo a continuous phase transition in which a giant connected component is formed. We calculate the position of the phase transition analytically using the Erdős-Rényi phenomenology.