BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Koji Hukushima

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.dis-nncond-mat.mtrl-sciInformation Theorymath.ITBiological Physicscond-mat.mes-hallcond-mat.str-elData Structures and AlgorithmsMachine Learningquant-ph

Catalog footprint

What is connected

27works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

27 published item(s)

preprint2022arXiv

Absence of Phase Transition in Random Language Model

The Random Language Model, proposed as a simple model of human languages, is defined by the averaged model of a probabilistic context-free grammar. This grammar expresses the process of sentence generation as a tree graph with nodes having symbols as variables. Previous studies proposed that a phase transition, which can be considered to represent the emergence of order in language, occurs in the random language model. We discuss theoretically that the analysis of the "order parameter" introduced in previous studies can be reduced to solving the maximum eigenvector of the transition probability matrix determined by a grammar. This helps analyze the distribution of a quantity determining the behavior of the "order parameter" and reveals that no phase transition occurs. Our results suggest the need to study a more complex model such as a probabilistic context-sensitive grammar, in order for phase transitions to occur.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Data-driven sensitivity analysis in a total-reflection high-energy positron diffraction (TRHEPD) experiment of the Si$_4$O$_5$N$_3$ / 6H-SiC (0001)-($\sqrt{3} \times \sqrt{3}$) R30$^\circ$

The present article proposes a data analysis method for experimentally-derived measurements, which consists of an auto-optimization procedure and a sensitivity analysis. The method was applied to the results of a total-reflection high-energy positron diffraction (TRHEPD) experiment, a novel technique of determining surface structures or the position of the atoms near the material surface. This method solves numerically the partial differential equation in the fully-dynamical quantum diffraction theory with many trial surface structures. In the sensitivity analysis, we focused on the experimental uncertainties and the variation over individual fitting parameters, which was analyzed by solving the eigenvalue problem of the variance-covariance matrix. A modern massively parallel supercomputer was used to complete the analysis within a moderate computational time. The sensitivity analysis provides a basis for the choice of variables in the data analysis for practical reliability. The effectiveness of the present analysis method was demonstrated in the structure determination of a Si$_4$O$_5$N$_3$ / 6H-SiC(0001)-($\sqrt{3} \times \sqrt{3}$) R30$^\circ$ surface. Furthermore, this analysis method is applicable to many experiments other than TRHEPD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-driven determination of the spin Hamiltonian parameters and their uncertainties: The case of the zigzag-chain compound KCu$_4$P$_3$O$_{12}$

We propose a data-driven technique to estimate the spin Hamiltonian, including uncertainty, from multiple physical quantities. Using our technique, an effective model of KCu$_4$P$_3$O$_{12}$ is determined from the experimentally observed magnetic susceptibility and magnetization curves with various temperatures under high magnetic fields. An effective model, which is the quantum Heisenberg model on a zigzag chain with eight spins having $J_1= -8.54 \pm 0.51 \{\rm meV}$, $J_2 = -2.67 \pm 1.13 \{\rm meV}$, $J_3 = -3.90 \pm 0.15 \{\rm meV}$, and $J_4 = 6.24 \pm 0.95 \{\rm meV}$, describes these measured results well. These uncertainties are successfully determined by the noise estimation. The relations among the estimated magnetic interactions or physical quantities are also discussed. The obtained effective model is useful to predict hard-to-measure properties such as spin gap, spin configuration at the ground state, magnetic specific heat, and magnetic entropy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Higher-order tensor renormalization group approach to lattice glass model

In this study, the higher-order tensor renormalization group (HOTRG) method is applied to a lattice glass model that has local constraints on the occupation number of neighboring particles represented by many-body interactions. This model exhibits first- and second-order transitions depending on a certain model parameter. The results obtained by using the HOTRG method for the model were confirmed to be consistent with those obtained by a Markov-chain Monte Carlo (MCMC) method for systems of relatively small sizes. The transition points are accurately estimated by the HOTRG calculation for the systems of large sizes, which is challenging to perform using the MCMC method. These results demonstrate that the HOTRG method can be an efficient method for studying systems with many-body interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lattice glass model in three spatial dimensions

The understanding of thermodynamic glass transition has been hindered by the lack of proper models beyond mean-field theories. Here, we propose a three-dimensional lattice glass model on a simple cubic lattice that exhibits the typical dynamics observed in fragile supercooled liquids such as two-step relaxation, super-Arrhenius growth in the relaxation time, and dynamical heterogeneity. Using advanced Monte Carlo methods, we compute the thermodynamic properties deep inside the glassy temperature regime, well below the onset temperature of the slow dynamics. The specific heat has a finite jump towards the thermodynamic limit with critical exponents close to those expected from the hyperscaling and the random first-order transition theory for the glass transition. We also study an effective free energy of glasses, the Franz--Parisi potential, as a function of the overlap between equilibrium and quenched configurations. The effective free energy indicates the existence of a first-order phase transition, consistent with the random first-order transition theory. These findings strongly suggest that the glassy dynamics of the model has its origin in thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Fidelity Susceptibility Approach to Quantum Annealing of NP-hard problems

The computational complexity conjecture of NP $\nsubseteq$ BQP implies that there should be an exponentially small energy gap for Quantum Annealing (QA) of NP-hard problems. We aim to verify how this computation originated gapless point could be understood based on physics, using the quantum Monte Carlo method. As a result, we found a phase transition detectable only by the divergence of fidelity susceptibility. The exponentially small gapless points of each instance are all located in the phase found in this study, which suggests that this phase transition is the physical cause of the failure of QA for NP-hard problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Compressed sensing in scanning tunneling microscopy/spectroscopy for observation of quasi-particle interference

We applied a method of compressed sensing to the observation of quasi-particle interference (QPI) by scanning tunneling microscopy/spectroscopy to improve efficiency and save measurement time. To solve an ill-posed problem owing to the scarcity of data, the compressed sensing utilizes the sparseness of QPI patterns in momentum space. We examined the performance of a sparsity-inducing algorithm called least absolute shrinkage and selection operator (LASSO), and demonstrated that LASSO enables us to recover a double-circle QPI pattern of the Ag(111) surface from a dataset whose size is less than that necessary for the conventional Fourier transformation method. In addition, the smallest number of data required for the recovery is discussed on the basis of cross validation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Irreversible simulated tempering

An extended ensemble Monte Carlo algorithm is proposed by introducing a violation of the detailed balance condition to the update scheme of the inverse temperature in simulated tempering. Our method, irreversible simulated tempering, is constructed based on the framework of the skew detailed balance condition. By applying this method to the ferromagnetic Ising model in two dimensions on a square lattice as a benchmark, the dynamical behavior of the inverse temperature and an autocorrelation function of the magnetization are studied numerically. It is found that the relaxation dynamics of the inverse temperature changes qualitatively from diffusive to ballistic by violating the detailed balance condition. Consequently, the autocorrelation time of magnetization is several times smaller than that for the conventional algorithm satisfying the detailed balance condition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase transitions and ordering structures of a model of chiral helimagnet in three dimensions

Phase transitions in a classical Heisenberg spin model of a chiral helimagnet with the Dzyaloshinskii--Moriya (DM) interaction in three dimensions are numerically studied. By using the event-chain Monte Carlo algorithm recently developed for particle and continuous spin systems, we perform equilibrium Monte Carlo simulations for large systems up to about $10^6$ spins. Without magnetic fields, the system undergoes a continuous phase transition with critical exponents of the three-dimensional \textit{XY} model, and a uniaxial periodic helical structure emerges in the low temperature region. In the presence of a magnetic field perpendicular to the axis of the helical structure, it is found that there exists a critical point on the temperature and magnetic-field phase diagram and that above the critical point the system exhibits a phase transition with strong divergence of the specific heat and the uniform magnetic susceptibility.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Statistical-mechanical Analysis of Linear Programming Relaxation for Combinatorial Optimization Problems

Typical behavior of the linear programming (LP) problem is studied as a relaxation of the minimum vertex cover, a type of integer programming (IP) problem. A lattice-gas model on the Erdös-Rényi random graphs of $α$-uniform hyperedges is proposed to express both the LP and IP problems of the min-VC in the common statistical-mechanical model with a one-parameter family. Statistical-mechanical analyses reveal for $α=2$ that the LP optimal solution is typically equal to that given by the IP below the critical average degree $c=e$ in the thermodynamic limit. The critical threshold for good accuracy of the relaxation extends the mathematical result $c=1$, and coincides with the replica symmetry-breaking threshold of the IP. The LP relaxation for the minimum hitting sets with $α\geq 3$, minimum vertex covers on $α$-uniform random graphs, is also studied. Analytic and numerical results strongly suggest that the LP relaxation fails to estimate optimal values above the critical average degree $c=e/(α-1)$ where the replica symmetry is broken.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Typical Performance of Approximation Algorithms for NP-hard Problems

Typical performance of approximation algorithms is studied for randomized minimum vertex cover problems. A wide class of random graph ensembles characterized by an arbitrary degree distribution is discussed with some theoretical frameworks. Here three approximation algorithms are examined; the linear-programming relaxation, the loopy-belief propagation, and the leaf-removal algorithm. The former two algorithms are analyzed using the statistical-mechanical technique while the average-case analysis of the last one is studied by the generating function method. These algorithms have a threshold in the typical performance with increasing the average degree of the random graph, below which they find true optimal solutions with high probability. Our study reveals that there exist only three cases determined by the order of the typical-performance thresholds. We provide some conditions for classifying the graph ensembles and demonstrate explicitly examples for the difference in the threshold.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eigenvalue analysis of an irreversible random walk with skew detailed balance conditions

An irreversible Markov-chain Monte Carlo (MCMC) algorithm with skew detailed balance conditions originally proposed by Turitsyn et al. is extended to general discrete systems on the basis of the Metropolis-Hastings scheme. To evaluate the efficiency of our proposed method, the relaxation dynamics of the slowest mode and the asymptotic variance are studied analytically in a random walk on one dimension. It is found that the performance in irreversible MCMC methods violating the detailed balance condition is improved by appropriately choosing parameters in the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Event-chain algorithm for the Heisenberg model: Evidence for $z \simeq 1$ dynamic scaling

We apply the event-chain Monte Carlo algorithm to the three-dimensional ferromagnetic Heisenberg model. The algorithm is rejection-free and also realizes an irreversible Markov chain that satisfies global balance. The autocorrelation functions of the magnetic susceptibility and the energy indicate a dynamical critical exponent $z \approx 1$ at the critical temperature, while that of the magnetization does not measure the performance of the algorithm. This seems to be the first report that the event-chain Monte Carlo algorithm substantially reduces the dynamical critical exponent from the conventional value of $z\simeq 2$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evidence of one-step replica symmetry breaking in a three-dimensional Potts glass model

We study a 7-state Potts glass model in three dimensions with first, second, and third neighbor interactions with a bimodal distribution of couplings by Monte Carlo simulations. Our results show the existence of a spin-glass transition at a finite temperature T_c, a discontinuous jump of an order parameter at T_c without latent heat, and a non-trivial structure of the order-parameter distribution below T_c. They are compatible with a one-step replica symmetry breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Detection of the Ergodicity Breaking in a Lattice Glass Model

We directly detect the ergodicity breaking in a lattice glass model by a numerical simulation. The obtained results nicely agree with those by the cavity method that the model on a regular random graph exhibits a dynamical transition with the ergodicity breaking at an occupation density. The present method invented for a numerical detection of the ergodicity breaking is applicable to glassy systems in finite dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A List Referring Monte Carlo Method for Lattice Glass Models

We present an effcient Monte-Carlo method for lattice glass models which are characterized by hard constraint conditions. The basic idea of the method is similar to that of the $N$-fold way method. By using a list of sites into which we can insert a particle, we avoid trying a useless transition which is forbidden by the constraint conditions. We applied the present method to a lattice glass model proposed by Biroli and M{é}zard. We first evaluated the efficiency of the method through measurements of the autocorrelation function of particle configurations. As a result, we found that the efficiency is much higher than that of the standard Monte-Carlo method. We also compared the efficiency of the present method with that of the $N$-fold way method in detail. We next examined how the efficiency of extended ensemble methods such as the replica exchange method and the Wang-Landau method is inflnuenced by the choice of the local update method. The results show that the efficiency is considerably improved by the use of efficient local update methods. For example, when the number of sites $N_{\rm site}$ is 1024, the ergodic time $τ_{\rm E}$ of the replica exchange method in the grand-canonical ensemble, which is the average round-trip time of a replica in chemical-potential space, with the present local update method is more than $10^2$ times shorter than that with the standard local update method. This result shows that the efficient local update method is quite important to make extended ensemble methods more effective.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of the One-Dimensional Ising Model without Detailed Balance Condition

We study an irreversible Markov chain Monte Carlo method based on a skew detailed balance condition for an one-dimensional Ising model. Dynamical behavior of the magnetization density is analyzed in order to understand the properties of this method. As a result, it is found theoretically that the relaxation time of the magnetization density is reduced by using some transition probabilities satisfying the skew detailed balance condition, in comparison to that with the corresponding transition probability with the detailed balance condition, and that one of the transition probabilities changes the dynamical critical exponent even with a local spin update.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Minimum vertex cover problems on random hypergraphs: replica symmetric solution and a leaf removal algorithm

We study minimum vertex cover problems on random α-uniform hypergraphs using two different approaches, a replica method in statistical mechanics of random systems and a leaf removal algorithm. It is found that there exists a phase transition at the critical average degree e/(α-1). Below the critical degree, a replica symmetric ansatz in the statistical-mechanical method holdsand the algorithm estimates a solution of the problem which coincide with that by the replica method. In contrast, above the critical degree, the replica symmetric solution becomes unstable and these methods fail to estimate the exact solution.These results strongly suggest a close relation between the replica symmetry and the performance of approximation algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multicanonical sampling of rare events in random matrices

A method based on multicanonical Monte Carlo is applied to the calculation of large deviations in the largest eigenvalue of random matrices. The method is successfully tested with the Gaussian orthogonal ensemble (GOE), sparse random matrices, and matrices whose components are subject to uniform density. Specifically, the probability that all eigenvalues of a matrix are negative is estimated in these cases down to the values of $\sim 10^{-200}$, a region where naive random sampling is ineffective. The method can be applied to any ensemble of matrices and used for sampling rare events characterized by any statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Response to a twist on systems with Zp symmetry

We study response to a twist in the two-dimensional p-state clock model, which has the discrete Zp symmetry. The response is measured in terms of helicity modulus, which is usually defined with respect to an infinitesimal twist. However, we demonstrate that such a definition is inappropriate for the clock model. The helicity modulus must be defined with respect to a finite, quantized twist which matches the discrete Zp symmetry of the model. Recent numerical results, which casted doubt on the standard picture that two Berezinskii-Kosterlitz-Thouless transitions occur for p > 4, are rather a consequence of use of the inappropriate quantity. Numerical calculation of the appropriately defined helicity modulus indeed supports the standard picture.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Typical behavior of the linear programming method for combinatorial optimization problems: From a statistical-mechanical perspective

Typical behavior of the linear programming problem (LP) is studied as a relaxation of the minimum vertex cover problem, which is a type of the integer programming problem (IP). To deal with the LP and IP by statistical mechanics, a lattice-gas model on the Erdös-Rényi random graphs is analyzed by a replica method. It is found that the LP optimal solution is typically equal to that of the IP below the critical average degree c*=e in the thermodynamic limit. The critical threshold for LP=IP is beyond a mathematical result, c=1, and coincides with the replica-symmetry-breaking threshold of the IP.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Replica symmetry breaking in an adiabatic spin-glass model of adaptive evolution

We study evolutionary canalization using a spin-glass model with replica theory, where spins and their interactions are dynamic variables whose configurations correspond to phenotypes and genotypes, respectively. The spins are updated under temperature T_S, and the genotypes evolve under temperature T_J, according to the evolutionary fitness. It is found that adaptation occurs at T_S < T_S^{RS}, and a replica symmetric phase emerges at T_S^{RSB} < T_S < T_S^{RS}. The replica symmetric phase implies canalization, and replica symmetry breaking at lower temperatures indicates loss of robustness.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Partial annealing of a coupled mean-field spin-glass model with an embedded pattern

A partially annealed mean-field spin-glass model with a locally embedded pattern is studied. The model consists of two dynamical variables, spins and interactions, that are in contact with thermal baths at temperatures T_S and T_J, respectively. Unlike the quenched system, characteristic correlations among the interactions are induced by the partial annealing. The model exhibits three phases, which are paramagnetic, ferromagnetic and spin-glass phases. In the ferromagnetic phase, the embedded pattern is stably realized. The phase diagram depends significantly on the ratio of two temperatures n=T_J/T_S. In particular, a reentrant transition from the embedded ferromagnetic to the spin-glass phases with T_S decreasing is found only below at a certain value of n. This indicates that above the critical value n_c the embedded pattern is supported by local field from a non-embedded region. Some equilibrium properties of the interactions in the partial annealing are also discussed in terms of frustration.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A statistical-mechanical study of evolution of robustness in noisy environment

In biological systems, expression dynamics that can provide fitted phenotype patterns with respect to a specific function have evolved through mutations. This has been observed in the evolution of proteins for realizing folding dynamics through which a target structure is shaped. We study this evolutionary process by introducing a statistical-mechanical model of interacting spins, where a configuration of spins and their interactions $\bm{J}$ represent a phenotype and genotype, respectively. The phenotype dynamics are given by a stochastic process with temperature $T_{S}$ under a Hamiltonian with $\bm{J}$. The evolution of $\bm{J}$ is also stochastic with temperature $T_{J}$ and follows mutations introduced into $\bm{J}$ and selection based on a fitness defined for a configuration of a given set of target spins. Below a certain temperature $T_{S}^{c2}$, the interactions $\bm{J}$ that achieve the target pattern evolve, whereas another phase transition is observed at $T_{S}^{c1}<T_{S}^{c2}$. At low temperatures $T_{S}<T_{S}^{c1}$, the Hamiltonian exhibits a spin-glass like phase, where the dynamics toward the target pattern require long time steps, and the fitness often decreases drastically as a result of a single mutation to $\bm{J}$. In the intermediate-temperature region, the dynamics to shape the target pattern proceed rapidly and are robust to mutations of $\bm{J}$. The interactions in this region have no frustration around the target pattern and results in funnel-type dynamics. We propose that the ubiquity of funnel-type dynamics, as observed in protein folding, is a consequence of evolution subjected to thermal noise beyond a certain level; this also leads to mutational robustness of the fitness.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Thermodynamic Construction of an One-Step Replica-Symmetry-Breaking Solution in Finite Connectivity Spin Glasses

An one-step replica-symmetry-breaking solution for finite connectivity spin-glass models with K body interaction is constructed at finite temperature using the replica method and thermodynamic constraints. In the absence of external fields, this construction provides a general extension of replica symmetric solution at finite replica number to one-step replica-symmetry-breaking solution. It is found that this result is formally equivalent to that of the one-step replica-symmetry-breaking cavity method. To confirm the validity of the obtained solution, Monte Carlo simulations are performed for K = 2 and 3. The thermodynamic quantities of the Monte Carlo results extrapolated to a large-size limit are consistent with those estimated by our solution for K = 2 at all simulated temperatures and for K = 3 except near the transition temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Computational Experiment on Glassy Dynamic Nature of the Field-Cooled Magnetization in an Ising Spin-Glass Model

The field-cooled (FC) process on an Ising spin-glass model is investigated by a standard Monte Calro (MC) method on one hand, and the equilibrium magnetization of the same system is evaluated by the exchang MC method on the other hand. The two types of simulation reveal intriguing glassy (nonequilibrium) dynamic properties of the FC magnetization (FCM) of the system. Particularly, the FCM decrease is observed when the FC process is halted at a low temperater, although its value is smaller than the corresponding equilibrium value. We present a comprehensive interpretation of such peculiar phenomena based on the scenario for the FCM process recently propsed by Jönsson and one of us (HT).

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Dynamical Critical Phenomena in three-dimensional Heisenberg Spin Glasses

Spin-glass (SG) and chiral-glass (CG) orderings in three dimensional (3D) Heisenberg spin glass with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo simulations. A characteristic time of relaxation, which diverges at a transition temperature in the thermodynamic limit, is obtained as a function of the temperature and the system size. Based on the finite-size scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs at a lower temperature, which is compatible with zero. Our results of the anisotropic case support the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses.

0 citations0 reviewsNo code linkedOpen access