Search

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.

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 kT/kN, 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 kT/kN and becomes indistinguishable from the translational band for large kT/kN.

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.

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.

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.

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.

A periodically modulated N-state model whose dynamics are 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 with numerical results, which are available for a wide range of modulation frequencies.

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.

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.

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.

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.

We extend a recent theory for the rheology of frictionless granular materials [1] 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.

Efficiency at maximum power (MP) output for an engine with a passive piston without mechanical controls between two reservoirs is studied theoretically. We enclose a hard core gas partitioned by a massive piston in a temperature-controlled container and analyze the efficiency at MP under a heating and cooling protocol without controlling the pressure acting on the piston from outside. We find the following three results: (i) The efficiency at MP for a dilute gas is close to the Chambadal–Novikov–Curzon–Ahlborn (CNCA) efficiency if we can ignore the sidewall friction and the loss of energy between a gas particle and the piston, while (ii) the efficiency for a moderately dense gas becomes smaller than the CNCA efficiency even when the temperature difference of the reservoirs is small. (iii) Introducing the Onsager matrix for an engine with a passive piston, we verify that the tight coupling condition for the matrix of the dilute gas is satisfied, while that of the moderately dense gas is not satisfied because of the inevitable heat leak. We confirm the validity of these results using the molecular dynamics simulation and introducing an effective mean-field-like model which we call the stochastic mean field model.

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.

We propose a non-singular representation for a non-Hermitian operator even if the parameter space contains exceptional points (EPs), at which the operator cannot be diagonalized and the usual spectral representation ceases to exist. Our representation has a generalized Jordan block form and is written in terms of extended pseudo-eigenstates. Our method is free from divergence in the spectral representation at EPs, at which multiple eigenvalues and eigenvectors coalesce and the eigenvectors cannot be normalized. Our representation improves the accuracy of numerical calculations of physical quantities near EPs. We also find that our method is applicable to various problems related to EPs in the parameter space of non-Hermitian operators. We demonstrate the usefulness of our representation by investigating Boltzmann's collision operator in a 1D quantum Lorentz gas in the weak-coupling approximation.

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.

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.

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.

The motion of an adiabatic piston under dry friction is investigated to clarify the roles of dry friction in nonequilibrium steady states. We clarify that dry friction can reverse the direction of the piston motion and causes a discontinuity or a cusplike singularity for velocity distribution functions of the piston. We also show that the heat fluctuation relation is modified under dry friction.

We derive some nonequilibrium identities such as the integral fluctuation theorem and the Jarzynski equality starting from a nonequilibrium state for dissipative classical systems. Thanks to the existence of the integral fluctuation theorem we can naturally introduce an entropy-like quantity for dissipative classical systems in far from equilibrium states. We also derive the generalized Green-Kubo formula as a nonlinear response theory for a steady dynamics around a nonequilibrium state. We numerically verify the validity of the derived formulas for sheared frictionless granular particles.

Using molecular dynamics simulations, we investigate transport properties of a classical two-dimensional electron system confined in a microchannel with a narrow constriction. As a function of the confinement strength of the constriction, the calculated conductance in the simulations exhibits steplike increases, as reported in a recent experiment [ Rees et al. Phys. Rev. Lett. 106 026803 (2011)]. It is confirmed that the number of the steps corresponds to the number of stream lines of electrons through the constriction. We verify that density fluctuation plays a major role in smoothing the steps in the conductance.

Recently, theorists have developed several theoretical results known as the fluctuation theorem (FT), which has been numerically and experimentally tested. FT estimates the probability to observe the negative entropy production in small systems subject to large thermal fluctuation.

We have derived the effect of steady heat flux on the rate of chemical reaction based on the line-of-centers model using the explicit velocity distribution function of the steady-state Boltzmann equation for hard-sphere molecules to second order. It is found that the second-order velocity distribution function plays an essential role for the calculation of it. We have also compared our result with those from the steady-state Bhatnagar–Gross–Krook (BGK) equation and information theory, and found no qualitative differences among them.