-
Quantum-centric Supercomputing for Materials Science: A Perspective on Challenges and Future Directions
Authors:
Yuri Alexeev,
Maximilian Amsler,
Paul Baity,
Marco Antonio Barroca,
Sanzio Bassini,
Torey Battelle,
Daan Camps,
David Casanova,
Young jai Choi,
Frederic T. Chong,
Charles Chung,
Chris Codella,
Antonio D. Corcoles,
James Cruise,
Alberto Di Meglio,
Jonathan Dubois,
Ivan Duran,
Thomas Eckl,
Sophia Economou,
Stephan Eidenbenz,
Bruce Elmegreen,
Clyde Fare,
Ismael Faro,
Cristina Sanz Fernández,
Rodrigo Neumann Barros Ferreira
, et al. (102 additional authors not shown)
Abstract:
Computational models are an essential tool for the design, characterization, and discovery of novel materials. Hard computational tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their simulation, analysis, and data resources. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of…
▽ More
Computational models are an essential tool for the design, characterization, and discovery of novel materials. Hard computational tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their simulation, analysis, and data resources. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of the computational tasks needed for materials science. In order to do that, the quantum technology must interact with conventional high-performance computing in several ways: approximate results validation, identification of hard problems, and synergies in quantum-centric supercomputing. In this paper, we provide a perspective on how quantum-centric supercomputing can help address critical computational problems in materials science, the challenges to face in order to solve representative use cases, and new suggested directions.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Update of $\mathcal{H}Φ$: Newly added functions and methods in versions 2 and 3
Authors:
Kota Ido,
Mitsuaki Kawamura,
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Youhei Yamaji,
Synge Todo,
Naoki Kawashima,
Takahiro Misawa
Abstract:
$\mathcal{H}Φ$ [$aitch$-$phi…
▽ More
$\mathcal{H}Φ$ [$aitch$-$phi$] is an open-source software package of numerically exact and stochastic calculations for a wide range of quantum many-body systems. In this paper, we present the newly added functions and the implemented methods in vers. 2 and 3. In ver. 2, we implement spectrum calculations by the shifted Krylov method, and low-energy excited state calculations by the locally optimal blocking preconditioned conjugate gradient (LOBPCG) method. In ver. 3, we implement the full diagonalization method using ScaLAPACK and GPGPU computing via MAGMA. We also implement a real-time evolution method and the canonical thermal pure quantum (cTPQ) state method for finite-temperature calculations. The Wannier90 format for specifying the Hamiltonians is also implemented. Using the Wannier90 format, it is possible to perform the calculations for the $ab$ $initio$ low-energy effective Hamiltonians of solids obtained by the open-source software RESPACK. We also update Standard mode $\unicode{x2014}$simplified input format in $\mathcal{H}Φ$$\unicode{x2014}$ to use these functions and methods. We explain the basics of the implemented methods and how to use them.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Absorbing Phase Transitions in Artificial Deep Neural Networks
Authors:
Keiichi Tamai,
Tsuyoshi Okubo,
Truong Vinh Truong Duy,
Naotake Natori,
Synge Todo
Abstract:
Theoretical understanding of the behavior of infinitely-wide neural networks has been rapidly developed for various architectures due to the celebrated mean-field theory. However, there is a lack of a clear, intuitive framework for extending our understanding to finite networks that are of more practical and realistic importance. In the present contribution, we demonstrate that the behavior of pro…
▽ More
Theoretical understanding of the behavior of infinitely-wide neural networks has been rapidly developed for various architectures due to the celebrated mean-field theory. However, there is a lack of a clear, intuitive framework for extending our understanding to finite networks that are of more practical and realistic importance. In the present contribution, we demonstrate that the behavior of properly initialized neural networks can be understood in terms of universal critical phenomena in absorbing phase transitions. More specifically, we study the order-to-chaos transition in the fully-connected feedforward neural networks and the convolutional ones to show that (i) there is a well-defined transition from the ordered state to the chaotics state even for the finite networks, and (ii) difference in architecture is reflected in that of the universality class of the transition. Remarkably, the finite-size scaling can also be successfully applied, indicating that intuitive phenomenological argument could lead us to semi-quantitative description of the signal propagation dynamics.
△ Less
Submitted 5 July, 2023;
originally announced July 2023.
-
Stochastic approximation analysis of dynamical quantum critical phenomena in long-range transverse-field Ising chain
Authors:
Sora Shiratani,
Synge Todo
Abstract:
The quantum phase transition of the long-range transverse-field Ising model is explored by combining a quantum Monte Carlo method with the optimal computational complexity scaling and stochastic parameter optimization that renders space and imaginary time isotropic, specifically achieved by tuning correlation lengths. Varying the decay rate of the long-range interaction, we exhaustively calculate…
▽ More
The quantum phase transition of the long-range transverse-field Ising model is explored by combining a quantum Monte Carlo method with the optimal computational complexity scaling and stochastic parameter optimization that renders space and imaginary time isotropic, specifically achieved by tuning correlation lengths. Varying the decay rate of the long-range interaction, we exhaustively calculate the dynamical critical exponent and the other exponents precisely in mean-field, nonuniversal, and Ising universality regimes. In our simulations, critical properties are extracted only from a set of simulations with different $L$, significantly improving computational cost compared to the standard finite-size scaling approach based on data collapse. We also perform a hypothesis test at the predicted universality boundary, which supports preceding reports arguing that conventional theoretical prediction fails to locate the universality boundary.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
A Noise-Robust Data Assimilation Method for Crystal Structure Prediction Using Powder Diffraction Intensity
Authors:
Seiji Yoshikawa,
Ryuhei Sato,
Ryosuke Akashi,
Synge Todo,
Shinji Tsuneyuki
Abstract:
Crystal structure prediction for a given chemical composition has long been a challenge in condensed-matter science. We have recently shown that experimental powder X-ray diffraction (XRD) data are helpful in a crystal structure search using simulated annealing, even when they are insufficient for structure determination by themselves (N. Tsujimoto et al., Phys. Rev. Materials 2, 053801 (2018)). I…
▽ More
Crystal structure prediction for a given chemical composition has long been a challenge in condensed-matter science. We have recently shown that experimental powder X-ray diffraction (XRD) data are helpful in a crystal structure search using simulated annealing, even when they are insufficient for structure determination by themselves (N. Tsujimoto et al., Phys. Rev. Materials 2, 053801 (2018)). In the method, the XRD data are assimilated into the simulation by adding a penalty function to the physical potential energy, where we used a crystallinity-type penalty function defined by the difference between experimental and simulated diffraction angles. To improve the success rate and noise robustness, we introduce a correlation-coefficient-type penalty function adaptable to XRD data with significant experimental noise. We apply the new penalty function to SiO$_2$ coesite and $ε$-Zn(OH)$_2$ to determine its effectiveness in the data assimilation method.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
MateriApps LIVE! and MateriApps Installer: Environment for starting and scaling up materials science simulations
Authors:
Yuichi Motoyama,
Kazuyoshi Yoshimi,
Takeo Kato,
Synge Todo
Abstract:
In our current era, numerical simulations have become indispensable theoretical and experimental tools for use in daily research activities, particularly in the materials science fields. However, the installation processes for such simulations frequently become problematic because they depend strongly on the device environment, and troubleshooting those processes is a challenging task for beginner…
▽ More
In our current era, numerical simulations have become indispensable theoretical and experimental tools for use in daily research activities, particularly in the materials science fields. However, the installation processes for such simulations frequently become problematic because they depend strongly on the device environment, and troubleshooting those processes is a challenging task for beginners. To minimize such difficulties, we created MateriApps LIVE! and MateriApps Installer, which can solve most of the related issues. Specifically, MateriApps LIVE! offers a virtual environment in which users can quickly try out computational materials science simulations on a personal computer while MateriApps Installer provides a comprehensive set of shell scripts for use when installing software on Unix, Linux, macOS, and supercomputer systems. Herein, we provide detailed descriptions of MateriApps LIVE! and MateriApps Installer together with illustrative examples of their use.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
Nested Iterative Shift-invert Diagonalization for Many-body Localization in the Random-field Heisenberg Chain
Authors:
Taito Kutsuzawa,
Synge Todo
Abstract:
We study the many-body localization of the random-field Heisenberg chain using the nested shift-invert Lanczos method with an iterative linear solver. We use the minimum residual method (MINRES) inside each Lanczos iteration. The memory consumption of the proposed method is only proportional to the dimension of the Hilbert space. We also introduce a preconditioner that takes into account the effec…
▽ More
We study the many-body localization of the random-field Heisenberg chain using the nested shift-invert Lanczos method with an iterative linear solver. We use the minimum residual method (MINRES) inside each Lanczos iteration. The memory consumption of the proposed method is only proportional to the dimension of the Hilbert space. We also introduce a preconditioner that takes into account the effects of disorder and interaction in the random-field Heisenberg chain. As a probe of many-body localization transition, we propose a unitary operator called the twist operator, which has a clear interpretation in the real space. We discuss its behavior for thermal and localized eigenstates. We demonstrate the efficiency of the nested iterative shift-invert diagonalization method with the proposed preconditioner for the many-body localization problem and estimate the transition point of the random-field Heisenberg chain more precisely based on the finite-size analysis of the expectation value of the twist operator.
△ Less
Submitted 30 March, 2022; v1 submitted 18 March, 2022;
originally announced March 2022.
-
Randomized-gauge test for machine learning of Ising model order parameter
Authors:
Tomoyuki Morishita,
Synge Todo
Abstract:
Recently, machine learning has been applied successfully for identifying phases and phase transitions of the Ising models. The continuous phase transition is characterized by spontaneous symmetry breaking, which can not be detected in general from a single spin configuration. To investigate if neural networks can extract correlations among spin snapshots, we propose a new test using the random-gau…
▽ More
Recently, machine learning has been applied successfully for identifying phases and phase transitions of the Ising models. The continuous phase transition is characterized by spontaneous symmetry breaking, which can not be detected in general from a single spin configuration. To investigate if neural networks can extract correlations among spin snapshots, we propose a new test using the random-gauge Ising model. We show that neural networks can extract the order parameter or the energy of the random-gauge model as in the ferromagnetic case. We also discuss how and where the information of random gauge is encoded in neural networks and attempt to reconstruct the gauge from the neural network parameters. We find that the fully connected network encodes the effect of random gauge to its weights naturally. In contrast, the convolutional network copes with the randomness by assigning different network parts to local gauge patterns. This observation indicates that although the latter demonstrates higher performance than the former for the present randomized-gauge test, the former is more effective and suitable for dealing with models with spatial randomness.
△ Less
Submitted 12 September, 2021;
originally announced September 2021.
-
Fast-update in self-learning algorithm for continuous-time quantum Monte Carlo
Authors:
Ruixiao Cao,
Synge Todo
Abstract:
We propose a novel technique for speeding up the self-learning Monte Carlo method applied to the single-site impurity model. For the case where the effective Hamiltonian is expressed by polynomial functions of differences of imaginary-time coordinate between vertices, we can remove the dependence of CPU time on the number of vertices, $n$, by saving and updating some coefficients for each insertio…
▽ More
We propose a novel technique for speeding up the self-learning Monte Carlo method applied to the single-site impurity model. For the case where the effective Hamiltonian is expressed by polynomial functions of differences of imaginary-time coordinate between vertices, we can remove the dependence of CPU time on the number of vertices, $n$, by saving and updating some coefficients for each insertion and deletion process. As a result, the total cost for a single-step update is drastically reduced from $O(nm)$ to $O(m^2)$ with $m$ being the order of polynomials in the effective Hamiltonian. Even for the existing algorithms, in which the absolute value is used instead of the difference as the variable of polynomial functions, we can limit the CPU time for a single step of Monte Carlo update to $O(m^2 + m \log n)$ with the help of balanced binary search trees. We demonstrate that our proposed algorithm with only logarithmic $n$-dependence achieves an exponential speedup from the existing methods, which suffer from severe performance issues at low temperatures.
△ Less
Submitted 22 June, 2021;
originally announced June 2021.
-
Bond-weighted Tensor Renormalization Group
Authors:
Daiki Adachi,
Tsuyoshi Okubo,
Synge Todo
Abstract:
We propose an improved tensor renormalization group (TRG) algorithm, the bond-weighted TRG (BTRG). In BTRG, we generalize the conventional TRG by introducing bond weights on the edges of the tensor network. We show that BTRG outperforms the conventional TRG and the higher-order tensor renormalization group with the same bond dimension, while its computation time is almost the same as that of TRG.…
▽ More
We propose an improved tensor renormalization group (TRG) algorithm, the bond-weighted TRG (BTRG). In BTRG, we generalize the conventional TRG by introducing bond weights on the edges of the tensor network. We show that BTRG outperforms the conventional TRG and the higher-order tensor renormalization group with the same bond dimension, while its computation time is almost the same as that of TRG. Furthermore, BTRG can have non-trivial fixed-point tensors at an optimal hyperparameter. We demonstrate that the singular value spectrum obtained by BTRG is invariant under the renormalization procedure in the case of the two-dimensional Ising model at the critical point. This property indicates that BTRG performs the tensor contraction with high accuracy while keeping the scale-invariant structure of tensors.
△ Less
Submitted 3 November, 2020;
originally announced November 2020.
-
Multithreaded event-chain Monte Carlo with local times
Authors:
Botao Li,
Synge Todo,
A. C. Maggs,
Werner Krauth
Abstract:
We present a multithreaded event-chain Monte Carlo algorithm (ECMC) for hard spheres. Threads synchronize at infrequent breakpoints and otherwise scan for local horizon violations. Using a mapping onto absorbing Markov chains, we rigorously prove the correctness of a sequential-consistency implementation for small test suites. On x86 and ARM processors, a C++ (OpenMP) implementation that uses comp…
▽ More
We present a multithreaded event-chain Monte Carlo algorithm (ECMC) for hard spheres. Threads synchronize at infrequent breakpoints and otherwise scan for local horizon violations. Using a mapping onto absorbing Markov chains, we rigorously prove the correctness of a sequential-consistency implementation for small test suites. On x86 and ARM processors, a C++ (OpenMP) implementation that uses compare-and-swap primitives for data access achieves considerable speed-up with respect to single-threaded code. The generalized birthday problem suggests that for the number of threads scaling as the square root of the number of spheres, the horizon-violation probability remains small for a fixed simulation time. We provide C++ and Python open-source code that reproduces all our results.
△ Less
Submitted 23 April, 2020;
originally announced April 2020.
-
K$ω$ -- Open-source library for the shifted Krylov subspace method of the form $(zI-H)x=b$
Authors:
Takeo Hoshi,
Mitsuaki Kawamura,
Kazuyoshi Yoshimi,
Yuichi Motoyama,
Takahiro Misawa,
Youhei Yamaji,
Synge Todo,
Naoki Kawashima,
Tomohiro Sogabe
Abstract:
We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundatio…
▽ More
We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundation of the shifted Krylov subspace methods. Applications in materials science are presented to demonstrate the advantages of the algorithm over the standard Krylov subspace methods such as the Lanczos method. We introduce benchmark calculations of (i) an excited (optical) spectrum and (ii) intermediate eigenvalues by the contour integral on the complex plane. In combination with the quantum lattice solver $\mathcal{H} Φ$, K$ω$ can realize parallel computation of excitation spectra and intermediate eigenvalues for various quantum lattice models.
△ Less
Submitted 28 May, 2020; v1 submitted 23 January, 2020;
originally announced January 2020.
-
Dynamical phase transition in Floquet optical bistable systems: An approach from finite-size quantum systems
Authors:
Tatsuhiko Shirai,
Synge Todo,
Seiji Miyashita
Abstract:
We study a dynamical phase transition in optical bistable systems subject to a time-periodic driving field. The phase transition occurs in the structure of limit cycle as a function of the frequency of the driving field. In the thermodynamic limit, a single limit cycle is divided into two separated limit cycles at the transition point. In finite-size systems, however, there is always a single limi…
▽ More
We study a dynamical phase transition in optical bistable systems subject to a time-periodic driving field. The phase transition occurs in the structure of limit cycle as a function of the frequency of the driving field. In the thermodynamic limit, a single limit cycle is divided into two separated limit cycles at the transition point. In finite-size systems, however, there is always a single limit cycle due to the quantum tunneling effect. We use a Floquet dissipative map, which is a time-evolution operator over one period in a dynamics given by a quantum master equation, and discuss the decay rate of relaxation dynamics into the limit cycle based on the dominant eigenvalue of the map. We found that the decay rate exhibits qualitatively different system-size dependence before and after the phase transition, and it shows a finite-size scaling of spinodal phenomena around the transition point. The present work provides a systematic way of studying dynamical phase transition observed in time-periodically driven open systems in terms of the Floquet dissipative map.
△ Less
Submitted 11 December, 2019; v1 submitted 23 October, 2019;
originally announced October 2019.
-
Anisotropic Tensor Renormalization Group
Authors:
Daiki Adachi,
Tsuyoshi Okubo,
Synge Todo
Abstract:
We propose a new tensor renormalization group algorithm, Anisotropic Tensor Renormalization Group (ATRG), for lattice models in arbitrary dimensions. The proposed method shares the same versatility with the Higher-Order Tensor Renormalization Group (HOTRG) algorithm, i.e., it preserves the lattice topology after the renormalization. In comparison with HOTRG, both of the computation cost and the me…
▽ More
We propose a new tensor renormalization group algorithm, Anisotropic Tensor Renormalization Group (ATRG), for lattice models in arbitrary dimensions. The proposed method shares the same versatility with the Higher-Order Tensor Renormalization Group (HOTRG) algorithm, i.e., it preserves the lattice topology after the renormalization. In comparison with HOTRG, both of the computation cost and the memory footprint of our method are drastically reduced, especially in higher dimensions, by renormalizing tensors in an anisotropic way after the singular value decomposition. We demonstrate the ability of ATRG for the square lattice and the simple cubic lattice Ising models. Although the accuracy of the present method degrades when compared with HOTRG of the same bond dimension, the accuracy with fixed computation time is improved greatly due to the drastic reduction of the computation cost.
△ Less
Submitted 13 September, 2019; v1 submitted 5 June, 2019;
originally announced June 2019.
-
Updated Core Libraries of the ALPS Project
Authors:
Markus Wallerberger,
Sergei Iskakov,
Alexander Gaenko,
Joseph Kleinhenz,
Igor Krivenko,
Ryan Levy,
Jia Li,
Hiroshi Shinaoka,
Synge Todo,
Tianran Chen,
Xi Chen,
James P. F. LeBlanc,
Joseph E. Paki,
Hanna Terletska,
Matthias Troyer,
Emanuel Gull
Abstract:
The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated electron systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation codes, and the applications strive to…
▽ More
The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated electron systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation codes, and the applications strive to make commonly used and proven computational algorithms available to a non-expert community. In this paper we present an update of the core ALPS libraries. We present in particular new Monte Carlo libraries and new Green's function libraries.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.
-
Parallel loop cluster quantum Monte Carlo simulation of quantum magnets based on global union-find graph algorithm
Authors:
Synge Todo,
Haruhiko Matsuo,
Hideyuki Shitara
Abstract:
A large-scale parallel loop cluster quantum Monte Carlo simulation is presented. On 24,576 nodes of the K computer, one loop cluster Monte Carlo update of the world-line configuration of the $S=1/2$ antiferromagnetic Heisenberg chain with $2.6 \times 10^6$ spins at inverse temperature $3.1 \times 10^5$ is executed in about 8.62 seconds, in which global union-find cluster identification on a graph…
▽ More
A large-scale parallel loop cluster quantum Monte Carlo simulation is presented. On 24,576 nodes of the K computer, one loop cluster Monte Carlo update of the world-line configuration of the $S=1/2$ antiferromagnetic Heisenberg chain with $2.6 \times 10^6$ spins at inverse temperature $3.1 \times 10^5$ is executed in about 8.62 seconds, in which global union-find cluster identification on a graph of about 1.1 trillion vertices and edges is performed. By combining the nonlocal global updates and the large-scale parallelization, we have virtually achieved about $10^{13}$-fold speed-up from the conventional local update Monte Carlo simulation performed on a single core. We have estimated successfully the antiferromagnetic correlation length and the magnitude of the first excitation gap of the $S=4$ antiferromagnetic Heisenberg chain for the first time as $ξ= 1.040(7) \times 10^4$ and $Δ= 7.99(5) \times 10^{-4}$, respectively.
△ Less
Submitted 17 October, 2018;
originally announced October 2018.
-
Search for Common Minima in Joint Optimization of Multiple Cost Functions
Authors:
Daiki Adachi,
Naoto Tsujimoto,
Ryosuke Akashi,
Synge Todo,
Shinji Tsuneyuki
Abstract:
We present a novel optimization method, named the Combined Optimization Method (COM), for the joint optimization of two or more cost functions. Unlike the conventional joint optimization schemes, which try to find minima in a weighted sum of cost functions, the COM explores search space for common minima shared by all the cost functions. Given a set of multiple cost functions that have qualitative…
▽ More
We present a novel optimization method, named the Combined Optimization Method (COM), for the joint optimization of two or more cost functions. Unlike the conventional joint optimization schemes, which try to find minima in a weighted sum of cost functions, the COM explores search space for common minima shared by all the cost functions. Given a set of multiple cost functions that have qualitatively different distributions of local minima with each other, the proposed method finds the common minima with a high success rate without the help of any metaheuristics. As a demonstration, we apply the COM to the crystal structure prediction in materials science. By introducing the concept of data assimilation, i.e., adopting the theoretical potential energy of the crystal and the crystallinity, which characterizes the agreement with the theoretical and experimental X-ray diffraction patterns, as cost functions, we show that the correct crystal structures of Si diamond, low quartz, and low cristobalite can be predicted with significantly higher success rates than the previous methods.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
Fast algorithm for generating random bit strings and multispin coding for directed percolation
Authors:
Hiroshi Watanabe,
Satoshi Morita,
Synge Todo,
Naoki Kawashima
Abstract:
We present efficient algorithms to generate a bit string in which each bit is set with arbitrary probability. By adopting a hybrid algorithm, i.e., a finite-bit density approximation with correction techniques, we achieve 3.8 times faster random bit generation than the simple algorithm for the 32-bit case and 6.8 times faster for the 64-bit case. Employing the developed algorithm, we apply the mul…
▽ More
We present efficient algorithms to generate a bit string in which each bit is set with arbitrary probability. By adopting a hybrid algorithm, i.e., a finite-bit density approximation with correction techniques, we achieve 3.8 times faster random bit generation than the simple algorithm for the 32-bit case and 6.8 times faster for the 64-bit case. Employing the developed algorithm, we apply the multispin coding technique to one-dimensional bond-directed percolation. The simulations are accelerated by up to a factor of 14 compared with an optimized scalar implementation. The random bit string generation algorithm proposed here is applicable to general Monte Carlo methods.
△ Less
Submitted 29 November, 2018; v1 submitted 15 August, 2018;
originally announced August 2018.
-
Localized Mode and Nonergodicity of a Harmonic Oscillator Chain
Authors:
Fumihiro Ishikawa,
Synge Todo
Abstract:
We present a simple and microscopic physical model that breaks the ergodicity. Our model consists of coupled classical harmonic oscillators, and the motion of the tagged particle obeys the generalized Langevin equation satisfying the second fluctuation dissipation theorem. It is found that although the nonergodicity strength, which is expected to detect the ergodicity breaking, for this model vani…
▽ More
We present a simple and microscopic physical model that breaks the ergodicity. Our model consists of coupled classical harmonic oscillators, and the motion of the tagged particle obeys the generalized Langevin equation satisfying the second fluctuation dissipation theorem. It is found that although the nonergodicity strength, which is expected to detect the ergodicity breaking, for this model vanishes, the velocity auto correlation function of the tagged particle asymptotically oscillates. We analyze the model by using the molecular dynamics and the exact diagonalization as well as the rigorous mapping to the generalized Langevin equation. Our analysis reveals that the asymptotic oscillation is caused by a localized mode with an isolated frequency from the continuous phonon spectrum.
△ Less
Submitted 27 December, 2018; v1 submitted 8 May, 2018;
originally announced May 2018.
-
Optical Bistability in a Low Photon-Density Regime
Authors:
Tatsuhiko Shirai,
Synge Todo,
Hans de Raedt,
Seiji Miyashita
Abstract:
We give a microscopic description of the optical bistability, where the transmission coefficient has two different values as a function of input light intensity, and the system exhibits a discontinuous jump with a hysteresis loop. We developed an efficient numerical algorithm to treat the quantum master equation for hybridized systems of many photons and a large number of two-level atoms. By using…
▽ More
We give a microscopic description of the optical bistability, where the transmission coefficient has two different values as a function of input light intensity, and the system exhibits a discontinuous jump with a hysteresis loop. We developed an efficient numerical algorithm to treat the quantum master equation for hybridized systems of many photons and a large number of two-level atoms. By using this method, we characterize the bistability from the viewpoint of eigenmodes and eigenvalues of the time evolution operator of the quantum master equation. We investigate the optical bistability within the low photon-density regime, where the hybridization of photon and atom degrees of freedom occurs and the resonance spectrum has a double peak structure. We compared it with the standard optical bistability between the low photon-density regime and the high photon-density regime, where the photons can be treated as a classical electromagnetic field and the resonance spectrum has a single peak structure. We discuss the steady-state properties of the optical bistability: dependencies of the photon number density on the intensity and the double peak structure of the photon number distribution inside the bistable region. As for the dynamical properties, we find that the relaxation timescale shows an exponential growth with the system size, and reveal how the hysteresis loop of the optical bistability depends on the size of the system and the sweeping rate of the driving amplitude. Finally, by investigating the effects of detuning frequency of the input field, we clarify the characteristic properties of the present optical bistability within the low photon-density regime, which are qualitatively different from the standard optical bistable phenomena.
△ Less
Submitted 25 April, 2018;
originally announced April 2018.
-
Magnetization process of the S = 1/2 two-leg organic spin-ladder compound BIP-BNO
Authors:
Kazuya Nomura,
Yasuhiro H. Matsuda,
Yasuo Narumi,
Koichi Kindo,
Shojiro Takeyama,
Yuko Hosokoshi,
Toshio Ono,
Naoya Hasegawa,
Hidemaro Suwa,
Synge Todo
Abstract:
We have measured the magnetization of the organic compound BIP-BNO (3,5'-bis(N-tert-butylaminoxyl)-3',5-dibromobiphenyl) up to 76 T where the magnetization is saturated. The S = 1/2 antiferromagnetic Heisenberg two-leg spin-ladder model accounts for the obtained experimental data regarding the magnetization curve, which is clarified using the quantum Monte Carlo method. The exchange constants on t…
▽ More
We have measured the magnetization of the organic compound BIP-BNO (3,5'-bis(N-tert-butylaminoxyl)-3',5-dibromobiphenyl) up to 76 T where the magnetization is saturated. The S = 1/2 antiferromagnetic Heisenberg two-leg spin-ladder model accounts for the obtained experimental data regarding the magnetization curve, which is clarified using the quantum Monte Carlo method. The exchange constants on the rung and the side rail of the ladder are estimated to be J(rung)/kB = 65.7 K and J(leg)/kB = 14.1 K, respectively, deeply in the strong coupling region: J(rung)/J(leg) > 1.
△ Less
Submitted 24 June, 2017;
originally announced June 2017.
-
Crystal Structure Prediction Supported by Incomplete Experimental Data
Authors:
Naoto Tsujimoto,
Daiki Adachi,
Ryosuke Akashi,
Synge Todo,
Shinji Tsuneyuki
Abstract:
The prediction of material structure from chemical composition has been a long-standing challenge in natural science. Although there have been various methodological developments and successes with computer simulations, the prediction of crystal structures comprising more than several tens of atoms in the unit cell still remains difficult due to the many degrees of freedom, which increase exponent…
▽ More
The prediction of material structure from chemical composition has been a long-standing challenge in natural science. Although there have been various methodological developments and successes with computer simulations, the prediction of crystal structures comprising more than several tens of atoms in the unit cell still remains difficult due to the many degrees of freedom, which increase exponentially with the number of atoms. Here we show that when some experimental data is available, even if it is totally insufficient for conventional structure analysis, it can be utilized to support and substantially accelerate structure simulation. In particular, we formulate a cost function based on a weighted sum of interatomic potential energies and a penalty function referred to as "crystallinity", which is defined using limited X-ray diffraction data. This method is applied to well-known polymorphs of $\rm{SiO_2}$ with up to 96 atoms in the simulation cell to find that it reproduces the correct structures efficiently with a very limited number of diffraction peaks. The penalty function is confirmed to destabilize the local minima of the potential energy surface, which facilitates finding the correct structure. This method opens a new avenue for determining and predicting structures that are difficult to determine by conventional methods, such as surface, interface, glass, and amorphous structures.
△ Less
Submitted 24 May, 2017;
originally announced May 2017.
-
Quantum Lattice Model Solver ${\mathcal H}Φ$
Authors:
Mitsuaki Kawamura,
Kazuyoshi Yoshimi,
Takahiro Misawa,
Youhei Yamaji,
Synge Todo,
Naoki Kawashima
Abstract:
${\mathcal H}Φ$ [$aitch$-$phi$] is a program package based on the Lanczos-type eigenvalue solution applicable to a broad range of quantum lattice models, i.e., arbitrary quantum lattice models with two-body interactions, including the Heisenberg model, the Kitaev model, the Hubbard model and the Kondo-lattice model. While it works well on PCs and PC-clusters, ${\mathcal H}Φ…
▽ More
${\mathcal H}Φ$ [$aitch$-$phi$] is a program package based on the Lanczos-type eigenvalue solution applicable to a broad range of quantum lattice models, i.e., arbitrary quantum lattice models with two-body interactions, including the Heisenberg model, the Kitaev model, the Hubbard model and the Kondo-lattice model. While it works well on PCs and PC-clusters, ${\mathcal H}Φ$ also runs efficiently on massively parallel computers, which considerably extends the tractable range of the system size. In addition, unlike most existing packages, ${\mathcal H}Φ$ supports finite-temperature calculations through the method of thermal pure quantum (TPQ) states. In this paper, we explain theoretical background and user-interface of ${\mathcal H}Φ$. We also show the benchmark results of ${\mathcal H}Φ$ on supercomputers such as the K computer at RIKEN Advanced Institute for Computational Science (AICS) and SGI ICE XA (Sekirei) at the Institute for the Solid State Physics (ISSP).
△ Less
Submitted 12 March, 2017; v1 submitted 10 March, 2017;
originally announced March 2017.
-
Updated Core Libraries of the ALPS Project
Authors:
A. Gaenko,
A. E. Antipov,
G. Carcassi,
T. Chen,
X. Chen,
Q. Dong,
L. Gamper,
J. Gukelberger,
R. Igarashi,
S. Iskakov,
M. Könz,
J. P. F. LeBlanc,
R. Levy,
P. N. Ma,
J. E. Paki,
H. Shinaoka,
S. Todo,
M. Troyer,
E. Gull
Abstract:
The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation code, and the applications strive to make commo…
▽ More
The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation code, and the applications strive to make commonly used and proven computational algorithms available to a non-expert community. In this paper we present an updated and refactored version of the core ALPS libraries geared at the computational physics software development community, rewritten with focus on documentation, ease of installation, and software maintainability.
△ Less
Submitted 15 December, 2016; v1 submitted 13 September, 2016;
originally announced September 2016.
-
Upper and Lower Critical Decay Exponents of Ising Ferromagnets with Long-range Interaction
Authors:
Toshiki Horita,
Hidemaro Suwa,
Synge Todo
Abstract:
We investigate the universality class of the finite-temperature phase transition of the two-dimensional Ising model with the algebraically decaying ferromagnetic long-range interaction, $J_{ij} = |\vec{r}_i -\vec{r}_j|^{-(d+σ)}$, where $d$ (=2) is the dimension of the system and $σ$ the decay exponent, by means of the order-$N$ cluster-algorithm Monte Carlo method. In particular, we focus on the u…
▽ More
We investigate the universality class of the finite-temperature phase transition of the two-dimensional Ising model with the algebraically decaying ferromagnetic long-range interaction, $J_{ij} = |\vec{r}_i -\vec{r}_j|^{-(d+σ)}$, where $d$ (=2) is the dimension of the system and $σ$ the decay exponent, by means of the order-$N$ cluster-algorithm Monte Carlo method. In particular, we focus on the upper and lower critical decay exponents, the boundaries between the mean-field-universality, intermediate, and short-range-universality regimes. At the critical decay exponents, it is found that the critical amplitude of the standard Binder ratio of magnetization exhibits the extremely slow convergence as a function of the system size. We propose more effective physical quantities, the combined Binder ratio and the self-combined Binder ratio, both of which cancel the leading finite-size corrections of the conventional Binder ratio. Utilizing these techniques, we clearly demonstrate that in two dimensions the lower and upper critical decay exponents are $σ= 1$ and 7/4, respectively, contrary to the recent Monte Carlo and the renormalization-group studies [M. Picco, arXiv:1207.1018; T. Blanchard, et al., Europhys. Lett. 101, 56003 (2013)].
△ Less
Submitted 27 September, 2017; v1 submitted 31 May, 2016;
originally announced May 2016.
-
Z_N Berry phase and symmetry protected topological phases of SU(N) antiferromagnetic Heisenberg chain
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
The local Z_N quantized Berry phase for the SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local Z_2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is effective as a tool to characterize and classify the various symmetry protected topological phases of one-di…
▽ More
The local Z_N quantized Berry phase for the SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local Z_2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is effective as a tool to characterize and classify the various symmetry protected topological phases of one-dimensional SU(N) spin systems. We extend the path-integral quantum Monte Carlo method for the Z_2 Berry phase in order to calculate the ZN Berry phase numerically. We demonstrate our method by calculating the Z_4 Berry phase for the bond-alternating SU(4) antiferromagnetic Heisenberg chain, which is represented by the Young diagram of four columns.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
-
Stochastic approximation of dynamical exponent at quantum critical point
Authors:
Shinya Yasuda,
Hidemaro Suwa,
Synge Todo
Abstract:
We have developed a unified finite-size scaling method for quantum phase transitions that requires no prior knowledge of the dynamical exponent $z$. During a quantum Monte Carlo simulation, the temperature is automatically tuned by the Robbins-Monro stochastic approximation method, being proportional to the lowest gap of the finite-size system. The dynamical exponent is estimated in a straightforw…
▽ More
We have developed a unified finite-size scaling method for quantum phase transitions that requires no prior knowledge of the dynamical exponent $z$. During a quantum Monte Carlo simulation, the temperature is automatically tuned by the Robbins-Monro stochastic approximation method, being proportional to the lowest gap of the finite-size system. The dynamical exponent is estimated in a straightforward way from the system-size dependence of the temperature. As a demonstration of our novel method, the two-dimensional $S=1/2$ quantum $XY$ model in uniform and staggered magnetic fields is investigated in the combination of the world-line quantum Monte Carlo worm algorithm. In the absence of the uniform magnetic field, we obtain the fully consistent result with the Lorentz invariance at the quantum critical point, $z=1$, i.e., the three-dimensional classical $XY$ universality class. Under a finite uniform magnetic field, on the other hand, the dynamical exponent becomes two, and the mean-field universality with effective dimension $(2+2)$ governs the quantum phase transition.
△ Less
Submitted 16 June, 2015;
originally announced June 2015.
-
Thermal Phase Transition of Generalized Heisenberg Models for SU(N) Spins on Square and Honeycomb Lattices
Authors:
Takafumi Suzuki,
Kenji Harada,
Haruhiko Matsuo,
Synge Todo,
Naoki Kawashima
Abstract:
We investigate thermal phase transitions to a valence-bond solid phase in SU(N) Heisenberg models with four- or six-body interactions on a square or honeycomb lattice, respectively. In both cases, a thermal phase transition occurs that is accompanied by rotational symmetry breaking of the lattice. We perform quantum Monte Carlo calculations in order to clarify the critical properties of the models…
▽ More
We investigate thermal phase transitions to a valence-bond solid phase in SU(N) Heisenberg models with four- or six-body interactions on a square or honeycomb lattice, respectively. In both cases, a thermal phase transition occurs that is accompanied by rotational symmetry breaking of the lattice. We perform quantum Monte Carlo calculations in order to clarify the critical properties of the models. The estimated critical exponents indicate that the universality classes of the square- and honeycomb-lattice cases are identical to those of the classical $XY$ model with a $Z_4$ symmetry-breaking field and the 3-state Potts model, respectively. In the square-lattice case, the thermal exponent, $ν$, monotonically increases as the system approaches the quantum critical point, while the values of the critical exponents, $η$ and $γ/ν$, remain constant. From a finite-size scaling analysis, we find that the system exhibits weak universality, because the $Z_4$ symmetry-breaking field is always marginal. In contrast, $ν$ in the honeycomb-lattice case exhibits a constant value, even in the vicinity of the quantum critical point, because the $Z_3$ field remains relevant in the SU(3) and SU(4) cases.
△ Less
Submitted 23 May, 2015;
originally announced May 2015.
-
Temperature-Dependence of Magnetically-Active Charge Excitations in Magnetite across the Verwey Transition
Authors:
M. Taguchi,
A. Chainani,
S. Ueda,
M. Matsunami,
Y. Ishida,
R. Eguchi,
S. Tsuda,
Y. Takata,
M. Yabashi,
K. Tamasaku,
Y. Nishino,
T. Ishikawa,
H. Daimon,
S. Todo,
H. Tanaka,
M. Oura,
Y. Senba,
H. Ohashi,
S. Shin
Abstract:
We have studied the electronic structure of bulk single crystals and epitaxial films of magnetite Fe$_3$O$_4$. Fe $2p$ core-level spectra show clear differences between hard x-ray (HAX-) and soft x-ray (SX-) photoemission spectroscopy (PES), indicative of surface effects. The bulk-sensitive spectra exhibit temperature ($T$)-dependent charge excitations across the Verwey transition at $T_V$=122 K,…
▽ More
We have studied the electronic structure of bulk single crystals and epitaxial films of magnetite Fe$_3$O$_4$. Fe $2p$ core-level spectra show clear differences between hard x-ray (HAX-) and soft x-ray (SX-) photoemission spectroscopy (PES), indicative of surface effects. The bulk-sensitive spectra exhibit temperature ($T$)-dependent charge excitations across the Verwey transition at $T_V$=122 K, which is missing in the surface-sensitive spectra. An extended impurity Anderson model full-multiplet analysis reveals roles of the three distinct Fe-species (A-Fe$^{3+}$, B-Fe$^{2+}$, B-Fe$^{3+}$) below $T_V$ for the Fe $2p$ spectra, and its $T-$dependent evolution. The Fe $2p$ HAXPES spectra show a clear magnetic circular dichroism (MCD) in the metallic phase of magnetized 100-nm-thick films. The model calculations also reproduce the MCD and identify the magnetically distinct sites associated with the charge excitations. Valence band HAXPES shows finite density of states at $E_F$ for the polaronic metal with remnant order above $T_V$, and a clear gap formation below $T_V$. The results indicate that the Verwey transition is driven by changes in the strongly correlated and magnetically active B-Fe$^{2+}$ and B-Fe$^{3+}$ electronic states, consistent with resistivity and bulk-sensitive optical spectra.
△ Less
Submitted 23 November, 2015; v1 submitted 21 April, 2015;
originally announced April 2015.
-
Generalized Moment Method for Gap Estimation and Quantum Monte Carlo Level Spectroscopy
Authors:
Hidemaro Suwa,
Synge Todo
Abstract:
We formulate a convergent sequence for the energy gap estimation in the worldline quantum Monte Carlo method. The ambiguity left in the conventional gap calculation for quantum systems is eliminated. Our estimation will be unbiased in the low-temperature limit and also the error bar is reliably estimated. The level spectroscopy from quantum Monte Carlo data is developed as an application of the un…
▽ More
We formulate a convergent sequence for the energy gap estimation in the worldline quantum Monte Carlo method. The ambiguity left in the conventional gap calculation for quantum systems is eliminated. Our estimation will be unbiased in the low-temperature limit and also the error bar is reliably estimated. The level spectroscopy from quantum Monte Carlo data is developed as an application of the unbiased gap estimation. From the spectral analysis, we precisely determine the Kosterlitz-Thouless quantum phase-transition point of the spin-Peierls model. It is stablished that the quantum phonon with a finite frequency is essential to the critical theory governed by the antiadiabatic limit, i.e., the $k=1$ SU(2) Wess-Zumino-Witten model.
△ Less
Submitted 30 May, 2015; v1 submitted 31 January, 2014;
originally announced February 2014.
-
Geometric Allocation Approaches in Markov Chain Monte Carlo
Authors:
Synge Todo,
Hidemaro Suwa
Abstract:
The Markov chain Monte Carlo method is a versatile tool in statistical physics to evaluate multi-dimensional integrals numerically. For the method to work effectively, we must consider the following key issues: the choice of ensemble, the selection of candidate states, the optimization of transition kernel, algorithm for choosing a configuration according to the transition probabilities. We show t…
▽ More
The Markov chain Monte Carlo method is a versatile tool in statistical physics to evaluate multi-dimensional integrals numerically. For the method to work effectively, we must consider the following key issues: the choice of ensemble, the selection of candidate states, the optimization of transition kernel, algorithm for choosing a configuration according to the transition probabilities. We show that the unconventional approaches based on the geometric allocation of probabilities or weights can improve the dynamics and scaling of the Monte Carlo simulation in several aspects. Particularly, the approach using the irreversible kernel can reduce or sometimes completely eliminate the rejection of trial move in the Markov chain. We also discuss how the space-time interchange technique together with Walker's method of aliases can reduce the computational time especially for the case where the number of candidates is large, such as models with long-range interactions.
△ Less
Submitted 24 October, 2013;
originally announced October 2013.
-
Path-integral Monte Carlo for the Gauge-fixed Berry Connection and the Local Z_2 Berry Phase
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
We present a general formula of the gauge-fixed Berry connection which can be evaluated by path- integral Monte Carlo method. We also propose that the gauge-fixed local Berry connection can be an- other effective tool to estimate precisely the quantum critical point. For a demonstration, we calculate the gauge-fixed Berry connection and the local Z2 Berry phase of the antiferromagnetic Heisenberg…
▽ More
We present a general formula of the gauge-fixed Berry connection which can be evaluated by path- integral Monte Carlo method. We also propose that the gauge-fixed local Berry connection can be an- other effective tool to estimate precisely the quantum critical point. For a demonstration, we calculate the gauge-fixed Berry connection and the local Z2 Berry phase of the antiferromagnetic Heisenberg model on a staggered bond-alternating ladder, and estimated quantum critical point is consistent with other methods.
△ Less
Submitted 1 October, 2013;
originally announced October 2013.
-
Monte Carlo simulation with aspect ratio optimization: Anomalous anisotropic scaling in dimerized antiferromagnet
Authors:
Shinya Yasuda,
Synge Todo
Abstract:
We present a method that optimizes the aspect ratio of a spatially anisotropic quantum lattice model during the quantum Monte Carlo simulation, and realizes the virtually isotropic lattice automatically. The anisotropy is removed by using the Robbins-Monro algorithm based on the correlation length in each direction. The method allows for comparing directly the value of critical amplitude among dif…
▽ More
We present a method that optimizes the aspect ratio of a spatially anisotropic quantum lattice model during the quantum Monte Carlo simulation, and realizes the virtually isotropic lattice automatically. The anisotropy is removed by using the Robbins-Monro algorithm based on the correlation length in each direction. The method allows for comparing directly the value of critical amplitude among different anisotropic models, and identifying the universality more precisely. We apply our method to the staggered dimer antiferromagnetic Heisenberg model and demonstrate that the apparent non-universal behavior is attributed mainly to the strong size correction of the effective aspect ratio due to the existence of the cubic interaction.
△ Less
Submitted 17 July, 2013;
originally announced July 2013.
-
Possibility of Deconfined Criticality in SU(N) Heisenberg Models at Small N
Authors:
Kenji Harada,
Takafumi Suzuki,
Tsuyoshi Okubo,
Haruhiko Matsuo,
Jie Lou,
Hiroshi Watanabe,
Synge Todo,
Naoki Kawashima
Abstract:
To examine the validity of the scenario of the deconfined critical phenomena, we carry out a quantum Monte Carlo simulation for the SU($N$) generalization of the Heisenberg model with four-body and six-body interactions. The quantum phase transition between the SU($N$) Néel and valence-bond solid phases is characterized for $N=2,3,$ and $4$ on the square and honeycomb lattices. While finite-size s…
▽ More
To examine the validity of the scenario of the deconfined critical phenomena, we carry out a quantum Monte Carlo simulation for the SU($N$) generalization of the Heisenberg model with four-body and six-body interactions. The quantum phase transition between the SU($N$) Néel and valence-bond solid phases is characterized for $N=2,3,$ and $4$ on the square and honeycomb lattices. While finite-size scaling analysis works well up to the maximum lattice size ($L=256$) and indicates the continuous nature of the phase transition, a clear systematic change towards the first-order transition is observed in the estimates of the critical exponent $y \equiv 1/ν$ as the system size increases. We also confirm the relevance of a squared valence-bond solid field $Ψ^2$ for the SU(3) model.
△ Less
Submitted 15 December, 2013; v1 submitted 1 July, 2013;
originally announced July 2013.
-
Parallelized Quantum Monte Carlo Algorithm with Nonlocal Worm Updates
Authors:
Akiko Masaki-Kato,
Takafumi Suzuki,
Kenji Harada,
Synge Todo,
Naoki Kawashima
Abstract:
Based on the worm algorithm in the path-integral representation, we propose a general quantum Monte Carlo algorithm suitable for parallelizing on a distributed-memory computer by domain decomposition. Of particular importance is its application to large lattice systems of bosons and spins. A large number of worms are introduced and its population is controlled by a fictitious transverse field. For…
▽ More
Based on the worm algorithm in the path-integral representation, we propose a general quantum Monte Carlo algorithm suitable for parallelizing on a distributed-memory computer by domain decomposition. Of particular importance is its application to large lattice systems of bosons and spins. A large number of worms are introduced and its population is controlled by a fictitious transverse field. For a benchmark, we study the size-dependence of the Bose-condensation order parameter of the hardcore Bose-Hubbard model with $L\times L\times βt = 10240\times 10240\times 16$, using 3200 computing cores, which shows good parallelization efficiency.
△ Less
Submitted 11 April, 2014; v1 submitted 1 July, 2013;
originally announced July 2013.
-
Long-Range Order of the Three-Sublattice Structure in the S = 1 Heisenberg Antiferromagnet on a Spatially Anisotropic Triangular Lattice
Authors:
Hiroki Nakano,
Synge Todo,
Toru Sakai
Abstract:
We study the S=1 Heisenberg antiferromagnet on a spatially anisotropic triangular lattice by the numerical diagonalization method. We examine the stability of the long-range order of a three-sublattice structure observed in the isotropic system between the isotropic case and the case of isolated one-dimensional chains. It is found that the long-range-ordered ground state with this structure exists…
▽ More
We study the S=1 Heisenberg antiferromagnet on a spatially anisotropic triangular lattice by the numerical diagonalization method. We examine the stability of the long-range order of a three-sublattice structure observed in the isotropic system between the isotropic case and the case of isolated one-dimensional chains. It is found that the long-range-ordered ground state with this structure exists in the range of 0.7 \simle J_2/J_1 \le 1, where J_1 is the interaction amplitude along the chains and J_2 is the amplitude of other interactions.
△ Less
Submitted 26 February, 2013;
originally announced February 2013.
-
Path-integral Monte Carlo method for the local Z_2 Berry phase
Authors:
Yuichi Motoyama,
Synge Todo
Abstract:
We present a loop cluster algorithm Monte Carlo method for calculating the local Z_2 Berry phase of the quantum spin models. The Berry connection, which is given as the inner product of two ground states with different local twist angles, is expressed as a Monte Carlo average on the worldlines with fixed spin configurations at the imaginary-time boundaries. The "complex weight problem" caused by t…
▽ More
We present a loop cluster algorithm Monte Carlo method for calculating the local Z_2 Berry phase of the quantum spin models. The Berry connection, which is given as the inner product of two ground states with different local twist angles, is expressed as a Monte Carlo average on the worldlines with fixed spin configurations at the imaginary-time boundaries. The "complex weight problem" caused by the local twist is solved by adopting the meron cluster algorithm. We present the results of simulation on the antiferromagnetic Heisenberg model on an out-phase bond-alternating ladder to demonstrate that our method successfully detects the change in the valence bond pattern at the quantum phase transition point. We also propose that the gauge-fixed local Berry connection can be an effective tool to estimate pricisely the quantum critical point.
△ Less
Submitted 7 February, 2013; v1 submitted 14 November, 2012;
originally announced November 2012.
-
General Construction of Irreversible Kernel in Markov Chain Monte Carlo
Authors:
Hidemaro Suwa,
Synge Todo
Abstract:
The Markov chain Monte Carlo update method to construct an irreversible kernel has been reviewed and extended to general state spaces. The several convergence conditions of the Markov chain were discussed. The alternative methods to the Gibbs sampler and the Metropolis-Hastings algorithm were proposed and assessed in some models. The distribution convergence and the sampling efficiency are signifi…
▽ More
The Markov chain Monte Carlo update method to construct an irreversible kernel has been reviewed and extended to general state spaces. The several convergence conditions of the Markov chain were discussed. The alternative methods to the Gibbs sampler and the Metropolis-Hastings algorithm were proposed and assessed in some models. The distribution convergence and the sampling efficiency are significantly improved in the Potts model, the bivariate Gaussian model, and so on. This approach using the irreversible kernel can be applied to any Markov chain Monte Carlo sampling and it is expected to improve the efficiency in general.
△ Less
Submitted 1 July, 2012;
originally announced July 2012.
-
Ordering and Excitation in Orbital Compass Model on a Checkerboard Lattice
Authors:
Joji Nasu,
Synge Todo,
Sumio Ishihara
Abstract:
We study an orbital compass model on a checkerboard lattice where orbital degree of freedom is represented by the pseudo-spin operator. Competition arises from an Ising interaction for the $z$ component of pseudo-spins along the vertical/horizontal bonds and an Ising interaction for the $x$ component along diagonal bonds. Classical and quantum compass models are analyzed by utilizing several analy…
▽ More
We study an orbital compass model on a checkerboard lattice where orbital degree of freedom is represented by the pseudo-spin operator. Competition arises from an Ising interaction for the $z$ component of pseudo-spins along the vertical/horizontal bonds and an Ising interaction for the $x$ component along diagonal bonds. Classical and quantum compass models are analyzed by utilizing several analytical methods and numerical simulations. At a fully frustrated point where the two Ising interactions compete with each other, a macroscopic number of orbital configurations are degenerate in a classical ground state. This degeneracy is lifted by thermal and quantum fluctuations, and a staggered long-range order of the $z$ component of the pseudo-spin is realized. A tricritical point for this order appears due to competition between the bond dependent Ising interactions. Roles of geometrical frustration on excitation dynamics are also examined.
△ Less
Submitted 16 March, 2012;
originally announced March 2012.
-
Critical temperature and correlation length of an elastic interaction model for spin-crossover materials
Authors:
Taro Nakada,
Takashi Mori,
Seiji Miyashita,
Masamichi Nishino,
Synge Todo,
William Nicolazzi,
Per Arne Rikvold
Abstract:
It has previously been pointed out that the coexistence of infinite-range and short-range interactions causes a system to have a phase transition of the mean-field universality class, in which the cluster size is finite even at the critical point. In the present paper, we study this property in a model of bistable molecules, whose size changes depending on the bistable states. The molecules can mo…
▽ More
It has previously been pointed out that the coexistence of infinite-range and short-range interactions causes a system to have a phase transition of the mean-field universality class, in which the cluster size is finite even at the critical point. In the present paper, we study this property in a model of bistable molecules, whose size changes depending on the bistable states. The molecules can move in space, interacting via an elastic interaction. It is known that due to the different sizes, an effective long-range interaction between the spins appears, and thus this model has a mean-field type of phase transition. It is found that the scaling properties of the shift of the critical temperature from the pure short-range limit in the model with infinite-range and short-range interactions hold also in the present model, regarding the ratio of the size of the two states as a control parameter for the strength of the long-range interaction. By studying the structure factor, it is shown that the dependence of the cluster size at the critical temperature also shows the same scaling properties as a previously studied model with both infinite-range and short-range interactions. We therefore conclude that these scaling relations hold universally in hybrid models with both short-range and weak long-range interactions.
△ Less
Submitted 9 February, 2012; v1 submitted 28 October, 2011;
originally announced October 2011.
-
Geometric Allocation Approach for Transition Kernel of Markov Chain
Authors:
Hidemaro Suwa,
Synge Todo
Abstract:
We introduce a new geometric approach that constructs a transition kernel of Markov chain. Our method always minimizes the average rejection rate and even reduce it to zero in many relevant cases, which cannot be achieved by conventional methods, such as the Metropolis-Hastings algorithm or the heat bath algorithm (Gibbs sampler). Moreover, the geometric approach makes it possible to find not only…
▽ More
We introduce a new geometric approach that constructs a transition kernel of Markov chain. Our method always minimizes the average rejection rate and even reduce it to zero in many relevant cases, which cannot be achieved by conventional methods, such as the Metropolis-Hastings algorithm or the heat bath algorithm (Gibbs sampler). Moreover, the geometric approach makes it possible to find not only a reversible but also an irreversible solution of rejection-free transition probabilities. This is the first versatile method that can construct an irreversible transition kernel in general cases. We demonstrate that the autocorrelation time (asymptotic variance) of the Potts model becomes more than 6 times as short as that by the conventional Metropolis-Hastings algorithm. Our algorithms are applicable to almost all kinds of Markov chain Monte Carlo methods and will improve the efficiency.
△ Less
Submitted 1 July, 2012; v1 submitted 17 June, 2011;
originally announced June 2011.
-
The ALPS project release 2.0: Open source software for strongly correlated systems
Authors:
B. Bauer,
L. D. Carr,
H. G. Evertz,
A. Feiguin,
J. Freire,
S. Fuchs,
L. Gamper,
J. Gukelberger,
E. Gull,
S. Guertler,
A. Hehn,
R. Igarashi,
S. V. Isakov,
D. Koop,
P. N. Ma,
P. Mates,
H. Matsuo,
O. Parcollet,
G. Pawlowski,
J. D. Picon,
L. Pollet,
E. Santos,
V. W. Scarola,
U. Schollwöck,
C. Silva
, et al. (7 additional authors not shown)
Abstract:
We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries…
▽ More
We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.
△ Less
Submitted 23 May, 2011; v1 submitted 13 January, 2011;
originally announced January 2011.
-
Markov Chain Monte Carlo Method without Detailed Balance
Authors:
Hidemaro Suwa,
Synge Todo
Abstract:
We present a specific algorithm that generally satisfies the balance condition without imposing the detailed balance in the Markov chain Monte Carlo. In our algorithm, the average rejection rate is minimized, and even reduced to zero in many relevant cases. The absence of the detailed balance also introduces a net stochastic flow in a configuration space, which further boosts up the convergence. W…
▽ More
We present a specific algorithm that generally satisfies the balance condition without imposing the detailed balance in the Markov chain Monte Carlo. In our algorithm, the average rejection rate is minimized, and even reduced to zero in many relevant cases. The absence of the detailed balance also introduces a net stochastic flow in a configuration space, which further boosts up the convergence. We demonstrate that the autocorrelation time of the Potts model becomes more than 6 times shorter than that by the conventional Metropolis algorithm. Based on the same concept, a bounce-free worm algorithm for generic quantum spin models is formulated as well.
△ Less
Submitted 13 October, 2010; v1 submitted 14 July, 2010;
originally announced July 2010.
-
Orbital ordering in e_g orbital systems: Ground states and thermodynamics of the 120 degree model
Authors:
Andre van Rynbach,
Synge Todo,
Simon Trebst
Abstract:
Orbital degrees of freedom shape many of the properties of a wide class of Mott insulating, transition metal oxides with partially filled 3d-shells. Here we study orbital ordering transitions in systems where a single electron occupies the e_g orbital doublet and the spatially highly anisotropic orbital interactions can be captured by an orbital-only model, often called the 120 degree model. Our a…
▽ More
Orbital degrees of freedom shape many of the properties of a wide class of Mott insulating, transition metal oxides with partially filled 3d-shells. Here we study orbital ordering transitions in systems where a single electron occupies the e_g orbital doublet and the spatially highly anisotropic orbital interactions can be captured by an orbital-only model, often called the 120 degree model. Our analysis of both the classical and quantum limits of this model in an extended parameter space shows that the 120 degree model is in close proximity to several T=0 phase transitions and various competing ordered phases. We characterize the orbital order of these nearby phases and their associated thermal phase transitions by extensive numerical simulations and perturbative arguments.
△ Less
Submitted 23 June, 2010;
originally announced June 2010.
-
Direct observation of localization in the minority-spin-band electrons of magnetite below the Verwey temperature
Authors:
H. Kobayashi,
T. Nagao,
M. Itou,
S. Todo,
B. Barbiellini,
P. E. Mijnarends,
A. Bansil,
N. Sakai
Abstract:
Two-dimensional spin-uncompensated momentum density distributions, $ρ_{\rm s}^{2D}({\bf p})$s, were reconstructed in magnetite at 12K and 300K from several measured directional magnetic Compton profiles. Mechanical de-twinning was used to overcome severe twinning in the single crystal sample below the Verwey transition. The reconstructed $ρ_{\rm s}^{2D}({\bf p})$ in the first Brillouin zone chan…
▽ More
Two-dimensional spin-uncompensated momentum density distributions, $ρ_{\rm s}^{2D}({\bf p})$s, were reconstructed in magnetite at 12K and 300K from several measured directional magnetic Compton profiles. Mechanical de-twinning was used to overcome severe twinning in the single crystal sample below the Verwey transition. The reconstructed $ρ_{\rm s}^{2D}({\bf p})$ in the first Brillouin zone changes from being negative at 300 K to positive at 12 K. This result provides the first clear evidence that electrons with low momenta in the minority spin bands in magnetite are localized below the Verwey transition temperature.
△ Less
Submitted 25 August, 2009;
originally announced August 2009.
-
Successive phase transitions at finite temperatures of the supersolid in the three-dimensional extended Bose-Hubbard model
Authors:
Keisuke Yamamoto,
Synge Todo,
Seiji Miyashita
Abstract:
We study the finite temperature properties of the extended Bose-Hubbard model on a cubic lattice. This model exhibits the so-called supersolid state. To start with, we investigate ordering processes by quantum Monte Carlo simulations, and find successive superfluid and solid phase transitions. There, we find that the two order parameters compete with each other. We obtain the finite temperature…
▽ More
We study the finite temperature properties of the extended Bose-Hubbard model on a cubic lattice. This model exhibits the so-called supersolid state. To start with, we investigate ordering processes by quantum Monte Carlo simulations, and find successive superfluid and solid phase transitions. There, we find that the two order parameters compete with each other. We obtain the finite temperature phase diagram, which contains the superfluid, the solid, the supersolid and the disordered phase. We develop a mean-field theory to analyze the ordering processes and compare the result with that obtained by simulations, and discuss the mechanism of the competition of these two orders. We also study how the supersolid region shrinks as the on-site repulsion becomes strong.
△ Less
Submitted 7 July, 2008;
originally announced July 2008.
-
Improved chain mean-field theory for quasi-one-dimensional quantum magnets
Authors:
Synge Todo,
Akira Shibasaki
Abstract:
A novel mean-field approximation for quasi-one-dimensional (Q1D) quantum magnets is formulated. Our new mean-field approach is based on the Bethe-type effective-field theory, where thermal and quantum fluctuations between the nearest-neighbor chains as well as those in each chain are taken into account exactly. The self-consistent equation for the critical temperature contains the boundary-field…
▽ More
A novel mean-field approximation for quasi-one-dimensional (Q1D) quantum magnets is formulated. Our new mean-field approach is based on the Bethe-type effective-field theory, where thermal and quantum fluctuations between the nearest-neighbor chains as well as those in each chain are taken into account exactly. The self-consistent equation for the critical temperature contains the boundary-field magnetic susceptibilities of a multichain cluster, which can be evaluated accurately by some analytic or numerical methods, such as the powerful quantum Monte Carlo method. We show that the accuracy of the critical temperature of Q1D magnets as a function of the strength of interchain coupling is significantly improved, compared with the conventional chain mean-field theory. It is also demonstrated that our new approximation can predict nontrivial dependence of critical temperature on the sign (i.e., ferromagnetic or antiferromagnetic) of interchain coupling as well as on the impurity concentration in randomly diluted Q1D Heisenberg antiferromagnets.
△ Less
Submitted 11 December, 2008; v1 submitted 20 May, 2008;
originally announced May 2008.
-
Green's function theory for spin-1/2 ferromagnets with an easy-plane exchange anisotropy
Authors:
Daisuke Yamamoto,
Synge Todo,
Susumu Kurihara
Abstract:
The many-body Green's function theory with the random-phase approximation is applied to the study of easy-plane spin-1/2 ferromagnets in an in-plane magnetic field. We demonstrate that the usual procedure, in which only the three Green's functions $<< S_i^μ;S_j^->>$ ($μ=+,-,z$) are used, yields unreasonable results in this case. Then the problem is discussed in more detail by considering all com…
▽ More
The many-body Green's function theory with the random-phase approximation is applied to the study of easy-plane spin-1/2 ferromagnets in an in-plane magnetic field. We demonstrate that the usual procedure, in which only the three Green's functions $<< S_i^μ;S_j^->>$ ($μ=+,-,z$) are used, yields unreasonable results in this case. Then the problem is discussed in more detail by considering all combinations of Green's functions. We can derive one more equation, which cannot be obtained by using only the set of the above three Green's functions, and point out that the two equations contradict each other if one demands that the identities of the spin operators are exactly satisfied. We discuss the cause of the contradiction and attempt to improve the method in a self consistent way. In our procedure, the effect of the anisotropy can be appropriately taken into account, and the results are in good agreement with the quantum Monte Carlo calculations.
△ Less
Submitted 3 September, 2008; v1 submitted 5 February, 2008;
originally announced February 2008.
-
Order-N Cluster Monte Carlo Method for Spin Systems with Long-range Interactions
Authors:
Kouki Fukui,
Synge Todo
Abstract:
An efficient O(N) cluster Monte Carlo method for Ising models with long-range interactions is presented. Our novel algorithm does not introduce any cutoff for interaction range and thus it strictly fulfills the detailed balance. The realized stochastic dynamics is equivalent to that of the conventional Swendsen-Wang algorithm, which requires O(N^2) operations per Monte Carlo sweep if applied to…
▽ More
An efficient O(N) cluster Monte Carlo method for Ising models with long-range interactions is presented. Our novel algorithm does not introduce any cutoff for interaction range and thus it strictly fulfills the detailed balance. The realized stochastic dynamics is equivalent to that of the conventional Swendsen-Wang algorithm, which requires O(N^2) operations per Monte Carlo sweep if applied to long-range interacting models. In addition, it is shown that the total energy and the specific heat can also be measured in O(N) time. We demonstrate the efficiency of our algorithm over the conventional method and the O(N log N) algorithm by Luijten and Bloete. We also apply our algorithm to the classical and quantum Ising chains with inverse-square ferromagnetic interactions, and confirm in a high accuracy that a Kosterlitz-Thouless phase transition, associated with a universal jump in the magnetization, occurs in both cases.
△ Less
Submitted 3 February, 2008;
originally announced February 2008.
-
The ALPS project release 1.3: open source software for strongly correlated systems
Authors:
A. F. Albuquerque,
F. Alet,
P. Corboz,
P. Dayal,
A. Feiguin,
S. Fuchs,
L. Gamper,
E. Gull,
S. Guertler,
A. Honecker,
R. Igarashi,
M. Koerner,
A. Kozhevnikov,
A. Laeuchli,
S. R. Manmana,
M. Matsumoto,
I. P. McCulloch,
F. Michel,
R. M. Noack,
G. Pawlowski,
L. Pollet,
T. Pruschke,
U. Schollwock,
S. Todo,
S. Trebst
, et al. (3 additional authors not shown)
Abstract:
We present release 1.3 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an international open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. Development is centered on common XML and binary data formats, o…
▽ More
We present release 1.3 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an international open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. Development is centered on common XML and binary data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs. The programs enable non-experts to start carrying out numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), as well as the density matrix renormalization group (DMRG). Changes in the new release include a DMRG program for interacting models, support for translation symmetries in the diagonalization programs, the ability to define custom measurement operators, and support for inhomogeneous systems, such as lattice models with traps. The software is available from our web server at http://alps.comp-phys.org/ .
△ Less
Submitted 11 January, 2008;
originally announced January 2008.