Source author record

Yutaka Okabe

Yutaka Okabe 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

17works
6topics
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

17 published item(s)

preprint2020arXiv

Machine-Learning Studies on Spin Models

With the recent developments in machine learning, Carrasquilla and Melko have proposed a paradigm that is complementary to the conventional approach for the study of spin models. As an alternative to investigating the thermal average of macroscopic physical quantities, they have used the spin configurations for the classification of the disordered and ordered phases of a phase transition through machine learning. We extend and generalize this method. We focus on the configuration of the long-range correlation function instead of the spin configuration itself, which enables us to provide the same treatment to multi-component systems and the systems with a vector order parameter. We analyze the Berezinskii-Kosterlitz-Thouless (BKT) transition with the same technique to classify three phases: the disordered, the BKT, and the ordered phases. We also present the classification of a model using the training data of a different model.

preprint2020arXiv

Machine-Learning Study using Improved Correlation Configuration and Application to Quantum Monte Carlo Simulation

We use the Fortuin-Kasteleyn representation based improved estimator of the correlation configuration as an alternative to the ordinary correlation configuration in the machine-learning study of the phase classification of spin models. The phases of classical spin models are classified using the improved estimators, and the method is also applied to the quantum Monte Carlo simulation using the loop algorithm. We analyze the Berezinskii-Kosterlitz-Thouless (BKT) transition of the spin 1/2 quantum XY model on the square lattice. We classify the BKT phase and the paramagnetic phase of the quantum XY model using the machine-learning approach. We show that the classification of the quantum XY model can be performed by using the training data of the classical XY model.

preprint2015arXiv

Probing phase transition order of $q$-state Potts models using Wang-Landau Algorithm

Phase transitions are ubiquitous phenomena, exemplified by the melting of ice and spontaneous magnetization of magnetic material. In general, a phase transition is associated with a symmetry breaking of a system; occurs due to the competition between coupling interaction and external fields such as thermal energy. If the phase transition occurs with no latent heat, the system experiences continuous transition, also known as second order phase transition. The ferromagnetic $q$-state Potts model with $r$ extra invisible states, introduced by Tamura, Tanaka, and Kawashima [Prog. Theor. Phys. 124, 381 (2010)], is studied by using the Wang-Landau method. The density of states difference (DOSD), $\ln g(E +ΔE) - \ln g(E)$, is used to investigate the order of the phase transition and examine the critical value of $r$ changing the second to the first order transition.

preprint2014arXiv

CUDA programs for GPU computing of Swendsen-Wang multi-cluster spin flip algorithm: 2D and 3D Ising, Potts, and XY models

We present sample CUDA programs for the GPU computing of the Swendsen-Wang multi-cluster spin flip algorithm. We deal with the classical spin models; the Ising model, the $q$-state Potts model, and the classical XY model. As for the lattice, both the 2D (square) lattice and the 3D (simple cubic) lattice are treated. We already reported the idea of the GPU implementation for 2D models [Comput. Phys. Commun. 183 (2012) 1155-1161]. We here explain the details of sample programs, and discuss the performance of the present GPU implementation for the 3D Ising and XY models. We also show the calculated results of the moment ratio for these models, and discuss phase transitions.

preprint2013arXiv

Dynamically Optimized Wang-Landau Sampling with Adaptive Trial Moves and Modification Factors

The density of states of continuous models is known to span many orders of magnitudes at different energies due to the small volume of phase space near the ground state. Consequently, the traditional Wang-Landau sampling which uses the same trial move for all energies faces difficulties sampling the low entropic states. We developed an adaptive variant of the Wang-Landau algorithm that very effectively samples the density of states of continuous models across the entire energy range. By extending the acceptance ratio method of Bouzida, Kumar, and Swendsen such that the step size of the trial move and acceptance rate are adapted in an energy-dependent fashion, the random walker efficiently adapts its sampling according to the local phase space structure. The Wang-Landau modification factor is also made energy-dependent in accordance with the step size, enhancing the accumulation of the density of states. Numerical simulations show that our proposed method performs much better than the traditional Wang-Landau sampling.

preprint2012arXiv

Difference of energy density of states in the Wang-Landau algorithm

Paying attention to the difference of density of states, Δln g(E) = ln g(E+ΔE) - ln g(E), we study the convergence of the Wang-Landau method. We show that this quantity is a good estimator to discuss the errors of convergence, and refer to the $1/t$ algorithm. We also examine the behavior of the 1st-order transition with this difference of density of states in connection with Maxwell's equal area rule. A general procedure to judge the order of transition is given.

preprint2012arXiv

GPU-based single-cluster algorithm for the simulation of the Ising model

We present the GPU calculation with the common unified device architecture (CUDA) for the Wolff single-cluster algorithm of the Ising model. Proposing an algorithm for a quasi-block synchronization, we realize the Wolff single-cluster Monte Carlo simulation with CUDA. We perform parallel computations for the newly added spins in the growing cluster. As a result, the GPU calculation speed for the two-dimensional Ising model at the critical temperature with the linear size L=4096 is 5.60 times as fast as the calculation speed on a current CPU core. For the three-dimensional Ising model with the linear size L=256, the GPU calculation speed is 7.90 times as fast as the CPU calculation speed. The idea of quasi-block synchronization can be used not only in the cluster algorithm but also in many fields where the synchronization of all threads is required.

preprint2012arXiv

GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional classical spin systems

We present the GPU calculation with the common unified device architecture (CUDA) for the Swendsen-Wang multi-cluster algorithm of two-dimensional classical spin systems. We adjust the two connected component labeling algorithms recently proposed with CUDA for the assignment of the cluster in the Swendsen-Wang algorithm. Starting with the q-state Potts model, we extend our implementation to the system of vector spins, the q-state clock model, with the idea of embedded cluster. We test the performance, and the calculation time on GTX580 is obtained as 2.51 nano sec per a spin flip for the q=2 Potts model (Ising model) and 2.42 nano sec per a spin flip for the q=6 clock model with the linear size L=4096 at the critical temperature, respectively. The computational speed for the q=2 Potts model on GTX580 is 12.4 times as fast as the calculation speed on a current CPU core. That for the q=6 clock model on GTX580 is 35.6 times as fast as the calculation speed on a current CPU core.

preprint2012arXiv

Large-scale Monte Carlo simulation of two-dimensional classical XY model using multiple GPUs

We study the two-dimensional classical XY model by the large-scale Monte Carlo simulation of the Swendsen-Wang multi-cluster algorithm using multiple GPUs on the open science supercomputer TSUBAME 2.0. Simulating systems up to the linear system size L=65536, we investigate the Kosterlitz-Thouless (KT) transition. Using the generalized version of the probability-changing cluster algorithm based on the helicity modulus, we locate the KT transition temperature in a self-adapted way. The obtained inverse KT temperature β_{KT} is 1.11996(6). We estimate the exponent to specify the multiplicative logarithmic correction, -2r, and precisely reproduce the theoretical prediction -2r=1/8.

preprint2012arXiv

Multi-GPU-based Swendsen-Wang multi-cluster algorithm for the simulation of two-dimensional q-state Potts model

We present the multiple GPU computing with the common unified device architecture (CUDA) for the Swendsen-Wang multi-cluster algorithm of two-dimensional (2D) q-state Potts model. Extending our algorithm for single GPU computing [Comp. Phys. Comm. 183 (2012) 1155], we realize the GPU computation of the Swendsen-Wang multi-cluster algorithm for multiple GPUs. We implement our code on the large-scale open science supercomputer TSUBAME 2.0, and test the performance and the scalability of the simulation of the 2D Potts model. The performance on Tesla M2050 using 256 GPUs is obtained as 37.3 spin flips per a nano second for the q=2 Potts model (Ising model) at the critical temperature with the linear system size L=65536.

preprint2010arXiv

Phase transition of two-dimensional generalized XY model

We study the two-dimensional generalized XY model that depends on an integer $q$ by the Monte Carlo method. This model was recently proposed by Romano and Zagrebnov. We find a single Kosterlitz-Thouless (KT) transition for all values of $q$, in contrast with the previous speculation that there may be two transitions, one a regular KT transition and another a first-order transition at a higher temperature. We show the universality of the KT transitions by comparing the universal finite-size scaling behaviors at different values of $q$ without assuming a specific universal form in terms of the KT transition temperature $T_{\rm KT}$.

preprint2010arXiv

Stability of critical bubble in stretched fluid of square-gradient density-functional model with triple-parabolic free energy

The square-gradient density-functional model with triple-parabolic free energy, that was used previously to study the homogeneous bubble nucleation [J. Chem. Phys. 129, 104508 (2008)], is used to study the stability of the critical bubble nucleated within the bulk under-saturated stretched fluid. The stability of the bubble is studied by solving the Schrödinger equation for the fluctuation. The negative eigenvalue corresponds to the unstable growing mode of the fluctuation. Our results show that there is only one negative eigenvalue whose eigenfunction represents the fluctuation that corresponds to the isotropically growing or shrinking nucleus. In particular, this negative eigenvalue survives up to the spinodal point. Therefore the critical bubble is not fractal or ramified near the spinodal.

preprint2004arXiv

Study of the Fully Frustrated Clock Model using the Wang-Landau Algorithm

Monte Carlo simulations using the newly proposed Wang-Landau algorithm together with the broad histogram relation are performed to study the antiferromagnetic six-state clock model on the triangular lattice, which is fully frustrated. We confirm the existence of the magnetic ordering belonging to the Kosterlitz-Thouless (KT) type phase transition followed by the chiral ordering which occurs at slightly higher temperature. We also observe the lower temperature phase transition of KT type due to the discrete symmetry of the clock model. By using finite-size scaling analysis, the higher KT temperature $T_2$ and the chiral critical temperature $T_c$ are respectively estimated as $T_2=0.5154(8)$ and $T_c=0.5194(4)$. The results are in favor of the double transition scenario. The lower KT temperature is estimated as $T_1=0.496(2)$. Two decay exponents of KT transitions corresponding to higher and lower temperatures are respectively estimated as $η_2=0.25(1)$ and $η_1=0.13(1)$, which suggests that the exponents associated with the KT transitions are universal even for the frustrated model.

preprint1999arXiv

Application of Monte Carlo Method to Phase Separation Dynamics of Complex Systems

We report the application of the Monte Carlo simulation to phase separation dynamics. First, we deal with the phase separation under shear flow. The thermal effect on the phase separation is discussed, and the anisotropic growth exponents in the late stage are estimated. Next, we study the effect of surfactants on the three-component solvents. We obtain the mixture of macrophase separation and microphase separation, and investigate the dynamics of both phase separations.

preprint1996arXiv

Statistical Dependence Analysis

We review our recent studies on the dynamical correlations in MC simulations from the view point of the statistical dependence. Attentions are paid to the reduction of the statistical degrees of freedom for correlated data. Possible biases on several cumulants, such as the susceptibility and the Binder number due to finite MC length are discussed. A new method for calculating the equilibrium relaxation time from the analysis of the statistical dependence is presented. We apply it to the critical dynamics of the Ising model to estimate the dynamical critical exponent accurately.

preprint1996arXiv

Universal Finite-Size-Scaling Functions

The idea of universal finite-size-scaling functions of the Ising model is tested by Monte Carlo simulations for various lattices. Not only regular lattices such as the square lattice but quasiperiodic lattices such as the Penrose lattice are treated. We show that the finite-size-scaling functions of the order parameter for various lattices are collapsed on a single curve by choosing two nonuniversal scaling metric factors. We extend the idea of the universal finite-size-scaling functions to the order-parameter distribution function. We pay attention to the effects of boundary conditions. Keywords: Universal Finite-Size-Scaling Function; Ising Model; Order-Parameter Probability Distribution Function.

preprint1993arXiv

Recycle of random sequences

The correlation between a random sequence and its transformed sequences is studied. In the case of a permutation operation or, in other word, the shuffling operation, it is shown that the correlation can be so small that the sequences can be regarded as independent random sequences. The applications to the Monte Carlo simulations are also given. This method is especially useful in the Ising Monte Carlo simulation.