Search

The Mpemba effect is a counterintuitive physical phenomenon where a hot system cools faster than a warm one. In recent years, theoretical analyses of the Mpemba effect have been developed for microscopic systems and experimentally verified. However, the conventional theory relies on a specific choice of distance measure to quantify relaxation speed, leading to several theoretical ambiguities. In this Letter, we derive a rigorous quantification of the Mpemba effect based on thermomajorization theory, referred to as the thermomajorization Mpemba effect. This approach resolves all existing ambiguities and provides a unification of the conventional Mpemba effect across all monotone measures. Furthermore, we demonstrate the universality of the thermomajorization Mpemba effect for Markovian dynamics, rigorously proving that it can occur in any temperature regime with fixed energy levels. Our findings offer a paradigm shift in understanding thermal relaxation processes and provide insights into designing systems that exhibit the Mpemba effect., Comment: 7+6 pages, 3+0 figures

Mpemba effects (MPEs), where a hotter system cools faster than a colder one, present intriguing anomalies in relaxation processes. Despite their universal observation and significant fundamental and practical implications, a comprehensive theoretical understanding based on microscopic properties remains elusive. In this Letter, we introduce two universal frameworks for classical systems to address this gap. Firstly, we reveal that MPEs, traditionally defined by macroscopic temperature comparisons, can be understood through microstate comparisons. This insight offers a straightforward and universal microscopic perspective on MPEs, relevant for experiments and numerical simulations to identify their microscopic origins. Secondly, we establish a "no-Mpemba theorem," a rigorous sufficient condition for the absence of MPEs, thereby identifying specific classes of systems devoid of these effects. Our findings are exemplified using ferromagnetic Ising models and one-dimensional multiparticle systems, demonstrating the practical applicability of our theoretical advancements., Comment: 7 pages, 2 figures (main text) + 10 pages, 4 figures (supplemental material). Minor updates

It has recently been established that open quantum systems may exhibit a strong spectral sensitivity to boundary conditions, known as the non-Hermitian/Liouvillian skin effect (NHSE/LSE), making the topological properties of the system boundary-condition sensitive. In this Letter, we ask the query: Can topological phase transitions of open quantum systems, captured by open boundary conditioned invariants, be observed in the dynamics of a system in a periodic boundary condition, even in the presence of NHSE/LSE? We affirmatively respond to this question, by considering the quench dynamics of entanglement spectrum in a periodic open quantum fermionic system. We demonstrate that the entanglement spectrum exhibits zero-crossings only when this periodic system is quenched from a topologically trivial to non-trivial phase, defined from the spectrum in open boundary conditions, even in systems featuring LSE. Our results reveal that non-Hermitian topological phases leave a distinctive imprint on the unconditional dynamics within a subsystem of fermionic systems., Comment: 5 pages, 3 figures, Supplemental Materials 3 pages, 1 supplemental figure

This study investigates the phenomenon of the early-stage dynamics of impact-induced hardening in dense suspensions, where materials undergo solidification upon impact. While Stokes flow theory traditionally applies to suspensions with micrometer-sized particles due to their low Reynolds numbers, suspensions containing larger particles defy such idealizations. Our work focuses on the early-stage impact-induced hardening of suspensions containing millimeter-sized particles through dynamic impact experiments. We are particularly interested in the maximum drag force $F_\mathrm{max}$ acting on the projectile as a function of the impact speed $u_0$. We successfully conducted experiments using these suspensions and confirmed the relation $F_\mathrm{max}\sim u_0^{3/2}$ for relatively large $u_0$ as observed in the previous studies suspensions of micrometer-sized particles. Our findings reveal that the early-stage behaviors of millimeter-sized particle suspensions align well with predictions from the floating model, typically applicable under Stokes flow conditions. This research sheds light on the complex dynamics of impact-induced hardening in dense suspensions, particularly with larger particles, advancing our understanding beyond conventional micrometer-sized systems., Comment: 9 pages, 8 Figures

We explore the role of exceptional points and complex eigenvalues on the occurrence of the quantum Mpemba effect. To this end, we study a two-level driven dissipative system subjected to an oscillatory electric field and dissipative coupling with the environment. We find that both exceptional points and complex eigenvalues can lead to $multiple$ quantum Mpemba effect. It occurs in an observable when time evolved copies corresponding to two different initial conditions, one initially having higher observable value compared to the other and both relaxing towards the same steady state, intersect each other more than once during their relaxation process. Each of the intersections denotes a quantum Mpemba effect and marks the reversal of identities between the two copies i.e. the copy with higher observable value before the intersection becomes the lower valued copy (and vice versa) after the intersection. Such multiple intersections originate from additional algebraic time dependence at the exceptional points and due to oscillatory relaxation in the case of complex eigenvalues. We provide analytical results for quantum Mpemba effect in the density matrix in presence of coherence. Depending on the control parameters (drive and dissipation), observables such as energy, von Neumann entropy, temperature etc. exhibit either single or multiple quantum Mpemba effect. However, the distance from steady state measured in terms of the Kullback-Leibler divergence shows only single quantum Mpemba effect although the corresponding speed gives rise to either single or multiple quantum Mpemba effect., Comment: 24 pages, 18 figures, published version

We demonstrate the quantum Mpemba effect in a quantum dot coupled to two reservoirs, described by the Anderson model. We show that the system temperatures starting from two different initial values (hot and cold), cross each other at finite time (and thereby reverse their identities i.e. hot becomes cold and vice versa) to generate thermal quantam Mpemba effect. The slowest relaxation mode believed to play the dominating role in Mpemba effect in Markovian systems, does not contribute to such anomalous relaxation in the present model. In this connection, our analytical result provides necessary condition for producing quantum Mpemba effect in the density matrix elements of the quantum dot, as a combined effect of the remaining relaxation modes., Comment: 14 pages, 9 figures

The stress-strain curve of two-dimensional frictional dispersed grains interacting with a harmonic potential without considering the dynamical slip under a finite strain is determined by using eigenvalue analysis of the Hessian matrix. After the configuration of grains is obtained, the stress-strain curve based on the eigenvalue analysis is in almost perfect agreement with that obtained by the simulation, even if there are plastic deformations caused by stress avalanches. Unlike the naive expectation, the eigenvalues in our model do not indicate any precursors to the stress-drop events., Comment: 22 pages, 16 figures. arXiv admin note: text overlap with arXiv:2207.06632

We propose a simple model comprising three particles to study the nonlinear mechanical response of jammed frictional granular materials under oscillatory shear. Owing to the introduction of the simple model, we obtain an exact analytical expression of the complex shear modulus for a system including many mono-dispersed disks, which satisfies a scaling law in the vicinity of the jamming point. These expressions perfectly reproduce the shear modulus of the many-body system with low strain amplitudes and friction coefficients. Even for disordered many-body systems, the model reproduces results by introducing a single fitting parameter.

We analyze the cluster formation in a non-ergodic stochastic system as a result of counter-flow, with the aid of an exactly solvable model. To illustrate the clustering, a two species asymmetric simple exclusion process with impurities on a periodic lattice is considered, where the impurity can activate flips between the two non-conserved species. Exact analytical results, supported by Monte Carlo simulations, show two distinct phases, free flowing phase and clustering phase. The clustering phase is characterized by constant density and vanishing current of the non-conserved species, whereas the free flowing phase is identified with non-monotonic density and non-monotonic finite current of the same. The $n$-point spatial correlation between $n$ consecutive vacancies grows with increasing $n$ in the clustering phase, indicating the formation of two macroscopic clusters in this phase, one of the vacancies and the other consisting of all the particles. We define a rearrangement parameter that permutes the ordering of particles in the initial configuration, keeping all the input parameters fixed. This rearrangement parameter reveals the significant effect of non-ergodicity on the onset of clustering. For a special choice of the microscopic dynamics, we connect the present model to a system of run and tumble particles used to model active matter, where the two species having opposite net bias manifest the two possible run directions of the run and tumble particles, and the impurities act as tumbling reagents that enable the tumbling process., Comment: 20 pages, 13 figures, published version

Using the Jacobian matrix, we obtain theoretical expression of rigidity and the density of states of two-dimensional amorphous solids consisting of frictional grains in the linear response to an infinitesimal strain, in which we ignore the dynamical friction caused by the slip processes of contact points. The theoretical rigidity agrees with that obtained by molecular dynamics simulations. We confirm that the rigidity is smoothly connected to the value in the frictionless limit. For the density of states, we find that there are two modes in the density of states for sufficiently small $k_{T}/k_{N}$, which is the ratio of the tangential to normal stiffness. Rotational modes exist at low frequencies or small eigenvalues, whereas translational modes exist at high frequencies or large eigenvalues. The location of the rotational band shifts to the high-frequency region with an increase in $k_{T}/k_{N}$ and becomes indistinguishable from the translational band for large $k_{T}/k_{N}$. The rigidity determined by the translational modes agrees with that obtained by the molecular dynamics simulations, whereas the contribution of the rotational modes is almost zero for small $k_{T}/k_{N}$., Comment: 28 pages, 16 figures

We demonstrate that discontinuous shear thickening (DST) can take place even in a moderately dense inertial suspension of hydrodynamically interacting frictionless soft particles. The results which demonstrate this fact are obtained using the Lubrication-Friction Discrete Element Method. Our simulation indicates that DST can be observed for lower densities if the inertia of suspended particles and their softness are both of a marked nature. We also confirm that the DST behavior is effectively approximated by the kinetic theory under these conditions without consideration of hydrodynamic interactions., Comment: 12 pages, 4 figures

We demonstrate that a discontinuous shear thickening (DST) can take place even in a moderately dense inertial suspension consisting of frictionless soft particles. This DST can be regarded as an ignited-quenched transition in the inertial suspension. An approximate kinetic theory well recovers the results of the Langevin simulation in the wide range of the volume fraction without any fitting parameters., Comment: 27 pages, 18 figures

We theoretically study the entropy production and the work extracted from a system connected to two reservoirs by periodic modulations of the electrochemical potentials of the reservoirs and the parameter of a system Hamiltonian under isothermal conditions. We find that the modulation of the parameters can drive a geometric state, which is away from a nonequilibrium steady state. Using this property, we construct a demon in which the entropy production during the first one-cycle is negative such that we can extract the work if we start from the nonequilibrium steady state without parameter modulations. We use the Anderson model to implement the demon in a realistic situation., Comment: 16 pages, 16 figures, v6: revised version

The elastic response of dense suspensions under an impact is studied using coupled Lattice Boltzmann Method and Discrete Element Method (LBM-DEM) and its reduced model. We succeed to extract the elastic force acting on the impactor in dense suspensions, which can exist even in the absence of percolating clusters of suspended particles. We then propose a reduced model to describe the motion of the impactor and demonstrate its relevancy through the comparison of the solution of the reduced model and that of LBM-DEM. Furthermore, we illustrate that the perturbation analysis of the reduced model captures the short-time behavior of the impactor motion quantitatively. We apply this reduced model to the impact of the foot-spring-body system on a dense suspension, which is the minimal model to realize walking on the suspension. Due to the spring force of the system and the stiffness of the suspension, the foot undergoes multiple bounces. We also study the parameter dependencies of the hopping motion and find that multiple bounces are suppressed as the spring stiffness increases.

We obtain an exact matrix product steady state for a class of multi species asymmetric simple exclusion process with impurities, under periodic boundary condition. Alongside the usual hopping dynamics, an additional flip dynamics is activated only in the presence of impurities. Although the microscopic dynamics renders the system to be non-ergodic, exact analytical results for observables are obtained in steady states for a specific class of initial configurations. Interesting physical features including negative differential mobility and transition of correlations from negative to positive with changing vacancy density, have been observed. We discuss plausible connections of this exactly solvable model with multi lane asymmetric simple exclusion processes as well as enzymatic chemical reactions., Comment: 43 pages, 15 figures, typos corrected, references added, new appendices added

We propose a geometrical engine undergoing an adiabatic (Thouless) pumping process for a small system connected to external isothermal reservoirs with the control of electrochemical potentials of the reservoirs and one parameter in the system Hamiltonian. Thanks to the geometrical nature of this process, the entropy production is characterized by the geometric metric tensor which is connected to the Fisher information and the Hessian of the density matrix in a nonequilibrium steady state. The existence of an inequality between the thermodynamic length and entropy production is established. We also establish that the work done on this system is characterized by a vector potential and is equivalent to the thermodynamic flux. To characterize the engine, we the introduce effective efficiency as the relation between the work and entropy production. Through the theoretical analysis of the quantum master equation for the Anderson model of a quantum dot within the wide-band approximation, we illustrate the explicit values of the work, thermodynamic length, and effective efficiency of the engine as functions of the phase difference of the externally controlled electrochemical potentials., Comment: 27 pages, 11 figures

The rheology of a dilute binary mixture of inertial suspension under simple shear flow is analyzed in the context of the Boltzmann kinetic equation. The effect of the surrounding viscous gas on the solid particles is accounted for by means of a deterministic viscous drag force plus a stochastic Langevin-like term defined in terms of the environmental temperature $T_\text{env}$. Grad's moment method is employed to determine the temperature ratio and the pressure tensor in terms of the coefficients of restitution, concentration, the masses and diameters of the components of the mixture, and the environmental temperature. Analytical results are compared against event-driven Langevin simulations for mixtures of hard spheres with the same mass density $m_1/m_2=(\sigma^{(1)}/\sigma^{(2)})^3$, $m_i$ and $\sigma^{(1)}$ being the mass and diameter, respectively, of the species $i$. It is confirmed that the theoretical predictions agree with simulations of various size ratios $\sigma^{(1)}/\sigma^{(2)}$ and for elastic and inelastic collisions in the wide range of parameters' space. It is remarkable that the temperature ratio $T_1/T_2$ and the viscosity ratio $\eta_1/\eta_2$ ($\eta_i$ being the partial contribution of the species $i$ to the total shear viscosity $\eta=\eta_1+\eta_2$) discontinuously change at a certain shear rate as the size ratio increases; this feature (which is expected to occur in the thermodynamic limit) cannot be completely captured by simulations due to small system size. In addition, a Bhatnagar--Gross--Krook (BGK)-type kinetic model adapted to mixtures of inelastic hard spheres is exactly solved when $T_\text{env}$ is much smaller than the kinetic temperature $T$. A comparison between the velocity distribution functions obtained from Grad's method, BGK model, and simulations is carried out., Comment: 39 pages, 22 figures

We numerically study impact processes on dense suspensions using the lattice Boltzmann method to elucidate the connection between the elastic rebound of an impactor and relations among the impact speed $u_0$, maximum force acting on the impactor $F_{\rm max}$, and elapsed time $t_{\rm max}$ to reach $F_{\rm max}$. We find that $t_{\rm max}$ emerges in the early stage of the impact, while the rebound process takes place in the late stage. We find a crossover of $F_{\rm max}$ from $u_0$ independent regime for low $u_0$ to a power law regime satisfying $F_{\rm max}\propto u_0^\alpha$ with $\alpha\approx 1.5$ for high $u_0$. Similarly, $t_{\rm max}$ satisfies $t_{\rm max}\propto u_0^{\beta}$ with $\beta\approx -0.5$ for high $u_0$. Both power-law relations for $F_{\rm max}$ and $t_{\rm max}$ versus $u_0$ for high $u_0$ are independent of the system size, but the rebound phenomenon strongly depends on the depth of the container for suspensions. Thus, we indicate that the rebound phenomenon is not directly related to the relations among $u_0$, $F_{\rm max}$ and $t_{\rm max}$. We propose a floating + force chain model, where the rebound process is caused by an elastic term that is proportional to the number of the connected force chains from the impactor to the bottom plate. On the other hand, there are no elastic contributions in the relations for $F_{\rm max}$ and $t_{\rm max}$ against $u_0$ because of the absence of percolated force chains in the early stage. This phenomenology predicts $F_{\rm max}\propto u_0^{3/2}$ and $t_{\rm max}\propto u_0^{-1/2}$ for high $u_0$ and also recovers the behavior of the impactor quantitatively even if there is the rebound process.

In this study, we numerically investigated the mechanical responses and trajectories of frictional granular particles under oscillatory shear in the reversible phase where particle trajectories form closed loops below the yielding point. When the friction coefficient is small, the storage modulus exhibits softening, and the loss modulus remains finite in the quasi-static limit. As the friction coefficient increases, the softening and residual loss modulus are suppressed. The storage and loss moduli satisfy scaling laws if they are plotted as functions of the areas of the loop trajectories divided by the strain amplitude and diameter of grains, at least for their small values., Comment: 7 pages, 11 figures

From a theoretical study of the mechanical response of jammed materials comprising frictionless and overdamped particles under oscillatory shear, we find that the material becomes soft, and the loss modulus remains finite even in an absorbing state where any irreversible plastic deformation does not exist. The trajectories of the particles in this region exhibit hysteresis loops. We succeed in clarifying the origin of the softening of the material and the residual loss modulus with the aid of Fourier analysis. We also clarify the roles of the yielding point in the softening to distinguish the plastic deformation from reversible deformation in the absorbing state.

We perform numerical simulations of a two-dimensional frictional granular system under oscillatory shear confined by constant pressure. We found that the system undergoes dilatancy as the strain increases. We confirmed that compaction also takes place at an intermediate strain amplitude for a small mutual friction coefficient between particles. We also found that compaction depends on the confinement pressure while dilatancy little depends on the pressure., Comment: 4 pages, 4 figures

The Mpemba effect (a counterintuitive thermal relaxation process where an initially hotter system may cool down to the steady state sooner than an initially colder system) is studied in terms of a model of inertial suspensions under shear. The relaxation to a common steady state of a suspension initially prepared in a quasi-equilibrium state is compared with that of a suspension initially prepared in a nonequilibrium sheared state. Two classes of Mpemba effect are identified, the normal and the anomalous one. The former is generic, in the sense that the kinetic temperature starting from a cold nonequilibrium sheared state is overtaken by the one starting from a hot quasi-equilibrium state, due to the absence of initial viscous heating in the latter, resulting in a faster initial cooling. The anomalous Mpemba effect is opposite to the normal one since, despite the initial slower cooling of the nonequilibrium sheared state, it can eventually overtake an initially colder quasi-equilibrium state. The theoretical results based on kinetic theory agree with those obtained from event-driven simulations for inelastic hard spheres. It is also confirmed the existence of the inverse Mpemba effect, which is a peculiar heating process, in these suspensions. More particularly, we find the existence of a mixed process in which both heating and cooling can be observed during relaxation., Comment: 28 pages, 20 figures

The Enskog kinetic theory for moderately dense inertial suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the background fluid on suspended particles is modeled via a viscous drag force plus a Langevin-like term defined in terms of the background temperature. In a previous paper [Hayakawa et al., Phys. Rev. E 96, 042903 (2017)], Grad's moment method with the aid of a linear shear-rate expansion was employed to obtain a theory which gave good agreement with the results of event-driven Langevin simulations of hard spheres for low densities and/or small shear rates. Nevertheless, the previous approach had a limitation of applicability to the high shear-rate and high density regime. Thus, in the present paper, we extend the previous work and develop Grad's theory including higher order terms in the shear rate. This improves significantly the theoretical predictions, a quantitative agreement between theory and simulation being found in the high-density region (volume fractions smaller than or equal to $0.4$)., Comment: 27 pages, 9 figures. arXiv admin note: text overlap with arXiv:1707.09694

We numerically study the impact-induced hardening in dense suspensions. We employ the lattice Boltzmann method and perform simulations of dense suspensions under impacts, which incorporate the contact between suspended particles with the free surface of the suspension. Our simulation for a free-falling impactor on a dense suspension reproduces experimental results, where rebound takes place for frictional particles at high-speed impact and high volume fraction shortly after the impact before subsequently sinking. We found that the shear stress of the suspension is not affected by the impact, which clearly distinguishes the impact-induced hardening from the discontinuous shear thickening. Instead, we found the existence of a localized region with distinctively high value of normal stress corresponding to the dynamically jammed region. Our simulation indicates that the frictional interaction between suspended particles is important for the impact-induced hardening to maintain the dynamically jammed region. Furthermore, persistent homology analysis successfully elucidates the topological structure of force chains.

A periodically modulated N-state model whose dynamics is governed by a time-convoluted generalized master equation is theoretically analyzed. It is shown that this non-Markovian master equation can be converted to a Markovian master equation having a larger transition matrix, which affords easier analysis. The behavior of this model is investigated by focusing on the cycle-averaged pumping current. In the adiabatic limit, the geometrical current is calculated analytically, and compared to numerical results which are available for a wide range of modulation frequencies., Comment: 12 pages, 14 figures, revised based on referee comments

We derive the fluctuation theorem for a stochastic and periodically driven system coupled to two reservoirs with the aid of a master equation. We write down the cumulant generating functions for both the current and entropy production in closed compact forms so as to treat the adiabatic and nonadiabatic contributions systematically. We derive the fluctuation theorem by taking into account the time reversal symmetry and the property that the instantaneous currents flowing into the left and the right reservoir are not equal. It is found that the fluctuation-dissipation relation derived from the fluctuation theorem involves an expansion with respect to the time derivative of the affinity., Comment: 22 pages, 7 figures, revised for publication

We investigate a heat engine under an adiabatic (Thouless) pumping process. In this process, the extracted work and lower bound on dissipated availability are characterized by a vector potential and a Riemannian metric tensor, respectively. We derive a trade-off relation between the power and effective efficiency. We also explicitly calculate the trade-off relation as well as the power and effective efficiency for a spin-boson model coupled to two reservoirs., Comment: 9 pages, 7 figures

We study nonadiabatic effects of geometric pumping. With arbitrary choices of periodic control parameters, we go beyond the adiabatic approximation to obtain the exact pumping current. We find that a geometrical interpretation for the nontrivial part of the current is possible even in the nonadiabatic regime. The exact result allows us to find a smooth connection between the adiabatic Berry phase theory at low frequencies and the Floquet theory at high frequencies. We also study how to control the geometric current. Using the method of shortcuts to adiabaticity with the aid of an assisting field, we illustrate that it enhances the current., Comment: 6+7 pages, revised for publication

We derive an extended fluctuation relation for an open system coupled with two reservoirs under adiabatic one-cycle modulation. We confirm that the geometric phase caused by the Berry-Sintisyn-Nemenman curvature in the parameter space generates non-Gaussian fluctuations. This non-Gaussianity is enhanced for the instantaneous fluctuation relation when the bias between the two reservoirs disappears.

Particle flows injected as beams and scattered by an intruder are numerically studied. We find a crossover of the drag force from Epstein's law to Newton's law, depending on the ratio of the speed to the thermal speed. These laws can be reproduced by a simple analysis of a collision model between the intruder and particle flows. The crossover from Epstein's law to Stokes' law is also found for the low-speed regime as the time evolution of the drag force caused by beam particles. We also show the existence of turbulent-like behavior of the particle flows behind the intruder with the aid of the second invariant of the velocity gradient tensor and the relative mean square displacement for the high-speed regime and a large intruder., Comment: 10 pages, 13 figures

Dense non-Brownian suspensions including both the hydrodynamic interactions and the frictional contacts between particles are numerically studied under simple and oscillatory shears in terms of the lattice Boltzmann method. We successfully reproduce the discontinuous shear thickening (DST) under a simple shear for bulk three-dimensional systems. For our simulation of an oscillatory shear in a quasi-two-dimensional system, we measure the mechanical response when we reduce the strain amplitude after the initial oscillations with a larger strain amplitude. Here, we find the existence of the shear-jammed state under this protocol in which the storage modulus $G^{\prime}$ is only finite for high initial strain amplitude $\gamma_0^{I}$. We also find the existence of the fragile state in which both fluid-like and solid-like responses can be detected for an identical area fraction and an initial strain amplitude $\gamma_0^{I}$ depending on the initial phase $\Theta$ (or the asymmetricity of the applied strain) of the oscillatory shear. We also observe the DST-like behavior under the oscillatory shear in the fragile state. Moreover, we find that the stress anisotropy becomes large in the fragile state. Finally, we confirm that the stress formula based on the angular distribution of the contact force recovers the contact contributions to the stress tensors for both simple and oscillatory shears with large strains.

Herein, we numerically study the rheology of a two-dimensional frictional granular system confined by constant pressure under oscillatory shear. Several scaling laws for the storage and loss moduli against the scaled strain amplitude have been found. The scaling laws in plastic regime for large strain amplitude can be understood by the angular distributions of the contact force. The scaling exponents are estimated by considering the physical mechanism., Comment: 13 pages, 18 figures

The drag acting on an intruder in a three-dimensional frictionless dry granular environment is numerically studied. It is found the followings: (i) There is no yield force for the motion of the intruder without the gravity. (ii) The drag is proportional to the cross section of the moving intruder. (iii) If the intruder is larger than surrounding grains, the drag is proportional to the moving speed $V$ of the intruder for dense systems, but it exhibits a crossover from quadratic to linear dependences of the moving speed when the volume fraction of the surrounding grains is much lower than the jamming point. (iv) There is a plateau regime where the drag is almost independent of $V$ if the size of the intruder is identical to those of the environmental grains and the volume fraction is near the jamming point., Comment: 6 pages, 7 figures, Granular Matter in press

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus $G'$ and the loss modulus $G''$ in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying $G'>0$) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed $G'$ after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus $G''$ exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

Modeling the dynamics of an individual active particle invariably involves an isotropic noisy self-propulsion component, in the form of run-and-tumble motion or variations around it. This nonequilibrium source of noise is neither white---there is persistence---nor Gaussian. While emerging collective behavior in active matter has hitherto been attributed to the persistent ingredient, we focus on the non-Gaussian ingredient of self-propulsion. We show that by itself, that is without invoking any memory effect, it is able to generate particle accumulation close to boundaries and effective attraction between otherwise repulsive particles, a mechanism which generically leads to motility-induced phase separation in active matter., Comment: 14 pages, 4 figures

Making use of the system of pulling a spherical intruder in static three-dimensional granular media, we numerically study the scaling law for the drag force $F_{\rm drag}$ acting on the moving intruder under the influence of the gravity. Suppose if the intruder of diameter $D$ immersed in a granular medium consisting of grains of average diameter $d$ is located at a depth $h$ and moves with a speed $V$, we find that $F_{\rm drag}$ can be scaled as $(D+d)^{\phi_\mu} h^{\alpha_\mu}$ with two exponents $\phi_\mu$ and $\alpha_\mu$, which depend on the friction coefficient $\mu$ and satisfy an approximate sum rule $\phi_\mu+\alpha_\mu\approx 3$. This scaling law is valid for the arbitrary Froude number (defined by $\mathrm{Fr}={2 V}\sqrt{{2D}/{g}}\big/(D+d)$), if $h$ is sufficiently deep. We also identify the existence of three regimes (quasistatic, linear, and quadratic) at least for frictional grains in the velocity dependence of drag force. The crossovers take place at $\mathrm{Fr}\approx 1$ between the quasistatic to the linear regimes and at $\mathrm{Fr}\approx 5$ between the linear to the quadratic regimes. We also observe that Froude numbers at which these crossovers between the regimes happen are independent of the depth $h$ and the diameter of the intruder $D$. We also report the numerical results on the average coordination number of the intruder and average contact angle as functions of intruder velocity., Comment: 15 pages, 20 figures

A systematic microscopic theory for the rheology of dense non-Brownian suspensions characterized by the volume fraction $\varphi$ is developed. The theory successfully derives the critical behavior in the vicinity of the jamming point (volume fraction $\varphi_{J}$), for both the pressure $P$ and the shear stress $\sigma_{xy}$, i.e. $P \sim \sigma_{xy} \sim \dot\gamma \eta_0 \delta\varphi^{-2}$, where $\dot\gamma$ is the shear rate, $\eta_0$ is the shear viscosity of the solvent, and $\delta\varphi = \varphi_J - \varphi > 0$ is the distance from the jamming point. It also successfully describes the behavior of the stress ratio $\mu = \sigma_{xy}/P$ with respect to the viscous number $J=\dot\gamma\eta_{0}/P$., Comment: 49 pages, 16 figures

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicate that while the theoretical predictions for the kinetic temperature agree well with simulations for densities $\varphi \lesssim 0.5$, the agreement for the other rheological quantities (the viscosity, the stress ratio and the normal stress differences) is limited to more moderate densities ($\varphi \lesssim 0.3$) if the inelasticity during collisions between particles is not large. [This paper has been published in Phys. Rev. E {\bf 96}, 42903 (2017) but we have realized that there are some typos and mistakes after its publication. So we add the Erratum which will be published in PRE in the top of this paper.], Comment: 30 pages, 15 figures +Erratum 2 pages, 2 figures

Rheology of a dilute cohesive granular gas is theoretically and numerically studied. The flow curve between the shear viscosity and the shear rate is derived from the inelastic Boltzmann equation for particles having square-well potentials in a simple shear flow. It is found that (i) the stable uniformly sheared state only exists above a critical shear rate and (ii) the viscosity in the uniformly sheared flow is almost identical to that for uniformly sheared flow of hard core granular particles. Below the critical shear rate, clusters grow with time, in which the viscosity can be approximated by that for the hard-core fluids if we replace the diameter of the particle by the mean diameter of clusters., Comment: 14 pages, 16 figures

We derive an extended fluctuation theorem for a geometric pumping in a spin-boson system under a periodic control of environmental temperatures by using a Markovian quantum master equation. We perform the Monte-Carlo simulation and obtain the current distribution, the average current and the fluctuation. Using the extended fluctuation theorem we try to explain the results of our simulation. The fluctuation theorem leads to the fluctuation dissipation relations but the absence of the conventional reciprocal relation., Comment: 15 pages, 8 figures

The shear modulus of jammed frictional granular materials with the harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude $\gamma_0$ discontinuously emerges at the jamming transition point. The storage modulus for small $\gamma_0$ differs from that of frictionless grains even in the zero friction limit, while they are almost identical with each other for sufficiently large $\gamma_0$, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for relatively larger strain. We propose a new scaling law to interpolate between the states of small and large $\gamma_0$., Comment: 10 pages, 13 figures

A kinetic theory for a dilute gas-solid suspension under a simple shear is developed. With the aid of the corresponding Boltzmann equation, it is found that the flow curve (stress-strain rate relation) has a S-shape as a crossover from the Newtonian to the Bagnoldian for a granular suspension or from the Newtonian to a fluid having a viscosity proportional to the square of the shear rate for a suspension consisting of elastic particles. The existence of the S-shape in the flow curve directly leads to a discontinuous shear thickening (DST). This DST corresponds to the discontinuous transition of the kinetic temperature between a quenched state and an ignited state. The results of the event-driven Langevin simulation of hard spheres perfectly agree with the theoretical results without any fitting parameter. The simulation confirms that the DST takes place in the linearly unstable region of the uniformly sheared state., Comment: 21 pages, 6 figures (to be published in Progress of Theoretical and Experimental Physics)

We extend a recent theory for the rheology of frictionless granular materials [K. Suzuki and H. Hayakawa, Phys. Rev. Lett. 2015, 115, 098001] to the case of frictional disks in two dimensions. Employing a frictional contact model for molecular dynamics simulations, we derive difference equations of the shear stress, the granular temperature, and the spin temperature from the generalized Green-Kubo formula, where all the terms are given by microscopic expressions. The numerical solutions of the difference equations not only describe the flow curve, but also reproduce the hysteresis of shear stress, which can be the signature of discontinuous shear thickening of frictional disks., Comment: 4 pages, 1 figure, the conference proceedings for Powders & Grains 2017

A simple kinetic theory to exhibit a discontinuous shear thickening (DST) is proposed. The model includes the collision integral and the friction from environment as well as a thermostat term characterized by $T_{\rm ex}$. The viscosity of this model is proportional to $\dot\gamma^2$ for large shear rate $\dot\gamma$, while it is Newtonian for low $\dot\gamma$. The emergence of the DST is enhanced for lower density and lower nonzero $T_{\rm ex}$., Comment: 4 pages, 2 figures, Powders and Grains 2017 (in press)

We study the non-Gaussian dynamics of a quasiclassical electronic circuit inductively coupled to a mesoscopic conductor. Non-Gaussian noise accompanying the nonequilibrium transport through the conductor significantly modifies the steady-state probability density function (PDF) of the flux in the dissipative circuit. We evaluate the quantum correction of the steady-state PDF by incorporating the quantum fluctuation of the circuit. The inverse formula to infer the statistical properties of the non-Gaussian noise from the steady-state PDF is extended to the classical-quantum crossover regime., Comment: 12 pages, 6 figures

We build up a phenomenological picture in terms of the effective dynamics of a tracer confined in a cage experiencing random hops to capture somec haracteristics of glassy systems. This minimal description exhibits scale invariance properties for the small-displacement distribution that echo experimental observations. We predict the existence of exponential tails as a cross-over between two Gaussian regimes. Moreover, we demonstrate that the onset of glassy behavior is controlled only by two dimensionless numbers: the number of hops occurring during the relaxation of the particle within a local cage, and the ratio of the hopping length to the cage size., Comment: 9 pages, 3 figures

This study numerically and analytically investigates the dynamics of a rotor under viscous or dry friction as a non-equilibrium probe of a granular gas. In order to demonstrate the role of the rotor as a probe for a non-equilibrium bath, the molecular dynamics (MD) simulation of the rotor is performed under viscous or dry friction surrounded by a steady granular gas under gravity. A one- to-one map between the velocity distribution function (VDF) of the granular gas and the angular distribution function for the rotor is theoretically derived. The MD simulation demonstrates that the one-to-one map accurately infers the local VDF of the granular gas from the angular VDF of the rotor, and vice versa., Comment: 14pages, 10 figures

We develop the kinetic theory of dilute cohesive granular gases in which the attractive part is described by a square well potential. We derive the hydrodynamic equations from the kinetic theory with the microscopic expressions for the dissipation rate and the transport coefficients. We check the validity of our theory by performing the direct simulation Monte Carlo., Comment: 22 pages, 11 figures

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter., Comment: 18 pages, 5 figures

We extend the dynamic van der Waals model introduced by A. Onuki [Phys. Rev. Lett. 94, 054501 (2005)] to the description of cohesive granular flows under a plane shear to study their hydrodynamic instabilities. Numerically solving the dynamic van der Waals model, we observe various heterogeneous structures of the density in steady states, where the viscous heating is balanced with the energy dissipation caused by inelastic collisions. Based on the linear stability analysis, we find that the spatial structures are determined by the mean volume fraction, the applied shear rate, and the inelasticity, where the instability is triggered if the system is thermodynamically unstable, i.e. the pressure, $p$, and the volume fraction, $\phi$, satisfy $\partial p/\partial\phi<0$., Comment: 15 pages, 7 figures. Submitted to Soft Matter