Source author record

Harukuni Ikeda

Harukuni Ikeda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Active Spherical Model

The spherical model is a popular solvable model and has been applied to describe several critical phenomena such as the ferromagnetic transition, Bose-Einstein condensation, spin-glass transition, glass transition, jamming transition, and so on. Motivated by recent developments of active matter, here we consider the spherical model driven by the Ornstein-Uhlenbeck type self-propulsion force with persistent time $τ_p$. We show that the model exhibits the Ising universality for finite $τ_p$. On the contrary, the model exhibits the random field Ising universality in the limit $τ_p\to\infty$.

preprint2022arXiv

Universality of classical and quantum SAT-UNSAT transitions of convex continuous satisfaction problems

Here we investigate the single-layer linearized perceptron near the SAT-UNSAT transition point as a prototypical model of the convex continuous satisfaction problems. The simplicity of the model allows us to take into account the effects of the quantum fluctuation, which have not been fully investigated before. We found that the classical and quantum models have different critical exponents and thus have different universality classes. We also briefly discuss the effects of the random field.

preprint2022arXiv

Vibrational density of states of jammed packing: mean-field theory

Several mean-field theories predict that Hessian matrices of amorphous solids can be written by using the random matrix in the limit of the large spatial dimensions $d\to\infty$. Motivated by these results, we here propose a way to map a Hessian of the amorphous solid to a random matrix. This is possible by determining the coefficients of a random matrix so that the trace of the random matrix coincides with the Hessian of the original system. We compare our result with that of previous numerical simulations of harmonic spheres in several spatial dimensions $d=3$, $5$, and $9$. For small pressure $p\ll 1$ (near jamming), we find a good agreement even in $d=3$, and obtain better agreements in larger $d$, suggesting that the approximation indeed becomes exact in the limit of large spatial dimensions.

preprint2021arXiv

Multiple glass transitions and higher order replica symmetry breaking of binary mixtures

We extend the replica liquid theory in order to describe the multiple glass transitions of binary mixtures with large size disparities, by taking into account the two-step replica symmetry breaking (2RSB). We determine the glass phase diagram of the mixture of large and small particles in the large-dimension limit where the mean-field theory becomes exact. When the size ratio of particles is beyond a critical value, the theory predicts three distinct glass phases; (i) the 1RSB double glass where both components vitrify simultaneously, (ii) the 1RSB single glass where only large particles are frozen while small particles remain mobile, and (iii) a new glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.

preprint2020arXiv

Force balance controls the relaxation time of the gradient descent algorithm in the satisfiable phase

We numerically study the relaxation dynamics of the single layer perceptron with the spherical constraint. This is the simplest model of neural networks and serves a prototypical mean-field model of both convex and non-convex optimization problems. The relaxation time of the gradient descent algorithm rapidly increases near the SAT-UNSAT transition point. We numerically confirm that the first non-zero eigenvalue of the Hessian controls the relaxation time. This first eigenvalue vanishes much faster upon approaching the SAT-UNSAT transition point than the prediction of Marchenko-Pastur law in random matrix theory derived under the assumption that the set of unsatisfied constraints are uncorrelated. This leads to a non-trivial critical exponent of the relaxation time in the SAT phase. Using a simple scaling analysis, we show that the isolation of this first eigenvalue from the bulk of spectrum is attributed to the force balance at the SAT-UNSAT transition point. Finally, we show that the estimated critical exponent of the relaxation time in the non-convex region agrees very well with that of frictionless spherical particles, which have been studied in the context of the jamming transition of granular materials.

preprint2020arXiv

Jamming and replica symmetry breaking of weakly disordered crystals

We discuss the physics of crystals with small polydispersity near the jamming transition point. For this purpose, we introduce an effective single-particle model taking into account the nearest neighbor structure of crystals. The model can be solved analytically by using the replica method in the limit of large dimensions. In the absence of polydispersity, the replica symmetric solution is stable until the jamming transition point, which leads to the standard scaling of perfect crystals. On the contrary, for finite polydispersity, the model undergoes the full replica symmetry breaking (RSB) transition before the jamming transition point. In the RSB phase, the model exhibits the same scaling as amorphous solids near the jamming transition point. These results are fully consistent with the recent numerical simulations of crystals with polydispersity. The simplicity of the model also allows us to derive the scaling behavior of the vibrational density of states that can be tested in future experiments and numerical simulations.

preprint2020arXiv

Jamming below upper critical dimension

Extensive numerical simulations in the past decades proved that the critical exponents of the jamming of frictionless spherical particles remain unchanged in two and three dimensions. This implies that the upper critical dimension is $d_u=2$ or lower. In this work, we study the jamming transition below the upper critical dimension. We investigate a quasi-one-dimensional system: disks confined in a narrow channel. We show that the system is isostatic at the jamming transition point as in the case of standard jamming transition of the bulk systems in two and three dimensions. Nevertheless, the scaling of the excess contact number shows the linear scaling. Furthermore, the gap distribution remains finite even at the jamming transition point. These results are qualitatively different from those of the bulk systems in two and three dimensions.

preprint2020arXiv

Jamming with tunable roughness

We introduce a new model to study the effect of surface roughness on the jamming transition. By performing numerical simulations, we show that for a smooth surface, the jamming transition density and the contact number at the transition point both increase upon increasing asphericity, as for ellipsoids and spherocylinders. Conversely, for a rough surface, both quantities decrease, in quantitative agreement with the behavior of frictional particles. Furthermore, in the limit corresponding to the Coulomb friction law, the model satisfies a generalized isostaticity criterion proposed in previous studies. We introduce a counting argument that justifies this criterion and interprets it geometrically. Finally, we propose a simple theory to predict the contact number at finite friction from the knowledge of the force distribution in the infinite friction limit.

preprint2020arXiv

Note: Relaxation time below jamming

Like other critical phenomena, the jamming transition accompanies the divergence of the relaxation time $τ$. A recent numerical study of frictionless spherical particles proves that $τ$ is inversely proportional to the lowest non-zero eigenvalue $λ_1$ of the dynamical matrix. In this note, we derive the scaling of $λ_1$ below the jamming transition point $φ_J$ by solving the linearized dynamical equation. The resultant critical exponent agrees with a previous theoretical result for sheared suspension obtained by applying the virtual work theorem to a simple shear, highlighting the universality of the relaxation dynamics below jamming.

preprint2020arXiv

Testing mean-field theory for jamming of non-spherical particles: Contact number, gap distribution, and vibrational density of states

We perform numerical simulations of the jamming transition of non-spherical particles in two dimensions. In particular, we systematically investigate how the physical quantities at the jamming transition point behave when the shapes of the particle deviate slightly from the perfect disks. For efficient numerical simulation, we first derive an analytical expression of the gap function, using the perturbation theory around the reference disks. Starting from disks, we observe the effects of the deformation of the shapes of particles by the $n$-th order term of the Fourier series $\sin(nθ)$. We show that the several physical quantities, such as the number of contacts, gap distribution, and characteristic frequencies of the vibrational density of states, show the power-law behaviors with respect to the linear deviation from the reference disks. The power-law behaviors do not depend on $n$ and are fully consistent with the mean-field theory of the jamming of non-spherical particles. This result suggests that the mean-field theory holds very generally for nearly spherical particles whose shape can be expressed by the Fourier series.

preprint2019arXiv

Infinitesimal asphericity changes the universality of the jamming transition

The jamming transition of non-spherical particles is fundamentally different from the spherical case. Non-spherical particles are hypostatic at their jamming points, while isostaticity is ensured in the case of the jamming of spherical particles. This structural difference implies that the presence of asphericity affects the critical exponents related to the contact number and the vibrational density of states. Moreover, while the force and gap distributions of isostatic jamming present power-law behaviors, even an infinitesimal asphericity is enough to smooth out these singularities. In a recent work [PNAS 115(46), 11736], we have used a combination of marginal stability arguments and the replica method to explain these observations. We argued that systems with internal degrees of freedom, like the rotations in ellipsoids, or the variation of the radii in the case of the \textit{breathing} particles fall in the same universality class. In this paper, we review comprehensively the results about the jamming with internal degrees of freedom in addition to the translational degrees of freedom. We use a variational argument to derive the critical exponents of the contact number, shear modulus, and the characteristic frequencies of the density of states. Moreover, we present additional numerical data supporting the theoretical results, which were not shown in the previous work.

preprint2016arXiv

The decoupling of the glass transitions in the two-component $p$-spin spherical model

Binary mixtures of large and small particles with disparate size ratio exhibit a rich phenomenology at their glass transition points. In order to gain insights on such systems, we introduce and study a two-component version of the $p$-spin spherical spin glass model. We employ the replica method to calculate the free energy and the phase diagram. We show that when the strengths of the interactions of each component are not widely separated, the model has only one glass phase characterized by the conventional one-step replica symmetry breaking. However when the strengths of the interactions are well separated, the model has three glass phases depending on temperature and component ratio. One is the "single" glass phase in which only the spins of one component are frozen while the spins of the other component remain mobile. This phase is characterized by the one-step replica symmetry breaking. The second is the "double" glass phase obtained by cooling further the single glass phase, in which the spins of the remaining mobile component are also frozen. This phase is characterized by the two-step replica symmetry breaking. The third is also the "double" glass phase, which however is formed by the simultaneous freezing of the spins of both components at the same temperatures and is characterized by the one-step replica symmetry breaking. We discuss the implications of these results for the glass transitions of binary mixtures.

preprint2015arXiv

Fredrickson-Andersen model on Bethe lattice with random pinning

We study the effects of random pinning on the Fredrickson-Andersen model on the Bethe lattice. We find that the nonergodic transition temperature rises as the fraction of the pinned spins increases and the transition line terminates at a critical point. The freezing behavior of the spins is analogous to that of a randomly pinned p-spin mean-field spin glass model which has been recently reported. The diverging behavior of correlation lengths in the vicinity of the terminal critical point is found to be identical to the prediction of the inhomogeneous mode-coupling theory at the A3 singularity point for the glass transition.

preprint2015arXiv

One-dimensional Kac model of dense amorphous hard spheres

We introduce a new model of hard spheres under confinement for the study of the glass and jamming transitions. The model is an one-dimensional chain of the $d$-dimensional boxes each of which contains the same number of hard spheres, and the particles in the boxes of the ends of the chain are quenched at their equilibrium positions. We focus on the infinite dimensional limit ($d \to \infty$) of the model and analytically compute the glass transition densities using the replica liquid theory. From the chain length dependence of the transition densities, we extract the characteristic length scales at the glass transition. The divergence of the lengths are characterized by the two exponents, $-1/4$ for the dynamical transition and $-1$ for the ideal glass transition, which are consistent with those of the $p$-spin mean-field spin glass model. We also show that the model is useful for the study of the growing length scale at the jamming transition.