Quantum Mpemba Effect in a Quantum Dot with Reservoirs.

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 quantum 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.

An exact expression of three-body system for the complex shear modulus of frictional granular materials.

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 monodispersed 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 the results by introducing a single fitting parameter.

Softening and Residual Loss Modulus of Jammed Grains under Oscillatory Shear in an Absorbing State.

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 nonzero 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.

Correction: Simulation of dense non-Brownian suspensions with the lattice Boltzmann method: shear jammed and fragile states.

Correction for 'Simulation of dense non-Brownian suspensions with the lattice Boltzmann method: shear jammed and fragile states' by Pradipto et al., Soft Matter, 2020, 16, 945-959, DOI: 10.1039/C9SM00850K.

Shear modulus and reversible particle trajectories of frictional granular materials under oscillatory shear.

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 small values of the areas.

Nonadiabatic Control of Geometric Pumping.

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.

Simulation of dense non-Brownian suspensions with the lattice Boltzmann method: shear jammed and fragile states.

Dense non-Brownian suspensions, including both hydrodynamic interactions and 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 after the reduction of the strain amplitude from the initial oscillations. Here, we find the existence of a shear-jammed state under this protocol in which the storage modulus G' is only finite for high initial strain amplitude γI0. We also find the existence of a fragile state in which both fluid-like and solid-like responses can be detected for an identical area fraction and an initial strain amplitude γI0 depending on the initial phase Θ (or the asymmetricity of the applied strain) of the oscillatory shear. We also observe a 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 a 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.

Size distribution of particles in Saturn's rings from aggregation and fragmentation.

Saturn's rings consist of a huge number of water ice particles, with a tiny addition of rocky material. They form a flat disk, as the result of an interplay of angular momentum conservation and the steady loss of energy in dissipative interparticle collisions. For particles in the size range from a few centimeters to a few meters, a power-law distribution of radii, ~r(-q) with q ≈ 3, has been inferred; for larger sizes, the distribution has a steep cutoff. It has been suggested that this size distribution may arise from a balance between aggregation and fragmentation of ring particles, yet neither the power-law dependence nor the upper size cutoff have been established on theoretical grounds. Here we propose a model for the particle size distribution that quantitatively explains the observations. In accordance with data, our model predicts the exponent q to be constrained to the interval 2.75 ≤ q ≤ 3.5. Also an exponential cutoff for larger particle sizes establishes naturally with the cutoff radius being set by the relative frequency of aggregating and disruptive collisions. This cutoff is much smaller than the typical scale of microstructures seen in Saturn's rings.

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach.

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.

Minimal model of stochastic athermal systems: origin of non-Gaussian noise.

For a wide class of stochastic athermal systems, we derive Langevin-like equations driven by non-Gaussian noise, starting from master equations and developing a new asymptotic expansion. We found an explicit condition whereby the non-Gaussian properties of the athermal noise become dominant for tracer particles associated with both thermal and athermal environments. Furthermore, we derive an inverse formula to infer microscopic properties of the athermal bath from the statistics of the tracer particle. We apply our formulation to a granular motor under viscous friction and analytically obtain the angular velocity distribution function. Our theory demonstrates that the non-Gaussian Langevin equation is the minimal model of athermal systems.

Hydrodynamic instabilities in shear flows of dry cohesive granular particles.

We extend the dynamic van der Waals model introduced by A. Onuki [Phys. Rev. Lett., 2005, 94, 054501] to the description of cohesive granular flows under a plane shear to study their hydrodynamic instabilities. By numerically solving the dynamic van der Waals model, we observed various heterogeneous structures of density fields in steady states, where the viscous heating is balanced with the energy dissipation caused by inelastic collisions. Based on the linear stability analysis, we found 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, ϕ, satisfy ∂p/∂ϕ < 0.

Geometric pumping induced by shear flow in dilute liquid crystalline polymer solutions.

We investigate nonlinear rheology of dilute liquid crystalline polymer solutions under time dependent two-directional shear flow. We analyze the Smoluchowski equation, which describes the dynamics of the orientation of a liquid crystalline polymer, by employing technique of the full counting statistics. In the adiabatic limit, we derive the expression for time integrated currents generated by a Berry-like curvature. Using this expression, it is shown that the expectation values of the time-integrated angular velocity of a liquid crystalline polymer and the time-integrated stress tensor are generally not zero even if the time average of the shear rate is zero. The validity of the theoretical calculations is confirmed by direct numerical simulations of the Smoluchowski equation. Nonadiabatic effects are also investigated by means of simulations and it is found that the time-integrated stress tensor depends on the speed of the modulation of the shear rate if we adopt the isotropic distribution as an initial state.

Effective viscosity and elasticity in dense suspensions under impact: Toward a modeling of walking on suspensions.

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 a 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.

Counterflow-induced clustering: Exact results.

We analyze the cluster formation in a nonergodic stochastic system as a result of counterflow, 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 nonconserved 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 nonconserved species, whereas the free-flowing phase is identified with nonmonotonic density and nonmonotonic 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 nonergodicity 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.

Theory of rigidity and numerical analysis of density of states of two-dimensional amorphous solids with dispersed frictional grains in the linear response regime.

Using the Jacobian matrix, we obtain a 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. 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}.

Eigenvalue analysis of stress-strain curve of two-dimensional amorphous solids of dispersed frictional grains with finite shear strain.

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.

Stochastic energetics for non-Gaussian processes.

By introducing a new stochastic integral, we investigate the energetics of classical stochastic systems driven by non-Gaussian white noises. In particular, we introduce a decomposition of the total energy difference into the work and the heat for each trajectory, and derive a formula to calculate the heat from experimental data on the dynamics. We apply our formulation and results to a Langevin system driven by a Poisson noise.

Mpemba effect in inertial suspensions.

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 quasiequilibrium 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 quasiequilibrium 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 quasiequilibrium 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.

Negative normal restitution coefficient found in simulation of nanocluster collisions.

The oblique impacts of nanoclusters are studied theoretically and by means of molecular dynamics. In simulations we explore two models--Lennard-Jones clusters and particles with covalently bonded atoms. In contrast with the case of macroscopic bodies, the standard definition of the normal restitution coefficient yields for this coefficient negative values for oblique collisions of nanoclusters. We explain this effect and propose a proper definition of the restitution coefficient which is always positive. We develop a theory of an oblique impact based on a continuum model of particles. A surprisingly good agreement between the macroscopic theory and simulations leads to the conclusion that macroscopic concepts of elasticity, bulk viscosity, and surface tension remain valid for nanoparticles of a few hundred atoms.