Search

In a sheared steady state, glasses reach a nonequilibrium criticality called yielding criticality. We report that the qualitative nature of this nonequilibrium critical phenomenon depends on the details of the system and that responses and fluctuations are governed by different critical correlation lengths in specific situations. This scale separation of critical lengths arises when the screening of elastic propagation of mechanical signals is not negligible. We also discuss the determinant of the impact of screening effects from the viewpoint of the microscopic dissipation mechanism.

The Herschel-Bulkley law, a universal constitutive relation, has been empirically known to be applicable to a vast range of soft materials, including sheared glasses. Although the Herschel-Bulkley law has attracted public attention, its structural origin has remained an open question. In this Letter, by means of atomistic simulation of binary Lennard-Jones glasses, we report that the instantaneous normal modes with negative eigenvalues, or so-called imaginary modes, serve as the structural signatures for the Herschel-Bulkley rheology in sheared glasses.

Plastic events in sheared glasses are considered an example of so-called avalanches, whose sizes obey a power-law probability distribution with the avalanche critical exponent τ. Although the so-called mean-field depinning (MFD) theory predicts a universal value of this exponent, τ_{MFD}=1.5, such a simplification is now known to connote qualitative disagreement with realistic systems. Numerically and experimentally, different values of τ have been reported depending on the literature. Moreover, in the elastic regime, it has been noted that the critical exponent can be different from that in the steady state, and even criticality itself is a matter of debate. Because these confusingly varying results have been reported under different setups, our knowledge of avalanche criticality in sheared glasses is greatly limited. To gain a unified understanding, in this work, we conduct a comprehensive numerical investigation of avalanches in Lennard-Jones glasses under athermal quasistatic shear. In particular, by excluding the ambiguity and arbitrariness that has crept into the conventional measurement schemes, we achieve high-precision measurement and demonstrate that the exponent τ in the steady state follows the prediction of MFD theory, τ_{MFD}=1.5. Our results also suggest that there are two qualitatively different avalanche events. This binariness leads to the nonuniversal behavior of the avalanche size distribution and is likely to be the cause of the varying values of τ reported thus far. To investigate the dependence of criticality and universality on applied shear, we further study the statistics of avalanches in the elastic regime and the ensemble of the first avalanche event in different samples, which provide information about the unperturbed system. We show that while the unperturbed system is indeed off-critical, criticality gradually develops as shear is applied. The degree of criticality is encoded in the fractal dimension of the avalanches, which starts from zero in the off-critical unperturbed state and saturates in the steady state. Moreover, the critical exponent τ is consistent with the prediction of the MFD τ_{MFD} universally, regardless of the amount of applied shear, once the system becomes critical.

Turbulence is ubiquitous in nonequilibrium systems, and it has been noted that even dense granular flows exhibit characteristics that are typical of turbulent flow, such as the power-law energy spectrum. However, studies on the turbulentlike behavior of granular flows are limited to two-dimensional (2D) flow. We demonstrate that the statistics in three-dimensional (3D) flow are qualitatively different from those in 2D flow. We also elucidate that avalanche dynamics can explain this dimensionality dependence. Moreover, we define clusters of collectively moving particles that are equivalent to vortex filaments. The clusters unveil complicated structures in 3D flows that are absent in 2D flows.

The collective dynamics of externally driven N_{p}-colloidal systems (1≤N_{p}≤4) in a confined viscous fluid have been investigated using three-dimensional direct numerical simulations with fully resolved hydrodynamics. The dynamical modes of collective particle motion are studied by changing the particle Reynolds number as determined by the strength of the external driving force and the confining wall distance. For a system with N_{p}=3, we found that at a critical Reynolds number a dynamical mode transition occurs from the doublet-singlet mode to the triplet mode, which has not been reported experimentally. The dynamical mode transition was analyzed in detail from the following two viewpoints: (1) spectrum analysis of the time evolution of a tagged particle velocity and (2) the relative acceleration of the doublet cluster with respect to the singlet particle. For a system with N_{p}=4, we found similar dynamical mode transitions from the doublet-singlet-singlet mode to the triplet-singlet mode and further to the quartet mode.

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the collective motion of model swimmers in bulk and confinement. Calculating the dynamic structure factor, we clarified that the swarming in bulk systems can be understood as a pseudoacoustic mode. Under confinement between flat parallel walls, this pseudoacoustic mode leads to a traveling wavelike motion. This swarming behavior is due purely to the hydrodynamic interactions between the swimmers and depends strongly on the type and strength of swimming (i.e., pusher or puller).

The dependence of the ion-temperature-gradient scale length on the hydrogen isotope mass was examined in conventional H-mode plasmas in JT-60U tokamak. While identical profiles for density and temperature were obtained for hydrogen and deuterium plasmas, the ion conductive heat flux necessary for hydrogen to sustain the same ion temperature profile was two times that required for deuterium, resulting in a clearly higher ion heat diffusivity for hydrogen at the same ion-temperature-gradient scale length. On the other hand, the ion-temperature-gradient scale length for deuterium is less than that for hydrogen at a given ion heat diffusivity.

A complex multistage transition of the edge radial electric field is observed in JT-60U H-mode phase without edge localized mode. An interesting feature is that the poloidal rotation velocity of the carbon impurity ions changes in the later H-phase without a comparable change in the main ion pressure gradient, indicating a change in the parallel momentum (and particle) balance channel.

The role of electron cyclotron resonance heating (ECRH) on the toroidal rotation velocity profile has been investigated in the JT-60U tokamak device by separating the effects of the change in momentum transport, the intrinsic rotation by pressure gradient, and the intrinsic rotation by ECRH. It is found that ECRH increases the toroidal momentum diffusivity and the convection velocity. It is also found that ECRH drives the codirection (co) intrinsic rotation inside the EC deposition radius and the counterdirection (ctr) intrinsic rotation outside the EC deposition radius. This ctr rotation starts from the EC deposition radius and propagates to the edge region.

We have observed a bursting mode in the high-beta plasmas above the ideal beta limit without a conducting wall. The mode frequency is chirping down as the mode amplitude increases, and its initial value is close to the precession frequency of the trapped energetic particle from the perpendicular neutral beams. The mode structure is radially extended with a peak around the q = 2 surface. This mode can finally trigger the resistive wall mode (RWM) despite enough plasma rotation for RWM stabilization. It is concluded that the mode is driven by trapped energetic particles. The mode is attributed to the interaction between the trapped energetic particles and a marginally stable mode in the wall-stabilized high-beta_{N} region.

A spontaneous transition phenomena between two states of a plasma with an internal transport barrier (ITB) is observed in the steady-state phase of the magnetic shear in the negative magnetic shear plasma in the JT-60U tokamak. These two ITB states are characterized by different profiles of the second radial derivative of the ion temperature inside the ITB region (one has a weak concave shape and the other has a strong convex shape) and by different degrees of sharpness of the interfaces between the L mode and the ITB region, which is determined by the turbulence penetration into the ITB region.

The toroidal plasma rotation generated by the external momentum input and by the plasma itself (intrinsic rotation) has been separated through a novel momentum transport analysis in the JT-60U tokamak device. The toroidal rotation, which is not determined by the momentum transport coefficients and the external momentum input, has been observed. It is found that this intrinsic rotation is locally determined by the local pressure gradient and increases with increasing pressure gradient. This trend is almost the same for various plasmas: low and high confinement mode, co and counterrotating plasmas.

Nondimensional parameter dependence of heat transport between edge localized modes (ELMs) is examined for H mode plasmas. The electron heat diffusivity between ELMs is reduced to the level of ion neoclassical transport in the plasma edge region which is affected by ELM burst. At lower edge collisionality, the heat flux assigned to the heat transport between ELMs is reduced and the ELM loss power is enhanced. During the inter-ELM phase, the energy confinement time becomes larger with decreasing the edge collisionality and poloidal Larmor radius.

The migration processes of the 30 degrees partial dislocation in silicon have been investigated using first-principles total-energy calculations. We have presented for the first time a reliable determination of the formation and migration energies of kinks, which are consistent with experiments. The obtained results indicate that the 30 degrees partial dominates the dislocation dynamics in Si and suggest the dislocation glide picture based on the no-collision model. We suggest that the interpretation of experimental studies on the dislocation dynamics should be reexamined in light of these results.